Next Article in Journal
Norm-Resolvent Convergence for Neumann Laplacians on Manifold Thinning to Graphs
Previous Article in Journal
Automatic Evaluation Method for Functional Movement Screening Based on a Dual-Stream Network and Feature Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation

1
Department of Finance, Chaoyang University of Technology, 168, Jifeng East Road, Wufeng District, Taichung City 41349, Taiwan
2
Ph.D. Program of Business Administration in Industrial Development, Department of Business Administration, Chaoyang University of Technology, 168, Jifeng East Road, Wufeng District, Taichung City 41349, Taiwan
*
Authors to whom correspondence should be addressed.
Mathematics 2024, 12(8), 1163; https://doi.org/10.3390/math12081163
Submission received: 18 March 2024 / Revised: 6 April 2024 / Accepted: 11 April 2024 / Published: 12 April 2024

Abstract

:
This work investigates a two-way communication retrial queue with synchronous working vacation and a constant retrial policy. During the idle time, a server makes an outgoing call after a random length. The service time of the incoming call and outgoing call obeys exponential distribution with different rates. If the incoming call finds all servers to be unavailable, it may or may not enter orbit. All servers immediately go on vacation simultaneously as soon as they find an empty system after the service finishes. During vacation, the servers can provide a service to those incoming calls, but this is at a lower-speed rate. The stationary probability distribution and the ergodic condition are obtained utilizing the matrix geometric technique. Some system characteristics are developed. Using MATLAB software, the variation in average orbit length, idle ratio, and the average number of servers in different server states is plotted for different values of the incoming/outgoing call rate and retrial rate. We further propose a multi-objective optimization model from which the optimal rate of outgoing calls and optimal vacation rate are explicitly obtained.

1. Introduction

This paper deals with investigations into retrial queue with two-way communication and synchronous working vacation. For more information about retrial queues, please refer to [1,2,3,4,5]. The main reason for this is that retrial queues provide a suitable model for the performance evaluation of computer networks, communications systems, and call centers. The characteristic of the retrial queue is that an arriving customer who cannot be serviced immediately joins an orbit and tries the request again after a period of time. Therefore, the analysis of a retrial queue is more difficult than that of its counterpart model without retrial.
In most of the literature on retrial queue, servers only serve incoming arrivals. Servers wait for the next arrival (from outside or orbit) after completing service. However, in some real-world situations, there is an opportunity for servers to make an outgoing phone call at their idle state. This is especially true for a business organization, such as a call center, where agents not only handle incoming calls but can also make calls to sell, advertise, and promote the business’s products and services. A characteristic of two-way communication is the fact that unoccupied servers can perform outgoing calls to the source. In such systems, the server’s utilization is always a critical issue; for an example, see [6,7,8,9]. Ayyappan and Gowthami [10] performed a stationary analysis of a feedback retrial queue with impatient customers, vacations, and two types of arrivals. Lee et al. [11] analyzed the waiting time distribution of a two-way communication retrial queue. Sztrik et al. [12] examined the unreliable operation of a finite-source two-way communication retrial queue and compared the different distribution of failure time on performance measures. Tóth and Sztrik [13] studied a finite-source retrial queue with two-way communication.
Since the vacation queue reflects the situation of servers utilizing this idle time for different purposes like doing supplementary jobs, equipment maintenance, testing, and so on, it has been extensively studied for modeling and analyzing practical problems, including communication networks, call centers, service centers, production lines, manufacturing systems, and more. A special class of vacation queues includes working vacation, in which the server provides a service at a lower rate rather than stopping working completely. A work on retrial queue with working vacation and starting failure was conducted by Yang and Wu [14]. Li et al. [15] analyzed the retrial queue with Bernoulli working vacation interruption. A retrial queue with general retrial times and single working vacation was studied by Li et al. [16] by utilizing the supplementary variable method. Gupta and Kumar [17] studied a retrial queue with a single waiting server subject to breakdown and repair under working vacation and interruption. In the work of Muthusamy et al. [18], they considered a working vacation retrial queue with three different classes of customers and optional reservice during working vacations in the Bernoulli schedule. Pazhani et al. [19] used the supplementary variable technique to obtain the probability generating function for the number of customers and the mean number of customers in the invisible waiting area of a retrial queue with working vacation and a single waiting server. Shanmugam and Saravanarajan [20] investigated an unreliable retrial queueing system with working vacation. The orbit and system lengths were derived through the supplementary variable method. Chen et al. [21] studied a single server retrial queue incorporating random working vacation and improved service efficiency during vacation policies and examined its optimal queueing strategies. Sundararaman et al. [22] examined a new type of queue in two waiting queues (original queue and orbit queue) with working vacations.
In these articles mentioned above, the authors mostly focus on a single server. However, real systems such as telecommunications or call centers are often multiserver rather than single servers. There is limited research on multiserver two-way communication retrial queues in the literature because the analysis of multiserver two-way communication retrial queues is more complex than single-server queues. In this paper, we consider a multiserver retrial queue with two-way communication with working vacation. Our advantages and contributions are as follows. We carry out an extensive analysis for a multiserver two-way communication retrial queue with synchronous working vacation. With the support of the matrix geometric technique, the steady-state probabilities and ergodic condition are derived. The effect of the parameters on the system characteristics is displayed numerically. We also construct a multi-objective optimization analysis, from which the optimal rate of outgoing calls and optimal vacation rate are explicitly obtained.
The paper structure is as follows. The description of the model and the balance equations governing the system’s behavior are given in Section 2. A detailed system analysis, including the expressions for the stationary probability distribution and the ergodicity condition, is presented in Section 3. In Section 4, the system characteristics are demonstrated. Section 5 provides some numerical illustrations followed by a concluding remark.

