Performance of an M / M / 1 Retrial Queue with Working Vacation Interruption and Classical Retrial Policy

An M/M/1 retrial queue with working vacation interruption is considered. Upon the arrival of a customer, if the server is busy, it would join the orbit of infinite size. The customers in the orbit will try for service one by one when the server is idle under the classical retrial policy with retrial rate nα, where n is the size of the orbit. During a working vacation period, if there are customers in the system at a service completion instant, the vacation will be interrupted. Under the stable condition, the probability generating functions of the number of customers in the orbit are obtained. Various system performance measures are also developed. Finally, some numerical examples and cost optimization analysis are presented.


Introduction
Queueing models with server vacations have been well studied in the past three decades and successfully applied in manufacturing and production systems, service systems, and communication systems.Some vacation queues can be found in Tian and Zhang [1].On the basis of ordinary vacation, Servi and Finn [2] first introduced a class of semivacation policy known as working vacation (WV), where the server provides service at a lower speed during the vacation period rather than stopping service completely.The motivation of analyzing M/M/1/WV queue is to model approximately a multiqueue system where each queue can be served at one of two service rates.From then onwards, several works [3][4][5] have appeared that analyzed the single server queue with working vacations.In order to utilize the server effectively, Li and Tian [6] introduced vacation interruption policy.During the working vacation period, if at least one customer is present in the system at a service completion epoch, the server will interrupt the vacation and resume regular service.Working vacation interruption has become an important aspect.Using the matrix analytic method, Li and Tian [7] analyzed a GI/Geo/1 queue with working vacations and vacation interruption and Li et al. [8] discussed the GI/M/1 queue.Using the method of a supplementary variable, Zhang and Hou [9] investigated an M/G/1 queue with working vacations and vacation interruption.
Retrial queueing systems are described by the feature that the arriving customers who find the server busy join the retrial orbit to try again for their requests.During the last two decades, retrial queues have been studied extensively because of their applications in telephone switching systems, telecommunication networks, and computer systems.Readers can refer to Artalejo and Gómez-Corral [10] and Choudhury et al. [11].Recently, the retrial queueing systems with working vacations have been investigated extensively.Do [12] first studied an M/M/1 retrial queue with working vacations.Li et al. [13] considered a discrete time Geo/Geo/1 retrial queue with working vacations and vacation interruption and Liu and Song [14] introduced nonpersistent customers into the Geo/Geo/1 retrial queue with working vacations.Tao et al. [15] discussed an M/M/1 retrial queue with collisions and working vacation interruption under N-policy.Using the matrix analytic method, Do [12], Li et al. [13], Liu and Song [14], and Tao et al. [15] obtained the stationary probability distribution and showed the conditional stochastic decomposition for the queue length.Using the method of a supplementary variable, Aissani et al. [16] and Jailaxmi et al. [17] both generalized the model of [12] to an M/G/1 queue.Gao and Wang [18] analyzed a Geo  /G/1 retrial queue with general retrial times and working vacation interruption, and the continuous-time M/G/1 queue was investigated by Gao et al. [19].Note that the retrial policy of the above papers is either constant retrial policy or general retrial policy, where only the customer at the head of the orbit can request a service.
Many of the queueing systems with repeated attempts operate under the classical retrial policy, where each block of customers generate a stream of repeated attempts independently of the rest of the customers in the orbit.In telephone systems, the calls which will be reattempted can be modeled by an infinite server queue where the time until reattempt follows an exponential distribution.Thus, it needs to discuss a retrial queue with working vacation under the classical retrial policy.Ayyappan et al. [20] first studied such a queueing system.In order to obtain the steady state probability vector, the orbit size is restricted to  such that any arriving customer finding the orbit full is considered lost.In this paper, we also analyze an M/M/1 retrial queue with multiple working vacations under the classical retrial policy, and vacation interruption policy is considered.Ayyappan et al. [20] mainly focus on numerical computation and approximation, and they apply the direct truncation method to find the steady state probability vector.We place more emphasis on the analytic solutions, and the orbit size is always infinite in our model.The novelty of this investigation is that we use probability generating function method to deal with the steady state joint distribution of the server state and the number of customers in the orbit, and the ergodicity condition is obtained by Foster's criterion.Finally, a cost minimization problem is discussed.
This paper is organized as follows.In Section 2, we establish the model.Using Foster's criterion, the stable condition is obtained.In Section 3, we get the balance equations and derive the probability generating function of the number of customers in the orbit.Various performance measures of this model are discussed in Section 4. Section 5 presents some numerical examples and cost optimization analysis.Finally, Section 6 concludes the paper.

