Analysis of General Input State Dependent Working Vacation Queue with Changeover Time

We consider a finite buffer GI/M(n)/1 queue with multiple working vacations and changeover time, where the server can keep on working but at a slower speed during the vacation period. Moreover, the amount of service demanded by a customer is conditioned by the queue length at the moment service is begun for that customer. We provide a recursive algorithm using the supplementary variable technique to numerically compute the stationary queue length distribution of the system. Finally, some numerical results of the model are presented to show the parameter effect on various performance measures.


Introduction
Queueing systems with vacations are considered to be effective tools in modeling and analyzing complex computer and communication networks and several other engineering systems where the server can utilize the idle time for different purposes; see Doshi [1] and Tian and Zhang [2].On the other hand, numerous situations exist where the server remains active during the vacation period.Servi and Finn [3] introduced a class of semivacation policies for the analysis of an //1 queue.In this the server renders service to the queue with a lower service rate, known as working vacation (WV).In multiple working vacations (MWV) policy when a vacation ends and the system is not empty, a service period begins with the normal service rate; otherwise, the server takes another vacation.
Baba [4] extended Servi and Finn's [3] work to GI//1/ MWV queue using the matrix-analytic approach and derived the expressions of distributions for queue length at an arrival epoch and the steady state distribution for the waiting time.Banik et al. [5] discussed the GI//1/ queue with MWV and obtained some important performance measures.Zhang and Hou [6] generalized GI//1/ queue with a variant of MWV.They employed the supplementary variable and embedded Markov chain methods to obtain the queue length distribution at different epochs.At a service completion epoch during regular service, the server remains in the system for some time instead of leaving for a WV immediately.This time is known as changeover time during which if a customer arrives, the server starts service with regular service rate; otherwise, it leaves for a WV.Li and Zhu [7] provided an explicit formula for the Laplace transform of the additional delay for //1 queues with delayed vacations and exhaustive service discipline.Dong and Doo [8] investigated the //1 queue with changeover time and multiple vacations.They derived the joint distribution of the queue length and the remaining service (or vacation) time at arbitrary time by using the supplementary variable method and calculated the virtual waiting time distribution.The steady state behavior of an / (,,) /1 queueing system with multiple vacations and delayed vacations has been investigated by Jain and Singh [9].Recently, Vijaya Laxmi and Seleshi [10] extended this study to a renewal input GI/ (,,) /1 queue with WV and changeover times.
Many standard queueing systems operate on the assumption that input and service parameters are dependent on the state of the system, which reduces the waiting time of customers.For example, in telecommunication systems at the packet switch (router) when its buffer size increases, a controller drops the arriving packets with an increasing probability.In human based service systems, it is known that there is a strong correlation between the volume of work demanded from a human and his/her productivity.A computational algorithm of GI/()/1/ queue with state dependent vacations and -policy with multiple vacations are analyzed by Chao and Rahman [11].Using the supplementary variable technique, a recursive algorithm has been explained by Goswami et al. [12] for obtaining the system length distributions at prearrival and arbitrary epochs for a finite buffer state dependent GI//1 queue with MWV.
The above literature survey indicates that so far no work has been focused on state dependent MWV queue with changeover time, which is applicable in production, manufacturing, traffic signal and telecommunication systems, and so forth.For example, in a packet switched network, the router is an interconnection device that attaches two or more networks.It takes charge of receiving packets and forwarding them to the next hop, according to some routing information in its routing table.If the routes are available in the routing table, router will serve the packets depending on the number of packets which is termed as state dependent routing.This procedure routes the packets to the least loaded disks.This offers significant throughput benefits over state independent routing.To collect the complete routing information, the router may exchange its routing information with the other routers.When there is no packet to be sent, a "sleep mode" is operated for the purpose of power saving.If a packet arrives during this sleep mode period, the router returns to the awakening state immediately and it transmits packets.Meanwhile, under the maintenance, the router can serve the packet at slower speed which can reduce waiting times and economize the cost.In queueing terminology, the router, routing at slower speed, state dependent routing, and sleep mode correspond to the server, working vacation, state dependent services, and changeover times, respectively.
Motivated by such situations of sleep mode operations, which conserve energy and further reduce the waiting times, this paper aims to focus on finite buffer state dependent queue with MWV and changeover time.We assumed the service times during service period, vacation period, and vacation times are exponentially distributed.We provide a recursive method using the supplementary variable technique and treating the remaining interarrival time as the supplementary variable, to develop the steady state system length distributions at prearrival and arbitrary epochs.A computational algorithm is presented to compute the stationary system length distribution.Some performance measures and numerical results have been illustrated in the form of tables and graphs.
The paper is structured as follows.Next section presents the description of the model.Steady state equations and relation between prearrival and arbitrary epochs are derived in Section 3. A computational algorithm is presented in Section 4 followed by various performance measures in Section 5. Section 5 also contains some numerical results to show the effectiveness of the model parameters and conclusions are given in Section 6.

