Scholarly article on topic 'An M/M/2 Queueing System with Heterogeneous Servers Including One with Working Vacation'

An M/M/2 Queueing System with Heterogeneous Servers Including One with Working Vacation Academic research paper on "Mathematics"

Share paper
OECD Field of science

Academic research paper on topic "An M/M/2 Queueing System with Heterogeneous Servers Including One with Working Vacation"

Hindawi Publishing Corporation International Journal of Stochastic Analysis Volume 2012, Article ID 145867,16 pages doi:10.1155/2012/145867

Research Article

An M/M/2 Queueing System with Heterogeneous Servers Including One with Working Vacation

A. Krishnamoorthy1 and C. Sreenivasan2

1 Department of Mathematics, Cochin University of Science and Technology, Cochin 682022, India

2 Department of Mathematics, Government College, Chittur, Palakkad 678104, India

Correspondence should be addressed to A. Krishnamoorthy, Received 24 February 2012; Accepted 13 May 2012 Academic Editor: Ho Lee

Copyright © 2012 A. Krishnamoorthy and C. Sreenivasan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper analyzes an M/M/2 queueing system with two heterogeneous servers, one of which is always available but the other goes on vacation in the absence of customers waiting for service. The vacationing server, however, returns to serve at a low rate as an arrival finds the other server busy. The system is analyzed in the steady state using matrix geometric method. Busy period of the system is analyzed and mean waiting time in the stationary regime computed. Conditional stochastic decomposition of stationary queue length is obtained. An illustrative example is also provided.

1. Introduction

Queueing models with vacation have gained significance in the last three decades due to their wide range of applications, especially in the communication and the manufacturing systems. Doshi [1] provides an excellent survey of related works prior to 1986. Takagi [2] and Tian and Zhang [3] provide a good account of developments in this field since then. The literature on the vacation queueing models is growing rapidly.

In multiserver queueing models, we come across two classes of vacation mechanisms: station vacation and server vacation. In the first case, all servers take vacation simultaneously whenever the system becomes empty and they return to the system all together. Thus, station vacation is group vacation for all servers. For example, when a system consists of a number of machines operated by a single individual this scenario occurs. In such a situation, the whole station needs to be treated as a single unit for vacation, when the system is utilized for a secondary task. In the second case, each server takes its own vacation whenever it completes a service and finds no customers waiting for service. This phenomenon also occurs in practice. For example, in a post office or bank, when a clerk completes a service and finds no customer

waiting, he or she might go to attend another task. This is what we refer to as the server vacation model. Analysis of a server vacation model is more complicated than that of a station vacation model. This is because at any time point, in the latter, we may have any number of servers between 0 and c on vacation. We need to track individual servers going on vacation and completing their vacation. Upon returning from a vacation some servers may find no customers waiting for service. These servers take another vacation. But if any server finds a waiting customer on returning from a vacation, it immediately starts service. For further details on queues with station and server vacations, we refer the reader to Chao and Zhao [4].

Most of the earlier work on multiserver queueing models deal with homogeneous servers; that is, the individual service rates are same for all servers in the system. But this assumption may be valid only when the service process is mechanically or electronically controlled. In a queueing system with human servers the above assumption is highly unrealistic. Often servers providing identical service, serves at different rates. This motivated the researchers to study multiserver queueing system with heterogeneous servers. Several authors have analyzed multiserver queues with vacation. Levy and Yechiali [5], Vinod [6], and Kao and Narayanan [7] discuss asynchronous multiple vacation models with exponentially distributed vacation times. They study the scenario, where any server goes on vacation whenever there are no customers waiting in the system at a service completion epoch. At a vacation termination instant, if there are no waiting customers, the server takes another vacation; if there is a customer waiting for service, the server resumes service. It should be remarked that Vinod [6] was the first to use matrix geometric solutions method to analyze multiserver vacation models. But none of the above models deal with heterogeneous servers.

Servi and Finn [8] introduced a working vacation model with the idea of offering services but at a lower rate, whenever the server is on vacation. Their model was generalized to the case of M/G/1 in [9,10] and to GI/M/1 model in [11]. A survey of working vacation models with emphasis on the use of matrix analytic methods is given in Tian et al. [12]. Working vacation models have a number of applications in practice. Two such examples are given in [12]. Recently, Li and Tian [13] studied an M/M/1 queue with working vacations in which vacationing server offers services at a lower rate, for the first customer arriving during a vacation. Very recently, Zhang and Hou [14] studied a MAP/G/1 queue with working vacations and vacation interruptions using supplementary variable method. In this model, the authors assume that the vacation times are exponentially distributed and that the server gets back to normal service mode, when at a service (offered during a vacation) completion, the system has at least one customer waiting in the queue. The server is allowed to take multiple vacations. However, no work on multiserver working vacation model has come to our notice.