2. System Description and Mathematical Model

This section lists the assumptions made in this paper. This allows for the formulation of a set of balance equations that gives a starting point for the system analysis.

2.1. System Description

We consider a multiserver retrial queue in which incoming calls follow a Poisson arrival process with rate λ. The service times of an incoming call are exponentially distributed with rate μ1. If an incoming call finds that the server is fully occupied, it may or may not enter orbit, and if it does, it retries to seek service after an exponentially distributed time at rate σ. Within the orbit, we apply a constant retrial policy, i.e., only the call at the head of orbit can request service. Suppose that an incoming call goes into orbit with probability b. Otherwise, the incoming call starts service immediately. During server’s idle time, it makes an outgoing call with rate α and provides the service for an exponentially distributed time with rate μ2. We denote s as the number of servers in the system. All the servers take the vacation simultaneously when they find the system empty, and vacation duration is also exponentially distributed with parameter θ. During the vacation period, the servers can serve incoming calls at a low service rate μ3, but they will not make outgoing calls. The work flow of the two-way communication retrial queue with synchronous working vacation is described in Figure 1.

2.2. Markov Chain and Balance Equations

We denote the number of servers busy serving incoming and outgoing calls at time t by n 1 t and n 2 t , respectively. Also, C t and M t are the server state and the number of customers in orbit at time t, respectively, where
C t = 0 ,   all   servers   are   on   vacation   at   time   t , 1 ,   no   servers   are   on   vacation   at   time   t .
Thus, the process n 1 t , n 2 t , C t , M t ; t 0 is a continuous time Markov chain on the state space Ω = i , j , k , n : 0 i s , j = 0 , k = 0 , n 0 i , j , k , n : 0 i + j s , k = 1 , n 0 .
Let π i , j , k n = lim t P n 1 t ,   n 2 t ,   C t ,   M t be the limiting probabilities for the steady-state distribution. The system of governing steady-state equations is framed as follows:
λ + θ π 0 , 0 , 0 0 = μ 3 π 1 , 0 , 0 0 + μ 1 π 1 , 0 , 1 0 + μ 2 π 0 , 1 , 1 0 ,
λ + s α π 0 , 0 , 1 0 = θ π 0 , 0 , 0 0 ,
λ + θ + σ π i , 0 , 1 n = μ 3 π 1 , 0 , 0 n ,   n 1 ,
λ + s α + σ π i , 0 , 1 n = μ 1 π 1 , 0 , 1 n + μ 2 π 0 , 1 , 1 n + θ π 0 , 0 , 0 n , n 1 ,
λ + i μ 3 + θ + 1 δ n , 0 σ π i , 0 , 0 n = λ π i 1 , 0 , 0 n + i + 1 μ 3 π i + 1 , 0 , 0 n + σ π i 1 , 0 , 0 n + 1 , 1 i s 1 ,   n 0 ,
λ + i μ 1 + s i α + 1 δ n , 0 σ π i , 0 , 1 n = λ π i 1 , 0 , 1 n + i + 1 μ 1 π i + 1 , 0 , 1 n + μ 2 π i , 1 , 1 n + θ π i , 0 , 0 n + σ π i 1 , 0 , 1 n + 1 , 1 i s 1 ,   n 0 ,
λ + i μ 1 + j μ 2 + s i j α + 1 δ n , 0 σ π i , j , 1 n = ( s i j + 1 ) α π i , j 1 , 1 n + λ π i 1 , j , 1 n + i + 1 μ 1 π i + 1 , j , 1 n + j + 1 μ 2 π i , j + 1 , 1 n + σ π i 1 , j , 1 n + 1 , 0 i s j 1 ,   1 j s 1 ,   n 0 ,
b λ + s μ 3 + θ π s , 0 , 0 n = b λ π s , 0 , 0 n 1 + λ π s 1 , 0 , 1 n + σ π s 1 , 0 , 0 n + 1 ,   n 0 ,
b λ + s j μ 1 + j μ 2 π s j , j , 1 n = b λ π s j , j , 1 n 1 + α π s j , j 1 , 1 n + λ π s j 1 , j , 1 n + θ δ j , 0 π s j , j , 0 n + σ π s j 1 , j , 1 n + 1 ,   0 j s ,   n 0 .
The normalizing condition is
n = 0 i = 0 s π i , 0 , 0 n + n = 0 i = 0 s j = 0 s i π i , j , 1 n = 1 ,
where δ n , 0 represents Kronecker’s delta and π i , j , k n = 0 if i ,   j ,   k ,   n Ω .