Description of the Model
Let us consider a GI/()/1//MWV queue with changeover time where  is the capacity of the system.We assume that the interarrival times of successive arrivals are independent and identically distributed (i.i.d.) random variables with cumulative distribution function (), probability density function (),  ≥ 0, Laplace-Stiltjes (LS) transform  * (), and mean interarrival time 1/ = − * (1) (0), where ℎ (1) (0) denotes the first derivative of ℎ() evaluated at  = 0.The service rate during a regular busy period is assumed to be an exponentially distributed random variable.When there are  customers present in the system before beginning a service, the server serves them with a rate   , 1 ≤  ≤ .The server takes WV whenever the system becomes empty.When a vacation ends and if there are customers in the queue, a regular busy period begins and server serves the queue with its usual service rate; otherwise, the server takes another vacation.The vacation rate and service rate during any WV period are also assumed to be exponentially distributed random variables with   ,   (1 ≤  ≤ ), respectively, where there are  customers present in the system during a vacation.Whenever the queue becomes empty, instead of leaving for a WV immediately the server will remain in the system for a certain period of time, called changeover time and which is assumed to be exponentially distributed with parameter .Let () and  be the mean service rate during regular busy period (WV) and mean vacation rate, respectively, and they are given by  = ∑  =1   /,  = ∑  =1   /,  = ∑  =1   /.The customers are served by a single server on first-come first-served (FCFS) discipline.The arrival times, service time, and changeover times are mutually independent.
Let us define the state of the server as if the server is in regular busy period or in changeover time. ( We define the joint probability densities of system length   (), state of the server (), and the remaining interarrival time for the customer (), respectively, as The above probabilities at steady state are denoted by  , ().

Relation between Steady State Distributions at Arbitrary and Prearrival Epochs.
To obtain the arbitrary epoch probabilities, we develop the relations between number of customers in the system at prearrival and arbitrary epochs as follows.Let  − , ,  = 0, 1, (0 ≤  ≤ ) represent the probability of  customers in the system at prearrival epoch when the server is in state .Applying Bayes theorem, for  = 0, 1, we have ,  = 0, 1, 0 ≤  ≤ .
Now, to obtain the steady state probabilities at arbitrary epochs, we develop a relation between prearrival and arbitrary epoch probabilities.Setting  = 0 in (11), (10), and ( 14)- (12) and using (28), we obtain Remark 1.When  → ∞, that is, there is no changeover time, the model reduces to GI/()/1/ queue with MWV and our results match numerically with Goswami et al. [12].
Remark 2. When  → ∞,   = 0, for all , that is, there is no changeover time or service during WV, the model reduces to GI/()/1/ queue with multiple vacations.In absence of -policy, the above results match analytically with the results given by Chao and Rahman [11].
Remark 3. When  → ∞,   = ,   = ,   = , for all , that is, there are state independent services and vacations, our model reduces to GI//1/ queue with MWV and results match with Banik et al. [5].
Remark 4. When  → ∞,  → 0, for all , there is no WV and the model becomes state independent services with vacations.The matching has been done numerically with those of Tian et al. [13].

Computational Algorithm
To demonstrate the computational aspect for obtaining the prearrival epoch probabilities  − , at steady state, the following algorithm is developed based on the analysis of Section 3. The computational complexity of the algorithm is ( 3 ), where  is the maximum capacity of the system.

Performance Measures and Numerical Results
In this section, we discuss some operating characteristics of the queueing system such as average number of customers in the queue (  ), average number of customers in the system (  ), and the blocking probability of the server ( loss ).They are given by The average waiting time of a customer in the queue (system)   (  ) using Little's rule is given by   =   / λ, (  =   / λ), where λ = (1 −  loss ) is the effective arrival rate.
Table 1 presents the sensitivity analysis for various models like state dependent and MWV with changeover time (state + MWV + chover), state dependent and MWV (state + MWV), and constant with MWV (constant + MWV) models.The results are obtained by considering exponential interarrival time distribution.Note that for exponential distribution, the prearrival and arbitrary epoch probabilities are the same due to the memoryless property.So, we have presented only arbitrary epoch probabilities in the table.Moreover, as  → ∞, we obtain models with state dependent and constant service rates.Further, the performance measures  loss ,   ,   are lower in state + MWV + chover queue.Finally, the model with changeover time gives better results and justifies our analysis.
Figure 1 compares the effect of arrival rate  on the expected queue length   among the models: (i) state + MWV + chover, (ii) constant + MWV + chover, (iii) state + MWV, and (iv) constant + MWV for the deterministic interarrival time distribution.It is clear that   increases as arrival rate increases in all the above models.Further, among the four models considered, the state dependent model with changeover time performs best.
Figure 2 illustrates the effect of mean service rate during WV () on the average waiting time in the system (  ) for different mean vacation rates  = 0. considering   = 0.1,   = 0.15,   = 0.2) for exponential interarrival time.Considering the parameters by  = 5,   = 0.5 (1.5) ∀, we observe that   decreases as  increases.Further, as the mean vacation rate  increases, the   decreases and meets at a point  =  = 1.5.Therefore, the WV queue utilizes the idle time effectively when  ≤ .
Figure 3 depicts the effect of arrival rate  on average system length   with different values of  for  3 distribution.It may be observed that   increases with arrival rate.Further, increase of  also results in the increase of queue lengths.The effect of buffer size  on blocking probability  loss for various interarrival time distributions is shown in Figure 4. We observe that as the buffer size increases the  loss decreases.Furthermore,  loss is the highest in case of HE 2 distribution, whereas deterministic distribution yields the lowest.

Conclusions
This paper presents a single server state dependent queue with MWV and changeover time.The interarrival time of customer arrival is arbitrarily distributed, while the service rates during regular busy period and during WV and vacation rate are exponentially distributed.A recursive method has been developed to obtain the steady state queue length distributions at prearrival and arbitrary epochs.Numerical results have also been discussed.The changeover period concepts can be added to analyze more complex models like GI//1/ and GI//1/ queues with various services and vacation policies.