Model Formulation and Stable Condition
Let us consider an M/M/1 retrial queue with working vacations and vacation interruption.Customers arrive according to a Poisson process with rate .Upon the arrival of customers, if the server is free, arriving customers get service immediately.If the server is occupied, customers are forced to wait in the orbit of infinite size.In this paper, we adopt the classical retrial policy with retrial rate , where  is the orbit size.Upon the arrival of retrials, if the server is busy, retrials will go back to the orbit.If the server is not occupied, arriving retrials get services immediately.The service rate is  when the system is not on vacation.The single server takes a working vacation when the system becomes empty.Vacation time follows an exponential distribution with parameter .During the vacation period, customers can be served with rate .At a service completion instant in the vacation period, if there are customers in the system at that moment, the vacation is interrupted and the server comes back to the normal working level.The working vacation will continue if the system is empty.At the end of each vacation, the server only takes another new vacation if the system is empty.
Next, we give an application example of this model.In our daily life, some social organizations can provide telephone psychological counseling.Here, we consider a telephone counseling system staffed with a psychological consultant (called main server) and an assistant (called substitute server).The assistant only provides service to people (called customers) when the consultant does not work.Generally speaking, there is a telephone operator who is responsible for establishing communications between servers and customers; the operator also needs to note down the calls' information in a registration form (corresponding to the "orbit").When a person makes a call and if the server is free, the operator notes down the information and the customer is served immediately by the consultant or the assistant.If the server is busy, on the other hand, the operator tells the customer to call again some time later (called retrial), and the customer's information should also be noted down.When a retrial customer makes a call again, the operator does not need to note down the information again.Moreover, if the customer in the registration form completes his service, the telephone operator will make a note; then he can know whether the customer completed his service or not.When the operator finds that all customers in the registration form have completed their services (i.e., there is no customer in the orbit and the system is empty), the consultant will rest from his work (begin a vacation).During the consultant's vacation period, the assistant will provide service.If a service is completed by the assistant and if there are customers in the registration form who have not received services (i.e., there are customers in the orbit and the system is nonempty), the consultant will come back from his vacation no matter whether his vacation ends or not (vacation interruption happens).Furthermore, when a vacation ends and if all customers in the registration form have completed their services, the consultant begins another vacation.Otherwise, the consultant takes over the assistant.In order to understand the customer's need, the consultant will restart service no matter how long the assistant has served.
Let () be the number of customers in the orbit at time  and let () be the state of server at time .There are four possible states of the server as follows: 0, the server is in a working vacation period at time  and the server is free, 1, the server is in a working vacation period at time  and the server is busy, 2, the server is during a normal service period at time  and the server is free, 3, the server is during a normal service period at time  and the server is busy. (1) Note that when the orbit is empty, it is impossible that the server is free in a normal service period.Thus, {(), ()} is a Markov process with state space Remark 1.Note that, in this model, we consider working vacation interruption policy, so the state {(, 0),  ≥ 1} does not exist.
Define ) , where In order to discuss the ergodicity of the chain {  ;  = 1, 2, . ..}, we use Foster's criterion with the mean drift where the test function () = .From the transition matrix, it can be easily proved that As  → ∞, there exists  = lim →∞   =  − 1. Applying Foster's criterion, we can get that the embedded Markov chain is ergodic if  < 1. Besides, since  +1 −   ≥ −1, applying the criterion from Sennott et al. [21], we can obtain that  < 1 is necessary for ergodicity.
Since the arrival process is Poisson, using PASTA property, it can be shown from Burke's theorem (see [22, pp. 187-188]) that the steady state probabilities of the Markov process {(), ()} exist if and only if the stable condition  < 1 holds.Now we define the limiting probability

Performance Measures
From ( 19), we have From ( 22), we get where From (24), using L'Hospital rule, we obtain Since taking ( 25), ( 27), and (29) into the above equation, we can derive Thus,  00 is found out.From ( 19), (22), and (24),  1 (),  2 (), and  3 () are completely computed.And the pgf of the number of customers in the orbit is given by The pgf of the number of customers in the system is given by Let [  ] denote the average number of customers in the orbit when the server's state is ,  = 1, 2, 3.
Differentiating (19) with respect to , we have From ( 18), we can get Differentiating (24) with respect to  and using L'Hospital rule twice, the expression for [ 3 ] is obtained as and differentiating (19) twice with respect to  and taking  = 1 yield   1 (1) as Therefore, the average orbit length ([]) is given by And the average system length ([ L]) is derived as The probability that the server is busy is