3. System Analysis

This section provides a system analysis. First, a quasi-birth-and-death process formulation is provided to identify all the components of the involved infinitesimal generator and block matrices. Next, the steady-state probabilities are derived in matrix form, and ergodic condition is exported.

3.1. Infinitesimal Generator and Matrices

According to the matrix geometric method, the infinitesimal generator of the continuous time Markov process n 1 t ,   n 2 t ,   C t ,   M t ; t 0 } is written in the form of a block matrix Q as
Q = G 0 G u G l G G u ,
where G0, Gu, Gl, and G are square matrices of order (s + 1)(s + 4)/2.
For the concerned model, we can represent the block structures of these matrices as
G u = G 0 u G 1 u G s u ,
G l = G 0 l G 1 l G s l ,
G = B 0 B 1 A 0 A 0 + A 1 A 1 A 1 +                                             A s 1 A s 1 A s 1 + A s A S ,
where G j u   0 j s , G j l   0 j s , B 0 , B 1 , A j   0 j s , A j +   0 j s 1 , and A j   1 j s are (sj + 1) × (sj + 1), (sj + 1) × (sj + 1), (s + 1) × (s + 1), (s + 1) × (s + 1), (sj + 1) × (sj + 1), (sj + 1) × (sj), and (sj + 1) × (sj + 2) matrices, respectively, with entries given by
G j u k , k = b λ ,   k = k = s j + 1 , 0 ,   in   other   cases ;
G j l k , k = σ ,   k = k + 1 ,   1 k s j , 0 ,   in   other   cases ;                                        
B 0 k , k = λ ,                                             k = k + 1 ,   1 k s ,                                   k 1 μ 3 ,                     k = k 1 ,   2 k s + 1 ,                           λ + k 1 μ 3 + θ + σ ,   k = k ,   1 k s ,     b λ + s μ 3 + θ ,         k = k = s + 1 ,                                             0 ,                                               in   other   cases ;                                                      
B 1 k , k = θ ,   k = k ,   1 k s + 1 , 0 ,   in   other   cases ;                          
A j + k , k = s j + 1 k α ,   k = k ,   1 k s j , 0 ,                                                   in   other   cases ;                          
A j k , k =       λ ,                                                                                                         k = k + 1 ,   1 k s j ,                                   k μ 1 ,                                                                                                   k = k 1 ,   2 k s j + 1 ,             λ + k 1 μ 1 + j μ 2 + s j + 1 k α + σ ,   k = k ,   1 k s , b λ + s j μ 1 + j μ 2 ,                                                                             k = k = s + 1 ,                             0 ,                                                                                                             in   other   cases ;                                                                  
A j k , k = j μ 2 ,   k = k ,   1 k s j + 1 ,         0 ,   in   other   cases .                          
The matrix G 0 is the same as G, but the element of G 0 s + 3 ,   s + 2 is shifted to the position of G 0 s + 3 ,   1 , the element of G 0 2 s + 3 ,   s + 2 is shifted to the position of G 0 2 s + 3 , 1 , and the retrial rate σ is ignored.

3.2. Stationary Distribution

