Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)2219210.22363/2658-4670-2019-27-1-5-20Research ArticleHeavy outgoing call asymptotics for retrial queue with two way communication and multiple types of outgoing callsNazarovAnatoly A<p>Professor, Doctor of Technical Sciences, Head of Department of Probability Theory and Mathematical Statistics, Institute of Applied Mathematics and Computer Science</p>nazarov.tsu@gmail.comPaulSvetlana V<p>Candidate of Physical and Mathematical Sciences, Assistant Professor of Department of Probability Theory and Mathematical Statistics, Institute of Applied Mathematics and Computer Science</p>paulsv82@mail.ruLizyuraOlga D<p>Master’s Degree Student of Institute of Applied Mathematics and Computer Science</p>oliztsu@mail.ruNational Research Tomsk State University1512201927152020112019Copyright © 2019, Nazarov A.A., Paul S.V., Lizyura O.D.2019<p>In this paper, we consider a single server queueing model M |M |1|N with two types of calls: incoming calls and outgoing calls, where incoming calls arrive at the server according to a Poisson process. Upon arrival, an incoming call immediately occupies the server if it is idle or joins an orbit if the server is busy. From the orbit, an incoming call retries to occupy the server and behaves the same as a fresh incoming call. The server makes an outgoing calls after an exponentially distributed idle time. It can be interpreted as that outgoing calls arrive at the server according to a Poisson process. There are N types of outgoing calls whose durations follow N distinct exponential distributions. Our contribution is to derive the asymptotics of the number of incoming calls in retrial queue under the conditions of high rates of making outgoing calls and low rates of service time of each type of outgoing calls. Based on the obtained asymptotics, we have built the approximations of the probability distribution of the number of incoming calls in the system.</p>retrial queueing systemincoming callsoutgoing callsasymptotic analysis methodGaussian approximationсистема очередей на повторный процесс, входящие вызовы, исходящие вызовы, метод асимптотического анализа, гауссовское приближение<p>Introduction Retrial queueing systems are characterized by the following distinctive fea- ture: a customer who cannot receive service remains in the system and tries to occupy the server after some random delay. The pool of unsatisfied customers is called the orbit. Retrial queues have applications in telecommunication, computer networks and in daily life [1, 2]. In retrial queues idle time of the server is the downtime and it should be reduced to increase the efficiency of the system. We consider systems where operator not only receives calls from outside but also makes outgoing calls in the idle time. In queueing theory a model with this feature have been considered previously [3]. However, the retrial behaviour of customers is not taken into account. In call centers operators could receive arriving calls but as soon as they have free time and are standby mode they could make outgoing calls [4-7]. Systems with this server behaviour are called retrial queues with two way communication. Retrial Queues with two way communication have been studied recently [8-11]. In these papers Markovian models with two way communication were considered. Model of retrial queue with two way communication and multiple types of outgoing calls was considered by Sakurai and Phung-Duc [12]. For this model numerical algorithm of calculating joint stationary distribution of system state was obtained. Multiserver retrial queue with two way communication was studied in [13]. Recently the two way communication retrial queues with finite source [14], with server-orbit interaction [15, 16], with finite orbit [17], with breakdowns [18] and with a constant retrial rate [19] were considered. Asymptotic analysis methods have applications in queueing theory. Nazarov, Paul and Gudkova propose an asymptotic analysis method to research M |M |1|1 retrial queue with two way communication under low rate of retri- als condition [20]. Nazarov, Paul and Phung-Duc extended this model to MMPP|M |1|1 retrial queues and derived asymptotics in heavy outgoing call conditions [21]. In this paper, we consider retrial queue with two way communication and multiple types of outgoing calls. We assume that each type of outgoing calls has different rate and service times follow distinct exponential distributions. The main aim of this paper is to derive asymptotics for the model in two limit conditions: i) high rate of outgoing calls and ii) low service rate of outgoing calls. In both cases, the number of incoming calls in the system increases. The rest of the current paper is organized as follows. In Section 2 and 3, we describe the model in detail and preliminaries for later asymptotic analysis. In Section 4 and 5, we present our main contribution to the model. In Section 6 we show the ranges of parameters under which our approximations are usable. Section 7 is devoted to concluding remarks. Model and preliminaries Model description Figure 1 shows the structure of the model. We consider a single server retrial queue with two way communication and multiple types of outgoing calls. Incoming calls arrive at the system according to a Poisson process with rate and try to occupy the server for an exponentially distributed time with rate 1. Incoming calls that find the server busy join the orbit and repeat their request for service after an exponentially distributed time with rate . When the server is idle it makes an outgoing call of type n in an exponentially distributed time with rate n. There are N types of outgoing calls whose durations follow N distinct exponential distributions. We assume that the durations of outgoing calls of type n follow the exponential distribution with rate n. ❄ . . . ❄ ✚ ✚ ✚ ✲ 1 ✚ ✚ ✚ ✚✚ n ✏ ✑ ✒✻ n Figure 1. Markovian retrial queue with two way communication and multiple types of outgoing calls Problem definition Let k(t) denote the state of the server at the time t � 0, 0, if the server is idle, k(t) = 1, if an incoming call is in service, n, if an outgoing call of type n is in service, n = 2, N + 1. Let i(t) denote the number of incoming calls in the system at the time t. It is easy to see that process {k(t), i(t)} forms a continuous time Markov chain. We assume that the Markov chain is ergodic and the stationary distribution of {k(t), i(t)} exists. Let P {k(t) = k, i(t) = i} = Pk (i) denote the stationary probability distri- bution of the system state which is the unique solution of Kolmogorov system of equations: - г + i + N +1 l n n=2 P0(i) + 1P1(i + 1) + N +1 n=2 nPn(i) = 0, (1) - ( + 1)P1(i) + P1(i - 1) + P0(i - 1) + iP0(i) = 0, - ( + n)Pn(i) + Pn(i - 1) + nP0(i) = 0, n = 2, N + 1. Let Hk (u) denote the partial characteristic functions Hk (u) = i=0 ejuiPk (i), k = 0, N + 1, where j = -1. Multiplying equations of system (1) by ejui and taking the sum over i yields - г N +1 l + n n=2 t - H0(u) + jH0(u) + 1e ju H1(u) + N +1 n=2 nHn(u) = 0, (2) ju - ( + 1)H1(u) + e H1(u) + eju 0 H0(u) - jHt (u) = 0, - ( + n)Hn(u) + ejuHn(u) + nH0(u) = 0, n = 2, N + 1. The characteristic function H(u) of the number of incoming calls in the retrial queue is expressed through partial characteristic functions Hk (u) by N +1 H(u) = Hk (u). The main content of this paper is the solution of system k=0 (2) by using an asymptotic analysis methods in two limit conditions: of the high rate of making outgoing calls and the low rate of service time of outgoing calls. Prelimit analysis In this section, we obtain expressions for the stationary distribution using the characteristic functions. First, we derive explicit expression for the characteristic function H(u) of the number of incoming calls in the system. Theorem 1. Explicit expression for the characteristic function H(u) of the number of incoming calls in M |M |1|N retrial queue is given as follows: H(u) = 1 ( N +1 1 + n \ 1 + 1 n=2 n + (1 - eju) n(n-) Г 1 - тт n l (1+2)+1 N +1 Г 1 - p l n , where 1 - eju N +1 n=2 N +1 1 - pneju k = , 1 = , 2 1 k k=2 = k , k k=2 pn = n + , n = + n - 1, n = 2, N + 1. Proof. From equations 2 and 3 of the system (2) we obtain expressions for partial characteristic functions: eju j t (3) H1(u) = 1 0 - ju 0 H (u) + (1 - eju) 1 n H (u), + (1 - e ) (4) Hn(u) = n H0(u), n = 2, N + 1. + (1 - eju) Substituting this equations into the first equation of the system (2), we find that г eju Ht N +1 1 + (1 - eju) neju l H (u). (5) 0(u) = j + 1 - eju 1 - eju n=2 0 n + (1 - eju) The solution of this differential equation is given by Г 1 - l (1+ ) N +1 Г 1 - p 2 тт n l n(n-) (n) (6) H0(u) = r0 1 - eju n=2 1 - pn eju , 0 where = , r = H0(0) = P {k(t) = 0}, 2 N +1 = n k , p = , 1 n = + n - 1, n = 2, N + 1. Substituting u = 0 into the system (2) yields: k=2 k n + - ( N +1 \ + n n=2 t r0 + j H0(u)|u=0 + N +1 k=1 krk = 0, (7) 0 - 1r1 + r0 - j Ht (u)| u= 0 = 0, - nrn + nr0 = 0, n = 2, N + 1, 0 where expression for Ht (u)| u=0 can be obtained substituting u = 0 into (5). It follows from equations 2 and 3 of the system (7) that г ( N +1 \l n r1 = + 1 1 - + 1 n n=2 r0, n r = n n r0, n = 2, N + 1. Furthermore, from the normalization condition: N +1 rk = 1, we obtain k=0 1 - n(1 - ) r0 = 1 , r1 = (1 + 1) 1 1 n , rn = , n = 2, N + 1, (1 + 1) where 1 = we obtain k N +1 . Substituting (6) into (3) and (4) and summing up results, k=2 k H(u) = 1 ( N +1 1 + n \ 1 + 1 n=2 n + (1 - eju) n(n-) Г 1 - тт n l (1+2)+1 N +1 Г 1 - p l n . 1 - eju n=2 1 - pneju Asymptotic analysis of the model under the high rate of making outgoing calls In this section, we will investigate system (2) by asymptotic analysis method under the high rate of making outgoing calls condition. In particular, we prove that asymptotic characteristic function of the number of incoming calls in the system corresponds to Gaussian distribution. Denoting n = n, we obtain - г + N +1 l n n=2 t - H0(u)+jH0(u) + 1e ju H1(u)+ N +1 n=2 nHn(u) = 0, (8) ju - ( + 1)H1(u) + e H1(u) + eju 0 H0(u) - jHt (u) = 0, - ( + n)Hn(u) + ejuHn(u) + nHn(u) = 0, n = 2, N + 1. First order asymptotic Theorem 2. Suppose i(t) is the number of incoming calls in the system of the stationary M |M |1|N retrial queue with outgoing calls, then the (9) holds where i(t) lim Eejw = ejw1 , (9) 11 N +1 n (10) 1 1 = ( , 1 = - ) . n n=2 Proof. We denote = 1/ in the system (8), and introduce the following notations u = w, H0(u) = F0(w, ), Hk (u) = Fk (w, ), k = 1, N + 1, in order to get the following system - ( + N +1 n=2 n)F0(w, ) + j F0(w, ) + e- w 1 N +1 jw F1(w, )+ + n=2 nFn(w, ) = 0, - ( + 1)F1(w, ) + e jw F1(w, ) + e jw F0(w, )- (11) jw - j F0(w, ) w = 0, - ( + n)Fn(w, ) + e Fn(w, )+ + nF0(w, ) = 0, n = 2, N + 1. Summing up equations of system (11), we obtain N +1 F0(w, ) + ( - 1e-jw)F1(w, ) + Fn(w, ) = 0. (12) n=2 Considering the limit as 0 in the system (11) and equation (12), then we will get N +1 - n=2 0 nF0(w) + jF t(w) + N +1 k=1 k Fk (w) = 0, 0 - 1F1(w) - jF t(w) = 0, - nFn(w) + nF0(w) = 0, n = 2, N + 1, (13) 1 1 - ( - )F (w) + N +1 n=2 Fn(w) = 0. We propose to get the solution of the system (13) in the form of Fk (w) = (w)rk, k = 0, N + 1. (14) Here rk , k = 1, N + 1 is the probability of the server state k; r0 has no sense of probability, since the probability that the server will be in the zero state as is zero: N +1 - n=2 nr0 + j t(w) (w) r0 + N +1 k=1 krk = 0, t(w) - 1r1 - j (w) r0 = 0, n n n 0 - r + r = 0, n = 2, N + 1, (15) - (1 - )r1 + N +1 n=2 rn = 0. As the relation j t(w) (w) does not depend on w, the function is obtained in the following form (w) = exp{jw1}, which coincides with (9). The value of the parameter 1 will be defined below. We rewrite the system (15) in the form N +1 - n=2 nr0 - 1r0 + N +1 k=1 krk = 0, - 1r1 + 1r0 = 0, - nrn + nr0 = 0, n = 2, N + 1, (16) - ( - )r + 1 1 N +1 n=2 rn = 0. The normalization condition for stationary server state probability distri- N +1 bution is rk = 1. We have k=1 - nrn + nr0 = 0, n = 2, N + 1, - (1 - )r1 + N +1 N +1 n=2 rn = 0, (17) k=1 rk = 1. The solution of the system (17) is given by 1 - n(1 - ) r0 = , r1 = 1 1 , rn = n11 , n = 2, N + 1, (18) where 1 = n N +1 . Substituting (18) into system (16), we obtain an equation n=2 n for 1, which coincides with (10). The first order asymptotic i.e. Theorem 2, only defines the mean asymptotic value 1 of a number of incoming calls in the system in prelimit situation of . For more detailed research of number i(t) of incoming calls in the system lets consider the second order asymptotic. D Second order asymptotic Theorem 3. In the context of Theorem 2 the following equation is true ( lim E exp jw i(t) - 1 (jw)2 = e , (19) 2 2 where 1(1 - )(2 + 1) + 21 N +1 n N +1 n (20) 2 = (1 , 1 = - )2 n=2 , 2 = n 2 . n=2 n Proof. We introduce the following notations in the system (8) and we get k Hk (u) = exp{ju1}H(2)(u), k = 0, N + 1, (21) ( - + N +1 n=2 \ n + 1 0 H(2)(u) + j N +1 dH(2)(u) 0 + du + 1e -ju H (2) 1 (u) + n=2 n nH(2)(u) = 0, (2) ju - ( + 1)H1 (u) + e 1 H(2)(u) + (eju 0 + 1)H(2)(u)- 0 dH(2)(u) (22) - ( + n)H(2)(u) + ejuH(2)(u)+ - j du = 0, n n 0 + nH(2)(u) = 0, n = 2, N + 1. Denoting = 1/2, and introducing the following notations u = w, H(2)(u) = 2F (2)(w, ), H(2) 0 0 (2) (23) we obtain k (u) = Fk (w, ), k = 1, N + 1, j 0 F (2)(w, ) w - ( 1 + 2 N +1 \ + n n=2 0 F (2)(w, )+ + 1e -jw F (2) 1 (w, ) + N +1 n=2 n nF (2)(w, ) = 0, (2) - ( + 1)F1 (w, ) + e jw 1 F (2)(w, )+ F (2)(w, ) (24) 0 + (ejw2 + 1)F (2)(w, ) - j 0 w = 0, - ( + n)F (2)(w, ) + ejwF (2)(w, )+ n n + nF (2) 0 (w, ) = 0, n = 2, N + 1. Summing up equations of the system (24), we obtain 2F (2) (2) N +1 n 0 (w, ) + ( - 1e-jw)F1 (w, ) + n=2 F (2)(w, ) = 0. (25) Our idea is to seek for a solution of the system (24) and equation (25) in the form F (2) 2 k (w, ) = 2(w){rk + jwfk } + o( ), k = 0, N + 1. (26) Substituting (26) to (24) and (25), laying out the exhibitors in tailor series and taking (16) into account, dividing these equations by and taking the limit as 0, we have ( N +1 \ - 1 + n n=2 N +1 f0 + kfk - 1r1 + k=1 2 t (w) w(w) r0 = 0, 1f0 - 1f1 + r1 - 2 t (w) w(w) r0 = 0, -nfn + rn + nf0 = 0, n = 2, N + 1, N +1 -(1 - )f1 + fn + 1r1 = 0. n=2 2(w) This equations imply that t w2(w) doesnt depend on w and thus the function 2(w) is given in the following form (jw)2 2(w) = exp which coincides with (19). We have t 2(w) w2(w) and then we obtain the system 2 , 2 = -2 ( - 1 + N +1 \ n n=2 f0 + N +1 k=1 kfk = 1r1 + 2r0, 1f0 - 1f1 = -r1 - 2r0, - nfn + nf0 = -rn, n = 2, N + 1, (27) 1 1 - ( - )f + N +1 n=2 fn = -1r1. Substituting values (18) into the system (27), we have f = n n f0 + (1 - )n 2 , n = 2, N + 1, n 1 n1 1 22 where f1 = 1 f0 + - 11 + , 1 - k N +1 1 = , 2 k N +1 = k . 2 k=2 k=2 k Substituting this expressions into equation 2 of the system (27), we obtain an equation for 1, which coincides with (20). Second order asymptotic i.e. Theorem 3, shows that the asymptotic prob- ability distribution of the number i(t) of incoming calls in the system is Gaussian with mean asymptotic 1 and variance 2. D Asymptotic analysis of the model under the low rate of service time of outgoing calls In this section, we will investigate system (2) by asymptotic analysis method under the low rate of service time of outgoing calls condition. Denoting n = n, we obtain - г N +1 l + n n=2 t - H0(u)+jH0(u)+1e ju H1(u)+ N +1 n=2 nHn(u) = 0, (28) ju - ( + 1)H1(u) + e H1(u) + eju 0 H0(u) - jHt (u) = 0, - ( + n)Hn(u) + ejuHn(u) + nHn(u) = 0, n = 2, N + 1. Theorem 4. Suppose i(t) is a number of incoming calls in a system of stationary M |M |1|N retrial queue with two way communication, then the following equation is true H(u) = lim Eejwi(t) = 1 N +1 n 1 N +1 тт 1 \ 1n - (1-) jw , (29) where 1 = 0 n N +1 . 1 n=2 n - jw n=2 - n n=2 n Proof. We denote = , lets substitute the following in the system (28) u = w, H0(u) = F0(w, ), Hk (u) = Fk (w, ), k = 1, N + 1. We will get the system ( - + N +1 n=2 \ n F0(w, ) + j -jw F0(w, )+ w N +1 + 1e F1(w, ) + F0(w, ) n=2 nFn(w, ) = 0, jw (30) - ( + 1)F1(w, ) - j w jw + e + ejw F1(w, )+ F0(w, ) = 0, - ( + n)Fn(w, ) + e Fn(w, )+ + nF0(w, ) = 0, n = 2, N + 1. Considering the limit as 0 in the system (30) then we will get -jF t(w) - 1F1(w) = 0, jF t(w) + 1F1(w) = 0. (31) 0 0 Summing up equations of the system (30) we have N +1 F0(w, ) + ( - 1e-jw)F1(w, ) + nFn(w, ) = 0. (32) n=2 Laying out the exhibitors in tailor series, dividing equations by and taking the limit as 0, taking (31) into account, we obtain ( N +1 \ - + n n=2 N +1 F0(w) - jw1F1(w) + nFn(w) = 0, n=2 0 -jF t(w) - 1F1(w) = 0, (jw - n)Fn(w) + nF0(w) = 0, n = 2, N + 1 N +1 -(1 - )F1(w) + Fn(w) = 0. n=2 From the last system of equations we have n (33) Fn(w) = n F0(w), - jw N +1 (34) F1(w) = 1 - n=2 Fn(w). Then N +1 n (35) F1(w) = 1 F0(w) - n=2 . n - jw Substituting (35) into (31), we obtain F t 1 N +1 n 1 0(w) = j ( F0(w) - ) n=2 . n - jw The solution of differential equation is given by N +1 1 F0(w) = C тт n=2 1 - jw \ 1n - (1-) n , (36) where C is an integration constant and its value will be obtained later. We N +1 denote asymptotic characteristic function (36) into (33) and (34), we obtain Fk (w) = (w). Substituting k=1 N +1 N +1 1 \- 1k k тт (1-) F1(w) = 1 - k=2 C k - jw k=2 1 - jw , k Fn(w) = n n C - jw 1 N +1 тт k=2 \- 1 - jw k 1k (1-) , n = 2, N + 1. Summing up equations, we have N +1 1n - N +1 (w) = 1 C тт 1 \ 1 - jw . (1-) n 1 - n n=2 n - jw n=2 Using condition (0) = 1, we obtain C = 1 - 11 , where 1 = N +1 n n . n=2 We obtain the characteristic function (29). D Approximation accuracy The accuracy of the approximation P (2)(i) is defined by using Kolmogorov i range ∆2 = max (P () - P (2)()) , which represents the difference be- =0 0�i�N tween distributions P (i) and P (2)(i), where P (i) is obtained by using inverse Fourier transform for the characteristic function of the M |M |1|N retrial queue and the approximation P (2)(i) is given by obtained asymptotics. We consider N = 3, = 0.2, 1 = 1 and = 1 for Tables 1 and 2. Kolmogorov range, 2 = 2, 3 = 3, 4 = 4, 2 = 1, 3 = 2, 4 = 3 Table 1 = 3 = 5 = 10 = 50 = 100 ∆2 0.066 0.043 0.023 0.01 0.007 Kolmogorov range, 2 = 2, 3 = 3, 4 = 4, 2 = 1, 3 = 2, 4 = 3 Table 2 = 0.05 = 0.035 = 0.02 = 0.01 ∆2 0.059 0.044 0.026 0.014 Conclusions In this paper, we have considered retrial queue with two way communication with multiple types of outgoing calls. We have found characteristic function of the number of incoming calls in the system. We have found the first and the second order asymptotics of the number of calls in the system under the condition of the high rate of making outgoing calls. Based on the obtained asymptotics we have built the Gaussian approximation of the probability distribution of the number of incoming calls in the system. We have found asymptotic characteristic function of the number of incoming calls in retrial queue under the condition of the low service rate of outgoing calls. In future we plan to consider this retrial queueing system under other asymptotic conditions.</p>[J. R. Artalejo, A. Gómez-Corral, Retrial queueing systems: a computational approach (2008). doi:10.1007/978-3-540-78725-9.][G. Falin, J. G. C. Templeton, Retrial queues, Vol. 75, CRC Press, 1997.][G. P. Basharin, K. E. Samuyjlov, On a single-phase queuing system with two types of requests and relative priority [Ob odnofaznoy sisteme massovogo obsluzhivaniya s dvumya tipami zayavok i otnositel’nym prioritetom], Journal of Academy of Sciences of the USSR. Technical Cybernetics [Izvestiya akademii nauk SSSR. Tekhnicheskaya kibernetika] 3 (1983) 48-56, in Russian][S. Bhulai, G. Koole, A queueing model for call blending in call centers, IEEE Transactions on Automatic Control 48 (8) (2003) 1434-1438][A. Deslauriers, P. L’Ecuyer, J. Pichitlamken, A. Ingolfsson, A. N. Avramidis, Markov chain models of a telephone call center with call blending, Computers & operations research 34 (6) (2007) 1616-1645][L. Brown, N. Gans, A. Mandelbaum, A. Sakov, H. Shen, S. Zeltyn, L. Zhao, Statistical analysis of a telephone call center: a queueing-science perspective, Journal of the American Statistical Association 100 (469) (2005) 36-50][S. Aguir, F. Karaesmen, O. Z. Akşin, F. Chauvet, The impact of retrials on call center performance, OR Spectrum 26 (3) (2004) 353-376. doi:10.1007/s00291-004-0165-7.][J. R. Artalejo, T. Phung-Duc, Markovian retrial queues with two way communication, Journal of Industrial and Management Optimization 8 (4) (2012) 781-806.][J. R. Artalejo, T. Phung-Duc, Single server retrial queues with two way communication, Applied Mathematical Modelling 37 (4) (2013) 1811-1822.][T. Phung-Duc, W. Rogiest, Two way communication retrial queues with balanced call blending, in: International Conference on Analytical and Stochastic Modeling Techniques and Applications, Springer, 2012, pp. 16-31. doi:10.1007/978-3-642-30782-9_2][H. Sakurai, T. Phung-Duc, Scaling limits for single server retrial queues with two-way communication, Annals of Operations Research 247 (1) (2016) 229-256. doi:10.1007/s10479-015-1874-9.][H. Sakurai, T. Phung-Duc, Two-way communication retrial queues with multiple types of outgoing calls, Top 23 (2) (2015) 466-492. doi:10.1007/s11750-014-0349-5.][T. Phung-Duc, K. Kawanishi, An efficient method for performance analysis of blended call centers with redial, Asia-Pacific Journal of Operational Research 31 (02) (2014) 1440008. doi:10.1142/s0217595914400089.][A. Kuki, J. Sztrick, Á. Tóth, T. Bérczes, A contribution to modeling two-way communication with retrial queueing systems, in: Information Technologies and Mathematical Modelling. Queueing Theory and Applications, Springer, Cham, 2018, pp. 236-247. doi:10.1007/978-3-319- 97595-5_19][V. Dragieva, T. Phung-Duc, Two-way communication M/M/1 retrial queue with server-orbit interaction, in: Proceedings of the 11th International Conference on Queueing Theory and Network Applications, ACM, 2016, p. 11][V. Dragieva, T. Phung-Duc, Two-way communication M/M/1/1 queue with server-orbit interaction and feedback of outgoing retrial calls, in: International Conference on Information Technologies and Mathematical Modelling, Springer, 2017, pp. 243-255. doi:10.1007/978-3-319-68069-9_- 20][S. Ouazine, K. Abbas, A functional approxymation for retrial queues with two way communication, Annals of Operations Research 247 (1) (2016) 211-227. doi:10.1007/s10479-015-2083-2.][M. S. Kumar, A. Dadlani, K. Kim, Performance analysis of an unreliable M/G/1 retrial queue with two-way communication, Springer, 2018, pp. 1-14. doi:10.1007/s12351-018-0430-1.][E. Morozov, T. Phung-Duc, Regenerative analysis of two-way communication orbit-queue with general service time, in: International Conference on Queueing Theory and Network Applications, Springer, 2018, pp. 22-32. doi:10.1007/978-3-319-93736-6_2.][A. A. Nazarov, S. Paul, I. Gudkova, Asymptotic analysis of Markovian retrial queue with two-way communication under low rate of retrials condition, in: Proceedings 31st European Conference on Modelling and Simulation, 2017, pp. 678-693][A. Nazarov, T. Phung-Duc, S. Paul, Heavy outgoing call asymptotics for MMPP/M/1/1 retrial queue with two-way communication, in: International Conference on Information Technologies and Mathematical Modelling, Springer, 2017, pp. 28-41. doi:10.1007/978-3-319-68069-9_3]