Neuts and Takahashi [15] observed that for queueing systems with more than two heterogeneous servers analytical results are intractable and only algorithmic approach could be used to study the steady state behavior of the system. Based on this observation, Krishna Kumar and Pavai Madheswari [16] analyzed M/M/2 queueing system with heterogeneous servers, where the servers go on vacation in the absence of customers waiting for service. In this paper, we discuss an M/M/2 queueing model with heterogeneous servers where one server remains idle but the other goes on working vacation in the absence of waiting customers.

This paper is organized as follows. In Section 2, model description is provided. In Section 3, the steady state analysis of the model is presented. In Section 4, we discuss an illustrative example.

International Journal of Stochastic Analysis 2. Mathematical Model

We consider an M/M/2 queueing model with heterogeneous servers, server 1 and server 2. Server 1 is always available, whereas server 2 goes on vacation whenever there are no customers waiting for service. Let the service rates of servers 1 and 2 be ¡i and ¡2, respectively, where ¡1 / ¡2. Customers arrive to the system according to a Poisson process with parameter X. The duration of vacation is exponentially distributed with parameter n At the end of a vacation, service commences if there is a customer waiting for service. Otherwise the server goes on another vacation. During vacation, if an arrival finds server 1 busy, server 2 returns to serve the customer but at a lower rate. To be precise, server 2 serves this customer at the rate d^2, 0 < d < 1. As this vacation gets over, server 2 instantaneously switches over to the normal service rate ¡2, if there is at least one customer waiting for service. Upon completion of a service at low rate, the server will (a) continue the current vacation if it is not finished and no customer is waiting for service; (b) continue the slow service if the vacation has not expired and if there is at least one customer waiting for service. For clarity, we make it clear that if an arriving customer finds a free server, he enters service immediately. Else he joins the queue.

2.1. The QBD Process

The model discussed in Section 2 can be studied as a quasi-birth-and-death (QBD) process. First, we set up the necessary notations.

At time t, let N(t) be the number of customers in the system and

J (t) =

0, if the server 2 is on vacation,

1, if the server 2 is working in vacation mode,

2, if the server 2 is working in normal mode,

Let X(t) = (N(t),J(t)). Then (X(t) : t > 0) is a continuous time Markov chain (CTMC) with states space

Q= {(0,0), (1,0), (1,1), (1,2)}U IM

l(i) = {(i,j) : i > 2, j = 1 or 2}.

Note that when N(t) = 0, the only possible value of J(t) is 0 and when N(t) = 1, J(t) has three possible values 0, 1, and 2.

The infinitesimal generator matrix Q of this Markov chain is given by

B00 B01 B10 B11 B12

B21 A1 A0

A2 A1 A0

where the block matrices appearing in Q are as follows:

B00 = -1, B01 = [1 0 0],

¡1 '-1 - ¡1 0 0 1 0

B10 = 0^2 , B11 = 0 -1 - 0^2 - n n , B12 = 1 0

.¡2 _ 0 0 -1 - ¡¡2 0 1

d^2 ¡1 0 ¡2 0 ¡1

-1 - ¡1 - 0^2 - n n

0 -1 - ¡1 - ¡2

1 0 0 1

¡1 + d^2 0 0 ¡1 + ¡2

3. Steady State Analysis

In this section, we will discuss the steady state analysis of the model under study.

3.1. Stability Condition

Theorem 3.1. The queueing system described above is stable if and only if p < 1 where p = 1/(^1 + ¡2)-

Proof. To establish the stability condition we use Pakes' lemma (see [17]). Let Ni be the number of customers in the system immediately after the departure of the ith customer. Then {Ni : i e N} satisfies the equation

Ni-1 - 1 + Vi if Ni-1 > 1

if Ni-1 = 0,

where Vi is the number of arrivals during the service of ith customer. Clearly [Ni : i e N} is an irreducible aperiodic Markov chain. Pakes' lemma asserts that an aperiodic irreducible Markov chain is ergodic, if there exists an e > 0 such that the mean drift

