Performance Analysis of the Sleep Mode in WiMAX 2 Networks with Multimedia Application

,


Introduction
Wireless mesh networks are undergoing rapid development.They are intended to deliver wireless services for a large variety of applications in personal, local, campus, and metropolitan areas.Some protocol designs for wireless mesh networks have focused on power efficiency mechanisms.IEEE 802.16 standard has been designed for fixed subscriber stations.As an enhancement of the IEEE 802.16 standard, IEEE 802.16e, called WiMAX [1], has improved the original standard supporting mobility so that the Mobile Station (MS) can move during services.Aiming at the next generation mobile WiMAX, called WiMAX 2, IEEE 802.16m [2] is currently being processed for standardization.
There have been many studies analyzing the performance of the sleep mode operations for Types I-III in WiMAX [3][4][5][6][7].In [3], the authors evaluated and compared the sleep mode operations for Type I and Type II power saving classes using the method of an embedded Markov chain.In order to avoid too frequent switching between the sleep state and the awake state, an enhanced power saving class Type III was provided in [4] and the system performances were analyzed for user initiated data packet arrivals.In [5], by increasing the unavailability interval, the authors proposed an enhanced power saving mechanism (ePSM), where both activated Type I and Type II power saving classes exist in an MS.The performance evaluation confirms that an ePSM can save more energy than conventional schemes.
Moreover, in [8], the authors used a two-phase Markovian Arrival Process to investigate the effect of the traffic parameters, such as the correlation parameter in the performance evaluation of a system for WiMAX.The optimal sleep mode strategy was obtained by minimizing the total average system cost.

2
Journal of Applied Mathematics However, an IEEE 802.16m amendment has been drafted to meet the requirements of WiMAX 2. IEEE 802.16m provides the enhanced performance required to support future advanced services and offers a novel sleep mode scheme which can reduce the power consumption and extend the lifetime of a battery-powered Mobile Station (MS) for a multimedia scenario.
WiMAX 2 has recently attracted a lot of research interest.In [9], the performance of the sleep mode in WiMAX 2 was mathematically analyzed, where the sleep cycle was supposedly fixed and the traffic indication (TRF-IND) was periodic.Moreover, the optimal TRF-IND interval was given to minimize the average power consumption of the MS while satisfying the Quality of Service (QoS) on the mean delay.In [10], the authors proposed a power-saving mechanism with binary exponential traffic indications in IEEE 802.16e/m.In [11], an efficient sleep mode operation was proposed by using a T AMS timer.Also, the authors analyzed the proposed scheme by using an embedded Markov chain.The optimal parameters were given to minimize the power consumption while satisfying the QoS requirement on the average message delay.
In [12], taking into account that the listening window can be extended and the sleep cycle length can be adjusted, the authors conducted an analytical study on the power consumption and the average data packet delay to minimize the power consumption while satisfying a user-specified packet delay constraint.In [13], a concise analytical model for the sleep mode operation of WiMAX 2 was proposed.Considering both downlink and uplink traffics, the performance measures such as mean waiting time were derived and the performance comparisons between the WiMAX 2 and WiMAX were conducted with simulations.In [14], considering both downlink traffic and uplink traffic, a queueing model of two servers sharing vacations with close-down time and multiple vacations was built for WiMAX 2. The arrival process was assumed to follow a Bernoulli arrival process.By employing a two-dimensional embedded Markov chain, the authors derived the performance measures.
However, in the research on the sleep mode of WiMAX 2 mentioned above, the authors assumed that a sleep cycle was constituted under the condition that the time length of the extended listening window was to be equal to or less than a sleep cycle.The authors also assumed that the average remainder time of the extended listening window was equal to the difference between the time length of the sleep cycle and the average time length of the first sleep window.Moreover, in order to simplify the analysis procedure, the authors assumed uncorrelated traffic and modeled the incoming traffic using a Poisson process or a Bernoulli process.We know that, in multimedia application, two subsequent slot's states are correlated and the assumption of a Poisson process or a Bernoulli process is not reasonable.
In this paper, we present a comprehensive performance analysis by releasing those assumptions introduced in past research mentioned above for actual systems by considering the sleep mode of WiMAX 2 in a scenario, where the real-time traffic includes the real-time traffic and the best-effort traffic mixed.Considering the correlation coefficient between the two subsequent time slots' states for the real-time traffic, we assume the arrival process to be a Discrete Time Markovian Arrival Process (D-MAP).The motivation for using D-MAP is that it can capture the digital nature of modern communication and model correlated traffic.Then, we build a multiple heterogeneous vacation queueing model to model the stochastic behavior of the sleep mode in the network system.Note that D-MAP can represent a variety of arrival processes, including the Bernoulli arrival process and the Markov Modulated Bernoulli process (MMBP).Therefore, the model presented in this paper is a generalization of the conventional models for analyzing the system performance of the sleep mode in WiMAX 2 with the real-time traffic including the real-time traffic and the best-effort traffic mixed.
By using the method of an embedded Markov chain, we present the probability distribution for the number of data packets.Considering the trade-off between the average response time of data packets and the energy saving ratio, we then develop a cost function to determine the optimal length of the sleep cycle in order to maximize the energy saving ratio while satisfying the QoS constraint on the average response time of data packets.With numerical results, we investigate the influence of the system parameters on the system performance.Finally, concerning the trade-off between different performance measures, we optimize the sleep cycle with the cost function.
The rest of this paper is organized as follows.In Section 2, we describe the system model for the multimedia application in which the real-time traffic includes the real-time and the best-effort traffic mixed.In Section 3, we analyze the system model in the steady state and give the probability distribution for the number of data packets.In Section 4, we obtain the formulas for the average response time of data packets, the energy saving ratio, and the higher moments as the standard deviation for the number of data packets.In addition, we construct a cost function to give the optimal design for the sleep cycle.In Section 5, we provide numerical results to illustrate the influence of the system parameters on the system performance.Finally, we draw our conclusions in Section 6.

