We a give deterministic (sample path) proof of a result that
extends the Pollaczek-Khintchine formula for a multiple vacation single-server queueing
model. We also give a conservation law for the same system with multiple classes. Our
results are completely rigorous and hold under weaker assumptions than those given in
the literature. We do not make stochastic assumptions, so the results hold almost surely
on every sample path of the stochastic process that describes the system evolution. The
article is self contained in that it gives a brief review of necessary background material.
1. Introduction
Consider a single-server queue with multiple vacations general arrival process and general-service times. The server works until all customers in the queue are served then takes a vacation; the server takes a second vacation if when he is back, there are no customers waiting, and so on, until he finds one or more waiting customers at which point he resumes service until all customers, including new arrivals, are served. A vacation is not initiated when there are customers in the system. Under the stochastic assumptions of stationarity and i.i.d. inter-arrival and service times, it is known that the mean customer delay in the queue is the sum of two components: mean queue delay and mean vacation time.
This model has applications in communication systems and repair systems among others. Choudhury [1] considers a batch arrival queue with a single vacation where a server takes exactly one vacation after the end of each busy period. This model has applications in manufacturing systems of job-shop type. The paper derives the steady-state distribution of queue length, busy period, and unfinished work using Laplace-Stieltjes transformation approach. Boxma and Groenendijk [2] give (pseudo) conservation laws for multiqueue single server systems with cyclic service. Their work is based on well-known stochastic decomposition results. Green and Stidham Jr. [3] use a sample-path approach (similar to ours) to address conservation laws and their applications to scheduling and fluid systems. They use the achievable region method to solve stochastic problems utilizing mathematical programming formulations.
References also include Bisdikian [4] who presents a simple method for obtaining conservation laws for single-server queues with batch arrivals and multiple classes using the ASTA (El-Taha [5]) property. Relevant references also include Green and Stidham [3], Boxma and Groenendijk [2], and Choudhury [1]. There is a vast literature on queues with vacations, and conservation laws. The reader may consult Takagi [6] and references therein.
Most analyses give the transform of the waiting time and queue length distribution and then invert these transforms to give explicit expressions. The reader is left with the task of understanding the details of the mathematical analysis. The contribution of this article is to give completely rigorous intuitive proofs that avoid transform methods. We accomplish this by using sample-path analysis. In other words, our proofs are intuitive, completely rigorous and hold pathwise in the sense that they are true on every realization of the stochastic process of interest. Our results turn out to be valid under weaker assumptions that those required in the literature.
This paper is organized as follows: in Section 2, we give preliminary background results that are used in proving main results. In particular, we review H=λG relation and give sample-path proofs for the residual service time and Pollaczeck-Khintchine formula for M/G/1 queues. In Section 3, we focus on systems with multiple vacations and give the main result which relates virtual delay at any given time with systems actual delay and vacation times. In Section 4, we give a sample path proof of the conservation law and consider a few special cases.
2. Preliminary and Background Results
In this section, we review a few preliminary results that are used in the proof of the main result. Our proof uses the sample path relation H=λG which is a generalization of the well-known Little's formula.
We are given a deterministic sequence of time points, {Tk,k≥1}, with 0≤Tk≤Tk+1<∞, k≥1, and we define N(t):=max{k≥1:Tk≤t}, t≥0, so that N(t) is the number of points in [0,t]. We assume that Tk→∞ as k→∞, so that there are only a finite number of events in any finite time interval (N(t)<∞ for all t≥0), and we note that N(t)→∞ as t→∞, since Tk<∞ for all k≥1. Associated with each time point Tk, there is a function fk:[0,∞)→[0,∞). The bivariate sequence {(Tk,fk(·)),k≥1} constitutes the basic data, in terms of which the behavior of the system is described. We assume that fk(t) is Lebesgue integrable on t∈[0,∞), for each k≥1.
With H(t) and Gk defined by H(t)=∑k=1∞fk(t) and Gk=∫0∞fk(t)dt, respectively, define the following limiting averages, when they exist:λ∶=limt→∞t-1N(t),H∶=limt→∞t-1∫0tH(s)ds,G∶=limn→∞n-1∑k=1nGk.
Following Stidham Jr. [7], Heyman and Stidham Jr. [8] and El-Taha and Stidham Jr. [9] suppose that the bivariate sequence {(Tk,fk(·)),k≥1} satisfies the following condition.
Condition A.
There exists a sequence {Wk,k≥1} such that
Wk/Tk→0 as k→∞; and
fk(t)=0 for t∉[Tk,Tk+Wk).
In economic terms, Condition A says that all the cost associated with the kth point (e.g., the kth customer) is incurred in a finite time interval beginning at the point (e.g., the arrival of the customer), and that the lengths of these intervals cannot grow at the same rate as the points themselves, as k→∞. This is a stronger-than-necessary condition for H=λG, but it is satisfied in most applications to queueing systems, in which the time points Tk and Tk+Wk correspond to customer arrivals and departures, respectively, and it is natural to assume that customers can only incur cost while they are physically present in the system. The following theorem is given by El-Taha and Stidham Jr. [9, Chapter 6].
Theorem 2.1.
Suppose t-1N(t)→λ as t→∞, where 0≤λ≤∞, and Condition A holds. Then if n-1∑k=1nGk→G as n→∞, where 0≤G≤∞, then t-1∫0tH(s)ds→H as t→∞, and H=λG, provided λG is well defined.
We next consider two applications of H=λG. The first deals with residual service times, and the second application is the well-known Pollaczek-Khintchine formula for M/G/1 queues.
Residual Service Times
Let Tk with T0=0, Tk≤Tk+1, k=0,… be a deterministic point process, and Ak=Tk+1-Tk be the kth interevent. Let R(s) be the residual time, that is, time until next event, that is,
R(s)=∑k=1∞(Tk+1-s)1{Tk≤s<Tk+1}.
Define the following limits when they exist:
EA=limn→∞n-1∑k=1nAk;EA2=limn→∞n-1∑k=1nAk2;R=limt→∞t-1∫0tR(s)ds.
We interpret EA as the asymptotic average time between events, EA2 as the asymptotic second moment of the time between events, and R as the asymptotic (long-run) time-average residual time of R(s). In a queueing setting, it is sometimes useful to think of Ak as the service time of kth arrival, and R(s) as the residual service time of the customer in service at time s, but we shall see other examples. Now, we state the following preliminary result.
Lemma 2.2.
The asymptotic average residual time is
R=EA22EA.
Proof.
This proof uses H=λG. Let fk(t)=(Tk+1-t)1{Tk≤t<Tk+1}. Now, H(t)=∑k=1∞fk(t) is the residual time at a randomly given time, and H is the asymptotic residual time. Moreover, Gk=∫0∞fk(t)dt=Ak2/2; thus, G=limn→∞n-1∑k=1nAk2/2=EA2/2. Therefore (with λ=1/EA)
H=λEA22=EA22EA,
which completes the proof.
An alternative proof of this result uses the simpler relation, Y=λX, of (El-Taha and Stidham Jr. [9]). Now, consider an M/G/1 queue with multiple vacations as defined in Section 3. Let Sk and Dk be the service requirement and queueing delay of kth arrival, respectively, and let fk(t)=(Sk-(t-Tk-Dk))1{Tk+Dk≤t<Tk+Dk+Sk}. Then we immediately see that the mean residual service time of the customer in service is given by R=λES2/2, noting that λ=1/EA in Lemma 2.2. Note that R may also be written as (ρ=λES)R=ρES22ES
which is the product of the probability that the server is busy and the mean residual service time. Similarly, thinking of Tk as time instants when a server takes a vacation, the asymptotic mean residual vacation time, VR, is also given by VR=EV22EV.
Pollaczek-Khintchine Formula
Consider any G/G/1-FIFO queue. Let Tk, k≥1 be the time instant of kth arrival; Sk and Dk be the kth arrival service requirement and delay (time in queue) respectively. Also let N(t)=max{k:Tk≤t} be the total arrivals during [0,t]. Define the following limits when they exist:
λ=limt→∞N(t)t,ES=limn→∞n-1∑k=1nSk,ES2=limn→∞n-1∑k=1nSk2,ED=limn→∞n-1∑k=1nDk,ED2=limn→∞n-1∑k=1nDk2,ESD=limn→∞n-1∑k=1nSkDk.
Note that we use the suggestive expectation notation even though the quantities are defined pathwise. Let
fk(t)={Sk,ifTk≤t<Tk+DkSk-(t-Tk-Dk),ifTk+Dk≤t<Tk+Dk+Sk0,otherwise.
Now,
H(t)=∑k=1∞fk(t)
is the total amount of work in the system at a randomly given time, and H is the asymptotic average amount of work in the system. Moreover,
Gk=∫0∞fk(t)dt=SkDk+Sk22,
so that
G=limn→∞n-1∑k=1n[SkDk+Sk22]=ESD+ES22.
Therefore,
H=λ[ESD+ES22],
The first term of the r.h.s. is the total amount of work in the system associated with customers waiting (excluding the one in service), and the second term is the residual service time of the customer in service.
M/G/1 Model
For an M/G/1-queue, we make the additional assumptions that S and D are independent, and arrivals are Poisson. We also assume that ρ=λES<1. Then by PASTA, we have ED=H which implies ED=λ(ES)(ED)+λES2/2. Therefore,
ED=λES22(1-ρ).
In the next section, we extend Pollaczek-Khintchine formula to systems with multiple vacations using H=λG and maintaining a pure sample path approach.
3. Systems with Multiple Vacations
Consider any G/G/1-FIFO stable queue with multiple vacations of length {Vk,k≥1}. The server takes a vacation as soon as the number of customers in the system drops to 0 and continues to take successive vacations until the state N, where N is the number of customers in the system, is ≥1 when the server returns from vacation. The process of taking vacations is repeated at the end of the following busy period. A vacation is not initiated when there are customers in the system. Let {vk,k=1,2,…} be the start time of kth vacation. Let M(t) be the number of vacations up to time t so that M(t)=max{k:vk≤t}. Define the following limits when they exist: γ=limt→∞M(t)t,EV=limn→∞n-1∑k=1nVk,EV2=limn→∞n-1∑k=1nVk2.
Let ρ:=limt→∞t-1∑kN(t)Sk, then the system is stable if ρ=λES<1. Note that the ρ represents the long-run fraction of time the server is busy and 1-ρ is the long-run fraction of time the server is idle (i.e., on vacation).
Definition 3.1.
The two sequences {Sk,k≥1} and {Dk,k≥1} are said to be asymptotically pathwise uncorrelated if ESD=ESED.
Note that the requirement thatESD=ESED is weaker than the corresponding stochastic assumption that the random variables S and D, representing customer service times and delays, respectively, are independent.
Theorem 3.2.
Consider the G/G/1-FIFO multivacation model. Suppose that ED<∞ and EV<∞, and that {Sk,k≥1} and {Dk,k≥1} are asymptotically pathwise uncorrelated. Then virtual delay in the system is given by
H=ρED+ρES22ES+(1-ρ)EV22EV.
Proof.
Let Tk and vk be time instants of customer arrivals, and vacation starts, respectively. Let A={Tk,k=1,2,…}, B={vk,k=1,2,…}. Note that A∩B=ϕ. Let fk1(t) and fk2(t) be defined such that
fk1(t)=Sk1{Tk≤t<Tk+Dk}+(Sk-(t-Tk-Dk))1{Tk+Dk≤t<Tk+Dk+Sk};fk2(t)=(Vk-(t-vk))1{vk≤t<vk+Vk}.
Note that fk1(t) is the same as the fk(t) given by (2.9). Here, fk1(t) is the work remaining to be done for the kth arrival at time t; fk2(t) is the vacation time remaining for the kth vacation at time t. Consider the bivariate sequence {Tk,fk1(t)} and let H1(t):=∑k=1∞fk1(t), so that H1(t) is the total amount of work in the system at time t. Using H=λG, we obtain (see (2.13))
H1∶=limt→∞t-1∫0tH1(t)dt=λ[ESD+ES22].
Similarly, consider the second bivariate sequence {vk,fk2(t)} and let H2(t):=∑k=1∞fk2(t), Gk2=∫0∞fk2(t)dt=Vk2/2. Here, H2(t) is the residual vacation time in the system at time t, and Gk2 is the remaining vacation time encountered by the kth arrival. Applying H=λG, we obtain
G2=limn→∞n-1∑k=1n(Vk22)=EV22,H2∶=limt→∞t-1∫0tH2(t)dt=γ[EV22],
where γ is the unconditional vacation rate. We need to compute γ.
Recall that M(t)=max{k:vk≤t} is the number of vacations up to time t. Let I(t) be the status of the server at time t, that is, I(t)=0 if the server is idle at time t and 1 otherwise. Note that the server can be idle (on vacation) while customers are waiting for a vacation to end. It follows that
EV∶=limn→∞n-1∑k=1nVk=limt→∞M(t)-1∑k=1M(t)Vk=limt→∞M(t)-1∫0t1{I(t)=0}dt=limt→∞tM(t)-1(1-ρ),
so that
γ:=limt→∞t-1M(t)=1-ρEV.
Therefore,
H2=(1-ρ)EV22EV.
Let H(t):=H1(t)+H2(t) be the virtual delay, that is, total amount of work and remaining vacation time in the system at time t, and H=limt→∞t-1∫0tH(s)ds is the asymptotic (long-run) average amount of work and residual vacation in the system, in other words, H is the asymptotic average virtual delay for a randomly arriving customer. Therefore,
H∶=H1+H2=λESD+λES22+(1-ρ)EV22EV.
Now using the condition that Sk and Dk are asymptotically pathwise uncorrelated and simplifying, we obtain the result.
Equation (3.9) gives virtual delay H under conditions weaker than those given in Theorem 3.2. Specifically, (3.9) holds without the assumption that Sk and Dk are asymptotically pathwise uncorrelated.
M/G/1 Queue with Vacations
Suppose that S and D are independent, arrivals are Poisson, and vacations Vk are i.i.d. having a general distribution function with mean EV. Also assume that the system is stable, that is, ρ=λES<1. Then, by PASTA, we have ED=H, so Theorem 3.2 implies
ED=ρED+ρES22ES+(1-ρ)EV22EV.
Therefore,
ED=λES22(1-ρ)+EV22EV.
For systems with exponential vacations
ED=λES22(1-ρ)+EV.
Additionally, if service times are exponential with ES=1/μ, we obtain
ED=ρμ(1-ρ)+EV.
4. Conservation Laws
Now we consider multiclass G/G/1 queue with multiple vacations, that is, a single-server system with general (not necessarily i.i.d.) interarrival and service times. Conservation laws hold for a wide variety of service disciplines like FIFO, LCFS, service in random order, and other priority rules. For nonvacation models, any scheduling rule has to be work-conserving, and the server is never idle when there is work in the system, but for vacation models this is not possible due to server vacation so the rule is adapted so that the server is never idle except when on vacation. We define a work-conserving scheduling rule similar to El-Taha and Stidham Jr. [9, pages 204–211], but adapted for a server that takes vacations.
Work-Conserving Scheduling Rules
Consider a bivariate sequence, {(Tn,Sn),n=1,2,…}, where Tn and Sn are the arrival instant and work requirement, respectively, of nth arrival. A scheduling rule is said to be nonanticipative if the decision about which job to process at time t depends only on {(Tn,Sn),n=1,2,…,N(t)}, where N(t)=max{n:Tn≤t}, and possibly on decisions taken before time t. It is nonidling if the server does not take a vacation when there is at least one job in the system. Definition 4.1.
A single-server model with nonpreemptive work-conserving scheduling (WCS) rules consists of
a single-server working at unit rate,
a set of non-anticipative and non-idling scheduling rules.
We also assume that the scheduling rules are service time independent, and nonregenerative in the sense that the decision to schedule a job does not use any information from previous busy cycles. It is immediate from the definition of a work-conserving scheduling system that U(t)=∑n=1∞fk1(t), the total work in the system at time t is invariant with respect to WCS rules. It follows that the limiting average total work in the system,
U=limt→∞t-1∫0tU(s)ds,
is also invariant. The residual vacation at time t is defined by V(t)=∑k=1∞fk2(t), and the asymptotic (long-run) time-average residual vacation is given by
VR=limt→∞t-1∫0tV(s)ds.
So VR is also invariant. Now let H be the limiting average virtual delay in the system as defined in Theorem 3.2. Note that H and U are different due to the fact that the server takes multiple vacations. It is straight forward to see that H(t)=U(t)+V(t), therefore,
H=U+VR,
so that H is invariant. Therefore, we have the following.
Corollary 4.2.
Consider the G/G/1 multivacation model with WCS rules that are nonpreemptive, regenerative and service time independent. Suppose that ED<∞ and EV<∞, and that {Sk,k≥1} and {Dk,k≥1} are asymptotically pathwise uncorrelated. Then virtual delay in the system is given by
H=ρED+ρES22ES+(1-ρ)EV22EV.
Basically, this Corollary extends Theorem 3.2 from FIFO to any nonpreemptive WCS rule. Next, we focus attention on multiclass systems.Multiclass Multiple Vacation Model
Consider the multi-class GI/GI/1 multiple vacation model as given in Section 3. Suppose that the discipline is work conserving when the server is not on vacation, when the server starts service, it continues until all work is cleared before taking vacation. The scheduling rule is nonanticipative and independent of previous arrival and service times, and within each class j, j=1,…,J. That is, assume WCS rules. We give the following conservation law for a multi-class single-server system with multiple vacations.
Theorem 4.3.
Consider the G/G/1 stable multi-class multiple vacation model with WCS rules that are non-preemptive, regenerative, and within each class, they are service time independent. Also, suppose ED<∞ and EV<∞, and that for each j{Skj,k≥1} and {Dkj,k≥1} are asymptotically pathwise uncorrelated. The vector (ED1,…,EDJ) of expected actual queue delays per customer satisfies the following conservation law:
∑j=1JρjEDj=H-∑j=1JρjESj22ESj-(1-ρ)EV22EV.
Proof.
Let Hj and Uj, j=1,…,J be the class j virtual delay and work in the system respectively, where H=∑j=1JHj and U=∑j=1JUj. Now (4.3) becomes
∑j=1JHj=∑j=1JUj+VR,
Recall that H, U, and VR are invariant with respect to WCS rules. For each j, j=1,…,J, let λj, Sj, and Dj be the class j mean arrival rate, service time, and queue delay, respectively. Also let ρj=λjESj, and suppose that ∑j=1Jρj<1. Then it follows from (3.4)
EUj=λj[ESjDj+ESj22].
Using the fact that the Sj and Dj are asymptotically pathwise uncorrelated and ρi=λjESj, we have
EUj=ρjEDj+ρjESj22ESj.
Let Bj=ESj2/2ESj. It follows by (4.6) and (4.8) that
H=∑j=1JρjEDj+∑j=1JρjBj+(1-ρ)EV22EV
is invariant over all WCS rules, since we have working conserving scheduling rules. Now let λ=∑j=1Jλj, and let S be a random variable with distribution function F(t)=∑j=1J(λj/λ)P(Sj≤t). Then
∑j=1JρjBj=λES22
is also invariant over all WCS rules. Thus ∑j=1JρjEDj is also invariant, and we have the following conservation law satisfied by the vector (ED1,…,EDJ) of expected queue delays
∑j=1JρjEDj=H-∑j=1JρjBj-(1-ρ)EV22EV
which leads to (4.5), thus proving the theorem.
In the M/G/1 multi-class multiple vacation model, we derive an explicit expression for H=ED, using the fact that H is invariant and, therefore, is equal to the delay in the M/G/1-FIFO discipline which is given by (3.11). Hence, the conservation equation (4.5) becomes, by substituting (3.11) and (4.10) in (4.5),
∑j=1JρjEDj=λES22(1-ρ)+EV22EV-λES22-(1-ρ)EV22EV,
which simplifies to
∑j=1JρjEDj=λρES22(1-ρ)+ρEV22EV.
Theorem 4.3 can be used to construct conservation laws for waiting time in the system Wj, number of customers in the system Lj, and number of customers in the queue Ljq using the fact that EWj=EDj+ESj and Little's formula.In this paper, we give two primary results. The first result presented in Theorem 3.2 gives a relation between virtual delay and actual delay for single-server multiple vacation model under conditions that are weaker than those given in the literature. This result is extended to multi-class models where a conservation law is given in Theorem 4.3 which is our second result. We use sample path analysis which allows us to give rigorous arguments by focusing on one realization of the stochastic process that describes the system evolution.
ChoudhuryG.A batch arrival queue with a vacation time under single vacation policy200229141941195510.1016/S0305-0548(01)00059-41920586ZBL1010.90010BoxmaO. J.GroenendijkW. P.Pseudo-conservation laws in cyclic-service systems198724494996491383410.2307/3214218GreenT. C.Stidham,S.Jr.Sample-path conservation laws, with applications to scheduling queues and fluid systems2000361–3175199180696610.1023/A:1019183220080ZBL0966.90022BisdikianC.A note on the conservation law for queues with batch arrivals199341683283510.1109/26.231904ZBL0800.90453El-TahaM.2010New York, NY, USAJohn Wiley & SonsEncyclopedia of Operations Research and
Management ScienceTakagiH.1991Amsterdam, The NetherlandsNorth-Holland Publishing1149382Stidham,S.Jr.A last word on L=λW1974222417421040803810.1287/opre.22.2.417ZBL0274.60080HeymanD. P.Stidham,S.Jr.The relation between customer and time averages in queues198028498399410.1287/opre.28.4.983584901ZBL0439.60089El-TahaM.Stidham,S.Jr.1999Boston, Mass, USAKluwer Academic PublishingInternational Series in Operations Research & Management Science1784065