Assume that Π = Π 0 , Π 1 , represents the steady-state probability vector of Q. Also, Π n = π 0 , 0 , 0 n π s , 0 , 0 n π 0 n π s n and π j n = π 0 , j , 1 n π s j , j , 1 n for 0 j s .
We can rewrite Equations (1)–(10) in the matrix form as ΠQ = 0 and n = 0 Π n e = 1 , where 0 and e represent a row and a column vector with an appropriate size with all zero and all one entries, respectively.
Note that from Neuts [23], there is a matrix R, such that Π n = Π n 1 R . This implies that Π n = Π 0 R n . Here, R is the minimal non-negative solution of
G u + R G + R 2 G l = 0 ,
with a spectral radius less than one. Since it is difficult for R to obtain the explicit expression by solving Equation (12), researchers could use different methods such as a successive substitution approach to evaluate R. In the numerical illustration, the following procedure will be implemented to approximate R.
We can define I as an identity matrix. The first equation of ΠQ = 0 can be written as the following by replacing Π 1 with Π 1 = Π 0 R :
Π 0 G 0 + Π 0 R G l = 0 .
And n = 0 Π n e = 1 can be re-written as
n = 0 Π n e = Π 0 I R 1 e = 1 ,
where I is an identity matrix. These two equations can be used to determine Π(0).

3.3. Ergodicity

Consider that the matrix H = G u + G + G l =. Let p = p 0 , 0 p 0 , 1 p s , 1 is the invariant probability vector of H, where p 0 , 0 = p 0 , 0 , 0 p 1 , 0 , 0 p s , 0 , 0 and p j , 1 = p 0 , j , 1 p 1 , j , 1 p s j , j , 1 0 ≤ j ≤ s. Therefore, p satisfies the following criteria: pH = 0 and pe = 1. The condition for the ergodicity of the Markov chain n 1 t , n 2 t , C t ; t 0 is p G u e < p G l e . The ergodicity condition can be transformed to
b λ + σ P B l o s s < σ ,
where P B l o s s = j = 0 s p s j , j , 1 denotes the blocking probability of the corresponding loss system.
The ergodicity condition proposed by p G u e < p G l e itself is unambiguous because the number of states for n 1 t ,   n 2 t ,   C t ; t 0 is finite. However, it does not seem to be easy to obtain a simple scalar form in terms of the given parameters. Below, a simple explicitness is obtained from (15) for the special case of s = 1. For the matrices, we have
G u = 0 0 0 0 b λ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 b λ 0 0 b λ ,
G l = 0 σ 0 0 0 0 0 0 0 0 0 0 0 σ 0 0 0 0 0 0 0 0 0 0 0 ,
G = λ + θ + σ λ μ 3 b λ + μ 3 + θ θ                                               0 0                                               θ 0 0 0                                           0 0                                           0 λ + α + σ λ μ 1 b λ + μ 1 α 0 0                                           0 μ 2                                           0 b λ + μ 2 ,
H = λ + θ + σ λ + σ θ 0 0 μ 3 μ 3 + θ 0 θ 0 0 0 λ + α + σ λ + σ α 0 0 μ 1 μ 1 0 0 0 μ 2 0 μ 2 .
We have p 0 , 0 , 0 = 0 , p 1 , 0 , 0 = 0 , p 0 , 1 , 1 = λ + σ p 0 , 0 , 1 / μ 1 , and p 1 , 0 , 1 = α p 0 , 0 , 1 / μ 2 . Hence, the ergodicity condition is
b λ λ + σ μ 1 + α μ 2 < σ .
If we set b = 1, then the ergodicity condition is consistent with Equation (15) in Phung-Duc and Rogiest [24].

4. System Characteristics and Cost Function

We compute system characteristics and average operating cost in terms of the steady-state probabilities in this section.

4.1. System Characteristics

If the system is in stability condition, we can obtain the expressions of the system characteristics listed below.
The average number of calls in orbit (average orbit length) is
E N = n = 0 n i = 0 s π i , 0 , 0 n + n = 0 n j = 0 s i = 0 s j π i , j , 1 n .
By Little’s Law, the average waiting time of a call in orbit is
E W = E N b λ = 1 b λ n = 0 n i = 0 s π i , 0 , 0 n + n = 0 n j = 0 s i = 0 s j π i , j , 1 n .
The average number of incoming calls during the vacation period is
E S 0 = n = 0 i = 0 s i × π i , 0 , 0 n .
The average number of incoming calls during the non-vacation period is
E S 1 = n = 0 j = 0 s i = 0 s j i × π i , j , 1 n .
The average number of outgoing calls during the non-vacation period is
E S 2 = n = 0 j = 0 s j i = 0 s j π i , j , 1 n .
The average number of servers being busy during the non-vacation period is defined as
E B = n = 0 j = 0 s i = 0 s j i + j π i , j , 1 n = E S 1 + E S 2 .
The average number of servers on vacation, whether busy or not, is calculated as
E V = n = 0 s i = 0 s π i , 0 , 0 n .
The average number of idle servers is calculated as
E I = s E B E V .
The server utilization at the steady state is
U = E S 0 + E S 1 + E S 2 s .