Advances in Operations Research
The probability that the server is free is The probability that the server is in a working vacation period is given by The probability that the server is in a normal service period is given by Let [] ([ W]) be the expected waiting (sojourn) time of a customer in the orbit (system), using Little's formula: . (44 The system busy period  is defined as the period that starts at an epoch when an arriving customer finds an empty system and ends at the departure epoch at which the system is empty.
Using the theory of regenerative process, we get where  00 is the time length that the system is in the state (0, 0).Since the interarrival time between two customers follows exponential distribution with parameter , we have [ 00 ] = 1/.Thus, [] =  −1 ( −1 00 − 1).

Numerical Results
In this section, under the stable condition, we present some numerical examples to illustrate the effect of the varying parameters on the mean orbit length [], the probability that the system is busy 1− 0,0 , and the probabilities of server's state (  ,   ,   , and   ).Moreover, a cost minimization problem is also discussed.When  < , we find that [] and 1 −  0,0 both decrease as  increases, and  has a noticeable effect on the performance measures which cannot be ignored.An especial case is  = ; that is, the service rate in the vacation period is equal to the service rate in the normal period, and the model we consider reduces to an M/M/1 retrial queue with classical retrial policy.
Clearly,  has no effect on [] and 1− 0,0 , which agrees with the intuitive expectation.The effect of  on [] is presented in Figure 4; as expected, [] decreases with increasing values of .When  is large, it can be observed that the effect of  on [] is not obvious as  increases; this is because we consider working vacation interruption policy.If the system is not empty at a service completion instant in the vacation period, the server will come back to the normal working level.Moreover, since constant retrial policy means that only the customer at the head of the orbit can request a service, we can find that [] of M/M/1 queue with classical retrial policy is smaller than that of M/M/1 queue with constant retrial policy.As shown in Figure 5,   decreases with increasing values of , and   has the opposite tendency.The reason is that as  increases, the mean orbit length [] and the mean service time 1/ both decrease.
Figure 6 indicates that [] decreases with increasing values of ; this is because the mean retrial time decreases as  increases.From the instant when the server becomes idle, an external potential primary customer and retrial customers compete to access the server, and the smaller the mean retrial time is, the bigger the probability that the server is busy is, which leads to decreasing of [].Furthermore, under the same condition, decreasing of [] results in the increase in the probability that the system is empty, and Figure 7 reveals that   increases with increasing values of , while   decreases as  increases.In Figure 6, compared with constant retrial policy, we can also find that classical retrial policy decreases the waiting jobs effectively.

Cost Analysis.
In this subsection, we establish a cost function to search for the optimal service rate , so as to minimize the expected operating cost function per unit time.
Define the following cost elements: : cost per unit time for each customer present in the orbit.
: cost per unit time for service during a normal service period.
: cost per unit time for service in a working vacation period.
: fixed cost per unit time during a vacation period.
Because the expected operating cost function is highly nonlinear and complex, it is not easy to get the derivative of it.We assume   = 6,   = 10,   = 8, and   = 4 and develop approximations by MATLAB program to find the optimum value of , say  * .
From Figure 8, we can see that there is an optimal service rate  to minimize the cost.Implementing the computer software MATLAB by the parabolic method and controlling the error by 10 −4 , we find the solution  * = 0.8584 with ( * ) = 69.0971.

Conclusion
This paper analyzes a single server retrial queue with working vacation interruption under the classical retrial policy.Using embedded Markov chain and Foster's criterion, we get the condition of stability.The pgf of the number of customers in the orbit is obtained, and some important performance measures are also discussed.Moreover, the effects of various parameters on the system performance measures are examined numerically.Under the stable condition, a cost minimization problem is considered.For future research, using the same method, one can investigate a similar model with batch arrival and without vacation interruption.
The various parameters of this model are arbitrarily chosen as  = 2,  = 1,  = 5,  = 3, and  = 0.2, unless they are considered as variables or their values are mentioned in Figures 2, 3, 4, and 5. 5.1.Sensitivity Analysis.Figures 2 and 3 illustrate the effect of  on [] and 1 −  0,0 for different values of , respectively.

Figure 2 :Figure 3 :
Figure 2: The effect of  on [] for different values of .

Figure 4 :Figure 5 :
Figure 4: The effect of  on [] for two different retrial policies.

Figure 6 :
Figure 6: The effect of  on [] for two different retrial policies.

Figure 7 :
Figure 7: The probability of server's state with the change of .

Figure 8 :
Figure 8: Effect of  on the expected operating cost per unit time.