System Model
In the sleep mode of WiMAX 2 networks, MS is provided with a set of sleep cycles no matter the system being at the awake period or the sleep period.Except for the initial sleep window, each subsequent sleep cycle consists of both a listening window and a sleep window.As specified in the standard of IEEE 802.16m, the real-time traffic includes the real-time traffic and the best-effort traffic mixed, and the length of the sleep cycle is fixed.Moreover, the value of the default listening window remains without change.Therefore, in this paper, we present the system model to describe the sleep mode in WiMAX 2 networks by setting a fixed time length for the sleep cycle [2].
During the listening window, data packets (if any) will be transmitted in the same way as in the awake period and the listening window will be extended if necessary.Therefore, in this paper, we regard the listening window with data packet  transmission as an extended listening window and the awake period following an extended listening window as a part of that extended listening window.Then the extended listening window is defined as a busy period .Here, the extended listening window is perhaps followed by an awake period.A listening window will terminate on reaching the end of the current listening window including any extensions of the listening window, or reaching the end of the sleep cycle.
We assume the initial sleep window following an extended window to be a vacation period  1 with the time length   1 and the subsequent sleep cycle without data packet transmission to be another vacation period  2 with the time length   2 , respectively.We note that the extended listening window will end at any instant during a sleep cycle; in other words, the initial sleep window will begin at any instant of an extended listening window.For this, the time length   1 of the vacation period  1 is supposed to be one half of the time length   2 of  2 as   2 = (1/2)  2 , where   2 is in fact the time length of the sleep cycle that will be set as required by the system.
We also define a busy cycle  as a time period between the ending instants of two consecutive busy periods.We further define a system vacation  as a time period from the instant where a busy period ends to the instant where the next busy period begins.We can see that the first vacation  1 and one or more subsequent vacations  2 (if any) in a busy cycle  combine to produce a system vacation .
With these assumptions stated above, we can model the system as a multiple heterogeneous vacation queueing model.The state transition of the system model is illustrated in Figure 1.
Figure 1 shows that: (1) if there is no data packet to be transmitted in the buffer of the serving Base Station (BS), the busy period  will be over and the vacation  1 will begin; (2) if there is at least one data packet arrival within the time length   1 of the vacation  1 , the MS will begin a new busy period  when  1 is over.Otherwise, the vacation  2 will begin; (3) if there is at least one data packet arrival within the time length   2 of the vacation  2 , the MS will enter a new busy period  when  2 is over.Otherwise, the MS will begin another vacation  2 again.
In this system model, we assume that the time axis is segmented into a sequence of equal intervals of unit duration, called slots.We also assume that data packets arrive only just before the end of a slot  =  − ,  = 1, 2, . . .and depart only just after the end of a slot  =  + ,  = 2, 3, . ... This is called a late arrival queueing system with delayed access.Various time epochs at which events occur are depicted in Figure 2.
The system consists of a single channel and a finite buffer size .Data packets are supposed to be transmitted according to a First-Come, First-Served (FCFS) strategy.
We define the embedded Markov point as the instant at which a data packet is completely transmitted.We order these embedded points by ( = 1, 2, . ..).
Moreover, we define the arrival phase by the stochastic state of the arrival process and define the state of the system by the number of data packets and the arrival phase at these embedded Markov points.Let   and   be the number of data packets and the arrival phase at the th embedded point.Thus, {(  ,   ),  = 1, 2, . ..} forms a discrete-time two-dimensional embedded Markov chain in a semi-Markov process with state space {(, ), 0 ≤  ≤ , 1 ≤  ≤ }, where  is the system capacity and  is the number of phases in the underlying Markov Chain (UMC). and  are the system parameters that we will give in the numerical results.
It should be noted that, in this paper, the embedded Markov points are chosen at the instant at which a data packet is completely transmitted and the condition of this system model being the stationary state is the number  of data packets in real-time traffic at the embedded Markov points which must be less than the system capacity .With the assumptions and the conditions, in this system model, data packets would not be lost at the embedded Markov points.
Considering the correlation of the real-time traffic, the arrival process is assumed to be a D-MAP with representation {D  ,  = 0, 1}, where D  is an  ×  matrix and  is the number of phases in the UMC, where  = 1, 2, . ... The element [D 0 ]  is the probability that there is a phase transition from  to  without any data packet arrivals, where 1 ≤  ≤ , 1 ≤  ≤ .The element [D 1 ]  is the probability that there is a phase transition from  to  with a data packet arrival.
The matrix D = D 0 + D 1 is the transition probability matrix of the UMC.The matrix D is assumed to be irreducible.
Let π be the stationary vector of size 1 ×  for the UMC; that is, π is a solution for the equations as follows: where e denotes a column one vector of size  × 1.
The average arrival rate of the D-MAP is given as follows: The matrix Probability Generating Function (PGF) of D is defined as follows: The transmission time of a data packet is denoted by . is supposed to follow a general distribution.The probability distribution   , PGF (), and the average value [] of  are given as follows:

