CC BY
23
Research article
UDC 519.872:519.217
PACS 07.05.Tp, 02.60.Pn, 02.70.Bf
DOI: 10.22363/2658-4670-2020-28-4-319-326
Stochastic analysis of a single server unreliable queue with balking and general retrial time
Mohamed Boualem
Research Unit LaMOS (Modeling and Optimization of Systems) Faculty of Technology, University of Bejaia, Algeria
(received: June 8, 2020; accepted: November 12, 2020)
In this investigation, we consider an M/G/l queue with general retrial times allowing balking and server subject to breakdowns and repairs. In addition, the customer whose service is interrupted can stay at the server waiting for repair or leave and return while the server is being repaired. The server is not allowed to begin service on other customers until the current customer has completed service, even if current customer is temporarily absent. This model has a potential application in various fields, such as in the cognitive radio network and the manufacturing systems, etc. The methodology is strongly based on the general theory of stochastic orders. Particularly, we derive insensitive bounds for the stationary distribution of the embedded Markov chain of the considered system.
Key words and phrases: Retrial queue, Markov chain, balking, breakdowns and repairs, stochastic orders, bounds, ageing classes
1. Introduction
The study on queueing models have become an indispensable area due to its wide applicability in real life situations. Retrial queues occupy an intermediate situation between an Erlang model with loss and classical model with wait, which constitute their limiting models in the case of low and high retrial rates. Retrial queueing systems are characterized by the requirement that customers finding the service area busy must join the retrial group and retry for service at random intervals. Queues in which customers are allowed to retry have been extensively used to model many problems in telephone switching systems, cognitive radio network, manufacturing systems, telecommunication networks and computer systems for competing to gain service from a central processor unit [1]—[3].
Retrial queueing systems with general service times and non-exponential retrial time distributions have been received little attention because of the complexity of the known results. Indeed, in most cases, we are faced with systems of equations whose resolution is complex, or having solutions not
© Boualem M., 2020
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
easily interpretable. For instance, Pollaczek-Khintchine formula requires a numerical inversion of the Laplace transform to compute the distribution of the waiting time. In many cases, even the Laplace transform or probability generating functions are not available in explicit forms. To overcome these difficulties, approximation methods are often used to obtain quantitative and/or qualitative estimates for certain performance measures. For all these reasons, in this study, a particular interest is devoted to the stochastic comparison method based on the general theory of stochastic orders [4], [5].
The stochastic comparison method is a mathematical tool used to study the performance of some systems modeled by continuous or discrete time Markov chains. The general idea of this method is to bound a complex system with a new system that is simpler to solve providing qualitative bounds for these performance measures. These methods represent one of the main research activities in various scientific fields, such as economy, biology, operation research, reliability theory, decision theory, retrial queues and queueing networks [4]-[15].
In the present study, stochastic comparison analysis of an M/G/l queue with server subject to breakdowns, general service times and non-exponential retrial time distributions by considering the both balking and reneging behavior of the customer is presented.
The rest of the paper is organized as follows. In the next Section, we describe the considered queueing system. In Section 3, we introduce an overview on stochastic orders. In Section 4, we present some lemmas that will be used in what follows. Section 5 focusses on monotonicity properties of the transition operator. Insensitive bounds for the stationary number of customers in the system at a departure epochs are discussed in Section 6.
2. Mathematical model description
Primary customers arrive in a Poisson process with rate A. If an arriving primary customer finds the server idle, it begins service immediately and leaves the system after service completion. If the server is found to be blocked, the arriving primary customer either enters a retrial queue (according to a FCFS discipline) with probability p or leaves the system with probability 1 — p. The service times follow a general distribution B(x) with corresponding Laplace-Stieltjes transform B*(s) and the first two moments ¡3X and .
The customers in the orbit try to require the service later and the inter retrial times have a general distribution A(x) with corresponding Laplace-Stieltjes transform A* (s). The retrial customer cancels its attempt for service if a primary customer arrives first and either returns to its position in the retrial queue with probability q or quits the system with probability 1 — q.
We assume that the server may fail (or break down), but only when a customer is in service. The time until failure is exponentially distributed with mean 1/fi, but failure of the server can occur only when a customer is being served. When the server fails, the repair starts immediately and the customer just under service waits for the server until repair completion in order to accomplish its remaining service. The repair time is a random variable with probability distribution function C(x), Laplace-Stieltjes transform C*(s) and the first two moments and 72.
Upon failure of the server, the customer in service either remains in the service position with probability a until the server is up or enters a service retrial orbit with probability 1—a and keeps returning at times exponentially distributed with mean 1/0, until the server is repaired. If the customer in service enters the service retrial orbit upon server failure, after repair the server must wait for the customer to return. We refer to this time as the reserved time.
The performance characteristics of our queueing system are available in Wu etal. [16]. Under the stability condition pX^1 (1+^( +li)) < 1—q+qA* (A), the one-step transition probabilities are given as:
Tj, if i = 0, j > 0, = { (1 — q + qA*(X))r0, if z>0, j = x-1, 13 | q(1 - Wy— + (1-q + qA* (X))r_+i, if i>0, 3>i — 1, () ,0, otherwise,
where
+oo
f (pXx)m
rm = I -:— exp {—pXx}dB• (x), m = 0,1,2,...
J ml
0
is the distribution of the number of primary customers which arrive during the generalized service time of the (n + 1)th customer. We define the generalised service time as the length of time from when a customer begins service until service completion (that includes: some eventual repair times (Xi) and reserved times (Yi)) with common distribution
^(?) = i it i a^(1 — aY_^e_^y(x-y)dB(y).
*=0 j=0 J W
3. Overview on stochastic orders and ageing notions
Stochastic ordering is useful for studying internal changes of performance due to parameter variations, to compare distinct systems, to approximate a system by a simpler one, and to obtain upper and lower bounds for the main performance measures of systems. For a comprehensive discussion on these stochastic orders, see [4], [5].
Let X and Y be two random variables non-negative with distribution functions F and G, respectively.
Laplace order (<L): X Y^ J e_sxdF(x)^ J e_sxdG(x),Vs > 0.
00
Stochastic order (<si): X Y ^ F(x) > G(x),Vx e R.
CO CO
Convex order ): X Y ^ I(1 — F(t))dt < I(1 — G(t))dt, Vx > 0.
Proposition 1. Let X be a random variable with distribution function F and finite mean m.
(a) F is New Better than Used in Expectation (NBUE) iff F F *,
(b) F is New Worse than Used in Expectation (NWUE) iff F * F,
(c) F is of class L iff F F *, where F * is the exponential distribution function with the same mean as F.
4. Comparison bounds for the probability rm
This subsection presents several useful lemmas which will be used later in establishing the main results. For this, we consider two M/G/l queue with general retrial times allowing balking and server subject to breakdowns and repairs with parameters A^, q(z), a(z), C(%), B(z), for i = 1,2 respectively.
Lemma 1. If X(1) < A(2), q(1) < q(2), a(1) > a(2), ¿r> < ¿2), C(1) C(2), B(1) B(2), then {rti}.
Lemma 2. If X(1) < A(2), q(1) < q(2), a(1) > a(2), /i(1) < /i(2), C(1) C(2), B(1) B(2), then }.
Proof. To show both lemmas, we have to establish the usual numerical inequalities
rtt = Y (for <s=^st),
n'm
^m = Y^) (for <s = ^v).
n'm
By definition, we have
= Y f = fY e^*dB*<i) (X), i = 1,2.
fm 0 fm
rm = Yr)) = I Y Y —1\—e p A 1 dB ^ l = 1,2-
j'm 0 j'm '
The rest of the proof is known in the more general setting of a random summation. n
5. Monotonicity properties of the transition operator
Let T be the transition operator which associates to every distribution u = {uj}' a distribution Tw = {u^}' such that = ^ .
i'O
Theorem 1. The transition operator T is monotone with respect to the orders stochastic t) and convex v).
Proof. It is well known that an operator is monotone with respect to <si if and only if Pij — Pi_1j > 0, and is monotone with respect to if and only
if Ti_lj + Tz+lj — 2% >0, Vi, J.
__CO — C
In our case, PZJ = £ pu and PZJ = ^pu. □
Now, we consider two M/G/l queue with general retrial times allowing balking and server subject to breakdowns and repairs with parameters X(z), q(z), a(%), , C, A(z), B(%), for i = 1,2 respectively. In the following theorem, we give comparability conditions of two transition operators T1 and T2.
Theorem 2. If X(1) < A(2), q(1) < q(2), a(1) > a(2), /i(1) < /i(2), C(1) <s C(2), B(1) <s B(2) and 1) <L 2), then T1 <s T2, where <s is one of the symbols <si or <v. I.e. for any distribution u, we have T1 u <s T2u.
Proof. It is well known that to prove T1 <s T2, we have to show the following numerical inequalities for the one-step transition probabilities P-1,
Pff (see Stoyan [5]): < , Vi, j, (for <s=<si), and P^ < ^, Vi,j, (for <s).
After algebraic manipulations, we obtain the result searched. □
6. Bounds for the stationary distribution
Now, we add the the corresponding stationary distributions of the number of customers in the system nn 2 and nn2.
Theorem 3. If X(1) < A(2), q(1) < q(2), a(1) > a(2), ¿r> < ¿2), C(1) <s
C(2), B(1) <s B(2) and AW <L A(2), then {^} <s n2)}, where <s is one of the symbols <si or <v.
Proof. Using Theorems 1 and 2 which state that T are monotone with respect to the order <s and T1 <s T2, we have by induction T1,nw <s T2,nu for any distribution u. Taking the limit, we obtain the stated result. □
Theorem 4. We consider an M/G/l queue with general retrial times allowing balking and server subject to breakdowns and repairs. If we have, the service time distribution B(x) and the repair time distribution C(x) are NBUE and the retrial time distribution is of class L, then {Kn} {k^ }, where {K*n} is the stationary distribution in the Markovian retrial queue with the same parameters.
Proof. Consider an auxiliary Markovian M/M/l retrial queue allowing balking and server subject to breakdowns and repairs, having the same parameters as the model considered in our study , but with exponentially distributed retrial time A*(x), service time B* (x) and repair time C* (x). We have also, if B(x) and C(x) are NBUE and A(x) is L, then B(x) <v B*(x), C(x) C* (x) and A(x) <L A*(x). Therefore, by using Theorem 3, we deduce the statement of this theorem. □
7. Conclusions
In this work, we investigated the monotonicity of the transition operator of the embedded Markov chain with respect to stochastic and convex orderings. In addition, we done comparability conditions of the transition operators of the considered systems, associated with two Markov chains. Further, we obtained comparability conditions for which the stationary distribution in the considered model is bounded by the stationary distribution of the Markovian retrial queue, if the service time and the repair time distributions are NBUE and the retrial time distribution is of class L.
Acknowledgments
The authors are pleased to thank the anonymous referees and the editor for their valuable comments and suggestions, which improved the content and the presentation of this paper.
References
[1] J. R. Artalejo and A. Gómez-Corral, Retrial queueing system: A computational approach. Berlin: Springer, Berlin, Heidelberg, 2008, 318 pp.
[2] A. A. Nazarov, S. V. Paul, and O. D. Lizyura, "Two-way communication retrial queue with unreliable server and multiple types of outgoing calls," Discrete and Continuous Models and Applied Computational Science, vol. 28, no. 1, pp. 49-61, 2020. DOI: 10.22363/2658-4670-2020-28-149-61.
[3] D. Zirem, M. Boualem, K. Adel-Aissanou, and D. Ai'ssani, "Analysis of a single server batch arrival unreliable queue with balking and general retrial time," Quality Technology & Quantitative Management, vol. 16, pp. 672-695, 2019. DOI: 10.1080/16843703.2018.1510359.
[4] M. Shaked and J. G. Shanthikumar, Stochastic Orders. New York: Springer-Verlag, 2007, 473 pp.
[5] D. Stoyan, Comparison methods for queues and other stochastic models. New York: Wiley, 1983, 217 pp.
[6] L. M. Alem, M. Boualem, and D. Ai'ssani, "Bounds of the stationary distribution in M/G/l retrial queue with two-way communication and n types of outgoing calls," Yugoslav Journal of Operations Research, vol. 29, no. 3, pp. 375-391, 2019. DOI: 10.2298/YJ0R180715012A.
[7] L. M. Alem, M. Boualem, and D. Ai'ssani, "Stochastic comparison bounds for an Mx ,M2/Gx ,G2/1 retrial queue with two way communication," Hacettepe Journal of Mathematics and Statistics, vol. 48, no. 4, pp. 11851200, 2019. DOI: 10.1572/HJMS.2018.629.
[8] M. Boualem, "Insensitive bounds for the stationary distribution of a single server retrial queue with server subject to active breakdowns," Advances in Operations Research, vol. 2014, no. 1, pp. 1-12, 2014. DOI: 10.1155/2014/985453.
[9] M. Boualem, A. Bareche, and M. Cherfaoui, "Approximate controllability of stochastic bounds of stationary distribution of an M/G/l queue with repeated attempts and two-phase service," International Journal of Management Science and Engineering Management, vol. 14, no. 2, pp. 79-85, 2019. DOI: 10.1080/17509653.2018.1488634.
[10] M. Boualem, M. Cherfaoui, and D. Aïssani, "Monotonicity properties for a single server queue with classical retrial policy and service interruptions," Proceedings of the Jangjeon Mathematical Society, vol. 19, no. 2, pp. 225-236, 2016.
[11] M. Boualem, M. Cherfaoui, N. Djellab, and D. Aïssani, "A stochastic version analysis of an M/G/l retrial queue with Bernoulli schedule," Bulletin of the Iranian Mathematical Society, vol. 43, no. 5, pp. 13771397, 2017.
[12] M. Boualem, M. Cherfaoui, N. Djellab, and D. Aïssani, "Inégalités stochastiques pour le modèle d'attente M/G/l/l avec rappels," French, Afrika Matematika, vol. 28, pp. 851-868, 2017. DOI: 10.1007/s13370-017-0492-x.
[13] M. Boualem, N. Djellab, and D. Aïssani, "Stochastic inequalities for M/G/l retrial queues with vacations and constant retrial policy," Mathematical and Computer Modelling, vol. 50, no. 1-2, pp. 207-212, 2009. DOI: 10.1016/j.mcm.2009.03.009.
[14] M. Boualem, N. Djellab, and D. Aïssani, "Stochastic approximations and monotonicity of a single server feedback retrial queue," Mathematical Problems in Engineering, vol. 2012, 12 pages, 2012. DOI: 10.1155/2012/ 536982.
[15] M. Boualem, N. Djellab, and D. Aïssani, "Stochastic bounds for a single server queue with general retrial times," Bulletin of the Iranian Mathematical Society, vol. 40, no. 1, pp. 183-198, 2014.
[16] X. Wu, P. Brill, M. Hlynka, and J. Wang, "An M/G/l retrial queue with balking and retrials during service," International Journal of Operational Research, vol. 1, pp. 30-51, 2005. DOI: 10.1504/IJOR.2005.007432.
For citation:
M. Boualem, Stochastic analysis of a single server unreliable queue with balking and general retrial time, Discrete and Continuous Models and Applied Computational Science 28 (4) (2020) 319-326. DOI: 10.22363/2658-46702020-28-4-319-326.
Information about the authors:
Mohamed Boualem — Full Professor, Professor of Applied Mathematics at the Department of Technology at the University of Bejaia, Algeria. (e-mail: [email protected], phone: +213792991121,
ORCID: https://orcid.org/https://orcid.org/0000-0001-9414-714X, ResearcherID: A-3176-2019, Scopus Author ID: 26633399700)
УДК 519.872:519.217
PACS 07.05.Tp, 02.60.Pn, 02.70.Bf
DOI: 10.22363/2658-4670-2020-28-4-319-326
Стохастический анализ системы типа «клиент—сервер» с ненадёжной очередью с блокировкой и общим временем обновления
Мохамед Буалем
Исследовательское подразделение LaMOS Факультет технологий, Университет Беджая, Алжир
В статье рассматривается система массового обслуживания типа M/G/1 с обобщённым временем обновления, допускающая блокировку, выход из строя и возобновление работы сервера. Кроме того, клиент, обслуживание которого прервано, может оставаться на сервере в ожидании восстановления его работы, а может покинуть систему и вернуться в период восстановления работы сервера. Серверу не разрешается начинать обслуживание других клиентов до тех пор, пока текущий клиент не завершит обслуживание, даже если он временно отсутствует. Эта модель имеет потенциальное применение в различных областях, таких как сеть когнитивного радио, производственные системы и т. д. Методология строго базируется на общей теории стохастических порядков. В частности, получены оценки стационарного распределения вложенной цепи Маркова рассматриваемой системы.
Ключевые слова: очередь с обновлением, цепь Маркова, блокировка, выход из строя и восстановление, стохастический порядок, границы, классы старения
