AN ALTERNATIVE APPROACH TO SOME CLASSICAL QUEUEING PROBLEMS
By
VATSLA VIR MANI
Department of Mathematics
THESIS SUBMITTED
IN FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF
DOCTOR OF PHILOSOPHY
to the
INDIAN INSTITUTE OF TECHNOLOGY, DELHI
NEW DELHI - 11 001 6
INDIA
NOVEMBER, 1993
CERTIFICATE
This is to certify that the thesis entitled
"AN ALTERNATIVE APPROACH TO SOME CLASSICAL QUEUEING PROBLEMS"
is a record of bonafide research work carried out by Miss Vatsla Virmani under my guidance and supervision. She has submitted her work to the Indian Institute of Technology, Delhi for the award of the degree of DOCTOR OF PHILOSOPHY.
I feel satisfied that the thesis has reached the standard, fulfilling the requirements of the regulations, relating to the degree. The results obtained in the thesis have not been submitted to any other University or Institute for the award of any degree or diploma.
( O.P. SHARMA ) Professor
Department of Mathematics
Indian Institute of Technology Hauz Khas, New Delhi- 110016
To my parents
ACKNOWLEDGEMENTS
My gratitude towards my revered Prof. O.P. Sharma of I.I.T. Delhi (Dept. of Maths) for his consistent encouragement, guidance and supervision throughout my research work is surely unbounded. The amount of time he spared, moreover his keen interest, care and valuable suggestions have given the present shape to the thesis.
I would fail if I do not express my grateful thanks to the Head, Prof. K.N. Mehta, and also all the worthy members of the Faculty of Maths for their generous help in fulfilling my ordained work for the last four years.
I have no words to lay at the feet of my beloved Parents- nay my guides and philosophers- who had laid down this object for me since my childhood. My sisters Leena and Swati equally deserve appreciation for their everlasting love, good wishes and moral support which have gone a long way to achieve my cherished goal.
Thanks are also due to Ms. Neelam Dhody for her excellent typing of the manuscript.
Last but not the least I thank my friends at I.I.T. Delhi for their sweet discourses and discussions pertaining to the topic of my research.
I.I.T. Delhi (VATSLA VIRMANI)
SYNOPSIS
The theory of queues deals with the investigation of stochastic behaviour of situations arising in connection with mass servicing involving random fluctuations. A queue is formed when demand for service is greater than the servicing capacity. These demands may emanate from finite or infinite sources. Generally a queueing system is studied under two states such as transient and steady state, and can be broadly classified into two categories namely Markovian and Non-Markovian. Mostly the researchers have investigated the classical queueing problems using a one-dimensional state model which has certain limitations. However recently a two-dimensional state model has been introduced to study the M/M/l/co queue which gives a better insight into the arrival and departure of units in the system. A part of this thesis (Chapters II, III & IV) is devoted to study certain Markovian queues using this model.
Another aspect of queueing system discussed in this thesis relates to the concept of 'two queues in parallel'. This system has received a lot of attention in the literature because it models many real life situations such as vehicles going through toll booths, jobs scheduled on a multiprocessor system etc.
Inspite of its practical applications there has yet not been found some simple solution even for its equilibrium probabilities. In this thesis an attempt is made to obtain the steady state probabilities for this model both analytically as well as
computationally after incorporating a modification which widens the scope of its applicability. The model is covered in chapters V and VI of the thesis.
A brief layout of the chapters is as follows:
Chapter I.
It contains the survey of the existing literature and explains some fundamental concepts of queueing theory emphasizing the aspects relevant to the thesis. This chapter also gives the basic definitions on which the later results in the thesis are developed.
Chapter 2.
It is a known fact that the derivation of transient results for M/M/l/m is through a complicated procedure. Transient solutions for this simple birth and death process reported in the literature involve either Bessel functions or integral forms of the Bessel functions. However, recently a two-dimensional state model for the standard M/M/l/w queue, which is initially empty, was introduced and the results obtained are independent of the Bessel functions and have been expressed as sum of two terms: one pertaining to steady state and other to the transient state. In this chapter this model has been used to obtain transient probabilities for M/M/l/c3 queue with arbitrary li' initial customers present in the system. The probability of n arrivals and k departures up to time t has been worked out in a closed form and
ii
p(r,t), the probability that there are r units present in the system at time t is obtained as
r - (A+g)t p r ; (At)k r+k+i , I 4
.1
"
m-1 p(r,t) = (1-p)p + e E (k-m) °
k! ml
k=0 m=0 + e-(X+g)t
E (fit) (µt) (4t)k
[ 1 1
k!(r+k-i)! (r+k)!(k-i)! ] k=0
where p = 0 A , whence by taking i=0, the results reduce to those obtained by Sharma (1990).
Chapter 3.
Saaty (1961) and Parthasaraket al. (1989) have obtained the time dependent state probabilities for two server Markovian queue in terms of Laplace transform variable or in complicated integrals, which are not amenable for further algebraic manipulations. In this chapter we use the two-dimensional state model developed in the second chapter to study the transient behaviour of M/M/2/m queue and obtain the probability of there being r units in the system at time t as
p(r,t) -
1 - p
e-(A+2g)t
E C(0,m,t), r = 0.
1 + p m=0 1 - p
= 2pr
- e-(A+2g)t r
p E C(r,m,t), r 1.
1 + p m=0
where p - A
2g , C(0,m,t) and C(r,m,t) are suitably defined. It is also shown that p(0,0) = 1.
iii
Chapter 4.
Busy period analysis forms an integral part of any queueing system. The literature available reports the busy period distribution in the transient state of a Poisson queue with infinite waiting space capacity for a single server queue only.
This chapter analysis the c-channel initial busy period distribution for M/M/c/m queue via two-dimensional state model. A closed form expression is obtained for the p.d.f. of the duration of a busy period initiated by an arbitrary initial unit as
b( t ) - i-c+1 e-(A+cg)t (At)k
(cgt)k+i-c+1
t
k=0 (k+i-c+1)! k!whence the moments of the length of the c-channel busy period is also obtained as follows:
E( ) _ (i-c+1) p-2
[ (i-c+2)p-1 + j
E 1
(-10 (cg)P (1-p)2P-1
(i-c+2+j)
p-1-j
=
P- lCi a-c+2-p), pi + (-1) p-1 1
x (i-c+2-p)p-1 pp
where, p = cg and (a)A
m = a(a+1)....(a+m-1).
Results for particular cases have also been worked out.
Chapter 5.
Conolly (1984) discussed the queueing model 'two queues in parallel' by taking the system with finite waiting room. In this chapter the above 'model is discussed along with the modification.
Here we assume that when one queue is empty and the other queue
i
vhas more than one customer waiting in it, then in that case one customer is transferred from this queue to the other queue. After deriving the difference equations satisfying the equilibrium probabilities, a simpler matrix method is proposed to solve them.
The results are also numerically worked out and a comparison is made between this model and Conolly's model for equal and unequal service rates. For the case when two servers are similar, working backward on the difference equations we obtain all the probabilities in terms of PIN •
Further taking qn = P(i + j = n), we obtain the probability that there are n customers in the system as
-2 0 = P00 ' qi = 2p1 p00, i=1,2,....,2N
A 1-p
where p = 24
and P0,0 — 1+p-2p2N+1
which shows that in this case our model works similar to the M/M/2/2N system. Certain other important results have also been worked out.
Chapter 6.
In this chapter we consider the system with two similar servers and the customer on arrival joins the shorter queue. If the difference between the longer and the shorter queue reaches some prescribed value N, then one customer leaves the longer queue and goes
to the shorter queue. The first part of this chapter deals with an analytic method to calculate the equilibrium probabilities and the second part discusses a numerical technique
to obtain these probabilities.
Defining the generating function as . . . xi'
Q3 = . +1 1=0p1
' 3
j= 0,1,....,N -1 we get a matrix equation for Q .(x) as
A(x)Q = b,
where A(x) is a NxN matrix, b and Q are column vectors and solving Q
det(A-(x)) it we obtain .(x) - 3
det(A(x))
and then differentiating Qi(x) with respect to x we get
pi,j+i V i and j.
Also introducing a power series expansion of p1.
,3 . as functions of egivenbypi,J =E9i+j
E19k
1b..
13(k), we obtain a recurrence k=0
relation and then using the fact that the sum of the probabilities is one we calculate bi,j(k) V k. It is remarked that the above described numerical method is economical in time, as compared to the method implemented by Gertsbakh (1984).
vi
CONTENTS
CHAPTER 1 : INTRODUCTION
Page No.
1.1 BASIC DEFINITIONS AND
FUNDAMENTAL CONCEPTS
2
1.2
QUEUES5
1.3
NOTATION8
1.4
QUEUEING PROCESS9
1.5
BIRTH AND DEATH PROCESS 111.6
BASIC RELATION15
1.7
LITERATURE REVIEW16
CHAPTER 2 : A TWO DIMENSIONAL STATE MODEL FOR WM/1/03 QUEUE
2.1
INTRODUCTION29
2.2
THE MODEL AND SOLUTION33
2.3
SUM IS UNITY46
2.4
NUMERICAL COMPUTATIONS49
CHAPTER 3 : A NEW APPROACH TO THE M/M/2 QUEUE
3.1
INTRODUCTION51
3.2
MODEL AND SOLUTION54
3 . 3
TIME DEPENDENT SOLUTIONFOR AN EMPTY SYSTEM
69
3 . 4
APPENDIX72
CHAPTER 4 : THE BUSY PERIOD ANALYSIS OF M/M/c/w QUEUEING SYSTEM
4.1
INTRODUCTION76
4.2
THE MODEL AND SOLUTION80
4.3
MOMENTS91
4.4
APPENDIX96
CHAPTER 5 : AUTOSTRADA QUEUEING PROBLEM - A MATRIX APPROACH
5.1
INTRODUCTION100
5.2
THE MODEL104
5.3
THE SOLUTION107
5.4
PARTICULAR CASE112
5.5
CONOLLY'S MODEL123
5.6
NUMERICAL INVESTIGATION125
CHAPTER 6 : THE MODIFIED SHORTER QUEUE PROBLEM
6.1
INTRODUCTION136
6.2
EQUILIBRIUM EQUATIONS137
6.3
ANALYTICAL METHOD FOR CALCULATINGEQUILIBRIUM PROBABILITIES
139
6.4
PARTICULAR CASES148
6.5
NUMERICAL TECHNIQUE FOR CALCULATINGSTATE PROBABILITIES