Performance Analysis
Moreover, let (), 0 ≤  ≤  be the probability vector of data packet numbers for the embedded Markov chain in a steady state.Therefore (), 0 ≤  ≤  can be given as follows: In order to compute (), 0 ≤  ≤ , some notations are introduced as follows.
Let A  (, ) be a matrix of size  × , whose element [A  (, )]  is the conditional probability that there are  data packet arrivals during the transmission time (lasting  slots) of a data packet, given that the transmission started with arrival phase  and ended with arrival phase .
Let A  () denote a matrix of size  ×  representing that there are  data packet arrivals within the transmission time  of a data packet.A  () can be given as follows: Let A  1 (  , ) be a matrix of size  × , whose element [A  1 (  , )]  is the conditional probability that there are  data packet arrivals during the vacation  1 (lasting   slots), given that the arrival process transfer is from phase  at the beginning instant of the vacation  1 to phase  at the end instant of the vacation  1 .
Let A  2 (  2 , ) be a matrix of size  × , whose element [A  2 (  2 , )]  is the conditional probability that there are  data packet arrivals during the vacation  2 (lasting   2 slots), given that the arrival process transfer is from phase  at the beginning instant of the vacation  2 to phase  at the end instant of the vacation  2 .
Considering the multiple vacations in this paper, let A  () be a matrix of size  ×  with element [A  ()]  representing the probability that there are  data packet arrivals within the system vacation, given that the system vacation  started with the arrival process of phase  and ended with the arrival process of phase .
We calculate A  () in three possible cases as follows.
(1) For the case of  >   2 ,if 1 ≤  ≤   , if if  > 2  , (2) For the case of if if  = , if  > , (3) For the case of   > , if  < , if  = , if  > , Among which, I and 0 are the identity and zero matrices of size  × , respectively.
It may be noted here that in order to know A  () and A  (), we need to compute A  (, ), A  1 (  , ), and A  2 (  2 , ) efficiently.
The transition probability matrix P of the twodimensional Markov process {((), ()),  ≥ 0} is given as follows: ) ) ) ) , where elements of the first row are obtained as follows: The last element for each row, except for the first row, is given as follows: Let  = ((0), (1), . . ., ()) be the stationary probability vector of size 1 × ( + 1) for the transition probability matrix P. Then we can obtain () by solving the equations as follows: Let  be the number of data packets at the embedded Markov points.The average value [] of  can be given as follows: 3.2.Busy Cycle.Let   be the time length of a busy period  and [  ] be the corresponding average value of   .By using the analysis results of the classical Geom/G/1 queueing model [16], [  ] can be given as follows: Let   be the number of sleep windows in a busy cycle .Then we can get the probability distribution of   as follows: where  = 1 −  = 1 − πD 1 e.
Let   be the time length in slots of the system vacation period .The average value [  ] of   is then given as follows: Let   be the time length of the busy cycle .The average value [  ] of   is then given as follows: Combining ( 25) and ( 28), the probability of the MS being awake is given as follows:

Performance Measures and Optimal Design
We define the response time  of a data packet as being the time period in slots that has elapsed from the arrival instant to the end of the transmission of that data packet.This is a very important performance measure for evaluating the user QoS for the system.By applying Little's law [17], the average response time  of data packets is obtained as follows: The energy saving ratio  is defined as the amount of energy conserved per slot, which is one of the most important performance measures for evaluating the energy saving efficiency of the sleep mode in IEEE 802.16m.The energy saving ratio  is given as follows: where  1 is the energy saved per slot for the system being in the sleep period and  2 is the additional energy consumption per slot due to the traffic indications introduced in IEEE 802.16m.We know that the design and development of the system require not only average performance measures but also higher moments, because the output stream from one system often forms the input stream to another.
It is well known that the diffusion approximation of the queueing length in the system buffer is important for evaluating the system performance [18].The diffusion degree of the queueing length can reflect the delay jitter of data packets.Here, we give the standard deviation  of the number of data packets in the system buffer as follows: where [] is given by ( 24) and [ 2 ] is the second moment for the number of data packets at the embedded Markov points given as follows: Obviously, a longer sleep cycle can enhance the energy saving ratio while increasing the average response time of data packets.On the other hand, a shorter sleep cycle can decrease the average response time of data packets while reducing the energy saving ratio.Therefore, it is important to determine an optimal value for the sleep cycle in order to maximize the energy saving ratio while satisfying the QoS constraint on the average response time of data packets.
To this end, we introduce a cost function (  2 ), which is proportional to the average response time  and is inversely proportional to the energy saving ratio .Thus, the cost function is given as follows: where  1 is the factor of the average response time of data packets to the cost and  2 is the factor of the energy saving ratio to the cost.

