ANALYSIS OF NON-MARKOVIAN BATCH ARRIVAL RETRIAL QUEUE WITH PRIORITY SERVICES, IMMEDIATE FEEDBACK, PUSH OUT, DIFFERENTIATED BREAKDOWNS, DELAYED REPAIR, RANDOMIZED VACATION
1G. Ayyappan, 2S. Nithya •
Department of Mathematics Puducherry Technological University Puducherry-605014, India [email protected], [email protected]
Abstract
Priority and ordinary customers arrive according to Poisson processes, and their service time based on the general distribution. The server constantly offers a single service for both priority and ordinary customers. We compute the Laplace transforms of the time-dependent probabilities of system states using the probability generating function and supplementary variable technique. Numerical results are obtained which are also examined to facilitate the sensitivity analysis of system descriptions.
Keywords: Batch Arrivals; Priority Queues; Immediate Feedback; Push out; Differentiated Breakdowns; Delayed Repair; Randomized Vacation. AMS Subject Classification (2010): 60K25, 68M30, 90B22.
1. Introduction
We see several queueing situations every day where customers must wait for service and there is a delay in providing it. Retrial queues in queueing theory have been the topic of a lot of exciting research over the last two decades. The concept of retrial queues has attracted the attention of numerous scholars and received important contributions from them. An M/M/1 retrial queueing system with Poisson arrival flows, impatient customers, breakdown, collisions was studied by Danilyuk et al. [9]. Nazarov et al. [16] investigated how, depending on whether the server is busy or idle, it is dependent on random failures and repairs in a retry queuing system with a finite number of sources and customer collision. For the aggregation of the customers and their group service, D'Arienzo et al. [10] created a single-server retrial queue with a MAP flow, PH service times, and a finite capacity. Ahuja et al. [1] explored the retrial queueing system with an optional service and finite population subject to balking. Pavai Madheswari et al. [17] analysed at an M/G/1 retrial queueing system with two service phases, the second of which is optional, and a server working on a Bernoulli vacation schedule. Innovative applications for performance study of various systems in telecommunications, data split networks, traffic management on high-speed networks, and production engineering make use of these queueing models.
The literature on retrial queueing has extensively researched retrial queues with various customer categories. An important component of priority discipline is preemptive and non-preemptive priorities. D'Apice et al. [11] considered a priority queueing model with many types of requests and restricted processor sharing. Ayyappan and Thilgavathy [3] determined priority queueing system with breakdown, repair, discouragement, single vacation, standby server,
negative arrival and impatient customers. Li et al. [15] investigate equilibrium queueing strategies in an unobservable non-preemptive priority queue with homogeneous customers. Ammar and Rajadurai [2] introduced preemptive priority retrial queueing system with disaster under working breakdown services.
Most of the time, while discussing queueing, it is generally thought that the server is always accessible. However, server failure made the significant impact in queueing system. Therefore, in order to establish retrial queueing models, it is essential to carry out investigation on the retrial queue with breakdowns. Choudhury and Kalita [8] described a non-Markovian queueing system with breakdown, delayed repair and two general heterogeneous service, optional service. Krishna Kumar et al. [14] examined a Markovian retrial queue where the server is subject to breakdowns and repairs. Gao et al. [12] studied an M/G/1 retrial queue with two types of breakdowns. Ayyappan and Gowthami [5] researched the single server classical queueing system MAP/PH/1 with breakdown, repair, Bernoulli vacation and setup time. Begum and Choudhry[7] investigated a M/(G1, G2)/1 queue with service interruption consisting of a definite repairability.
Customers may be serviced more than once for particular reasons in several queueing situations. Customers have to re-join the queue and wait in queue after the service is completed. Optional re-service is a concept that can be considered as immediate feedback in this regard. The customer receives their service in the first step, and if there is any issue with it or they need it again, they will receive it immediately without having to wait in queue. Re-service has several practical applications in places like bank desks, functioning ATMs, large supermarkets, and medical facilities etc., The idea of immediate feedback (re-service) has been addressed by some authors, including Azhagappan and Deepa [6], Ayyappan and Deepa [4] and Jose and Deepthi [13]. According to the previously mentioned literature, a customer who wants to receive more service must visit the server once more at that moment.
The interesting parameter in this chapter is the randomized vacation policy. It is described as follows: After the vacation completion, if there is at least one unit present in the system, then the server immediately commence the service. Otherwise, the server will decide either remains idle or go for another vacation, if no units present in the system. The concept of variant vacation policy was proposed by Takagi [18], which is a generalization of the single and multiple vacation for the M/G/1 queueing system. Ke et al. [19] studied an MX/G/1 queueing system with a randomized vacation policy and at most J vacations. Geo and Yao [20] developed this vacation policy for an MX/G/1 queueing system, in which the server takes randomized vacation policy and at most J working vacations.
There are many papers dealing with unit's abandoned behaviour. Recently, Gao et al. [21] studied an M/G/1 retrial queue with abandoned customers and multi-optional vacations. Kr-ishnamoorthy et al. [22] presented an M1, M2/PH/1 retrial queueing system with pre-emptive priority service, orbital search and abandoned units in which the retrial is failed, then the failed units abandoned the system with certain probabilities. In this model, the arriving ordinary unit may remove the ordinary unit, who is getting service from the system. Here, the interrupting ordinary unit is referred as the abandoned unit.
We consider non-Markovian batch arrival retrial queue with priority services, discouragement, re-service, differentiate breakdown, restoration and delayed vacation. Priority customers and ordinary customers arrive according to Poisson processes, and their service time based on the general distribution. The server consistently provides a single service for both priority customers and ordinary customers. In the event of the server being unavailable, ordinary customers may choose to balk the system. Server failure may happen at any time during normal engaged period. The two types of system failure are hard and soft failures. Hard failure can be characterised as an equipment breakdown which demands the availability of a skilled repair person, which is an extensive process. Soft failure is described as breakdown based on by circumstances as instead of
mechanical components, and it can be generally resolved by restarting the system. Customer may re-enter the system as a feedback customer for receiving normal service due to inadequate quality of service after every priority service is completed. The server goes on vacation after priority services are completed; the time it takes the server to go on vacation is known as delay time.
The article's remaining content is formatted as follows: In Section 2, the mathematical model is presented. and the distribution of queue sizes is analysed in Section 3. Section 4 contains the exact expression for the governing equation. Section 5 of this article discusses steady state analysis. Section 6 lists stability condition. Section 7 provides an illustration of how system performance measures have an impact. Section 8 exhibits particular cases. In Sections 9 and 10, conclusions are drawn after deriving numerical and graphical results.
1.1. Integrating the Model into Real-life Situations
In the online food delivery network, independent contractors are responsible for delivering food orders to customers in their designated service areas. These contractors are self-employed individuals who have chosen to work with the food delivery platform, offering their transportation services to ensure that customers receive their meals promptly and efficiently. The nature of the food delivery business can sometimes lead to fluctuations in the availability of delivery orders within a specific area. This could be due to various factors, such as changes in customer demand, local events, or even the time of day. When independent contractors find that there are no orders available in their assigned service area, they may face a few challenges and considerations. To maintain their income and professional engagement during periods of low order availability, these contractors might choose to take a break from their work or explore other temporary job opportunities. By doing so, they can keep themselves occupied and ensure a steady income while waiting for orders to become available in their area again. This approach allows them to manage their workload and personal commitments more effectively.
Once the independent contractors have completed their temporary work or vacation, they can return to the online food delivery network and resume their services when orders are available in their area. This flexibility is crucial for contractors, as it enables them to balance their professional and personal lives while staying connected to the platform and being ready to serve customers when needed. In conclusion, the online food delivery networks independent contractors operate within separate service areas, and when delivery orders are scarce in their region, they may opt for vacations or other work to maintain their income and professional engagement. Upon the return of available orders, these contractors can resume their food delivery services, ensuring a balance between their work and personal lives. This flexibility helps them adapt to the dynamic nature of the food delivery business and maintain a sustainable work arrangement.
2. Mathematical Description
• Arrival Process :
Two distinct customer arrive in batches through separate Poisson compound processes. Ap, Ao > 0 are used to indicate, for PC and OC, the respective arrival rates. Assume the initial order probability for both priority and ordinary customers Apcidt (i = 1,2,3,...) and Aocjdt (j = 1,2,3, ...Respectively. The system has iand j batch size customers enters within a short period of time (t, t + dt). Here, 0 < c, cj < 1, Ei/j=1 ci/j = 1.
• Retrial Service Process :
Customers who are on retrial are treated the same as ordinary customers. Customers on a retrial are regular consumers. These customers will eventually return to orbit and seek their services again if the server is engaged or unavailable. Retrial service time is characterised by
Figure 1: Schematic representation
a rate of ¡3(u), which is defined by the probability density function m(u) and the probability distribution function M(u). This rate follows a general distribution.
• Regular Service Process :
Customers with different queues are served in batches, including priority and ordinary customers. The server offers a single service at rates of ^(u), i = 1,2 and service rates are distributed generally, characterized by probability distribution function B{(s) and probability density function bi(s), where i = 1,2. If the priority queue is empty, ordinary customers can begin receiving service.
• Immediate Feedback:
Only customers with a priority are offered the feedback service. Once every customer in the priority queue have received their services, if customers are dissatisfied with the service they received. Customers could either abandon the system with probability (1 — r) or they could get a re-service with probability 'r' without joining in queue.
• Push Out:
When the server is attending to an ordinary customers, a newly arriving ordinary customers has the potential to disrupt the ongoing service. It can either immediately take over the service area with probability 'q' or joins the orbit with probability 'q' (= 1 — q).
• Differentiate Breakdown:
The service channel is susceptible to failure at any moment, even when the server is operating at its usual engaged pace during any phase of service. Consequently, the server will be inaccessible for a brief duration. Both hard and soft failure rates follow exponential distributions, with rates denoted by a\ & a2 respectively.
• Delayed repair and Repair :
The server is subject to hard failure. The breakdown server does not send for repair instantly. There exists a delay period before the repair process initiates. Following this delay, the repair process commences to restore functionality. The probability distribution function D(s) and R(2) (s), along with the probability density function d(s) and r(2) (s) are employed
to characterize the delay time and the duration of hard failure repairs, respectively. Let £(u) and n2(u) be the completion rate for delay repair time and hard failure repair time. Soft failure repair time distributed exponentially with rate n1.
• Randomized Vacation: When the server determines that the system is empty upon completing a service, it opts for a vacation. Upon concluding the vacation, if the server observes an empty system again it decides whether to embark on another vacation with prob. 'p' or to remain idle with a probability of '(1-p)'. The duration of the vacation adheres to a general distribution characterized by a rate of y(u). The time spent on vacation follows a probability distribution function V(s), accompanied by its corresponding probability density function v(s).
3. Analysis of queue size distribution
The formation of governing equations is the main focus of this section. This model has been solved using the probability generating function and supplementary variable technique with respect to the non-Markovian queueing system.
Assuming that M(0) = 0, M(to) = 1, B{(0) = 0, B{(to) = 1, V(0) = 0, V(to) = 1, D(0) = 0, D(to) = 1, and R(2)(0) = 0, R(2)(to) = 1 are continuous at u = 0 for i = 1,2. So that the function ß(u), (u), *2(u), j(u), £(u) and y2(u) are the conditional completion rates for retrial, priority and ordinary customers service rate, vacation and delay time to hard failure repair, hard failure repair respectively.
Also ß(u) = 1-Mg), *(u) = 1-gö, Y(u) = £(u) = 1-D& and
n2(u) = ^RSg); i = 1,2. are the hazard rate functions of M(.), Bi(.) for i = 1,2, V(.), D(.) and R(2)(.) respectively.
Markov process for the given model is {Np(t), No(t), Y(t), M0(t), B0(t), B°(t), (V)0(t), (D)0(t), (R(2))0(t)}, where Np(t) and No(t) denote the number of customers in the priority queue and ordinary queue respectively. M0(t), B0(t), B°(t), (V)0(t), (D)0(t), (R(2))0(t)} are the elapsed retrial, service, vacation, delay time to hard failure repair and hard failure repair time of the server at time 't'.
Y(t) represents the server state. Here Y(t) = (0,1,2,3,4,5,6), denotes: 0, the server is idle; 1, retrial state; 2, engaged with PC; 3, engaged with OC; 4, on vacation ; 5, delayed to hard failure repair, and 6, hard failure repair.
Let's represent the probability as i0,„2 (t), indicating the probability that at time t, i0,„2 (t) equals the event that Np (t) = 0, No (t) = 0, and Y(t) = 0, where t > 0. We consider probability densities for this scenario.
I0,n2(u, t)du = Pr{Np(t) = 0,No(t) = n2, Y(t) = 1; u < I0(t) < u + du},n2 > 1 P^L(u, t)du = Pr{Np(t) = n1,No(t) = n2, Y(t) = 2; u < B0(t) < u + du}, P^nz(u, t)du = Pr{Np(t) = n1,No(t) = n2, Y(t) = 3; u < B^(t) < u + du}, Vn1,n2(u, t)du = Pr{Np(t) = n1,No(t) = n2, Y(t) = 4; u < V0(t) < u + du}, Dn1,n2(u, t)du = Pr{Np(t) = n1,No(t) = n2, Y(t) = 5; u < D0(t) < u + du},
42k (u, t)du = Pr {Np (t) = n1, No (t) = n2, Y(t) = 6; u < R0 (t) < u + du},
for u > 0, t > 0, n1 > 0 andn2 > 0.
G. Ayyappan, S. Nithya
ANALYSIS OF NON-MARKOVIAN BATCH ARRIVAL RETRIAL RT&A, No 3 (79)
QUEUE WITH PRIORITY SERVICES IMMEDIATE FEEDBACK, PUSH OUT..........Volume 19, September 2024
4. Equation Governing the System
d iM
—Io,o(t) = -(\p + \o)I0,0(t) + (1 - p) J Vo,o(u, t)j(u)du, (1)
d d
— Io,n(u, t) + duI0,n(u, i) = -(Ap + Ac + ß(u))Io,n(u, t), (2)
diP^h(u, t) + du^n2(u, t) = -(Ap + Ac + ai + U2 + U1 (u))P^n2(u, t)
ni (1) n2 (1)
+ Ap (1 - V ) E c^n-i,n2 (u t) + Ac (1 - hn2 ) E CjPZ2-j K t) (3)
i=1 j=1
l/iïn,(u, t) + duP"2^(u, t) = -(Ap + bAc + a1 + a2 + m(u))P^2n2(u, t)
n1 n2
+ Ap (1 - Ô0n1 ) E CiPn2-i,n2 (u, t) + Ac (1 - S0n2 ) E CjPni- (u, t), (4)
i=1 j=1
d d n1
jtR{nln2 (u, t) + dûR^ (u, t) = -(A p + bAc + m MÜ^ (t) + A p (1 - Son1 ) E cR{n!-hn2 (t)
i=1
+ a1 f"(p"1"(u, t)+ P{n\]n2(u, t))du + Ap(1 - 00n2) E CjR^n2-i(t),
JO j=1
(5)
d d
d~tDni,n2 (u, t) + (u, 0 = -(A p + bXo + £ (u))Dnlrn2 (u, t)
ni n2
+ Ap(1 - 00nl ) E CiDni-i,n2 (u, t)+ Ao (1 - ¿0n2) E cjDni,n2-i(u, t), (6)
i = i j = i
YtR(ni,n2 (u, t) + duR^ (u, t) = -(A p + bAo + n2 (u^R^ (u, t)
ni (2) n2 (2) + Ap(1 - ¿0ni ) E CiR%-i,n2 (u, t) + Ao (1 - ^) E cRZ-(u, t), (7)
i = i j = i
d d
dt Vni,n2(u, t) + du Vni,n2(u, t) = -(Ap + bAo + Y(u))Vni,n2(u, t)
ni n2 + Ap(1 - 5oni ) E CiVni-i,n2(u, t) + Ao(1 - Ô0n2) E CjVni,n2-j(u, t). (8)
i=i j=i
The preceding set of equations must be solved under the following boundary conditions at u = 0, P^k(0, t) = (1 - r) J P,(l+l,n2 (u, t)pi(u)du + rj pg^(u, t)pi(u)du
+ J0 Pnf+ i,n2 (u t)H2(u)du + R(n+,n2 (t)ni + Apcni+ik,n2 (t)
+ f Vm+i,n2(u, t)y(u)du + ^ Rn2)+i,n2(u, t)n(u)du, (9)
P02n)2(0,t) = AoCn2+ih,0(t) + f0 hn2+i(u,t)p(u)du + EAoCi(u,t)
^ Î0,n2+i-i(u, t)du + Aoqj P0$2 (u, t)^(u)du, ( 1 0)
R(ni,n2 (0, t) = J Dn1,n2 (u, t)du, (11)
D„1,„2 (0, t) = *2 J™ 1,n2 (u, t)du + 0-2 J™ p£(2 (u, t)du, (12)
V0,0(0, t) = (1 - r) J P((,10)(u,(u)du + y pOo*(u)(u, t)du
+ Ç rO2O (u, t)n2(u)du + rO1(2(t)m + P J™ Vo,o(u, t)7(u)du, (13)
Vn1,n2 (0, t)= 0, (14) 1o,n2(0, t) = (1 - r) ^ P((,1n)2(u, t)^1 (u)du + ^ P$2p-2(u)(u, t)du
rO2(2 (u, t)n2(u)du + rO!(2 (t)m f™ Vo ,n2 (u t)Y(u)du. (15)
70
!• TO
H
P^^ (0) = P®* (0) = Dn1,n2 (0) = Vn1,n2 (0) = r(11)(2 (0) = Rgk (0) = 0,
10,0 = 1,10,n2 (0) = 0, for n2 > 1 are the initial conditions. (16)
The PGF defined as,
TO TO
I(u, t, Zo) = ^ Z^10,n2 (u, t); A(u, t, Zp, Zo)= ^ ^ Z^ ZrpAn1,n2 (u, t);
«2=1 n1=0 n2=0
TO TO
A(u, t,Zp)= £ z«An1 (u, t); A(u, t,Zp)= £ z«2A«2(u, t); (17)
«1=0 n2=0
here A = P(1), P(2), D, V, R(1), R(2).
We derive the following equations by applying Laplace transforms to equations (1) to (15) along with (16) and (17).
70 (u, s, Zo) = 70 (0, s, Zo )e-(s+Ap)u-l0 p(t)dt, (18)
P(1)(u, s, Zp, Zo) = P(1)(0, s, Zp, Zo )e-^1 (s,Z)u-/0u MW, (19)
P(2)(u, s, Zp, Zo) = P(2)(0, s, Zp, Zo )e-<fc (s,Z)u-Lu MW, (20)
V (u, s, Zp, Zo) = V (0, s, Zp, Zo )e-<fc(s,Z)u-J0u Y(t)dt, (21)
D(u, s, Zp, Zo) = D(0, s, Zp, Zo )e-^2 (s,Z)u-!0 £ (t)dt, (22)
R(2)(u, s, Zp, Zo) = R(2)(0, s, Zp, Zo )e-<fc (s,Z)u-Lu n2(t)dt. (23)
where,
^(s,z) = s + AP + A0(1 - C(zo)) + a1 + «2,
^(s,z) = s + Ap + Ao(1 - C(Zo)),
^(s,z) = s + Ap + Ao(1 - C(Zo)) + m,
<1 (s,z) = s + Ap(1 - C(zp)) + Ao(1 - C(zo)) + a1 + «2,
<£2(s,z) = s + Ap(1 - C(zp)''0 + Ao(1 - C(zo)),
$3(s,z) = s + Ap(1 - C(zp)) + Ao(1 - C(zo)) + n1,
a!(s,z) = s + Ap(1 - C(g(zo))) + Ao(1 - C(zo)) + «1 + «2,
02(s,z) = s + Ap(1 - C(g(zo))) + Ao(1 - C(zo)),
0-3(s,z) = s + Ap(1 - C(g(zo))) + Ao(1 - C(zo)) + n1.
TO
70(0, s, Zo ) =
[(1 - C(g(Zo)))Ap
1 - M(s + Ap + Ao)
s + Ap + Ao
P(2)(0, s, Zo )
AoC(Zo)70,0(s)1 - ApC(g(Zo)) - [70,0 (s)(s + Ap + Ao) - : [ M(s + A p + Ao) + C(Zo )A,
1 - M(s + Ap + Ao) i. s + A p + Ao -I - [70,0(s)(s + Ap + Ao) - 1 + [(1 - V(02(s,z)))V0,0]]
1 - M(s + A p + Ao)
s + A p + Ao
1 - B2(01 (s,z))
(Zo - A pqZo h 01 (s, z) 1 - M(s + A p + Ao)
][(1 - C(g(zo)))Ap
s + A p + Ao 1 - M(s + A p + Ao)
_
s + A p + A,
] - [M(s + Ap + Ao) + C(Zo)Ao ][B2(01(s, z))U R(2)(02(s, z))D(02(s, z))
+
1 - B2(01(s,z))
03 (s, z)J Ap[C(zp) - C(g(Zo))]
01 (s, z)
1 - M(s + A p + Ao)
s + A p + A,
- V(02(s,z)))V0,0 + P(2)(0,s,Zo) + [B2(01 (s,z)) - B2(01(s,z))
P(1)(0, s, Zp, Zo )
+
«2 R(2)(02 (s, z))D(02(s, z)) + *2~R(2\02 (s, Z))D(02 (s, Z)) + ^^
70(0, s, Zo )(V(02 (s, Z))
s,z)) - B2(01(s, 1 - B2(01 (s,z))
01 (s, z) 1 - B2(01(s, z))
03 (s, z)n 01 (s, z)
i [Zp - ((rzp + (1 - r))B1(01(s,z)) + (Zpa2R( )(02(s,z)) |d(02 (s, z))
+
1 - B1 (01 (s,z))
«171
03(s,z) L 01 (s, z)
)]
(24)
(25)
Theorem.1 When the system is operating normally, experiencing a breakdown, going on a randomized vacation, delay time to repair, or being repair, the probability generating function of the number of customers in the relevant queue will be provided using Laplace transforms.
]
]
70 (s, Zo ) = 70(0, s, Zo )
1 - M(s + Ap + Ao) s + A p + Ao
P(1)(s, Zp, Zo ) = P(1)(0, s, Zp, Zo )
1 - Bj(0j(s,z))
P(2)(s, Zp, Zo ) = P(2)(0, s, Zo )
01 (s, z) 1 - B2(01 (s,z))
01 (s, z)
(26)
(27)
(28)
V(s, zp, Zo ) = V (0, s, Zo ) D(s, zp, Zo) = D(0, s, Zo)
1 - V(fc (s, z))
fc2 (s, Z) 1 - D(fc2 (s, Z))
$2 (s, Z)
R(2\s, Zp, Zo) = R(2)(0, s, Zp, Zo)
1 - R(2)($2(s,Z))
fc(s, Z)
(29)
(30)
(31)
Proof: The following result is reached by applying the renewal theory's solution and solving the previous equations (26) to (31) with respect to u through integration.
[1 - H(u)]e-sudu
1 - h(s)
(32)
The LST of the H(u) distribution function of a random variable is expressed by h(s) . The follow-
15(1),
ing states' precise probability generating function results are as follows: Io(s,zo ), P(1)(s, Zp, Zo ),
P(2)( (31).
_(2) _ _ _(2)
,Zp,Zo), D(s,Zp,Zo), V(s,Zp,Zo) and R (s,Zp,Zo) are obtained by using equation (26) to
s
5. Steady State Analysis
Steady state analysis refers to the examination of a system's behavior once it has reached a stable condition where its key parameters remain relatively constant over time. In this state, the system's inputs and outputs balance out, resulting in a consistent and unchanging pattern of behavior. Steady state analysis is often used in various fields such as engineering, economics and physics to understand long-term behavior and performance characteristics of systems.
According to Tauberian property,
lim,sf (s)
s^r 0
lim f (t).
t—>■ M
The queue size's PGF is as follows, in spite of the system's current state:
Wq (Zp, Zo )
Nr(Zp, Zo ) Dr(Zp, Zo)'
(33)
Nr(Zp, Zo ) = Nt(Z)D2 (Z)D3(Z)fc1 (z)$2(z)$3 (z)
1 - M(Ai + A2)
Ax + A2
+ N2(z)D1 (z)D3(z)
f1 (z)(1 - B2h(z)) +
1 - M(s + A1 + A2)
s + A1 + A2
N3(z)D1 D2(1 - B2fc1 (z))f1(z),
Dr(zp, Zo) = D1 (z)D2(z)D3 (z)D4(z)fr (z)fc (z)fo(z),
G. Ayyappan, S. Nithya
ANALYSIS OF NON-MARKOVIAN BATCH ARRIVAL RETRIAL RT&A, No 3 (79) QUEUE WITH PRIORITY SERVICES IMMEDIATE FEEDBACK, PUSH OUT..... Volume 19, September 2024
where, Ni(z) Di(z) N2(z)
VoolVoi(z) - 1] - (Ap + Aq)100 + [B2n(z) + («2R(2)0(z)Dai(z) + )
«1
1 - B2oi (z)
03 (z)yL oi (z)
1 - C(g(zo))Ap
1M(A p + Aq )
Ap + Aq
AqC(zq))loo[1 - C(g(zo))Ap
1 - M(Ap + Aq )
A p + Aq
Di(z) = [zq - Aoqz,
Voo[Voi(z) - 1]] [M(Ap + Aq) - AqC(zq) 1 - B201 (z)
] + [1 - (Ap + Aq)loo 1M(Ap + Aq )-
A p + Aq
],
01 (z)
+ («2 R(2) 02 (z)D02 (z) +
][1 - C(g(zo))Ap
1 M(A p + Aq
A p + Aq
][[Bi 01 (z)
«1 ) 03 (z))
1M(Ap
1 - B201 (z)
[ M(A p + Aq ) - AqC(zq )
-1 - M(Ap + Aq )
01 (z) Aq )
A p + Aq
N3(z) = Ap[C(zp) - C(g(zo))]
A p + Aq
],
Io(o,zq) + Voo[V0i(z) - V02(z)]
[Bi01(z) - B201 (z)(«2R(2)02(z)D02(z)
+
«1
1 - B201 (z)
03(zH 01(z)
- («2R(2)02(z)D02(z) + )
1 - B201 (z)
03 (z)yL 01 (z)
]P(2)(o, zo),
D3(z) =
zp - [((1 - r) + rzp)B101 (z) + («2R(2)02(z)D02(z) + )
1 - B101 (z)
03(z) L 01(z)
01 (z 02 (z 03 (z 01(z 02 (z 03 (z
01 (z
02 (z
03 (z
: Ap + Aq (1 - C(zo)) + «1 + «2, --Ap + Aq (1 - C(zo)), : Ap + Aq (1 - C(zo)) + 71, : A p (1 - C(zp)) + Aq (1 - C(zo )) + «1 + «2, ■-Ap(1 - C(zp)) + Aq(1 - C(zo)), ■-Ap(1 - C(zp)) + Aq(1 - C(zo))+ n1, A p (1 - C(g(zo))) + Aq (1 - C(zo)) + «1 + «2, Ap(1 - C(g(zo))) + Aq(1 - C(zo)), Ap(1 - C(g(zo))) + Aq(1 - C(zo)) +
6. Stability condition
The stability requirement is a criterion that establishes whether a QS can manage incoming traffic without increasing indefinitely over time. A stable QS maintains consistent queue length and performance measurements over time, even with changing arrival rates.
We apply the normalising condition to determine Io,o.
Io,o + Io I + P(1)(1,1) + P11 (1,1) + V (1,1) + R1 (1,1) + D(1,1) + R(2) (1,1) = 1. (34)
]
[Di (1,1)D2 (1,1)d3 (1,1)(ai + p + Ao )(-E(X))]
- [
1 - M(Ap + Ao)
A p + Ao
N1 (1,1)D3 (1,1)D2(1,1)
(«1 + a2)(A p + Ao )(-E(X))m (A p + Ao )(-E(X))+ N2(1,1)d3 (1,1)D1 (1,1) f '1 (1,1) (1 - B2(«1 + «2)) + N3(1,1)D1 (1,1)D2(1,1)
h = y f1 (1,1)D3(1,1)(1 - B1 («1 + a2))]_
, {D1 (1,1)D2 (U)D3(U)(a1 + *2 )m (Ap + Ao )(-E(X
and the utilization factor is given by -1 - M(A p + Ao)
N1 (1,1)D3 (1,1)D2(1,1)
P
Ap + Ao
(«1 + a2)(A p + Ao )(-E(X))m (A p + Ao )(-E(X)) + N2(U)D3(U)D1 (1,1)f '1 (1,1) (1 - B2(«1 + a2)) + N3(1,1)D1 (1,1)D2(1,1) f1 (1,1)D3(1,1)(1 - B1 («1 + a2))]
N3 (1) = A F
Ap + Ao
1o(0,1)[1 - E(X1)]E(X),
D1 (1,1)D2 (1,1)d3 (1,1)(a1 + a2)m(A p + Ao )(-E(X))} The steady state stability requirement for the model is p < 1.
where,
-1 M(A p + Ao )-
d3(1) = 1 - [rB2(A1 + A2) + (U2(E(R(2)) + E(D)) - ^ - 1)
m
f1 (1) = [m + «1 - *2m(E(R(2)) + E(D))][-(Ap + Ao)E(X)],
1 - B1(«1 + a2)
(35)
(36)
(tt1 + ^2)
],
7. Performance Assessments
Performance measures in QS are metrics used to evaluate and quantify various aspects of system behavior, efficiency and effectiveness. These measures help assess how well a QS is performing and provide insights into its operational characteristics.
The following is the expected queue size for PC and orbit size for OC
Lq1 = ipW (zp,1)lzp=1, (37)
L12 = dkWq (1, Zo ^=1. (38)
Where,
Dr[ (1)Nr'1 (1) - Dr'1 (1)Nr'1 (1)
3(Dr'1 (1 ))2 Dr2 (1)Nr'2 (1) - Dr'' (1)Nr2 (1) 3(Dr'2 (1))2
Lq1 = 1W ^ A 1W, (39)
Lq = ^ •2\±J1'"2\A-J ^ '2W1V'2W (4o)
The following is the expected waiting time for priority queue:
Lq
% = s (4i)
The following is the expected waiting time for orbit
L
A,
Wq2 = (42)
8. Particular Cases
Case 1: In the absence of priority queue, ordinary customers arrive individually without retrials, breakdowns, or push-out mechanisms. In this scenario, the model can be simplified of queue type M/G/1 with a general randomized vacation policy.
n .. Io,o[1 - B2(Aq (1 - zo ))]{V(Aq (1 - zo )) - (1 - zo)(1 - p)V (Aq ) - 1} PII(z) =-
V (Zo )
(1 - Zo )(1 - p)V (io H Z - B2(io (1 - Zo ))} Io,o[1 - V(io (1 - Zo ))]
(1 - p)(1 - Zo)V(io)
The aforementioned outcome bears resemblance to the findings of Chen et al. [23] albeit without the inclusion of a second optional service.
Case 2: In the absence of priority queue, ordinary customers arrive individually without any breakdowns or vacations. In such a scenario, this model can be simplified a RQ of type M/G/1 with abandoned customers.
. 1o,oZ(1 - Zo )[1 - B2(ip (1 - Zo )+ ioqZo )][1 - M(io )] MZ) =- r
Pu (Z) =
D(Zo )
IojoM(io)(1 - Zo)[1 - B2(io(1 - Zo) + ioqZo)]
D(zo) where,
D(z) = B2(Aq(1 - zo) + Aoqzo){(1 - zq + qzo)(zq + (1 - zq)M(Aq)) - zQq} - zq(1 - zq). These findings align with the results reported by Krishna Kumar et al. [24].
9. Numerical Results
The numerical and graphical analyses of this model are covered in this section. We assumed that the distribution of service period, failure, repair and vacation period are all exponential.
Table 1: The impact of the priority arrival rate (Ap)
ip Io P Lq1 Wq1 Lq2 Wq2
o.5 o.9o96 o.o9o4 3.3151 4.2339 1.3o73 o.6535
o.6 o.8464 o.1536 3.6657 4.4oo2 1.4388 o.7194
o.7 o.7722 o.2278 4.o216 4.5624 1.5698 o.7849
o.8 o.6847 o.3153 4.3764 4.7233 1.6996 o.8498
o.9 o.58o9 o.4191 4.7243 4.8869 1.8275 o.9138
1.o o.4573 o.5427 5.o592 5.o592 1.9526 o.9763
1.1 o.3o9o o.691o 5.3756 5.2492 2.o737 1.o368
1.2 o.1299 o.87o1 5.6679 5.47o5 2.189o 1.o945
Table 1 exhibit that when an arrival rate (Ap) for PQ escalates, then (Lq1 /Lq2) and (Wq1 / Wq2) also rises at Ac = 2, a1 = 2, a2 = 3, y = 4, y1 = 5, n2 = 8, j = 20, p = 0.3, t = 15, & = 8, q = 0.4, r = 0.1 and Ap = 0.5 to 1.3.
Table 2: Impact of hard failure repair rate (^2)
^2 I0 P Lqi Wqi Lq2 2 £
5.0 0.0435 0.9565 3.2910 3.2910 1.6708 0.8354
5.5 0.1807 0.8193 3.0328 3.0328 1.6612 0.8306
6.0 0.2742 0.7258 2.8290 2.8290 1.6507 0.8253
6.5 0.3421 0.6579 2.6642 2.6642 1.6392 0.8196
7.0 0.3937 0.6063 2.5283 2.5283 1.6267 0.8133
7.5 0.4344 0.5656 2.4145 2.4145 1.6128 0.8064
8.0 0.4673 0.5327 2.3177 2.3177 1.5975 0.7987
8.5 0.4944 0.5056 2.2346 2.2346 1.5804 0.7902
9.0 0.5173 0.4827 2.1623 2.1623 1.5612 0.7806
Table 2 exhibit that when an hard failure repair rate (n2) escalates, then (Lq1/Lq2) and (Wq1/Wq2) also decreases at Ap = 1, Ac = 2, a1 = 1, a2 = 5, y = 4, n1 = 4, 7 = 15, p = 0.4, t = 10, & = 12, q = 0.3, r = 0.1 and n2 = 5.0 to 10.0.
Table 3: Impact of ordinary arrival rate (\o)
Ao I0 P Lqi Wq1 Lq2 2 sr
1.0 0.8423 0.1577 1.6669 3.3339 0.1541 0.1541
1.1 0.8402 0.1598 1.7448 3.4896 0.1954 0.1777
1.2 0.8371 0.1629 1.8228 3.6455 0.2440 0.2034
1.3 0.8329 0.1671 1.9008 3.8016 0.3005 0.2311
1.4 0.8275 0.1725 1.9790 3.9579 0.3655 0.2610
1.5 0.8204 0.1796 2.0572 4.1144 0.4396 0.2931
1.6 0.8116 0.1884 2.1355 4.2710 0.5237 0.3273
1.7 0.8004 0.1996 2.2139 4.4277 0.6182 0.3636
Table 3 exhibit that when an ordinary arrival rate (Ac) for PQ escalates, then the (Lq1 /Lq2) and the (Wq1/Wq2) also rises at Ap = 0.5, = 0.8, a2 = 4, y = 3, m = 3, 7 = 17, p = 0.3, t = 10, & = 12, q = 0.2, r = 0.2, n2 = 7 and Ac = 1.0 to 2.0.
We obviously follow the exponential distribution for service time, breakdown, repair and vacation time in graphical representations. Figures 2-4 illustrate the 2D graphs. (Lq1, Lq2) increases when the priority arrival rate (Ap) increases, as demonstrates in Figure 2. Figure 3 demonstrates the behaviour of the queue sizes (Lq1, Lq2), which depends on the hard failure repair rate (n2). The length of the queue grows as the soft failure rate improves. Figure 4 depicts the behaviour of the queue sizes (Lq1, Lq2), which is affected by the average customer arrival rate (Ac).
6
g
n le
e u e
u
q
d te ct
e p
x
w
o.6 o.8 1 1.2
Priority arrival rate (Ap)
Lq1 Lq2
Figure 2: Lq Vs Ap
4
2
Figure 3: Lq Vs ni
g
n le
e u e
u
q
d te ct
e p
x
w
Lq1 Lq2
1.2 1.4 1.6 1 ordinary customers arrival rate (Ao)
Figure 4: Lq Vs Ao
2
1
o
1o. Conclusion
In this research, we examined a single server retrial queueing system with non-preemptive priority service, immediate feedback, push out, differentiated breakdowns, delayed repair, randomized vacation. The analytical findings that are supported by numerical examples can be used to a wide
range of real-world situations to produce results. The supplementary variable technique is used to determine the PGFs for the number of users in the system when it is free, busy, and under repair. The system's and orbit's average queue lengths contain many expressions. The mean busy period and other significant system performance measures are obtained. The conditional decomposition law is finally demonstrated to be effective for this retrial queueing system.In real-world queueing scenarios, our queueing system is more flexible when dealing with real-time systems used by many industries.
References
[1] Ahuja, A., Anamika Jain and Madhu Jain. (2018). Finite population multi-server retrial
queueing system with an optional service and balking, International Journal cf Computers and Applications, 25(3):1-8.
[2] Ammar, S.I. and Rajadurai, P (2019). Performance Analysis of Preemptive Priority Retrial Queueing System with Disaster under Working Breakdown Services, Symmetry, 11(3): 419.
[3] Ayyappan, G and Thilagavathy, K. (2021). Analysis of MAP(1),MAP(2)/PH/1 Non-preemptive Priority Queueing Model Under Classical Retrial Policy with Breakdown, Repair, Discouragement, Single Vacation, Standby Server, Negative Arrival and Impatient Customers, International Journal of Applied and Computational Mathematics, 7:174-184
[4] Ayyappan, G. and Deepa, T. (2021). Analysis of an Admission Control M[X1 /G(a, b) /1 Queue with Multiple Vacation, Restricted Re-service, Closedown and Setup Times, Proceedings of the National Academy of Sciences, India Section A: Physical Sciences,91:35-47.
[5] Ayyappan, G. and Gowthami, R. (2021). A MAP/PH/1 queue with Setup time, Bernoulli vacation, Reneging, Balking, Bernoulli feedback, Breakdown and repair, Reliability: Theory & Applications, 16(2):191- 221.
[6] Azhagappan, A. and Deepa, T. (2023). Variant Impatient Behavior Of A Markovian Queue With Balking Reserved Idle Time And Working Vacation, RAIRO Operations Research, 54:783793.
[7] Begum, A. and Choudhury, G. (2022). Analysis of an M/(G1, G2)/1 Queue with Bernoulli Vacation and Server Breakdown, International Journal of Applied and Computational Mathematics, 9(6):1-22.
[8] Choudhury, G. and Kalita, C.R. (2018). An M/G/1 queue with two types of general heterogeneous service and optional repeated service subject to server's breakdown and delayed repair, Quality Technology & Quantitative Management,15(5):622-654.
[9] Danilyuk, E., Plekhanov, A., Moiseeva, S.P and Janos Sztrik. (2022). Asymptotic Diffusion Analysis of Retrial Queueing System M/M/1 with Impatient Customers, Collisions and Unreliable Servers ,MDPI Journals, 11(12):699.
[10] D'Arienzo, M.P., Dudin, A.N., Dudin, S.A and Manzo, A. (2020). Analysis of a retrial queue with group service of impatient customers, Journal of Ambient Intelligence and Humanized Computing, 11:2591-2599.
[11] D'Apice, C., Dudin, A., Dudin, S. and Manzo, R. (2022). Priority queueing system with many types of requests and restricted processor sharing, Journal of Ambient Intelligence and Humanized Computing, 14:12651-12662.
[12] Gao, S., Zhang, J., and Wang, X. (2020). Analysis of a Retrial Queue With Two-Type Breakdowns and Delayed Repairs, IEEE Access, 8:172428-172442.
[13] Jose, J.K. and Deepthi, V. (2021). M/PH/1 queueing model with re-servicing, Communications in Statistics - Simulation and Computation, 52(7):2865-2885.
[14] Krishna Kumar, B., Rukmani, R., Thanikachalam, A. and Kanakasabapathi, V. (2018). Performance analysis of retrial queue with server subject to two types of breakdowns and repairs, Operational Research, 18:521-559.
[15] Li, Q., Guo, P. and Yulan Wang (2020). Equilibrium analysis of unobservable M/M/n priority queues with balking and homogeneous customers, Operations Research Letters, 48(5):674-681.
[16] Nazarov, A., Janos Sztrik, Anna Kvach and Adam Toth. (2o22). Asymptotic Analysis of Finite-Source M/GI/1 Retrial Queueing Systems with Collisions and Server Subject to Breakdowns and Repairs,Methodology and Computing in Applied Probability, 24:15o3-1518.
[17] Pavai Madheswari, S., Krishna Kumar, B. and Suganthi, P. (2o19). Analysis Of M/G/1 Retrial Queues With Second Optional Service And Customer Balking Under Two Types Of Bernoulli Vacation Schedule, RAIRO Operations Research, 53:415-443.
[18] Takagi, H. Queueing Analysis: A Foundation of Performance Evaluation, Vol. 1: Vacation and Priority Systems. Elsevier Science Ltd, North Holland, Amsterdam, New York, (1991).
[19] Ke, J. C., Wu C. H. and Zhang, Z. G. (2o1o). Recent developments in vacation queueing models : A short survey. International Journal of Operations Research, 7(4):3-8.
[20] Gao, S. and Yao, Y. (2o14). An MX/G/1 queue with randomized working vacations and at most j vacations. International Journal of Computer Mathematics, 91(3): 368-383.
[21] Gao, S., Liu B. and Wang, J. (2o18). On a retrial queue with abandoned customers and multi-optional vacations. International Journal of Computer Mathematics: Computer Systems Theory, 3(3):177-195.
[22] Krishnamoorthy, A., Joshua V. C. and Mathew, Ambily P. (2o17). A retrial queueing system with abandonment and search for priority customers, In Vishnevskiy, V. M., Samouylov K. E. and Kozyrev, D. V., editors, Distributed Computer and Communication Networks. Springer International Publishing. Cham. pages 98-1o7.
[23] Chen, Y., Lin X. W. Wei C. M. and Fan, Z. (2o18). An M/G/1 queue with second optional service and general randomized vacation policy, In Cao, Bing-Yuan, editor, Fuzzy Information and Engineering and Decision. Springer International Publishing. Cham. pages 297-3o7.
[24] Krishnakumar, B., Vijayakumar A. and Arivudainambi, D. (2oo2). An M/G/1 retrial queueing system with two-phase service and preemptive resume. Annals of Operations Research, 113(1):61-79.