CC BY
34
NUMBER OF CALLS IN A CYCLIC-WAITING SYSTEM
L. Lakatos, S. V. Serebriakova •
Eotvos Lorand University, Budapest, Hungary, V. M. Glushkov Institute of Cybernetics of NAS of Ukraine, Kyiv, Ukraine e-mail: [email protected], svitlana.pustova @ gmail.com
Abstract
A single-server cyclic-waiting system, i.e. a queuing system, in which a call gets service either at the time of arrival or after some cycle time T, has been studied in the article. The distribution functions of the intervals between arrivals and service times have exponential form; the service discipline is FCFS. Analytical results for the mean number of calls in the system are indicated. Numerical and graphical results have been given.
Keywords: cyclic-waiting system, queuing system, FCFS-service.
1 INTRODUCTION
Queuing systems (Gnedenko 1989) have several applications in different technical and biological systems, especially in telecommunication (Pustova 2010) and network systems. In order to adequately describe a real queuing system, one should consider different aspects of the real system, such as number of servers, availability of buffer and its type, order and duration of service, etc. Retrial queuing systems (Artalejo & Gomez-Corral 2008, Kovalenko & Koba 2010) are one of important tools in research and optimization of characteristics of different systems. Such systems take into account not only the initial flow of calls, but also the repeated ones. Retrial queuing systems have the following behavior. If all servers are busy upon call's arrival the call doesn't get lost, but goes to the virtual queue (an orbit), and after some time (a cycle) repeats its attempt to get service. Ignoring the retrial effect can cause significant accuracy errors in research (Aguir et al. 2004, Aguir et. al. 2008).
Retrial queuing systems may also take into account the service discipline: random service, FCFS (first come first served), LCFS (last come first served), etc. The analysis of such disciplines might increase the complexity of research, so it becomes hard to obtain analytical results.
We consider a so-called cyclic-waiting system (Lakatos 2010) having a special feature comparing to the common retrial systems, namely, calls get service in the order of arrival: if there is any call in the orbit, it will be served earlier than the arriving new ones. The calls from the orbit are served after some time, multiple to a given cycle time T.
The process of airplane landing initiated research of cyclic-waiting systems (Lakatos 1994 & 1998). An airplane lands on a given route. If a runway is available the airplane starts the landing procedure immediately. If there is a threat to the non-compliance of the separation standards or
there already are some planes in the waiting area (in the circle), waiting for permission to land, the aircraft is sent to the queue (to the circle). While on circle, the aircraft waits for permission to land. After the permission is given, the airplane has to get to the corresponding geometric point on the circle route, and start the landing procedure from there. The time period from the moment of landing clearance to the moment until the plane reaches the corresponding geometric position is called "idle": the system is ready to take a plane, but the landing is possible only on reaching the corresponding point on the circle. Thereby, the airplane landing process can be described as a cyclic-waiting system.
Kovalenko (2002) considered the multichannel retrial queuing system with a constant cycle orbit, and gave estimation of the loss probability in light traffic. Rogiest et al. (2006) applied cyclic-waiting system to the research of optical buffers. Lakatos & Zbaganu (2007) investigated different generalizations of the initial model. Karasz (2008) studied some cyclic-models with two different types of calls. Kovalenko (2015) started to study a new class of multiple-cyclic waiting systems, namely, considered a two-cyclic queuing system.
2 ANALYTICAL MODELING
2.1 Problem statement
We consider a single server retrial queuing system with exponentially distributed intervals between arrivals and service times, with corresponding parameters 2 and ¡i (Fig. 1). If the server is free, the arrived call gets service immediately. Otherwise, if the server is occupied, an arrived call goes to the orbit (i.e. a virtual queue), in which it retries its requests with a determined cycle time T > 0. The service is defined by FCFS (first come first served) discipline. It means that if there's any call cycling in the orbit, the call from the initial flow doesn't get service immediately, but goes to the orbit instead and after some time, the multiple of T, repeats its attempt to get service.
Figure 1. A cyclic-waiting system structure.
2.2 Ergodicity conditions
We will use results, obtained in (Lakatos 1994 & 2010) to find some characteristics of the system.
Let tn, n > 0 denote the moment of the service beginning of the nth call. Let's introduce an embedded Markov chain, the states of which are defined as the number of calls in the system at moments (tn - 0). The matrix of the chain's transition probabilities is defined as
ao ai a2 a3
ao ai a2
0 bo \ K
0 0 K \
the elements of which are defined by generating functions
J _
$ i ß Xz (1-е-^)е-Я(1-2)Г
(2)
¡=0
(where
X + p
= \e~XxiJe-^dx
is the probability of calls' non-arrival, and multiplier z in (2) represents the mandatory occurrence of a call), and
" 1
_ e-U2-z)T jjj _ e-[X(\-2)+ß\T j _
Го ' {\-e-lT)[\-e-W-^T] [2-z
Я ^ _ e-[A(2-z)+M]r _ g-A(l-z)r j
(3)
X{2-z) + fi
Let pi, i = 0,1,... be the ergodic probabilities of the chain. Hence, the generating function for the number of calls in the system is
= = />0-jp^j-'
1=0
where A(z) and B(z) are defined by (3) and (4) respectively, and p0 is defined by
/>o = l—
Then the ergodicity condition is
1-е
-XT
ре-ХГ(\-е-КГ)
<1.
(4)
(5)
(6)
2.3 Mean number of calls in the system
Let us determine mean number of calls in the system, in other words P'(1). Let us denote
N[z) = B{z)-z, S(z) = XzB{z) + pB(z) - XzA(z) - yzA(z),
then (4) can be written as
If z = 1 then
Let us obtain P'(z):
P\z)-
m=g
S(z) N{z)
N(z)
N( 1) = 1-1 = 0, S(l) = X + ß-X-in = 0.
S'(z)N(z)-S(z)N'(z) n2(z)
(7)
If z = 1 in (7) the indeterminate form of 0 / 0 occurs. By applying L'Hopital's rule to (7) we get
(S'(z)N(z)-S(z)N\z)y _ S"(z)N(z) + S'(z)N'(z)-S'(z)N'(z)-S(z)N"(z)
(N2(z)Y
2 N{z)N\z)
S"(z)N{z)-S{z)N'\z) 2N(z)N'{z) .
39
Similarly, if z = 1 in (8) the indeterminate form of 0 / 0 occurs. Analogously, for (8) we obtain (S"(z)N(z) - S(z)N"(z)Y S"\z)N{z) + S"{z)N\z) - S\z)N"{z) - S{z)N"\z)
(2N(z)N'(z))' 2[{N'{z)) + N(z)N"(z)]
If z = 1 (9) is transformed to
S"(1)7V"(1)-S'(1)7V"(1)
(9)
P'(l)= w """ w. (10)
2[N\l)f
Let us find the unknown derivatives in (10):
N\z) = B\z)-1, (11)
N"{z) = B"(z), (12)
S\z) = XB(z)+XzB\z) + fiB'(z) - XA(z) - XzA\z) - nA(z) - fizA\z), (13)
S"(z) = AB'(z) + XB'(z) + XzB"(z) + fiB"(z) -
-AA'(z) - XA\z) - AzA"(z) - nA\z) - nA\z) - f.izA"(z). (14)
If z = 1 the (11) - (14) are transformed to
AT'(1) = 5'(1)-1, (15)
= (16) S'(l) - B\\){1 + A\ 1)(A + & - n, (17)
S"(l) = 2XB\\) + XB"{ 1) + ¡iB"(\) - 2XA'(\) - AA"( 1) - 2juA"(l) =
= 2AB'(1)+B"( 1)(A + M)~ 2(A + M)A'( 1)-(A+¡i)A"{\). (18)
After some computations the derivatives for A(z) and B(z) that are part of the (15) - (18) can be written as
Ad+ W (19)
(X + Ml-e-P)
A"( 1) + +
(A + Ju)(l-e^)2 ,
B'(l) =-^-—U + H + X2T + X/nTe-(X+»)T - (A + /i)(e"A:r + e^) +
+(X + fi)e~^)T-{X +¡¿We'*1), (21)
+-2XTe~^ l{l_e-Xr){x_e-ßT
(1 -e-^f {l-e-XT){\-e-P)2\
A + /1
s
-2lTe~XT{\ - 2A7e_'"r(l -e~XT)~ X2T2eXT(\ - e~ßT)+
+ll2T2e-(l+ß)T _x2T2e-ßT(1 _e-xr) + 2A3r(l-g-^>r) _
(22)
Therefore, formulas (11) - (22) define expression (10) entirely. Due to cumbersome form, the final expression of P'(z) (as well as P'(1) ) isn't given here, though it can be easily obtained with the help of formulas (10), (11) - (22).
2.3 Analysis of the modelling results
Let us denote the left side of (6) as H (we will use this notation for graphical results):
U_X l-e-f Ve-^il-e-^)'
According to formula (6) we have obtained graphical representation for egrodicity conditions of the system under consideration (see Figs 2-3): everything below the line H = 1 qualifies for ergodicity condition (6).
_X_
--T= 1 ..... T- 2---T= 3 • T- 5 -T- 7
- - T= 10 -H=1
Figure 2. Graphical representation for ergodicity conditions, where T = {1, 2, 3, 5, 7, 10}, 2 = [0.1;0.9], u = 1.0.
' 2 ' 3 ' 4 ' 5 6 7 8 9 10
r-0.1 - -- - r-0.5 ---r=i.o o T= 1.5
T- 2.0 - - T-2.3 -H= 1
Figure 3. Graphical representation for ergodicity conditions, where T = {0.1, 0.5, 1.0, 1.5, 2.0, 2.3}, 2=1.0, u= [1.11;10.0].
The values of 2 and u for the results on Figures 2-3 have been chosen so that the ratio 2/ u = [0.1; 0.9]. From Figure 2 one may see that with the increase of the value of cycle time T the left part of the (6) increases sharply. The same situation persists on Figure 3.
Using results from Figures 2-3, we have also obtained a graphical dependence for the mean number of calls P'(1) from the mean service time 1 / ¡i and cycle time T (see Fig. 4), as well as from the intensity of arrivals 2 and cycle time T (see Fig. 5).
Figure 4. Graph for the mean number of calls (2 = 1.0).
Figure 5. Graph for the mean number of calls (/ = 1.0).
From Figures 4-5 one may see that the mean number of calls rises with the increasing service time as well as with the increasing intensity of arrivals. On both graphs (see Figs 4-5) P'(1) increases with the increasing cycle time T.
3 CONCLUSIONS
In the article we have obtained analytical and graphical results for the mean number of calls in the single-channel cyclic-waiting system. Ergodicity conditions have been shown. The analytical result for P'(1) can be used, for example, to obtain the average time a call spends in a system: from Little's law W = P'(1) / X.
The cyclic-waiting system under consideration may have different practical applications (Koba & Pustova 2012, Pustova 2016), such as an airplane landing process, a simple call center, a telephone answering machine, an optical fiber buffer (Rogiest et al. 2006, Langenhorst et al. 1996), an optical delay line, a microring resonator (Bogaerts et al. 2012), a compact optical buffer with ring resonators (IBM, 2006). All these systems have a common feature: a new call, which can't get service until there is at least another one in the system, goes to the orbit and cycles there.
References
[1] Aguir, M. S., Karaesmen, F., Aksin, Z., Chauvet F. 2004. The impact of retrials on call center performance. Operations Research 26: 353-376.
[2] Aguir, M. S., Aksin, O. Z., Karaesmen, F., Dallery Y. 2008. On the interaction between retrials and sizing of call centers. European Journal of Operational Research 191: 398-408.
[3] Artalejo, J.R., Gómez-Corral, A. 2008. Retrial queueing systems: a computational approach. Berlin: Springer-Verlag.
[4] Bogaerts, W., De Heyn, P., Van Vaerenbergh, T., De Vos, K., Kumar Selvaraja, S., Claes, T., Dumon, P., Bienstman, P., Van Thourhout, D., Baets, R. 2012. Silicon microring resonators. Laser Photonics 6(1): 47-73.
[5] Gnedenko, B.V., Kovalenko I.N. 1989. Introduction to queuing theory. Boston: Birkhauser.
[6] IBM. 2006. Compact optical buffer with ring resonators : http://domino.research.ibm.com/ comm/research_projects.nsf/pages/photonics.ringbuffer.html
[7] Langenhorst, R., Eiselt, M., Pieper, W., Grosskopf, G., Ludwig, R., Kuller, L., Dietrich, E., Weber, H.G. 1996. Fiber loop optical buffer Journal of Lightwave Technology 14(3): 324-335.
[8] Karasz, P. 2008. A special discrete cyclic-waiting queuing system. Central European Jornal of Operational Research 16(4): 391-406.
[9] Koba, E.V., Pustova, S.V. 2012. Lakatos queuing systems, their generalization and application. Cybernetics and Systems Analysis 48(3): 387-396.
[10] Kovalenko, I.N. 2002. Loss probability in the M/G/m queuing system with T-retrial calls in light traffic (in Russian) Dopovidi NANU 5: 77-80.
[11] Kovalenko, I.N. 2015. A two-cyclic queuing system. Cybernetics and Systems Analysis 51(1): 51-55.
[12] Kovalenko, I.N., Koba, E.V. 2010. On the classification of retrial queuing systems. Cybernetics and Systems Analysis 46(3): 420-425.
[13] Lakatos, L. 1994. On a simple continuous cyclic-waiting problem Annales Universitatis Scientarium Budapestinensis (Sectio Computatorica) 14: 105-113.
[14] Lakatos, L. 1998. On a simple discrete cyclic-waiting queueing problem Journal of Mathematical Sciences 92: 4031-4034.
[15] Lakatos, L. 2010. Cyclic-waiting systems. Cybernetics and Systems Analysis 46(3): 477-484.
[16] Lakatos, L, Zbaganu, G. 2007. Waiting time in cyclic-waiting systems. Annales Universitatis Scientarium Budapestinensis (Sectio Mathematica) 27: 217-228.
[17] Pustova, S.V. 2010. Investigation of call centers as retrial queuing systems. Cybernetics and Systems Analysis 46(3): 494-499.
[18] Pustova S.V. 2016. Applications of the cyclic-waiting systems (in Ukrainian). Dopovidi NANU - submitted.
[19] Rogiest, W., Laevens, K., Fiems, D., Bruneel, H. 2006. Analysis of a Lakatos type queuing system with general service times Proc. of ORBEL 20 (Quantitative Methods for Decision Making): 95-97.