4.2. Cost Function

Assume that the structure of the average operating cost is linear cost structure based on the cost elements associated with different system characteristics. First, define the cost elements involved in the average operating cost as follows:
r1 ≡ the cost of a retrial customer;
r2 ≡ the cost of an idle server during the non-vacation period;
r3 ≡ the cost of a busy server during the non-vacation period;
r4 ≡ the cost of a vacation server;
r5 ≡ the cost of providing a vacation.
Then, the average operating cost is
O C = r 1 E N + r 2 E I + r 3 E B + r 4 E V + r 5 θ .

5. Numerical Illustration

If the ergodicity condition is satisfied, this section presents some numerical results on the system characteristics and optimization results. We perform sensitivity analysis to illustrate the effect of system descriptors on the system characteristics.

5.1. Sensitivity Analysis

For the sensitivity analysis, the default parameters are fixed as s = 5, λ = 4, α = 3, μ1 = 1.5, μ2 = 1.5, μ3 = 0.8, σ = 10, θ = 0.3, and b = 0.8, unless their values are mentioned in tables and figures.
Table 1 and Figure 2 clearly show that the average orbit length increases whenever the incoming call rate increases. At the same time, the server idle time is reduced, and the average number of vacation servers decreases.
Table 2 and Figure 3 demonstrate that the increasing nature of the average orbit length first decreases and then increases as the outgoing call rate increases. For the larger outgoing call rate, if the server makes outgoing calls more frequently, it will have fewer opportunities to service incoming calls, resulting in a longer orbit length. On the other hand, the smaller the outgoing call rate, the longer the server takes to make outgoing calls, which results in the fact that the longer the server is idle, the more incoming calls it can serve. This threshold is related to the arrival rate of incoming calls. And the average number of vacation servers decreases. Simultaneously, the server idle time is reduced. Table 3 and Figure 4 display that the server idle time decreases slightly whenever the retrial rate increases. At the same time, the average number of vacation servers increases, and the average orbit length of incoming calls decreases.
Table 4 and Figure 5 clearly show that the non-balking rate has an increasing trend for the average orbit length, while the non-balking rate has an increasing tendency for the server idle time and the average number of vacation servers. This is reasonable since when the non-balking rate is increasing, this means that there are more incoming calls joining the orbit.

5.2. Optimization

The aim of optimization is to determine the optimal outgoing calls rate and optimal vacation rate. Except for operating cost consideration, average orbit length and server utilization are also important system characteristics in the modeling of any queue. Therefore, we will design a tri-objective optimization problem as follows:
Minimize: y1(α, θ) = OC
Minimize: y2(α, θ) = 1-U
Minimize: y3(α, θ) = E[W]
Subject to: p G u e < p G l e ,
α > 0,
θ > 0.
As can be seen from the expressions of the average operating cost and loss function, the analytical solution is difficult to derive. Computational software is applied to solve the optimization problem numerically. Here, the computer software MATLAB R2022a is used to implement a multi-objective genetic algorithm to solve the above tri-objective optimization problem. Readers can refer to Konak et al. [25] for a detailed overall flow of the general multi-objective genetic algorithm. The multi-objective genetic algorithm pseudocode is shown below (Algorithm 1).
Algorithm 1. Multi-Objective Genetic Algorithm
Begin
  • Initialize population
  • Evaluate fitness of each individual in the population
While not termination criteria do
  • Select the best individuals to reproduce
  • Apply crossover and mutation operations
  • Evaluate population
  • Update population
End While
  • Output the best solution
End
The optimization analysis has been performed for the default parameters, fixed as follows:
s = 5, λ = 4, μ1 = 1.5, μ2 = 1.5, μ3 = 0.8, σ = 10, b = 0.8, r1 = 20, r2 = 25, r3 = 100, r4 = 50, r5 = 5.
Figure 6 reveals that the higher the average operating cost, the lower the expected loss is. Table 5 presents the Pareto optimal solution sets. At this time, the operation manager needs to decide which set of solutions to implement. There is no right or wrong choice; it depends on the decision maker’s own experience and judgment. To help the decision maker make decisions based on the Pareto optimal solution sets, Solution #3 and Solution #4 are selected from Table 5 as reference schemes. Solution #3 is the solution that achieves the minimum average waiting time in orbit. This is a good solution if the decision maker cares about service quality. Solution #4 is chosen based on the assigned weight of these three objectives (0.5, 0.2, 0.3). The steps are as follows:
-
Step 1: Normalize the three objective values for each solution.
-
Step 2: Multiply their respective weights and add these three values.
-
Step 3: Sort the values obtained in step 2, and the first solution is the optimal solution.