fy = E [(Ni+1 - Ni) | Ni = j\

is finite for all j e N and tj < -e for all j e N except perhaps for a finite number. In the present model, value of the mean drift is

t, = f-1 + P if j 21 (3.3)

VP if, = 0

Thus, if p < 1, the Markov chain {Ni : i e N} is ergodic and hence the condition is sufficient.

To prove the necessity of the condition, assume that p > 1. We use Theorem 3.1 in Sennot et al. [18], which states that {Ni : i e N} is nonergodic if it satisfies Kaplan's condition, tj < cc, for j > 0 and there is a j0 such that tj > 0, for j > j0. When p > 1 Kaplan's condition is readily satisfied. Hence, the Markov chain is not ergodic. □

3.2. Steady State Probability Vector

Let x, partitioned as x = (x0, xi, x2,...), be the steady state probability vector of Q. Note that x0 is a scalar, x1 = (x10,x11,x12) and x* = (xi1 ,xi2) for i > 2. The vector x satisfies the condition xQ = 0 and xe = 1, where e is a column vector of 1's of appropriate dimension. Apparently when the stability condition is satisfied, the subvectors of x, corresponding to the different levels, are given by the equation xj = x2Rj-2, j > 3, where R is the minimal nonnegative solution of the matrix quadratic equation (see [20])

R2A2 + RA1 + A0 = 0. (3.4)

Knowing the matrix R, x0, x1 and x2 are obtained by solving the equations

xoBoo + X1B10 = 0, (3.5)

x0B0i + xiBn + x2B21 = 0, (3.6)

X1B12 + X2(A1 + RA2) = 0, (3.7)

subject to the normalizing condition

X0 + X1e + X2(I - R)-1e = 1, (3.8)

Theorem 3.2. The matrix R of (3.4) is given by R = [R0U R^ j, where R11 = (X + ¡1 + + n -\J(X + ¡1 + + n)2 - 4X(^1 + ))/2(^1 + d^2), R12 = p- (¡1 + )R11/(¡1 + ¡2) and R22 = p.

Proof. Since A0, A1, and A2 are upper triangular, R is essentially an upper triangular matrix. The value of R11, follows from the assertion that R is the minimal nonnegative solution of (3.4). The rest of the proof is an easy consequence of the condition RA2e = A0e. □

Remark 3.3. Though R has a nice structure which enables us to make use of the properties

like Rk =

R11 R12 2 k=0R11Rk2j 1

Rk r22

, for k > 1, due to the form of the expression for Rn it may not

be easy to carry out the computations required in the forthcoming discussions. Hence, we explore the possibility of algorithmic computation of R. The computation of R matrix can be carried out using a number of well-known methods such as logarithmic reduction. We list here only the main steps involved in logarithmic reduction algorithm for computation of R. For full details on the logarithmic reduction algorithm we refer the reader to [19].

Logarithmic Reduction Algorithm for R:

Step 0. H ^ (-Ai)-1Ao, L ^ (-Ai)-1A2, G = L, and T = H. Step 1.

U = HL + LH M = H2 H (I - U)-1M

M L2 (3.9)

L (I - U)-1M G G + TL T TH,

continue Step 1 until ye - Ge||«» < e. Step 2. R = -A0(A1 + A0G)-1.

3.3. Busy Period Analysis

For the system under study, busy period is the interval between arrival of a customer to the empty system and the first epoch thereafter when the system becomes empty again. Thus, it is precisely the first passage time from the state (1,0) to the state (0,0). For the working vacation model, busy cycle for the system is the time interval between two successive departures, which leave the system empty. Thus, the busy cycle is the first return time to state (0,0) with at least one visit to any other state. But before analyzing the busy period structure we need to introduce the notion of fundamental period. For the QBD process under consideration, it is the first passage time from level i, where i > 3, to the level i - 1. The cases i = 2, i = 1, and i = 0, corresponding to the boundary states, need to be discussed separately. It should be noted that due to the structure of the QBD process, the distribution of the first passage time is invariant in i.

Let Gjj' (k, x) denote the conditional probability that a QBD process starting in the state (i,j) at time t = 0 reaches the level i - 1 for the first time no later than time x, after exactly k transitions to the left and does so by entering the state (i-1,;'). For convenience we introduce the joint transform

Gjf (z,s) = VzM e~sxdGjf (k,x); |z| < 1, Re(s) > 0, (3.10)

k=1 ->0

International Journal of Stochastic Analysis and the matrix

(5 (z,s) = (Gjr (z,s)). (3.11)

The matrix (5(z,s) is the unique solution to the equation (see [20])

G(z, s) = z(sI - A1)-1 A2 + (sI - A1)-1Aq(G2(z/s). (3.12)

The matrix G = G(1,0) takes care of the first passage times, except for the boundary states. If we know the R matrix then G matrix can be computed using the result (see [19])

G = -(A1 + RA2 )-1A2. (3.13)

Otherwise, we may use logarithmic reduction method to compute G. For the boundary level states 2, 1, and 0, let Gj'^^x), Gjj,/0)(k, x), and Gj./0) (k, x) be the conditional probability discussed above for the first passage times from level 2 to level 1, level 1 to level 0, and the first return time to the level 0, respectively. Then as in (3.12) we get

G(2,1) (z, s) = z(sI - A1 )-1B21 + (si - A1 )-1A(]G(z, s)G(2,1) (z, s), (3.14)

G(1'0)(z,s) = z(sI - Bn)-1B10 + (si - Bn)-1B12G(24)(z,s)G(1'0)(z,s), (3.15)

GG(0'0)(z,s)= [-jL 0 0 GG(1-0)(z,s). (3.16)

.s + 1

Note that G(1,0)(z, s) is a 3 x 1 matrix. Thus, the Laplace Stieltjes transform (LST) of the busy period is the first element of G(1,0) (1, s). For convenience, we use the notations

G21 = G (24)(1,0), G10 = G (1-0)(1,0), G00 = Gm(1,0). (3.17)

Due to the positive recurrence of the QBD process, matrices G, G21, G10, and Goo are all stochastic. If we let

C0 = (-A1)-1A2, C2 = (-A1)-1A0, (3.18)

then G is the minimal nonnegative solution (see [20]) to the matrix equation

G = C0 + C2G2. (3.19)

From (3.14), (3.15), and (3.16), we get

G21 = -(Ai + AqG)-1B2i, G10 = -(B11 + B12G21) 1Bio, Goo = [10 0] G10


respectively. Equation (3.12) is equivalent to

ZA2 - (sI - A1 )G(z, s) + A0G2(z, s) = 0.


dG (z,s)

dG(z, s)




Differentiation of (3.21) with respect to s and z followed by setting z = 1 and s = 0 leads to (see [20])

M = -A-1G + C2 (GM + MG),

M = C0 + C2( GM + ÎMG),


with 0 as starting value for M and M, successive substitutions in the above equations yield the values of M and A4. Applying an exactly similar reasoning to (3.14), (3.15), and (3.16), we get

M21 = -(A1 + A0G)-1(I + A0MG21, M10 = (B11 + B12G21)-1 (I + B12M21)G10,


G10 + [1 0 0 M10,

M21 = -

Mio = -

Moo = -

dG (21)(zrs)

dG(1/0) (z, s)

dG (0'0)(z,s)





Note that M10 is a 3 x 1 matrix and M00 is a scalar. The first element of the matrix M10 and M00 are mean lengths of a busy period and a busy cycle, respectively. The second and third elements of the matrix Mi0 are the first passage times to the state (0,0) from (1,1) and (1,2), respectively. With the notations

d(5(2,1) (z, s)

dG(1,0)(z, s)




it follows from (3.14) and (3.15) that

M21 = -(A1 + A0G)-1(B21 + A0MG21), M10 = -(B11 + B12G21) 1(B10 + B12 M21G10).


The first component of the vector M10 is the mean number of service completions in a busy period.

3.4. Stationary Waiting Time Distribution in the Queue

Let W(t) be the distribution function for the waiting time in the queue of an arriving (tagged) customer. Note that if there is no customer in the system, the arrival receives service immediately. If either of the two servers is not busy then also there would be no delay in getting service. Thus, the probability that the customer gets his service without waiting is x0 + x10 + x11 + x12. Hence, with probability 1 - x0 -x10 - x11 - x12, the customer has to wait before getting the service. The waiting time may be viewed as the time until absorption in a Markov chain with state space

Q1 = {*} U(2,3,...}. (3.28)

Here * is the absorbing state, which corresponds to taking the tagged customer into service and is obtained by lumping together the level states 0 = {(0,0)} and 1 = {(1,0), (1,1), (1,2)}.

For i > 2, the level i is given by i = {(i,j), j = 1 or 2}. The states other than the absorbing state correspond to the number of customers present in the system as the tagged customer arrives. Once the tagged customer joins the queue, the subsequent arrivals will not affect his waiting time in the queue. Hence the parameter 1 does not show up in the generator matrix Q of this Markov process, given by

Q = 2 A2e D

where D =

- d^2 - n n 0 -^1 - ^2.


Define vector

y(t) = (y. (t), yi(t), ys(i),...), (3.30)

yi(t) = (yi1(t),yi2(t)), for i > 2. (3.31)

The components of the yi(t) are the probabilities that at time t, the CTMC with generator Q is in the respective states of level i. Note that the scalar y.(t) is the probability that the process is in the absorbing state at time t. By the PASTA property, we get

y(0) = (X0 + X11 + X10 + X12, X2, X3,...). (3.32)


W(t) = y. (t), for t > 0. (3.33)

The LST of W(t) is given by (see [20])

W(s) = ¿yi(0)[(sI - D)-1AJi-2(sI - D)-1A2e. (3.34)

The mean waiting time can be obtained from W(s) as

c i-1 c

E(W) = -W'(0) ^X2+i£tW(-D)-1Ui-jUe + £x2+iUi(-D)-2A2e, (3.35)

i=1 j=0 i=0

where U = (-D)-1A2 is a stochastic matrix. Hence, (3.34) can be simplified as

E(W) = -W'(0) = £X2+iXUj (-D)-1e + £X2+iUi(-D)-1e.

i=1 ;=0


H = £ X2+iUi, i=0


Since U is stochastic, we get

He = x2(I - R) 1e = 1 - x0 - x10 - x11 - x12


This result can be used to find an approximate value of H and hence that of the second term in (3.36) to any desired degree of accuracy. Thus, only the first term in (3.36) demands serious computation. For this we make use of the ideas in [15,16, 21]. Now consider the matrix

0 1 0 1


which has the property that

UU2 = U2U = U2.


Then we get

(I - U + U2) = I - Ui + iU2, for i > 1.


By the classical theorem on finite Markov chains, the matrix (I - U + U2) is nonsingular (see [22]). In view of the last equation, the first term in (3.36) becomes

£x2+i(I - Ui + U2) i=1

(I - U + U2)-1(-D)-1e


With this simplification, we get

E(W) = [x2 (R(I - R)-1 +1 + R(I - R)-2U^ - h] (I - U + U2)-1 (-D)



3.5. Conditional Stochastic Decomposition of Queue Length

In this section, we provide a stochastic decomposition of queue length in the stationary regime, subject to the condition that both servers are busy. Note that from (3.5)-(3.8) we get x0, xio, xn, X12, x2i, and x22. Let Qv be the queue length of the vacation model under study, subject to the condition that both servers are busy. Then we have the following.

Theorem 3.4. If p < 1, then Qv = Q0 + Qd, where Q0 and Qd are two independent random variables. Q0 is the queue length of the M/M/2 queueing model with heterogeneous servers without vacation and Qd can be interpreted as the additional queue length due to vacation and slow service, subject to the condition that both servers are busy.

Proof. Let Pb denote the Probability that both servers are busy. Then

to to to n-3

Pb = 2 xn2 = 2 x22pn-2 + Z x2lRl2% Rl1j p^-


so that

n=2 n=2 n=3 j=0

= X22^pk + X21R12^Rnk%pk; k = n - 3

k=0 k=0 k=0

= (1 - P) ~\X22P + X2lRl2(1 - Rll)-^,

P- = (1 - p) (x22P + X2iRi2(1 - Rii)-1)-1 = (1 - P)6, (3.45)

6 = (X22P + X21^12(1 - R11)-1)(3-46)

Qv(z), the generating function of the queue length subject to the condition that both servers are busy, is given by

-1 TO -1 TO

Qv (z) = Xn1 zn-2 = -Ë X22pn-2zn-2

Pbn=2 P n=2

Pr t (

J Pn-j-3 \zn-3



X21R12XR11jpn J " lz"

n=3 \ j=0

By following a computational procedure similar to that of Pb, we arrive at

r. ( ) 1 - P \c X21R12 M

Qv (z) = --\ S( X22 pz + -—-— ) \

1 - pz[ \ r 1 - Rnz/J (3.48)

= Q0(z)Qd(z),



From (3.49) it follows that Q0(z) is the generating function of an M/M/2 heterogeneous queuing model without vacations, which is precisely the case ¡5 = 1 in [23]. Equation (3.50) suggests that Qd has a geometric distribution with parameter 1 - R11. □

Remark 3.5. Due to the algorithmic approach used in the derivation stationary waiting time distribution, a similar decomposition result for the waiting time distribution is far from reality.

3.6. Key System Performance Measures

In this section, we list a number of key system performance measures along with their formulae in addition to the busy period structure and the mean waiting time discussed above.

(1) The probability that the system is empty: PEMP = x0.

(2) The probability that the server 1 is idle: PIDL = x0 + x11 + x12.

(3) The probability that the server 2 is on vacation: PVAC = x0 + x10.

(4) The probability that the server 2 is working in vacation mode: PSLOW =

x;1 = x11 + x21 /(1 - R11).

(5) The probability that the server 2 is working in normal mode: PNORM = 1 - x0 -PSLOW.

(6) The mean number of customers in the system: ¡NS = X;c=1 jx;e = x10 + x11 + x12 + x2(I - R)-2R-1e - x2R-1e.

4. Numerical Results

4.1. Illustrative Example

We analyze the effect of the parameters X, n and d on the key performance measures. Table 1 analyzes the effect of X, Table 2 explains the effect of n, and Table 3 examines the effect of d on the performance measures. To this end, we use the following abbreviations in addition to the notations used in Section 3.6.

¡wtq: Mean waiting time in the queue. ¡LBP: Mean length of a busy period. ¡LBC: Mean length of a busy cycle.

¡NSBP: Mean number of service completions in a busy period.

(i) Since ¡1 and ¡2 are fixed, the traffic intensity p increases with X. Due to this Pn0rm, ¡NS and ¡wtq increase and PVAC and PIDL decrease as X increases. Note that the busy period starts with the Markov chain in the state (1,0); that is, with server 2

Q0(z) = P

Qd(z) = öix22pz +

1 - pz' X21R12

1 - Ruz

Table 1: Case A: = 10, fi2 = 5, n = 1 and 0 = 0.6. Case B: ^ = 5, fi2 = 10, n = 1 and 0 = 0.6.

X A/B Pldl Pvac pslow pnorm ^ns ^wtq ^lbp ^lbc ^nsbp

2 A 0.829 0.913 0.065 0.022 1.676 0.002 0.080 0.496 1.167

B 0.707 0.916 0.070 0.014 2.176 0.005 0.125 0.482 1.17

4 A 0.687 0.745 0.168 0.087 1.934 0.009 0.067 0.246 1.128

B 0.520 0.755 0.188 0.057 2.301 0.016 0.097 0.236 1.156

6 A 0.550 0.568 0.246 0.186 2.189 0.017 0.062 0.166 1.121

B 0.374 0.574 0.296 0.130 2.60 0.031 0.088 0.164 1.203

8 A 0.414 0.405 0.285 0.311 2.547 0.027 0.062 0.132 1.172

B 0.251 0.399 0.363 0.238 3.178 0.047 0.094 0.143 1.362

10 A 0.282 0.262 0.278 0.460 3.193 0.036 0.072 0.122 1.340

B 0.150 0.245 0.365 0.390 4.282 0.062 0.120 0.153 1.755

12 A 0.158 0.141 0.218 0.641 4.702 0.043 0.106 0.144 1.849

B 0.074 0.123 0.282 0.595 6.626 0.069 0.191 0.216 2.818

14 A 0.049 0.042 0.092 0.866 11.885 0.047 0.294 0.324 4.659

B 0.020 0.034 0.113 0.853 16.043 0.064 0.566 0.585 8.418

Table 2: Case A: X = 12, = 10, fi2 = 5 and 0 = 0.6. Case B: X = 12, = 5, fi2 = 10 and 0 = 0.6.

n A/B Pldl Pvac pslow pnorm ^ns ^wtq ^lbp ^lbc ^nsbp

0.1 A 0.112 0.092 0.333 0.575 5.575 0.059 0.137 0.174 2.20

B 0.028 0.044 0.451 0.506 14.857 0.209 0.536 0.56 6.952

0.2 A 0.125 0.104 0.306 0.589 5.314 0.054 0127 0.165 2.086

B 0.04 0.064 0.408 0.528 10.582 0.137 0.365 0.39 4.907

0.3 A 0.133 0.113 0.287 0.600 5.169 0.051 0.121 0.159 2.019

B 0.048 0.078 0.379 0.543 9.044 0.111 0.301 0.326 4.137

0.4 A 0.140 0.119 0.272 0.608 5.066 0.049 0.117 0.155 1.974

B 0.054 0.088 0.357 0.555 8.232 0.097 0.266 0.290 3.714

0.5 A 0.144 0.125 0.26 0.615 4.983 0.047 0.114 0.152 1.941

B 0.059 0.097 0.340 0.564 7.724 0.088 0.243 0.268 3.442

0.6 A 0.148 0.129 0.250 0.622 4.914 0.046 0.112 0.150 1.915

B 0.063 0.103 0.325 0.572 7.374 0.082 0.227 0.252 3.248

on vacation. Hence, initially Pslow increases with 1. For this reason ¡LBP, ¡LBC, and ¡NSBP show an early downward trend. But as 1 further increases PSLOW declines as expected due to the high traffic intensity. Hence, ¡LBP and ¡NSBP reverse the direction of change. Due to the effect of PVAC and PIDL, this reversal occurs only at a later stage for ¡LBC. It is worth comparing the values of the measures in cases A and B. Even though the net service rate ¡1 + ¡i2 = 15 in both cases, the effect of the vacation parameter n becomes more predominant when ¡1 < ¡2. Due to this the measures PVAC and PIDL take smaller values and the measures ¡NS, ¡wtq, ¡lbp, and ¡NSBP take larger values in case B, compared to their values in case A.

(ii) As n increases, the mean duration of vacation 1/n decreases. But as the mean duration of vacation decreases, the probability of the expiry of the vacation without initiating the slow service increases. The chance of an early expiry of vacation always results in an increase in Pnorm and PVAC and a decrease in PSLOW. Note that PVAC + PSLOW decreases as n increases and PVAC + PSLOW < PNORM for any value of n in the given range. So PIDL increases with n. Thus, the proportion of time in which

Table 3: Case A: X = 12, fi1 = 10, = 5 and n = 1. Case B: X = 12, fi1 = 5, = 10 and n = 1.

e A/B PlDL Pvac Pslow Pnorm ^ns ^wtq ^lbp ^lbc ^nsbp

0.1 A 0.112 0.073 0.167 0.761 4.402 0.044 0.158 0.195 2.439

B 0.034 0.042 0.188 0.770 8.45 0.105 0.507 0.531 6.601

0.2 A 0.121 0.085 0.178 0.737 4.504 0.044 0.144 0.182 2.286

B 0.039 0.053 0.207 0.740 8.126 0.098 0.420 0.444 5.556

0.3 A 0.13 0.098 0.190 0.713 4.581 0.044 0.133 0.171 2.152

B 0.046 0.066 0.226 0.708 7.762 0.091 0.345 0.370 4.664

0.4 A 0.139 0.112 0.200 0.689 4.638 0.044 0.123 0.161 2.036

B 0.054 0.082 0.246 0.673 7.379 0.084 0.283 0.308 3.918

0.5 A 0.149 0.126 0.209 0.664 4.678 0.044 0.114 0.152 1.936

B 0.063 0.101 0.265 0.635 6.995 0.076 0.232 0.257 3.307

0.6 A 0.158 0.141 0.218 0.641 4.702 0.043 0.106 0.144 1.849

B 0.074 0.123 0.282 0.595 6.626 0.069 0.191 0.216 2.818

both servers work at the normal rate increases as n increases. Hence, the measures ¡NS, ¡wtq, ¡LBP, ¡LBC, and ¡NSBP decrease as n increases. The argument given the last paragraph holds good for the difference in magnitude of the measures in cases A and B.

(iii) As 0 increases, the service rate 0^2 of the second server during vacation mode of service increases. As a result, server 2 clears out customers at an increased rate in slow service mode. This produces an increase in PVAC, PSLOW, and PIDL and a decrease in Pn0rm as expected. Consequently, ¡LBP, ¡LBC and ¡NSBP decrease as 0 increases. The huge difference in the value of net service rate ¡1 + 0^2 between cases A and B, during vacation mode of service, is the reason for the pattern of behavior of ¡ NS in these two cases. Increase in 0 does not affect ¡wtq significantly in case A but it affects the measure in case B. This is because the effect of 0 becomes significant only when ¡2 is large compared to ¡1.

5. Concluding Remarks

In this paper, we studied an M/M/2 queueing model with heterogeneous servers. One server follows multiple vacation policy. But this server offers service at a lower rate during vacation if customers arrive. The other server remains in the system even when it is empty. The busy period of the system was analyzed in the stationary regime. Mean waiting time of a customer has been computed. Conditional stochastic decomposition of queue length has been derived. An illustrative numerical example to bring out the qualitative nature of the model has been presented.


The authors would like to thank the anonymous referees for their valuable suggestions and comments which greatly improved the presentation of the paper. C. Sreenivasan's research is supported by the University Grants Commission, Government of India, under Faculty Development Programme (Grant no. F.FIP/11th Plan/KLCA042TF02).


B. T. Doshi, "Queueing systems with vacations—a survey," Queueing Systems, vol. 1, no. 1, pp. 29-66, 1986.

H. Takagi, Queueing Analysis Volume 1: Vacations and Priority System, North-Holland, Amsterdam, The Netherlands, 1991.

N. Tian and Z. G. Zhang, Vacation Queueing Models, International Series in Operations Research & Management Science no. 93, Springer, New York, NY, USA, 2006.

X. Chao and Y. Q. Zhao, "Analysis of multi-server queues with station and server vacations," European Journal of Operational Research, vol. 110, no. 2, pp. 392-406,1998.

Y. Levy and U. Yechiali, "An M/M/s queue with servers' vacations," Information Systems and Operational Research, vol. 14, no. 2, pp. 153-163,1976.

B. Vinod, "Exponential queues with servers' vacations," Journal of the Operational Research Society, vol. 37, no. 10, pp. 1007-1014,1986.

E. P. C. Kao and K. S. Narayanan, "Analyses of an M/M/N/ queue with server's vacations," European Journal of Operational Research, vol. 54, no. 2, pp. 256-266,1991.

L. Servi and S. Finn, "M/M/1 queue with working vacations (M/M/1/WV)," Performance, vol. 50, pp. 41-52, 2002.

J. Kim, D. Choi, and K. Chae, "Analysis of queue-length distribution of the M/G/1with working vacations," in Proceedings of the International Conference on Statistics and Related, Honolulu, Hawaii, USA, 2003.

101 D. A. Wu and H. Takagi, "M/G/1 queue with multiple working vacations," Performance Evaluation,

vol. 63, no. 7, pp. 654-681, 2006. 111 Y. Baba, "Analysis of a GI/M/1 queue with multiple working vacations," Operations Research Letters,

vol. 33, no. 2, pp. 201-209, 2005. 121 N.-S. Tian, J.-H. Li, and Z. G. Zhang, "Matrix analytic method and working vacation queues—a survey," International Journal of Information and Management Sciences, vol. 20, no. 4, pp. 603-633, 2009. 131 J. Li and N. Tian, "The M/M/1 queue with working vacations and vacation interruptions," Journal of

Systems Science and Systems Engineering, vol. 16, no. 1, pp. 121-127, 2007. 141 M. Zhang and Z. Hou, "Performance analysis of MAP/G/1 queue with working vacations and

vacation interruption," Applied Mathematical Modelling, vol. 35, no. 4, pp. 1551-1560, 2011. 151 M. F. Neuts and Y. Takahashi, "Asymptotic behavior of the stationary distributions in the GI/PH/c queue with heterogeneous servers," Zeitschrift fur Wahrscheinlichkeitstheorie und Verwandte Gebiete, vol. 57, no. 4, pp. 441-452,1981. 161 B. Krishna Kumar and S. Pavai Madheswari, "An M/M/2 queueing system with heterogeneous servers and multiple vacations," Mathematical and Computer Modelling, vol. 41, no. 13, pp. 1415-1429, 2005.

171 W. Szpankowski, "Stability conditions for multidimensional queueing systems with computer applications," Operations Research, vol. 36, no. 6, pp. 944-957, 1988. 181 L.I. Sennott, P. A. Humblet, and R. L. Tweedie, "Mean drifts and the nonergodicity of Markov chains,"

Operations Research, vol. 31, no. 4, pp. 783-789,1983. 191 G. Latouche and V. Ramaswami, Introduction to Matrix Analytic Methods in Stochastic Modeling, ASA-SIAM Series on Statistics and Applied Probability, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, Pa, USA, 1999. 201 M. F. Neuts, Matrix-Geometric Solutions in Stochastic Models, vol. 2 of Johns Hopkins Series in the Mathematical Sciences, Johns Hopkins University Press, Baltimore, Md, USA, 1981, [1994 version is Dover Edition1 .

211 M. F. Neuts and D. M. Lucantoni, "A Markovian queue with N servers subject to breakdowns and

repairs," Management Science, vol. 25, no. 9, pp. 849-861,1979. 221 J. G. Kemeny and J. L. Snell, Finite Markov Chains, The University Series in Undergraduate Mathematics, D. Van Nostrand, Princeton, NJ, USA, 1960. 231 V. P. Singh, "Two-server Markovian queues with balking: heterogeneous versus homogeneous servers," Operations Research, vol. 18, no. 1, pp. 145-159,1970.

Copyright of International Journal of Stochastic Analysis is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.