Numerical Results
In this section, we provide numerical results from analysis and simulation to illustrate the system performance for the sleep mode in WiMAX 2. The simulation results match well with the results from the analysis.As an example, we set system parameters as follows:  = 15, [] = 1, and  = 0.2.Since D-MAP can represent a variety of arrival processes, we decide on a special case of D-MAP, namely, a Bernoulli modulated ON-OFF traffic model, as an example of the system input.The ON-OFF traffic model refers to the arrival process modulated by a two-state Markov chain.When the system is in the ON state, the arrival process follows a Bernoulli process with the parameter .If the system is in the OFF state, no arrival occurs.
We assume that the system is in the ON state with the probability .The correlation coefficient between the two subsequent slot's states is 1 − 1/, where  is the correlation parameter for the average time lengths for the system being in the ON and OFF states, respectively.The mean time length for the system being in the ON state is /(1 − ), and correspondingly / is the average time length for the system being in the OFF state.
Therefore, the matrix PGF D() for the traffic model in this paper can be given as follows: Letting  = 0 in (35), we can derive the formula of D 0 as follows: ) .
(36) Letting  = 1 in (35), we can obtain the formula of D as follows: Combining ( 36) and (37), we can obtain the formula of D 1 as follows: Meanwhile, we assume that the transmission time  of a data packet follows a geometric distribution with the parameter  (0 ≤  ≤ 1).Then the PGF () and the average value [] of  are given as follows: Now, we will discuss the influence of the sleep cycle on the system performance measures for different correlation parameters  with numerical results.Obviously, the smaller the value of correlation parameters  is, the less the correlation degree of the real-time traffic is.Particularly,  = 1 means that there is only the best-effort traffic in the system and  = ∞ means that there is only the real-time traffic when the traffic is mixed.The ratio of the best-effort traffic to the real-time traffic in the mixed traffic is given as 1/.When we let   denote the quantity of the best-effort traffic,   denotes the quantity of the real-time traffic, and   denotes the quantity of the mixed traffic, we can have By referencing [9,10], we set the system parameters for numerical results in Table 1.
The average response time  of data packets versus the sleep cycle   2 with different correlation parameters  is depicted in Figure 3.
It is illustrated that for the same correlation parameter , the longer the time length   2 of the sleep cycle is, the greater the average response time  of data packets is.The reason is that when the time length of the sleep cycle is longer, the system is more likely to be in the sleep period, so the average   response time of data packets is longer.On the other hand, for all the time lengths   2 of the sleep cycle, the larger the correlation parameter  is, the longer the average response time  of data packets will be.That is because the greater the correlation parameter is, the denser the data packet arrivals are and the longer the queueing length is, so the longer the average response time of data packets will be.Moreover, we observe that if the best-effort traffic only is considered in multimedia application, the average response time of data packets will be underevaluated.
The energy saving ratio  versus the time length   2 of the sleep cycle with different correlation parameters  is demonstrated in Figure 4.
It can be observed that for the same correlation parameter , the energy saving ratio  will increase along with the  time length   2 of the sleep cycle.The reason is that when the time length of the sleep cycle is longer, the system is more likely to be in a sleep period.Energy will be saved in the sleep period, so the energy saving ratio will be greater.It can also be found that for all the time lengths   2 of the sleep cycle, the energy saving ratio  will decrease as the correlation parameter  decreases.This is because the bigger the correlation parameter  is, the denser the data packet arrivals are, the busier the MS is, and the less energy will be saved.Moreover, we can find that if the best-effort traffic is only considered for a multimedia scenario, the energy saving effect will be overevaluated.
In Figures 5 and 6, we show how the average response time  of data packets and the energy saving ratio  change versus the correlation parameter  from 1 to 25 with different time lengths   1 of the initial sleep window from 5 to 35.
Moreover, we observe that for the same correlation parameter , with a smaller value for the time length   1 of the initial sleep window, such as   1 = 5, the average response time  and the energy saving ratio  have smaller results.But when the time length   1 of the initial sleep window becomes a higher value, such as   1 = 35, results of the average response time  of data packets and the energy saving ratio  are greater.The reason is that when the time length   1 of the initial sleep window is longer, the system is more likely to be in the sleep period, so the longer the response time of data packets is and the greater the energy saving ratio will be.
We also find that the differences among values of the energy saving ratio  are bigger when the time length   1 is smaller.This is because for a smaller value of the time length   1 , there will be more sleep windows in a busy cycle; therefore a slight increase in the time length   1 will result in a huge increase in the value of the energy saving ratio .On the other hand, when the time length   1 is greater, the system is more likely to transfer from a sleep state to an awake state and the fewer sleep windows will be in a busy cycle.Therefore, an Correlation parameter M Analysis Simulation Energy saving ratio  (mW) Standard deviation  of the number of data packets Arrival ratio  of data packets increase in the time length   1 will result in a gradual increase in the energy saving ratio .
Considering the special case for  = 1, namely, there is only the best-effort traffic, for all the time length   1 , the average response times  of data packets are shorter and the energy saving ratios  are bigger.That is to say, if the correlation of the real-time traffic is neglected, the system performance of the sleep mode with multimedia application in WiMAX 2 networks will be overoptimistically.
The change trend of the standard deviation  of the number of data packets versus the arrival ratio  of data packets is depicted in Figure 7.
In Figure 7, we find that for all the time lengths   1 of the initial sleep cycle, the average transmission time  [𝑆] of data packets and the correlation parameter , the higher the arrival ratio  of data packets is, the larger the standard deviation  of the number of data packets will be.This is because the higher the arrival ratio  of data packets is, the more data packets will arrive within a certain time period.Note that the minimal number of data packets is 1; a greater average number of data packets in the buffer during the extended listening period means a maximal number of data packets and thus the standard deviation  of the number of data packets will be greater.
We notice that, for the same time length   1 of the initial sleep cycle, the correlation parameter , and the arrival ratio  of data packets, the standard deviation  of the number of data packets will increase with an increase in the average transmission time [] of data packets.The reason is that the longer the average transmission time [] of data packets is, the maximal value of the queuing length in the system will be during the extended listening period.So the standard deviation  of the number of data packets will be greater.
In Figure 7, we also find that, for the same time length   1 of the initial sleep cycle, the same average transmission time [] of data packets, and the same arrival ratio  of data packets, the larger the correlation parameter  is, the greater the standard deviation  of the number of data packets will be.For example, the standard deviation  of the number of data packets for the case of  = 20 is much greater than that for the case of  = 16, for each case of   1 .This is because in the  = 20 case, the quantity of the real-time traffic is much greater than that of the best-effort traffic in the mixed traffic, so it is more likely that some of the real-time traffic happens to be waiting in the buffer.This results in the standard deviation  of the number of data packets being greater.
In Figure 7, we also observe that, for the same correlation parameter , the same average transmission time [] of data packets, and the same arrival ratio  of data packets, the longer the time length   1 is, the greater the standard deviation  of the number of data packets will be.This is because the longer the time length   1 is, the greater the maximal number of data packets arriving within the time length   1 will be, so that the number range of the data packets will become wide at the beginning instant of the data packets' transmission.It causes the standard deviation  of the data packets' number to be greater.
Figure 7 clearly shows that the standard deviation  of the number of data packets is heavily dependent on the correlation of the real-time traffic.If we underestimate the correlation degree of the real-time traffic in multimedia application, as in the system models presented in the available literature, the standard deviation  of the number of data packets will be underevaluated.It means that the results from analyzing the system model will be considerably separated from actual systems.
To maximize the energy saving ratio  while satisfying the QoS constraint on the average response time  of data packets, we investigate the changing trend for the cost function (  2 ) of the system given by (34) versus the time length   2 of the sleep cycle with different correlation parameters .
Taking  1 = 0.3 and  2 = 0.5 as an example, the results are shown in Figure 8.  From Figure 8, we see that two stages are experienced by the cost function (  2 ).In the first stage, the cost function (  2 ) will decrease along with an increase in the time length   2 of the sleep cycle.During this stage, the longer the time length of the sleep cycle is, the longer the MS will stay in the sleep period and the more energy will be saved, so the less the cost function will be.In the second stage, the cost function (  2 ) will increase along with the time length   2 of the sleep cycle.During this period, the longer the time length of the sleep cycle is, the longer the MS will stay in the sleep period, the longer the average response time  of data packets is, and thus the greater the cost function will be.
Conclusively, there is a minimal cost (  2 ) when the time length   2 of the sleep cycle is set to an optimal value.The optimal time length of the sleep cycle and the minimal cost for the different correlation parameters are summarized in Table 2.

Conclusions
How to improve the energy saving mechanism in battery powered mobile stations is one of the most pressing questions facing researchers of communication networks.As an enhancement of WiMAX, WiMAX 2 can save more energy.In this paper, concerning the real-time traffic including

Figure 1 :
Figure 1: The state transition of the system model.

Figure 2 :
Figure 2: Various time epochs in this late arrival system with delayed access.
length T V 2 of the sleep cycle (ms) Average response time  of data packets (ms)

Figure 3 :
Figure 3: Average response time  of data packets versus sleep cycle   2 .

Figure 5 :
Figure 5: Average response time  of data packets versus correlation parameter .

Figure 7 :
Figure 7: Standard deviation  of the number of data packets versus arrival ratio  of data packets.

Table 2 :
Optimal time length   2 of the sleep cycle.