6. Conclusions

This work investigated a synchronous working vacation model of a retrial queue, in which there are two types of arrivals, incoming calls made by customers and outgoing calls made by idle servers. The servers have no information about how many calls exist in orbit. So even with many calls in orbit, outgoing call activity is still possible. We utilized the matrix geometric technique to derive the stationary probabilities, ergodic condition, and system characteristics. At last, we implemented numerical examples to observe the trends of system characteristics for different parametric values and proposed a related optimization issue. The sensitivity analysis in this study can provide useful insights to decision makers to improve service quality. Through numerical examples, it is found that the server idle time is reduced with an increase in the incoming/outgoing call rate and retrial rate. The orbit length is found to decrease with an increase in the incoming/outgoing call rate and retrial rate. The average number of vacation servers decreases with an increase in the incoming call rate and outgoing call rate. The graph behavior was found to be consistent with theoretical expectations. Via a multi-objective genetic algorithm, the entire set of Pareto optimal solutions is determined. We also provide some reference solutions to illustrate how the decision maker can make decisions based on the Pareto optimal solution set.
One limitation of this paper is to assume that all the involved random variables in the model introductions are exponentially distributed. Extending our analysis to consider the general distribution is a future research direction. Additionally, similar models incorporating practical concepts such as working breakdown, batch arrival, immediate feedback, optional type service, negative customers, priority service, and many others can be studied.

Author Contributions

