ANALYSIS OF SINGLE SERVER FEEDBACK RETRIAL QUEUE WITH BERNOULLI WORKING VACATION AND STARTING FAILURE
The suggested queueing model describes a single-server feedback retrial queueing system with starting failure, Bernoulli working vacation and vacation interruptions. The server departs on a working vacation as soon as orbit is empty. During the working vacation period, the server provides a slower level ofservice. The supplementary variable method was utilized to determine the steady-state probability-generating functions for the system and its orbit. If there are consumers in the system at the end of each vacation, the server becomes idle and ready to serve new customers. The average busy time and the average busy cycle are presented as important system performance indicators. Additionally, the adaptive neuro-fuzzy interface system has compared the numerical results with the neuro-fuzzy results. Finally, particle swarm optimization (PSO) were utilized to obtain the best (optimal) cost for the system in this study. We have examined the convergence of these optimization strategies.
Keywords: Retrial queues, Feedback, Supplementary variable technique, Starting Failure and Working Vacation,ANFIS.
In a queueing system (QS), queues involving continuous tries occur when a consumer comes and identify the server is occupied. The client is instructed to leave the service region and join a virtual area referred to as the 'orbit'. Subsequently, the customer within the orbit can make a service request after a period of time. In a vacation periods, the server halts its service entirely, becoming unavailable to the primary clients for a short duration, which is termed a "vacation." However, during the working vacation (WV) period, the server provides services to consumers, albeit at a reduced service rate. Also, the server's vacation may be ignored if customers arrive during the vacation period, and the server may resume operation in its regularly scheduled manner. It is known as the vacation interruption(VI) strategy. Major uses for this QS include delivering network services, online services, file transfer services, mail services and so on. A more realistic RQ with feedback happens in many real-world scenarios; for instance, in multiple-access telecommunications systems, where data returned as failures is forwarded again, it may be treated as a retrial queue with feedback.
In an M/G/1 retrial queue (RQ) with general retrial times, consumers who find the server busy join the orbit according to the first-comefiirst-served (FCFS) principle as studied by Gomez-Corral [1]. Such an instance occurs in certain communication protocols, in production lines at
Keerthiga S1 and Indhira K*
-^Department of Mathematics, School of Advanced Sciences, Vellore Institute of Technology, Vellore - 632 014, Tamil Nadu, India. keerthiga.2020@vitstudent.ac.in, kindhira@vit.ac.in.
Abstract
1. Introduction
1.1. Survey of Literature
stores, etc. The RQ has been extensively studied by Falin and Templeton [2], Artalejo and Corral [3], Artalejo [4], etc. Many authors have investigated a single server retrial queue (SSRQ) with WVs and VIs,including Zhang and Hou [16], Gao and Liu [6], Gao et al. [5], Zhang and Liu [22], and Rajadurai et al. ([7], [8], [9]). Mokaddis et al.[10] explored the M/G/1 retrial queue with Bernoulli feedback,Starting Failure (SF) and a single vacation (SV). Clients in orbit connect to the server via FCFS discipline, and an arbitrary distribution is assumed for the retry time. The server goes on vacation when there are no clients on the system. If the server comes back from vacation and there are no consumers, it waits for the first client to arrive on the system from the outside.
Krishna Kumar et al.[11] researched a RQ with feedback and a server exposed to SF, as well as a general stochastic decomposition rule for M/G/1 vacation models. Rajadurai [12] investigated a single server preemptive priority RQ based on Bernoulli working vacation (BWV) and VIs. Rajadurai et al.[13] explored a SSRQ system with BWV and VI. Performance indicators and analytical illustrations are provided. Jain and Kumar [26] analyzed bulk arrival general service RQ subject to balking, feedback and vacation interruption under multiple WV policy. Pazhani Bala Murugan and Keerthana [27] investigated an M/G/1 feedback RQ with WV and a waiting server. Keerthiga and Indhira [28] examined SSRQ with two phases of service for retrial customers. Agarwal et al. [29] discussed detection of optimal WV service rate for retrial priority G-queue with immediate Bernoulli feedback. Rachita Sethi et al.[17] researched a threshold-based repair facility for machining systems with a WV approach. WV was established to allow repairmen to offer service at a reduced rate as opposed to entirely discontinuing operations. The idea of F-policy is used to govern its arrival in the system. The implementation of a threshold N-policy to start the repair reduces the system's cost. Performance measurements are computed utilizing the 4th-order Runge-Kutta approach, and the numerical findings obtained are compared to the adaptive neuro fuzzy inference system (ANFIS).
Charu Bhargava and Madhu Jain [18] studied the "Modelling and Analysis of a Markovian multi server queue with an (e,d) SV procedure, server failures and repairs". Some stationary performance indicators are established after service completion using the matrix geometric technique. Additionally,the direct search method is utilised to estimate the best no. of idle, vacationing, and total no. of servers at the most affordable price. Also, the acquired numerical outcomes have been compared using a soft computing technique (SCT) based on an ANFIS. Radhika Agarwal et al.[19] analyzed the performance metrics that are used in improving service standards using the SVT and compared the analytical outcomes to the neuro fuzzy outcomes via the ANFIS (SCT). In addition, single and bi-objective minimization issues are explored with minimum attained via "PSO and a multi-objective GA" respectively.
In this research, we have extended the work of Rajadurai et al. [13] by including the ideas of feedback and SF. By using PSO, we have also performed a cost analysis of the model under consideration. Because the suggested solution improves repeatedly and the system gives us the best option that is feasible, this approach has gained a lot of reputation in recent years. When it comes to queueing analysis, this technique may be used to get productive outcomes, whether the goal is to save overall costs or maximize performance metrics. This framework aids in our analysis of various real-world queuing scenarios, allowing us to enhance the customer experience. To that extent, this article contributes. In areas with heavy traffic and congestion, this kind of project is highly pertinent and beneficial.
To the best of author's knowledge, there has been no previous research that has examined in this work. Therefore, to fill up this gap, in this article, we consider the feedback RQ with WV and VI subjected to server breakdown and repair. SVT has been used, and for some of the variables, a 3D graphical representation has also been provided. To attain optimal operating conditions, minimize expected costs, and maximize economic performance PSO a well-known meta-heuristic technique are applied. The aforementioned framework may be used in a wide variety of situations, including but not limited to: telephone switching, telecommunications, computer networks, online ticket booking centers, aviation traffic control, quality control procedures, and inspection testing of items. The purpose of this investigation is to estimate the queue length and orbit size dist., which will be implemented to calculate the system's performance metrics.
The following is an overview of our article: Section 2 provides a detailed discussion of the queueing paradigm. Section 3 specifically determines the system's steady state (SS) behavior and the queue length's PGF at a random epoch. Section 4 includes various substantial indicators of system behavior. Section 5 discussed particular cases. Sections 6 and 7 provide numerical outcomes and cost optimization. Finally, Section 8 provides a conclusion and overview of the study.
2. Model description A comprehensive explanation of this framework is given below:
The arrival process: New customers join the system from the outside, according to a Poisson process (PP), at a rate of p.
The retrial process: According to FCFS discipline, when a customer visits while the server is occupied or unavailable, the consumer departs the service area and joins a group of blocked clients known as "orbit." The server appears to be accessible to the clients at the front of the orbit queue. Every customer's successive inter-retrial duration's are determined by an arbitrary probability distribution function (PDF) B(x), with an associated density function (df) b(x) and the "Laplace-Stieltjes transform" (LST) ft* (0).
The service process: When a new or repeated consumers enters at the server while it's idle, the server promptly begins its regular service for the incoming customers. The service time follows a general dist., and PDF C(x), a df c(x), a LST a* (0) & the 1st and 2nd moments are C- and C2.
The Bernoulli working vacation process: The server goes on a WV whenever the orbit is
empty, and the duration of this vacation follows an exponential dist., with a specified parameter 7. If a consumer visits during vacation time, the server will continue to operate at a slower service rate. The WV time is a slower-paced operating period. In the event that any clients in the orbit reach the instant of service completion during the vacation time, the server will end the vacation and return to its normally busy state, which is known as VI. On the other hand, if there are no clients in the system at the completion of the vacation, the server rejoins the system and waits to serve a new client with prob.ri (SWV) or makes for another WV with prob. r2 = 1 — r1 (MWV). When a vacation is over and there are still consumers in its orbit, the server resumes usual operation. During the WV period, the service time is determined by a general random variable Hv with a dist., function Hv(t), LST H* (q>) & the 1st and 2nd moments are h1 and h2, respectively.
Feedback Procedure: After getting their normal services, dissatisfied customers have two options: they may either exit the system with probability w = (1 — w) or they can return to the orbit as unsatisfied clients and get a service again with probability w .
Starting Failure: The customer will almost definitely start receiving service right away if the server is successfully activated. If the server is unable to start the service, the consumer exits the service area, enters the orbit, and repeats the request for the service after some time. The server is instantly repaired if a failure occurs. SF happens with prob. A and successful service begins with prob.A.
Repair Process: If the server fails to start, the repair process begins immediately. During the repair process, the server refuses to serve external or repeat consumers. Repair times have a distribution function I(x) and a corresponding density function i(x), and the first two moments are I1 and I2, respectively.
The system's stochastic processes are considered to be independent of each other.
2.1. Practical justification of the recommended paradigm
The suggested scenario has beneficial applications in a telecommunications. For example, we investigate a communication system designed for making reservations at restaurants. Let us consider a scenario in which a restaurant uses a phone system to accept reservations and provide a range of other services. client can use this system to reserve a table for themselves. The supervisor who answers all calls is in charge of this phone system. The consumer is able to pick up the phone (leave the system) or inquire regarding event reservations, purchase tickets for an upcoming musical performance, etc. after reserving a table. A caller must reaffirm their reservations if there is a possibility of a misinterpretation stemming from an unclear network or other related difficulties (feedback). When the manager is occupied overseeing other areas of the restaurant, he is unable to answer calls (vacation mode). In these circumstances, the junior manager often serves, albeit somewhat more slowly (WV).In this stage, the supervisor returns right away (i.e., a vacation interruption happens) if there are any calls in the system after the phone call is over (at service completion). However, the supervisor continues to take care of other restaurant-related matters if no calls come in after completing his secondary work (vacation mode). It is likely that when a consumer calls, the line is busy and that the client will call back after some period of time (retrial). It is possible that during a phone conversation, a bad signal, inadequate network coverage, or a virus attack (SF) might occur, causing the client to lose service. Once the communication system's signal is repaired, it functions flawlessly.
In steady state (SS), we presume that B(0) = 0, B(œ) = 1, C(0) = 0, C(œ) = 1, and Hv(0) = 0, Hv(œ) = 1,1(0) = 0,I(œ) = 1, are continuous at (p = 0. So that the function q((p), Z(Pp), K(Pp), and u( (p), are the hazard rates of the conditions (retrial, normal service/vacation and repair) are
Thus, the state of the system B0(£),C0(£), HV(£), and I0(£) are required to construct a bivariate Markov process {N(£); £ > 0 }, where L(£) belongs to the server stage (0,1,2,3) based on if the server is idle, typical operative period,slow service and repair time.
3.1. Ergodicity Condition
Let ; & = 1,2,...} represent a series of epochs in which either a service time is reduced or completed. U& = {L(£&+), X(£&+)} is a random vector sequence. The embedded Markov chain generated by the RQ system. Its state space is S={0,1,2,3} x N.
3. Scrutiny of the steady state probabilities
0, if the server is free
1, if the server is active period
2, if the server is operative mode on WV period
3, if the server is on repair
3.2. Theorem
The embedded Markov chain {Ua; aeN} is ergodic iff p < B(p) for our system to be stable, where p = ApC1 + A (1 + pl1) + Aw.
3.3. System of governing equations
For the procedure {N(£),£ > 0}, we specify the prob., <0(£) = P{L(£) = 0,X(£) = 0} and X0(£) = P{L(£) = 1, X(£) = 0} the probability densities, Xa(<p,£)dcp = P{L(£) = 1, X(£) = a, cp < 8°(£) < <p + d<p},
for £ > 0, p > 0 and a > 1. Ya(<p, £)dcp = P{L(£) = 2, X(£) = a, p <C0(£) < p + dcp},
for £ > 0, cp > 0, a > 0. Av,a(cp, £)dcp = P{L(£) = 3, X(£) = a, cp <H°v(£) < cp + dcp},
for £ > 0, cp > 0 and a > 0. na(cp, £)dcp = P{L(£) = 4, X(£) = a, cp <I0(£) < p + dcp}, for £ > 0, cp > 0, a > 1.
In subsequent parts, the following probabilities are applied:
1. The prob., of the server being idle and on WV at time £ is denoted by <0 (£).
2. The prob., of the server being idle and on typical active period at time £ is denoted by x0(£)■
3. If there are accurately a clients in the orbit at time £ and the elapsed retrial time of the test clients undergoing retrial is between cp and cp + dcp, then the prob., that this is the case is Xa(^ £).
4. When there are a consumers in the orbit, the prob., of the test customer's elapsed regular service time ranging between cp and cp + dcp is Ya (cp, £).
5. Av,a(cp, £)dcp and na(cp, £)dcp is the prob., that there are precisely a patrons in the orbit, with the elapsed (reduced service time and repair time) of the test patron being between cp and cp + d cp at time £.
Suppose that the sequel fulfills the stability condition, thus we can provide x0 = lim^^^x0(£)
and limiting densities are
Xa (cp) = lim£^Xa (cp, £) for cp > 0 and a > 1.
Ya(cp) = lim£(cp, £) for cp > 0 and a > 0.
Av,a(cp) = lim£^^Av,a(cp, £) for cp > 0 and a > 0.
na(cp) = lim^^n(cp, £) for cp > 0 and a > 1.
Applying the SVT, we create the following system of equations.
PX0 = 7^0 (1)
p TO p TO
(p + 7)^0 = Yr2$0 J Av,0(p)k(pp)dpp + J Y&(p)Z(pp)dp (2)
dd"(p) + (p + n(p))x&(p) = 0& > 1 (3)
dp Y0 (p ) + (p + Z (p ))Y0( p ) = 0, & = 0. (4)
dpY&(p) + (p + Z(p))Y&(p) = pY&- (p), & > 1 (5)
dp A0,v (pp ) + (p + 7 + k( pp ))A0,v (pp) = 0, & = 0. (6)
dpA&,v(p)) + (p + 7 + k(pp))A&,v(pp) = pAv,&-1(pp), & > 1. (7)
dp n (pp ) + (p + u( pp ))n0 (p> ) = 0, & = 0. (8)
dpn&(p) + (p + u(p))n&(pp) = pn&_x(p),& > 1. (9) At pp = 0 the steady state boundary conditions are as follows:
/* TO !• TO /* TO
(0) = co Y& (p )Z (pp )d pp + o/ Y&_1 (p )Z (pp )d pp + co Av,& (pp )k( pp )d pp 0 0 0 p TO /• TO
+ oy Av,&-1(pp)k(pp)dpp + y n&(p)u(pp)dpp (10)
/• TO /• TO
Y0(0) = Ay X1 (p)n(pp)dpp + ApX0 + 7 Jo A0,v(pp)dpp, & = 0 (11)
/»TO /»TO /»TO
Y&(0) = A/ X&+1 (p)n(pp)dpp + AW Y&(p)dp + 7/ A&,v(p)dp, & > 1 (12)
70 J 0 J 0
(■» = r " = 0 (13)
[0, a > 1
/" œ
n1 (0) = Ay X1 (P)n(P)dp( + Âpx0 (14)
/• œ /• œ
n(0) = A J^ Xa(p)n(pP)dpP + Ap J^ Xn-t(pP)dpP, a > 2 (15)
The normalizing condition is
œ ,, œ œ ,, œ œ ,, œ
X0 + ^0 + E / Xa(pp)dpp + E / Ya(pp)dpp + E / Aa,v(pp)dp
a=1 a=0 0 a=0 0
œ œ
+ E/ na ( pp )d pp = 1
a=1 0
(16)
3.4. The steady state solution
The PGF is used to compute the steady state solution for the RQ model. To solve the aforementioned equations, the generating functions for |V| < 1 are described as below:
TO TO
X(cp, V) = E Xa(cp)Va; X(0, V) = E Xa(0)Va;
a=1 a=1
TO TO
a
Y(cp, V) = E Ya(cp)Va; Y(0, V) = E Y0(0)Va; i = 1,2
a=0 n=0
TOTO
Av (cp, V) = E Ava(cp)Va; Av (0, V) = E Av,a(0)Va;
a=0 a=0
TOTO
n(cp, V) = E n(cp)Va; n(0, V) = E na(0)Va
a=1 a=1
Next multiply the SS eqn. and SS boundary conditions from (3) to (15) by Va and adding over a, (a = 0,1,2,...)
—x(cp, V) + [p + n(cp)]x(cp, V) = 0 (17)
—Y(cp, V) + [p(1 — V)+ Z(cp)]Y(cp, V) = 0 (18)
—Av(cp, V) + [y + p(1 — V) + k(cp)] Av(cp, V) = 0 (19)
—n( cp, V) + [p(1 — V)+ v( cp )]n( cp, V) = 0 (20)
Solving the partial differential eqns. (17) to (20), we obtain
X(cp, V) = x(0, V)[1 — B(cp)]e—p* (21)
Y( <p, V) = Y(0, V) [1 — C (cp )]e—F (V) * (22)
Av(cp, V) = Av(0, V)[1 — Hv(cp)]e—Fv(V)* (23)
n( cp, V) = n(0, V) [1 — I (p )]e—F (V) * (24)
where F(V) = p(1 — V), Fv(V) = y + p(1 — V)
Multiplying equation (10) and (12,13,15) by appropriate powers of V, adding over n with few mathematical manipulations, we obtain
!• TO !• TO
x(0, V) = (w + wV) J Y( cp, V)Z (cp )d cp +(w + wV) J Av (p, V)K(p)dp (25)
p TO
+ J n( (p, V)v( <p )d <p — (p + Yn)<0
A f TO f TO f TO
Y(0, V) = V Jo x(cp, V)n(p)d<p + XpJ^ x(cp, V)d<p + yJ Av(cp, V)d<p + Apx0 (26) Av (0, V) = p<0 (27)
TOTO
n(0, V) = XVp j^ x(cp, V)dcp + XJ x(cp, V)n(<p)d<p + VpXx0 (28)
Using eqn (21,23 and 27) in eqn (26)
" V + (1 — V)B (p)'
Y(0, V) = XX(0, V)
+ p<0V (V) + XYH <0 (29)
Keerthiga S,Indhira K
ANALYSIS OF SINGLE SERVER FEEDBACK RETRIAL QUEUE WITH RT&A, No 2 (78) BERNOULLI WORKING VACATION AND STARTING FAILURE_Volume 19 June, 2024
Similarly using equation (21) in (28)
n(0, 0) = 0jnAxo + Ax(0, 0) [0 + (1 - 0)B(p)] (30)
Substituting equations (22),(23) and (24) in (25), we obtain
x(0,0) = (a + a0)Y(0,0)C(F(0)) + (a + wd)Av(0, 0)H(Fv(0)) + n(0,0)I(F(0)) - p$0 - YH$0
(31)
Using equations (27),(29) and (30) in equation (31)
(a + a0)p$0[Hv (Fv (0)) + V (0)C (F (0))] + (a + a0)Ajr1 faC (F (0))
x(0 0) = 0 j _+ 0A7n$01 (F(0)) - p$0 - Yr1 $0 y (32)
^) * 0 - [(a + a0)AC(F(0)) + 0A1 (F(0))][0 + (1 - 0)B(p)] < (3)
substituting equation (32) in (29) and (30),we obtain
A[(a + a0)p$0[H v (Fv (0)) + V (0)C(F (0))] + (a + a0)Ajn $nC(F (0)) Y(0 0) _ , + 0Am$0I (F(0)) - u$0 - rn$0] [0 + (1 - 0)B(p)] + p$0V(0) + Ayn $0
д - [(а + wd)AC(F(д)) + дАХ(F(д))] [д + (1 - д)В(ц)]
А[(а + ад)цфо[Н (Fv (д)) + V (д)С (F (д))] + (а + ад)А7Гх фо C (F (д)) |
(33)
п(0 д) i_+ Щупфо!(F(д)) - цфо - 7Пфо] + дА7Пфо (34)
(0' д) ^ д - [(о + wd)AC"(F(д)) + дАX(F(д))][д +(1 - д)В(ц)] ' (3)
Substituting equations (27) and (32) to (34) in (21) to (24)
( (а + ад)цфо\Hv(Fv(д)) + V(&)C"(F(д))] + (а + ад)А7гхft,C(F(д))
х(ф д) = д J _+ дАтлфрХ(F(д)) - цфо - 7г1фо \ (35)
Л(^ ) | д - [(а + ад)А€^(д)) + д^X(F(д))][д + (1 - д)В(ц)] ' ( )
х [1 -B(ф)]е-цф
Г А[(а + ад)цфо[Hv(Fv(д)) + V(&)C"(F(д))] + (а + ад)А7тфC(F(д)) Y( _ д) = J + дА7пфо1 (F(д)) - цфо - тпфо] [д + (1 - д)В(ц)] + цфоV(д) + Атпфо (^ ) | д - [(а + ад)AC"(F(д)) + дАХ (F(д))][д + (1 - д)В(ц)]
(36)
х [1 -C(ф)]e-F
Лv (ф, д) = цфо[1 -Hv (ф )]e-Fv (37)
J А[(о + ад)цфо [Hv(Fv(д))+ V(&)C(F(д))] + (а + ад)А7п(д)) I
п(гй д) = I _+ дА7Г\фоX(F(д)) - цфо - фо] + дА7г1фо I (38)
() | д - [(а + ад)АС(Г(д)) + дАХ (F(д))][д + (1 - д)В(ц)]( )
х [1 -I(ф)]e-F(д)*
3.5. Theorem
The stationary dist. of the no. of clients in the orbit while the server is free,normal operative service, slow service and the prob. that the server is idle is described by p < B(p) under the stability condition
x(û) = û{
(œ + aû)fo [H (Fv (û)) + V (û)C(F (0))] + (a + aû) frn foC(F (û))
_+ û^foI(F(û)) - fo - ^
û - [{a + œû)AC(F(û)) + ûAI(F(û))] [û + (1 - û)B(p)]
x [1 -B(p)]
(39)
Y(û)
A[(œ + aû)pfo [Hv (Fv (û)) + V (û)C(F (û))] + (a + œû)A7n foC(F (û)) + ûX7ri$oï(F(û)) - pfo - Ynfo] [û + (1 - û)B(p)] + pfoV(û) + A7ri$o
û - [(œ + œû)AC(F(û)) + ûAI(F(û))] [û + (1 - û)B(p)]
(40)
[1 -C (F (û))] p(1 - û)
Av (û)
Mo V (û) Y
(41)
n(û)
A[(a + aû)pfo [Hv (Fv (û)) + V (û)C(F (û))] + (a + aû)A7r1po C (F (û)) ï
_+ ûAjn fol (F (û)) - Mo - Yn fo ] + ûAjn fo I
û - [(a + aû)AC (F (û)) + ûAI (F (û))] [û + (1 - û)B (p)]
(42)
v [1 -I(F(û))]
X P(1 - û)
Proof. Taking the equations. (35) - (38) and integrating them with regard to p and obtain the partial PGF's x(û) = fo°° x( p, û)d p, Y(û) = /o°°Y( p, û)d p, Av (û) = /"Av ( p, û)d p,
n(û) = f°n( p, û)d p,.
We can find the prob. that the server is free by using the normalisation condition (xo) and (fo) by establishing functions as, when there is no consumer in the orbit û = 1 in (3.39) - (3.42) and using the "L'Hospital rule" if it is required, we examine xo + fo + x(1) + Y(1) + Av(1) + n(1) = 1. ■
x
3.6. Theorem
The stability constraint p < B(p) used to determine the PGF of the no. of clients in the system and the orbit size dist. at a stationary point in time is given by
Hs w==Ne($ (43)
h o=Ne$! (44)
Proof. The "PGF of the no.of consumer in the system (Hs(0)) and in the orbit (Ho(&))" is calculated by applying Hs(0) = xo + po + x(#) + 0{Y(0) + Av(o)} + n(0). and H0(0) = Xo + Po + X(0) + {Y(0) + Av(0)} + n(0). Insert the eqns. (39) - (42) in the earlier results,then the eqns. (43) and (44) may be computed immediately. ■
4. Measures of system performance
This part calculates many appropriate system prob., system efficiency metrics, and signifies the mean busy period and cycle that occur while the system is in various phases.
4.1. System state probabilities
By putting 0 ^ 1 in equations. (39) - (42) and applying "L Hospital's rule" wherever possible. we obtain the following findings.
(i)Pr(The server being available for the duration of the retrial)
x(1) = $o [1 -B(y)]
[[%[1 - Hv(7)] - FCi[1 - Hv(7)]] - -bri[œ + FCi]
_+ fa [1 - yh ]
B(y) + Xœ - XyC1 - X(1 + yZ1 )
(ii)Pr(The server is operative on usual service period)
Y(1) =
$0XC1[[[Y [1 -Hv(7)] - yC [1 -Hv(7)]]] - X7T1 [œ + yC1
_+ Â7T1 [1 - yh] + y[yH + Y [1 - Hv(y)]]]
B(y) + Xœ - XyC1 - X(1 + yI1 )
(45)
(46)
(iii)Pr(The server is on WV)
Av (1)
$oy[1 -Hv (7)]
7
(iv)Pr(The server is under repair time during usual active period)
n = n(1)
$0XI1[[[7 [1 - Hv(7)] - yC1 [1 - Hv(7)]]] - X7T1 [œ + yC1]
_+ X7T1 [1 - yh ] + X7T1 ]
B(y) + Xœ - XyC1 - X(1 + yI1 )
(47)
(48)
4.2. Average system size and its orbit
In a steady state, the system,
(i) Differentiating the equation (44) and the predicted no. of clients in the orbit (Lq) is established with regard to 0 and 0 = 1.
Lq = H0(1) = lim ÍH0(d) 7 ^ 1 ad
Ne¡(1)De¡(1) - De¡(1)Ne¡(1) 3(D< (1))2
(49)
(ii) The predicted no. of clients in the system (Ls) is determined by differentiating the eqn. (43) with regard to 0 and giving 0 = 1 yields.
Ls = Hs (1) = lim ddHs(d) = $0
NeS'(1)Dej(1) - De'¡(1)Ne¡(1) 3(Deq (1))2
(50)
(iii) The mean waiting time of consumers in the system and queue [Ws and Wq ] are computed utilizing "Little's method" Ws = L and Wq = respectively.
4.3. Mean busy period and the busy cycle
Let A(Tb) and A(Tc) be the predicted sizes of the busy period and cycle, respectively under steady state conditions. The outcomes are derived directly from the justification of a different renewal procedure [5] , which concludes in
$0
A(To)
A(Tb)+ A(To)'
; A(Tb) = \ ($0 - l); A(Tc) = ¿0 = A(To) + A(Tb).
(51)
where To is the period of time spent in the system's null state. Because there is an exponential difference in time between the arrivals of two customers and A(T0) = (1/p) with the parameter p.
5. Particular Cases
We examine a few real-world examples of our technique that are consistent with the existing
research in this area.
Case (i): No feedback,No VI and No SF
If w = 1,A = 1, and y = 0. The model may be lowered to a M/G/l RQ with WV and the findings
match those of Arivudainambi et.al.[14]
Case (ii): No retrial, No feedback and No starting failure.
Let r2 = 0,w = 1, A = 1 and BB(p) ^ 1. our framework has been simplified to an "M/G/l queue with WVs and VI". Our results agree with Zhang and Hou [15].
6. Numerical Analysis
This section will demonstrate the different settings for system performance measures by using MATLAB. We investigate exponentially distributed retrial, service, slower pace service, vacation and repair periods. Numerical measurements are selected at random in order to fulfil the stability criteria. Tables 1 to 3 provides assessed outcomes of the idle prob., X0, $0 the "mean queue size (Lq), mean waiting time in the queue (Wq)" in our QM.
Table 1 shows that the retrial rate (n) escalates,X0 escalates, but Lq, Wq decreases for the value of w = 0.19, p = 0.9, A = 0.8, y = 3, HHv(y) = 0.9, n = 0.5.
Table 2 demonstrates that the vacation rate (y) mounts,$0 increases,Lq, Wq subsides for the value of w = 0.19, p = 1.5, A = 0.19, HHv(y) = 0.9, r1 = 0.9.
Table 3 clearly displays that feedback rate (w) mounts, ^0, Lq, Wq diminshes for the value of p = 0.9, A = 0.10, HHv(y) = 0.9, r1 = 0.19, y = 0.9.
Table 1: The impact of Retrial rate (n) on X0, Lq, Wq
Retrial rate (n) X0 Lq Wq
2.0 2.0718 0.0883 0.0982
2.5 2.2132 0.0885 0.0984
3.0 2.3172 0.0821 0.0912
3.5 2.3969 0.0728 0.0809
4.0 2.4599 0.0621 0.0690
4.5 2.5110 0.0507 0.0564
5.0 2.5532 0.0391 0.0434
The Figure 1 (a) indicates that retrial rate (n) escalates, (Lq) and (Wq) increases. The Figure 1 (b) displays that vacation rate (y) escalates, (Lq) and (Wq) decreases. The Figure 1 (c)
Table 2: The impact of Vacation rate (y) on <p0, Lq, Wq
Vacation rate (y) p0 Lq Wq
0.31 0.3724 0.6878 1.3756
0.32 0.3756 0.4931 0.9863
0.33 0.3785 0.3460 0.6921
0.34 0.3811 0.2890 0.4780
0.35 0.3835 0.1658 0.3316
0.36 0.3857 0.1211 0.2422
0.37 0.3878 0.1005 0.2010
Table 3: The impact of Feedback rate (w) on xo, Lq, Wq
Feedback rate (w) X0 Lq Wq
0.10 0.8306 0.2211 0.2457
0.20 0.8269 0.2008 0.2231
0.30 0.8233 0.1812 0.2013
0.40 0.8197 0.1623 0.1804
0.50 0.8161 0.1443 0.1603
0.60 0.8125 0.1270 0.1411
0.70 0.8090 0.1105 0.1227
(a) Retrial rate n Vs Lq, Wq
(b) Vacation rate 7 Vs Lq, Wq.
(c) Feedback rate w Vs Lq, Wq. Figure 1: Effects of a few parameters on 3D representation.
demonstrates that feedback rate (w) increases, (Lq) and (Wq) diminshes.
We may use the numerical findings above to determine the influence of features on the
system's assessment criteria with certainty, the outcomes correspond to actual circumstances.
6.1. ANFIS Computing
Using the fuzzy toolbox of MATLAB software, an ANFIS network can be executed to compare outcomes from analyses. ANFIS, a soft computing approach, is an effective tool for identifying important results that are useful in busy everyday environments. This approach aids in the identification of approximate solutions for measurements whose definite outcomes would otherwise be difficult to determine.
In our framework, a neuro-fuzzy technique is used to compute the expected no. of consumers in the queue (Wq) by changing the retrial rate (n), vacation rate (y) and feedback rate (w)), as shown in the numerical results in Figure 2 (a — c). We consider the parameters (n, y and w) as linguistic variables (LV) that are performed for four epochs each. The analytic (ANFIS) outcomes are exhibited by the solid (dashed) lines. In the context of fuzzy systems, these factors are regarded as LV and are used as input variables in ANFIS networks.
The Gaussian function provides the membership functions for each of these input variables. The following are the linguistic values for each parameter: low, average, high, and excessive. The diagrams demonstrate agreement between the analytical findings for the paradigm and the neuro fuzzy results achieved through the ANFIS approach.
(a) Retrial rate q Vs Wq
(b) vacation rate 7 Vs Wq
(c) Feedback rate w Vs Wq Figure 2: Effects of a few parameters on 2D representation.(ANFIS)
7. Cost Optimization
Our research aims to maintain system accessibility while optimizing system costs. Consequently, we establish the predicted cost function for system performance metrics and then accomplish a numerical analysis of the machining system under study. In order to calculate the best average cost per unit of time (TC), the parameters must be determined. This section discusses the best cost construct for the suggested approach using the standard cost notation form, and it provides the estimated total cost per unit of time as follows:
ETC = ShLs + Sb Y + Sv A + Sr n + S1Z + S2k (52)
where,
• Sh = Holding cost per unit consumer.
• Sb = Cost per unit time while the server provides service during a usual busy period.
• Sv = Cost per unit time in the system when the server is on vacation.
• Sr = Cost per unit time for providing repair to the failed server.
• S1 = Cost per unit time consumer served by the mean service rate Z.
• S2 = Cost per unit time consumer served by the mean vacation rate k.
Equation (52) has an estimated total cost function that is multivariate and nonlinear. Therefore, developing an analytical solution for optimal parameter values say Z* and k* is problematic. In order to determine the most suitable numerical value for the decision parameters, the well-recognised meta heuristic technique: The optimization approach used is called Particle Swarm Optimization (PSO).
Choosing at random the present values for the cost element and the parameters are : p = 0.05, B(p) = 4.5, X = 5.9, Hv(7) = 0.75, (Z* = 0.3670,k* = 0.3900).
Our goal is to identify the best values that will allow us to minimise the cost function. The five sets of cost factors that we have chosen are listed below.
Table 4: Cost Sets values for different cost aspects
Cost sets Sh Sb Sv Si S2 Sr
Set 1 $15 $75 $20 $19 $12 $10
Set 2 $20 $85 $25 $17 $19 $23
Set 3 $25 $95 $30 $15 $15 $29
Applying the PSO algorithm using MATLAB software to the previously specified cost factors. In this research, we have 100 candidates, 500 iterations throughout, and a range of parameters between 0.006 and 0.65 for the lower and upper bounds. Tables 5 demonstrates that the effects of p, w, y on TEC* using PSO.
Table 5: The PSO approach is executed by changing p, w and y to determine the minimal cost for different cost sets.
Parameters (TEC* )
Cost set 1 Cost set 2 Cost set 3
0.20 $75.1471 $59.7276 $66.7365
¥ 0.25 $115.1554 $87.6434 $92.1116
0.30 $151.7660 $114.6657 $119.3708
1.00 $73.5533 $58.6272 $65.7761
œ 1.15 $74.5054 $59.2872 $66.3535
1.20 $74.8255 $59.5074 $66.5452
1.25 $78.6575 $61.4278 $67.9418
Y 1.30 $76.0953 $60.0313 $66.8629
1.35 $73.5533 $58.6272 $65.7761
By changing a few of the variables, we were able to determine the entire system's cost, and we found that for p = 1.00, y = 1.35 and (Z* = 0.3670, k* = 0.3900) the lowest cost was $58.6272.
7.1. Particle Swarm Optimization
A precise evaluation of the QM is highly essential to offer adequate service and decrease congestion with the increasing expansion of computer networking and communications. If the greatest number of consumers can access an affordable system, then this is feasible. As a result, solutions that incorporate cost optimisation are very beneficial and advised. Jain et al. [25] and Jain and Meena [24] presented a cost investigation of a queueing model that includes an unreliable server and vacation periods. Utilizing the particle swarm optimization (PSO) approach, we have attempted to deal with the cost constraints in networking systems. This approach may sort through a very vast number of possible solutions to identify the most appropriate one. PSO has an additional benefit over a variety of optimisation strategies in that it does not require the objective function to be possible to differentiate.
Using this procedure, we first initiate a given population consisting of several particles or candidates. These candidates are then forced to travel inside the search space while adhering to the specified parameters and the goal function over their location and velocity. Every particle's fitness value is computed and to assess the values of the global best (gbest) and personal best (pbest) in more detail. The new gbest value is the particle whose pbest value is greater than gbest. This technique keeps on going until the predetermined number of iterations is reached. The algorithmic rule for PSO was first proposed by Kennedy and Eberhart [20]. The price optimization of a discrete-time RQ with SF utilizing this method has been examined by Upadhyaya [21]. Zhang et al.[22] examined set up cost and numerical answers for a single server recurrent model with state-dependent service using the PSO algorithmic approach. We have cited Malik et al.[23] investigation as it pertains to the operation of the PSO and GA algorithm.
Iteration (a) Convergence in PSO
(b) Expected Total Cost
Figure 3: (a)2D and 3D visualization of PSO optimization
7.2. Convergence in PSO
Convergence holds significant importance within meta-heuristic optimization algorithms, representing the gradual improvement of potential solutions towards an optimal or near-optimal solution. The convergence pattern of an algorithm reflects its efficacy in exploring the solution space adeptly and moving closer to the global optimum. The findings from these figures suggest that employing the concept of a working vacation enhances the system's stability and reliability. This is due to the consistent availability of the server during this period. When a machine fails, the server responds quickly to the problem, but it provides a slower rate of service than when the machine is operating normally. Within PSO, particles converge toward the most optimal solution they are aware of, coming together as their movements become restricted and the finest solution steadies. Fig 3(a) Demonstrates that PSO achieves convergence towards the optimal cost. Fig 3(b)Displays the convexity and optimality of the cost function concerning the cost sets utilized in the optimization analysis.
8. Conclusion
In this research, a single server feedback retrial queueing system with starting failures, Bernoulli working vacations, and vacation interruptions was investigated. The number of customers in the system and its orbit are used to find the PGFs. This is done by using the "supplementary variable method". The average orbital queue length and the system average queue length have precise expressions. Numerical examples are used to verify the analytical conclusions. The mean busy period as well as other significant system performance indicators are determined. Also,numerical outcomes are compared to ANFIS. We have demonstrated how to optimize the functioning of a real-world service system using the PSO meta-heuristic algorithm. This suggested paradigm may be used in communication networks, supermarkets, management and production industries, etc. Basically, it is nearly impossible to construct a paradigm in which the server never defects or deactivates in any of these enormous sectors. As a result, this analysis is pertinent and in favour of scenarios in which a server can remain idle to maximise the consumption of resources. This model's construction helps to prevent the regular overcrowding issues that networking and communication systems suffer. The suggested model may be expanded in the future to incorporate other factors, such as modified vacation policy, randomized policy, consumer impatience , priorities, and setup times.
• Acknowledgments : Not applicable.
• Competing interest :The authors declare that they have no competing interests.
• Funding :The authors declare that there is no funding source for this research.
• Author's contributions : All the authors made substantial contributions to the conception or design of the work.
[1] Gomez-Corral,A.(1999).Stochastic analysis of a single server retrial queue with general retrial times.Naval Res. Logist., Vol.46 (5):561-581.
[2] Falin,G.I. and Templeton, J.G.C.(1997).Retrial queues.London:Chapman and Hall,CRC Press,75.
[3] Artalejo, J., and Corral,A.G.(2008).Retrial Queueing Systems, Springer, Berlin, Germany.
[4] Artalejo,J.R.(2010).Accessible bibliography on retrial queues,progress in 2000-2009.Mathematical and computer modelling. 51:1071-81 .
[5] Gao, Wang, Li.(2014). An M/G/1 retrial queue with general retrial times, working vacations and vacation interruption.Asia-Paci J. Oper. Res.,.31 (2):1440006. https://doi.org/10.1142/S0217595914400065.
[6] Gao.S., and Liu,Z.(2013). An M/G/1 queue with single working vacation and vacation interruption under Bernoulli schedule, Applied Mathematical Modelling.37(3):1564-1579.
[7] Rajadurai,P., Saravanarajan,M.C. and Chandrasekaran,V.M.(2016). Analysis of an M/G/1 feedback retrial queue with unreliable server, non-persistent customers, single working vacation and vacation interruption.Int. J. Services and Operations Management. 24(2):235-266.
[8] Rajadurai,P., Chandrasekaran,V.M. and Saravanarajan,M.C.(2016). Analysis of an M[X1 /G/1 unreliable retrial G-queue with orbital search and feedback under Bernoulli vacation schedule.OPSEARCH.53(1):197-223.
[9] Rajadurai,P.,Saravanarajan,M.C. and Chandrasekaran,V.M.(2018). Analysis of an unreliable retrial G-queue with working vacations and vacation interruption under Bernoulli schedule.Ain Shams Engineering Journal.9(4):567-580. http://D0I.10.1016/jj.asej.2016.03.008.
[10] Mokaddis,G.S.,Metwally,S.A. and Zaki,B.M.(2007).A Feedback Retrial Queuing System with Starting Failures and Single Vacation.Tamkang Journal of Science and Engineering.10(3): 183-192.
Declarations
References
Krishna Kumar,B.,Pavai Madheswari, S. and Vijayakumar,A.(2002).The M/G/1 retrial queue with feedback and starting failures..Applied Mathematical Modelling26:1057-
1075.https://doi.org/10.1016/S0307-904X(02)00061-6.
P.Rajadurai.(2019).A study on M/G/1 preemptive priority retrial queue with Bernoulli working vacations and vacation interruption.Int. J. Process Management and Benchmarking.9 (2):193-215. https://doi.org/10.1504/IJPMB.2019.099331.
Rajadurai,P.,Saravanarajan,M. C. and Chandrasekaran, V. M.(2016).single server retrial queue with Bernoulli working vacation and vacation interruption.International Journal of Applied Engineering Research.11 (1).
Arivudainambi,D.,Godhandaraman,P. and Rajadurai,P.(2013).Performance analysis of a single server retrial queue with working vacation, OPSEARCH.51(3):434-462. https://doi.org/10.1007/s12597-013-0154-1.
Zhang,M., and Hou,Z.(2010).Performance analysis of M/G/1 queue with working vacations and vacation interruption.J. Comput. Appl. Math.234(10):2977-2985.https://doi.org/10.1016/j.cam.2010.04.010.
Zhang, M., and Hou,Z.(2012). M/G/1 queue with single working vacation.Journal of Appl. Math.Comput..39:221-234.
Rachita Sethi,Amita Bhagat and Deepika Garg.(2019). ANFIS based Machine Repair Model with Control Policies and Working Vacation. International Journal of Mathematical, Engineering and Management Sciences.4(6):1522.https://dx.doi.org/10.33889/IJMEMS.2019.4.6-120. Charu Bhargava and Madhu Jain.(2014). Unreliable multiserver queueing system with modified vacation policy, OPSEARCH, 2(51):159-182. DOI 10.1007/s12597-013-0138-1. Radhika Agarwal,Divya Agarwal,Shweta Upadhyaya and Izhar Ahmad.(2023).Üptimization of a stochastic model having erratic server with immediate or delayed repair,Annals of Operations Research.331(2):605-628. https://doi.org/10.1007/s10479-022-04804-2. Kennedy, J.,and Eberhart, R.(1995). Particle swarm optimization. In: Proceedings of the IEEE International Conference on Neural Networks (Perth, Australia): pp. 1942-1948. Upadhyaya, S.(2020). Cost optimization of a discrete-time retrial queue with Bernoulli feedback and starting failure. Int. J. Ind. Syst. Eng.36:165-196.
Zhou, M.,Liu, L.,Chai,X. and Wang,Z.(2018).Equilibrium strategies in a constant retrial queue with setup time and the N-policy. Commun. Stat. Theory Methods.49(7):1695-1711. Malik, G.,Upadhyaya,S.and Sharma,S.(2021). Cost inspection of a Geo/G/1 retrial model using particle swarm optimization and Genetic algorithm. Ain Shams Eng. J. 12(2):2241-2254. Jain, M., and Meena, R.K.(2020).Availability analysis and cost optimization of M/G/1 fault-tolerant machining system with imperfect fault coverage. Arabian Journal for Science and Engineering. 45(3):2281-2295.
Jain, M.,Kumar,P.and Meena,R.K.(2020).Fuzzy metrics and cost optimization of a fault tolerant system with vacationing and unreliable server. Journal of Ambient Intelligence and Humanized Computing. https://doi.org/10.1007/s12652-020-01951-x.
Jain, M.,and Kumar,A.(2023). Unreliable Server M[X1 /G/1 Retrial Feedback Queue with Balking, Working Vacation and Vacation Interruption. Proc. Natl. Acad. Sci., India, Sect. A Phys. Sci. 93: 57-73. https://doi.org/10.1007/s40010-022-00777-w.
Pazhani Bala Murugan,S., and Keerthana,R.(2023). An M/G/1 Feedback retrial queue with working vacation and a waiting server. Journal of computational analysis and applications.31(1). Keerthiga, S.,and Indhira, K.(2023).Two phase of service in M/G/1 queueing system with retrial customers. Journal of Analysis .https://doi.org/10.1007/s41478-023-00635-x. Agarwal, D., Agarwal,R. and Upadhyaya,S.(2024).Detection of optimal working vacation service rate for retrial priority G-queue with immediate Bernoulli feedback. Results in Control and Optimization.14:100397.https://doi.org/10.1016/j.rico.2024.100397