Conceptualization, T.-H.L., K.-C.C., C.-M.C. and F.-M.C.; methodology, T.-H.L. and K.-C.C.; software, T.-H.L. and C.-M.C.; validation, T.-H.L., K.-C.C. and F.-M.C.; writing—original draft preparation, T.-H.L. and C.-M.C.; writing—review and editing, K.-C.C.; supervision, F.-M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by National Science of Technology Council under grants 112-2221-E-324-016-MY2.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kim, J.; Kim, B. A survey of retrial queueing systems. Ann. Oper. Res. 2016, 247, 3–36. [Google Scholar] [CrossRef]
  2. Phung-Duc, T. Retrial Queueing Models: A Survey on Theory and Applications. arXiv 2019, arXiv:1906.09560. [Google Scholar]
  3. Ke, J.C.; Chang, F.M.; Liu, T.H. M/M/c balking retrial queue with vacation. Qual. Technol. Quant. Manag. 2019, 16, 54–66. [Google Scholar] [CrossRef]
  4. Zhang, Y. Strategic behavior in the constant retrial queue with a single vacation. RAIRO Oper. Res. 2020, 54, 569–583. [Google Scholar] [CrossRef]
  5. Liu, T.H.; Chang, F.M.; Ke, J.C.; Sheu, S.H. Optimization of retrial queue with unreliable servers subject to imperfect coverage and reboot delay. Qual. Technol. Quant. Manag. 2022, 19, 428–453. [Google Scholar] [CrossRef]
  6. Govindan, A.; Jayaraj, U. Analysis of mixed priority retrial queueing system with two way communication and working breakdown. J. Math. Model. 2018, 6, 195–212. [Google Scholar]
  7. Nazarov, A.; Phung-Duc, T.; Paul, S. Slow retrial asymptotics for a single server queue with two-way communication and Markov modulated Poisson Input. J. Syst. Sci. Syst. Eng. 2019, 28, 181–193. [Google Scholar] [CrossRef]
  8. Ayyappan, G.; Udayageetha, J.; Somasundaram, B. Analysis of non-pre-emptive priority retrial queueing system with two-way communication, Bernoulli vacation, collisions, working breakdown, immediate feedback and reneging. Int. J. Math. Oper. Res. 2020, 16, 480–498. [Google Scholar] [CrossRef]
  9. Kumar, M.S.; Dadlani, A.; Kim, K. Performance analysis of an unreliable M/G/1 retrial queue with two-way communication. Oper. Res. 2020, 20, 2267–2280. [Google Scholar]
  10. Ayyappan, G.; Gowthami, R. Analysis of MAP, PH2OA/PH1I, PH2O/1 retrial queue with vacation, feedback, two-way communication and impatient customers. Soft Comput. 2021, 25, 9811–9838. [Google Scholar] [CrossRef]
  11. Lee, S.W.; Kim, B.; Kim, J. Analysis of the waiting time distribution in M/G/1 retrial queues with two way communication. Ann. Oper. Res. 2022, 310, 505–518. [Google Scholar] [CrossRef]
  12. Sztrik, J.; Tóth, Á.; Pintér, Á. The Effect of Operation Time of the Server on the Performance of Finite-Source Retrial Queues with Two-Way Communications to the Orbit. J. Math. Sci. 2022, 267, 196–204. [Google Scholar] [CrossRef]
  13. Tóth, Á.; Sztrik, J. Analysis of retrial queueing systems with two-way communication and impatient customers using simulation. Ann. Math. Et Inform. 2023, 58, 160–169. [Google Scholar] [CrossRef]
  14. Yang, D.Y.; Wu, C.H. Performance analysis and optimization of a retrial queue with working vacations and starting failures. Math. Comput. Model. Dyn. Syst. 2019, 25, 463–481. [Google Scholar] [CrossRef]
  15. Li, T.; Zhang, L.; Gao, S. An M/G/1 retrial queue with balking customers and Bernoulli working vacation interruption. Qual. Technol. Quant. Manag. 2019, 16, 511–530. [Google Scholar] [CrossRef]
  16. Li, T.; Zhang, L.; Gao, S. An M/G/1 retrial queue with single working vacation under Bernoulli schedule. RAIRO-Oper. Res. 2020, 54, 471–488. [Google Scholar] [CrossRef]
  17. Gupta, P.; Kumar, N. Performance analysis of retrial queueing Model with working vacation, interruption, waiting server, breakdown and repair. Int. J. Sci. Res. 2021, 13, 833–844. [Google Scholar] [CrossRef]
  18. Muthusamy, S.; Devadoss, N.; Ammar, S.I. Reliability and optimization measures of retrial queue with different classes of customers under a working vacation schedule. Discret. Dyn. Nat. Soc. 2022, 2022, 6806104. [Google Scholar] [CrossRef]
  19. Pazhani, S.; Murugan, B.; Keerthana, R. An M/G/1 feedback retrial queue with working vacation and a waiting server. J. Comput. 2023, 31, 61–79. [Google Scholar]
  20. Shanmugam, B.; Saravanarajan, M.C. Unreliable retrial queueing system with working vacation. AIMS Math. 2023, 8, 24196–24224. [Google Scholar] [CrossRef]
  21. Chen, Z.; Xu, H.; Huo, H. Optimal queuing strategies for an M/G/1 retrial queue system with RWV and ISEV policies. ANZIAM J. 2024, 1–27. [Google Scholar] [CrossRef]
  22. Sundararaman, M.; Narasimhan, D.; Rajadurai, P. Performance Analysis of Two Different Types of Waiting Queues with Working Vacations. J. Math. 2024, 2024, 5829171. [Google Scholar] [CrossRef]
  23. Neut, M.F. Matrix-Geometric Solutions in Stochastic Models. Master’s Thesis, John Hopkins University, Baltimore, MD, USA, 1981. [Google Scholar]
  24. Phung-Duc, T.; Rogiest, W. Two Way Communication Retrial Queues with Balanced Call Blending. In Proceedings of the Analytical and Stochastic Modeling Techniques and Applications. ASMTA 2012, Grenoble, France, 4–6 June 2012. [Google Scholar]
  25. Konak, A.; Coit, D.W.; Smith, A.E. Multi-objective optimization using genetic algorithms: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 992–1007. [Google Scholar] [CrossRef]
Figure 1. Two-way communication retrial queue with synchronous working vacation.
Figure 1. Two-way communication retrial queue with synchronous working vacation.
Mathematics 12 01163 g001
Figure 2. Average orbit length and idle ratio vs. λ.
Figure 2. Average orbit length and idle ratio vs. λ.
Mathematics 12 01163 g002
Figure 3. Average orbit length and idle ratio vs. α.
Figure 3. Average orbit length and idle ratio vs. α.
Mathematics 12 01163 g003
Figure 4. Average orbit length and idle ratio vs. σ.
Figure 4. Average orbit length and idle ratio vs. σ.
Mathematics 12 01163 g004
Figure 5. Average orbit length and idle ratio vs. b.
Figure 5. Average orbit length and idle ratio vs. b.
Mathematics 12 01163 g005
Figure 6. Pareto front solutions found by multi-objective genetic algorithm.
Figure 6. Pareto front solutions found by multi-objective genetic algorithm.
Mathematics 12 01163 g006
Table 1. Effect of incoming call rate.
Table 1. Effect of incoming call rate.
λE[N]1-UE[B]E[V]
30.2240.2243.7200.308
3.20.4940.2143.7770.283
3.40.5880.2043.8320.259
3.60.7000.1953.8870.236
3.80.8330.1863.9400.214
40.9910.1773.9930.193
4.21.1820.1684.0450.173
4.41.4150.1604.0970.153
4.61.7040.1524.1470.135
4.82.0700.1444.1970.117
52.5460.1364.2460.100
Table 2. Effect of outgoing call rate.
Table 2. Effect of outgoing call rate.
αE[N]1-UE[B]E[V]
00.6590.3771.8861.947
0.20.6350.3612.1261.691
0.40.6190.3442.3591.458
0.60.6130.3272.5771.249
0.80.6160.3102.7791.065
10.6260.2932.9620.907
1.20.6430.2773.1270.771
1.40.6670.2623.2740.656
1.60.6950.2483.4050.559
1.80.7280.2363.5210.477
20.7640.2243.6240.407
Table 3. Effect of retrial rate.
Table 3. Effect of retrial rate.
σE[N]1-UE[B]E[V]
38.6020.1734.1070.049
43.3140.1744.0650.102
52.1700.1754.0410.134
61.6690.1764.0250.154
71.3890.1764.0130.168
81.2090.1774.0050.179
91.0830.1773.9980.187
100.9910.1773.9930.193
Table 4. Effect of non-balking rate.
Table 4. Effect of non-balking rate.
bE[N]1-UE[B]E[V]
0.00.0000.2253.6970.302
0.10.0450.2213.7250.292
0.20.1010.2163.7550.281
0.30.1720.2103.7870.270
0.40.2620.2053.8220.257
0.50.3760.1983.8600.243
0.60.5260.1923.9000.228
0.70.7230.1853.9450.211
0.80.9910.1773.9930.193
0.91.3650.1694.0460.173
1.01.9080.1604.1040.151
Table 5. Optimal solutions for the multi-objective optimization model.
Table 5. Optimal solutions for the multi-objective optimization model.
#α*θ*y2(α*, θ*)y3(α*, θ*)y1(α*, θ*)#α*θ*y2(α*, θ*)y3(α*, θ*)y1(α*, θ*)
10.00290.17190.33030.3469328.0764211.82146.18860.23020.2028458.0735
211.99839.36150.06075.6859887.95842211.96587.92050.06085.5383871.2552
30.00376.90890.47700.0577362.55942310.44857.65450.06822.3350662.1581
40.02250.10940.29090.5078329.48542411.27088.79650.06393.4762742.4771
511.69418.16590.06204.5257807.2304250.00344.61100.47530.0581352.2698
69.10948.12640.07641.4475604.6404260.85237.84450.31770.1256429.1866
79.80668.64420.07181.8235632.98242711.39639.29690.06343.7387762.0079
81.33397.33810.26700.1639447.73122811.63047.60430.06234.3340792.0429
98.35307.75550.08191.1600582.2932292.79205.00380.18080.2834475.5974
101.41715.01610.26080.1703438.9278309.15108.34050.07611.4663607.0263
113.03158.27270.17120.3045496.8050314.54921.84340.13100.4534489.1407
1211.92888.77670.06105.3783865.23803210.20037.39720.06952.1116646.0564
1311.56638.13760.06264.1553783.16703311.34128.54990.06363.6194750.5374
1410.93167.16980.06562.9075697.32453411.10957.74310.06473.1834718.1778
150.02250.10940.29090.5078329.48543511.74528.97430.06184.6904821.8985
160.41437.19250.38450.0905399.44453611.83158.68440.06144.9937839.9988
171.90768.65920.22410.2100473.06133711.86638.80010.06125.1255849.0706
1810.72581.96360.06672.6367653.55733811.50967.45070.06284.0075770.1729
197.76848.32200.08690.9913572.4912399.86918.69110.07151.8647635.9991
200.00290.17190.33030.3469328.0764408.61378.17240.07991.2486590.8052
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, T.-H.; Chiou, K.-C.; Chen, C.-M.; Chang, F.-M. Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation. Mathematics 2024, 12, 1163. https://doi.org/10.3390/math12081163

AMA Style

Liu T-H, Chiou K-C, Chen C-M, Chang F-M. Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation. Mathematics. 2024; 12(8):1163. https://doi.org/10.3390/math12081163

Chicago/Turabian Style

Liu, Tzu-Hsin, Kuo-Ching Chiou, Chih-Ming Chen, and Fu-Min Chang. 2024. "Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation" Mathematics 12, no. 8: 1163. https://doi.org/10.3390/math12081163

APA Style

Liu, T. -H., Chiou, K. -C., Chen, C. -M., & Chang, F. -M. (2024). Multiserver Retrial Queue with Two-Way Communication and Synchronous Working Vacation. Mathematics, 12(8), 1163. https://doi.org/10.3390/math12081163

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop