Abstract
A queueing-inventory system with server’s vacations and impatient customers is investigated. Continuous review policy is considered for the queueing-inventory system. The server leaves for a vacation once the inventory is zero. The customer arrived waits for service with a probability or is lost forever with a probability when the server is off due to a vacation. Customers are impatient during period when the customer is waiting for service in the queue or when the server is off due to a vacation. The matrix geometric solution of the steady-state probability is obtained using the method of truncated approximation. Furthermore, some performance measures are obtained using the steady-state probability. Finally, the effect of system parameters on the performance measures, the optimal policy and the optimal cost are presented by numerical examples.
1. Introduction
Over the last decades, many researchers have captured much attention on queueing systems with an attached inventory. Such a queueing system is called queueing-inventory system (QIS). A QIS is different from both the traditional inventory system and the traditional queuing system (see Zhao and Lian [1]). The QIS has a wide range of applications in fields such as integrated supply chain management [2], vehicle maintenance [3] and medical service [4].
The QIS was first introduced by Sigman and Simchi-Levi [5] and Melikov and Molchanov [6]. Sigman and Simchi-Levi [5] studied an M/G/1 QIS, where customers arrived when the inventory is zero were backlogged. A light traffic heuristic was developed to obtain the performance measures. Melikov and Molchanov [6] introduced QIS in investigating transportation/storage system (TSS), in which user requests were patient that they waited in the queue until the inventory was replenished. Besides, both exact and approximate solution methods were developed to derive the performance measures. Subsequently, more and more scholars are interested in QISs. Schwarz et al. [2] and Schwarz and Daduna [7] studied the M/M/1 QISs for cases of lost sales and backordering, respectively. Yue et al. [8] analyzed a QIS with batch demands for two cases, namely partial-lost sales and full-lost sales.
To utilize the server’s idle time, the QIS with servers’ vacation were studied by some scholar. Lakshmanan et al. [9] considered the multiple working vacations in QIS model, where any customer enters a single finite retrial orbit when the demanding item is stocked out or the server is busy. They obtained the joint steady-state distributions of the system, and the performance measures and long run expected total cost. Manikandan and Nair [10] studied an M/M/1 QIS model with working vacations, vacation interruptions and lost sales. They obtained various performance measures and the stationary waiting time distribution. Zhang et al. [11] investigated an M/M/1 QIS with server’s multiple vacations and a random order size policy, where lost sales was considered. The explicit product form solution for the steady-state distribution was obtained. The effect of parameters on performance measures was studied by using numerical examples and simulation method. For more results on this topic, we refer to the two survey papers given by Krishnamoorthy et al. [12, 13].
Vacation queueing models with customers’ impatience have been greatly studied. Altman and Yechiali [14] studied the M/M/1, M/G/1 and M/M/c queues with server’s vacation and customers’ impatience. They assumed customers become impatient when the server goes on vacation. Furthermore, they compared the proportion of customer abandonments under the single-vacation and multiple-vacation. And then, Altman and Yechiali [15] extended the impatience phenomenon studied by Altman and Yechiali [14] in an infinite-server queue system. Furthermore, more results on queue system with servers’ vacation and customer’s impatience can be refer to Bouchentouf et al. [16, 17], Bouchentouf and Guendouzi [18], Kadi et al. [19], Padmavathy et al. [20], Yue et al. [21–23].
The early analytical results of server’s vacations on QIS with impatient customers appeared in Viswanath et al. [24] who studied a very general model, where the customer demands followed the Markovian arrival process and the lead time followed a correlated process similar to it. Besides, the service times and the vacation time followed the phase type distributions. They assumed that whenever there were no on-hand inventory items or no customers in the system, the server took a vacation. Melikov et al. [25] studied a perishable QIS with early and delayed server vacations and they considered a finite queue of impatient claims. An asymptotic method was developed to obtain the performance measures. Koroliuk et al. [26] discussed a perishable QIS with impatient customers and a server’s vacation, where the restocking follows the two-level policy. The server took a vacation if either there was no on-hand inventory item or no customer or both none. They developed an asymptotic method to obtain the performance measures.
In this paper, we consider an M/M/1 QIS with server’s vacation and impatient customers. The server’s vacation and the sales policy of this model differ from the models that were studied in [24–26]. They all assumed that the server went on vacation if either there were no on-hand inventory items or no customers in the system or both none, and also assumed that the customers became impatient during the time of waiting for service. In our model, when the on-hand inventory is zero, the server leaves for a vacation. If the server finds that the inventory is not empty at the end of a vacation, the server returns back to the system and serves a customer waiting for service. Otherwise, the server takes another vacation immediately. Besides, we consider a mixed sales policy, i.e., the customer arrived waits for service with a probability or is lost forever with a probability when the server is off due to a vacation. Such mixed sales policy can be seen as a mixed of the lost sales and the backordering. If , this mixed sales becomes the lost sales, and if , it becomes the backordering case. To the best of our knowledge, all the research work on QIS with server’s vacation and impatient customer considers the lost sales and backordering separately. In presented model, we consider a mixed sales policy which integrates the lost sales and the backordering in a QIS with server’s vacation and impatient customers. How does the mixed sales policy affect the performance measures? In contrast with the lost sales policy () and the backordering policy (), is there a better mixed sales policy ()? This motivates us to the study of the present model.
A representative example of an application of our model is online sales of electronic products. The upgrading of electronic products such as mobile phone, smart TV, CPU, Graphics card is very fast because of their own particularity. Online retailers often sell their new products by online reservation. Some customers like to make reservation to buy a new generation of products, while some customers do not like to make reservation buy new products. This corresponds to the mixed sales policy of our model. The period of reservation can be regards as servers’ vacation. In this period of time, the servers do not provide service for customers until the deadline of reservation. Besides, some customers who make reservation may cancel their reservation. This customer is impatient customer of our model.
The main contributions of this paper are summarized as follows: (a) We derive the matrix geometric solution of the steady-state probability using the method of truncated approximation. (b) We investigate the effect of the probability and impatient rate on some performance measures. (c) The optimal policy and the optimal cost are computed by using the genetic algorithm. (d) We derive the queuing performance measures including the mean number of customers in system, the mean number of customers in queue, the mean sojourn time and the mean waiting time.
The rest of this paper is organized as follows. Section 2 describes the system model. In Section 3, we obtain the matrix geometric solution of the steady-state probability by using truncated approximation method. Furthermore, we derive some performance measures. Section 4 investigates the impact of the probability and impatient rate on some important performance measures, the optimal policy and the optimal cost. Section 5 concludes the paper.
2. System Model
Consider a QIS with impatient customers and multiple vacations. The arrival process of customers is Poisson process with rate . Customers are served one by one by a server under a First-Come, First-Served (FCFS) discipline. The service time follows exponential distribution with parameter .
Multiple vacations policy is considered. When the on-hand inventory is zero, the server leaves for a vacation. If the server finds that the inventory is not empty at the end of a vacation, the server returns back to the system and serves a customer waiting for service. Otherwise, the server takes another vacation immediately. The vacation time follows an exponential distribution with parameter .
The new arriving customer waits for service with probability or is lost forever with probability when the server is off due to a vacation. Customers are impatient during period when the customer is waiting for service in the queue or when the server is off due to a vacation. The impatient times follows an exponential distribution with parameter .
An ordering policy is adopted. When on-hand inventory level drops to , an order is placed such that the on-hand inventory level is restocked to level . The lead time follows an exponential distribution with parameter .
3. Steady-State Analysis
In this section, we develop a three-dimensional Markov process and derive the steady-state probability using method of matrix analysis. Then, we calculate the steady-state probability by using truncated approximation method and further derive some performance measures.
3.1. Steady-State Distribution
Let be the state process of the system, where denotes the number of customers at time , denotes the inventory level at time , denotes the inventory level at time , which is defined as follows: if the server is off due to a vacation, and if the server is on. Then, the process is a Markov process with state space:is the collection of states with , called the level . The state-transition diagram of the QIS is presented in Figure 1.

The state process is a level-dependent quasi-birth-and-death (LDQBD) process, and the infinitesimal generator of the process is as follows:where , , and are all square matrices of the order , and they are given bywhereandwhere is an identity matrix of the order .
The process is a level dependent quasi-birth-and-birth process, the steady-state probability distribution is defined as follows:
Let be the steady-state probability vector of the process , where
Then, satisfies the following set of equations:
The balance equations of the system is given by
The common method to solve the above balance equations is the probability generating function (PGF) method. However, these balance equations are so complex that they are difficult to obtain the explicit solution of the system by using PGF method to solve these equations. Alternatively, we solve these equations by using the matrix geometric solution (MGS) method. Comparing the PGF method, the advantage of MGS method is that we can obtain the steady-state probabilities of the system in matrix form, and then we can computer these steady-state probabilities in algorithmic approach.
Theorem 1. The components of the steady-state probability vector is given bywhere satisfies the set of equations:The matrix is the minimal nonnegative solution of the following equations:
Proof of Theorem 1. From equation (8), we have
Letwhere the matrix satisfies equation (13). Substituting equation (16) into the left side of equation (14), we have
Equation (14) is verified by using equation (11).
Substituting equation (16) into the left side of equation (15), we have
Thus, equation (15) is verified by using equation (13).
Applying the normalizing condition , we obtained equation (12). Thus, we complete the Proof of Theorem 1.
It is difficult to obtain the exact expressions of the matrix , from (13). However, based on Theorem 1, the steady-state probability can be computed approximately by using the method of truncated approximation given by Bright and Taylor [27]. The approximate steady-state probability obtained is denoted by . For more details on the procedure of the approximate algorithm for computing the steady-state probability, we refer the reader to Bright and Taylor [27]. The specific steps of method of truncated approximation are as follows:
If the condition of the process is normal, the is given bywhere Step 1: Solving . First, let , and then given such that . Step 2: Determining the value of . is the sum of the first terms in the sum of equation (19). We take where is chosen such that , where denotes the largest element in matrix . Step 3: we can obtain by equation (19). Solve subject to . Then, we obtain . Compute . Normalise , ,…, , s.t. until . Set . Step 4: solving . From equation (13), we have
Thus, we can obtain , , , by equation (21).
3.2. Performance Measures
In this subsection, we derive the following performance measures of the system by using the approximate steady-state probabilities .(1)The mean inventory level is given by where .(2)The mean number of replenishment per unit of time is given by where .(3)The mean number of order size per unit of time is given by(4)The probability that the server is off due to a vacation is given by(5)The mean number of lost customers due to customer impatient is given by(6)The mean number of lost customers due to server’s vacation is given by(7)The mean total number of lost customers per unit of time is given by(8)The mean arrival rate of customers who are admitted to the system per unit of time is given by(9)The mean number of customers in the system is given by where and .(10)The mean number of customers in the queue is given by(11)The mean sojourn time is given by(12)The mean waiting time is given by
4. Numerical Analysis
In this section, we conduct a numerical investigation to explore how the probability and the impatient rate affect some performance measures. Besides, we develop a cost model and compute the optimal policy and the optimal cost by using a genetic algorithm. The effects of system parameters on the optimal value and the optimal cost are also investigated.
4.1. Effects of the Mixed Probability and the Impatient Rate
We conduct a sensitivity analysis to show the effect of the two system parameters and on some performance measures. We plot the curves for , , , , and by varying the probability and the impatient rate , respectively. The results are presented by Figures 2–6. The system parameters are fixed as , , , , and .





The main observations can be summarized as follows:(i)Figure 2 shows that the mean inventory level first decreases slightly with an increase in the probability , and attains a minimum, and then increases significantly. This means there exists a mixed sales policy that minimizes the mean inventory level. Also, we observe that when is small, e.g. , the effect of on is much less than the case of large . Besides, Figure 2 shows that the effect of on for the case of small value is much higher than the cases of large values .This shows the sensitivity of on for the case of small value of is much less than the case of large value of , and that the sensitivity of to is much higher for the case of small value of than for the case of large value of .(ii)From Figure 3, we find that the mean number of order size increases with and decreases with . This can be explained as follows: The probability of customers willing to wait when the inventory level is zero increases with the probability , and the mean impatient time of customers increases with increasing of . This leads to the demands for inventory items which results in the increasing the mean number of order size . So, ,i.e., the lost sales policy is best policy among the mixed sales policy from the point of minimizing the mean number of order size .(iii)It is observed from Figures 4 and 5 tha decreases with either the probability or the impatient rate . However, increases with either the probability or the impatient rate . This agrees with our expectation. Since the probability of customers’ willing to wait increases with , and the mean impatient time increases with . Thus, the probability that the server go for vacations decreases with and . This explains why increase with and , and why decreases with and . This observation shows that the lost sales policy minimizes , and backordering policy minimizes .(iv)From Figure 6, we observe that the mean waiting time of customer increases with and decreases with . The reason for this observation is the same as we explained for our observation from Figure 3.
In summary, we have the following conclusion from the observations above:(i)There exists a mixed sales policy that minimizes the mean inventory level.(ii)The sensitivity of on for the case of small value of is much less than the case of large value of , and the sensitivity of on is much higher for the case of small value of than for the case of large value of .(iii)The lost sales policy is best policy among the mixed sales policy from the point of minimizing the mean number of order size .(iv)The lost sales policy minimizes , and backordering policy minimizes .
4.2. Optimization
In this section, we develop the following expected cost function by performance measures given in Subsection 3.2. The expected cost function per unit time can be given bywhere is a holding cost of inventory per unit time, is a fixed cost for placing an order, is a cost of per unit item for each order per unit time, is a cost incurred due to lost sales and impatience of customers per unit time, is a waiting cost per unit of time for each waiting customer in queue.
Unfortunately, we failed to obtain the explicit expressions for the expected cost function, and it is difficult to show the convexity of this cost function. Hence, we compute the optimal policy and the optimal cost by using the genetic algorithm (see [28]). The results are given in Tables 1–6. We assume the system parameters as , , , , , , unless their values are mentioned in the respective tables. The cost parameters are set as , , , , .
The main observations can be summarized as follows:(i)Both the optimal policy and the optimal cost increase with an increase in the arrival rate .(ii)The optimal value and the optimal cost decrease with an increase in the replenishment rate . The optimal value firstly increases and then decreases with the replenishment rate .(iii)The optimal policy and the optimal cost decrease with an increase in the vacation rate .(iv)The optimal value decreases slightly with an increase in the service rate , but the optimal value is independent of the service rate . We also observe that the optimal cost decreases with the service rate .(v)The optimal policy decreases with an increase in the impatient rate , and the optimal cost increases with the impatient rate .(vi)The optimal value and the optimal cost decrease with an increase in the probability , and the optimal value increases slightly with the probability when varies from 0.2 to 1.0.
Next, we explore the concavity and convexity of the cost function and the mean number of replenishment about the service rate. The curves are given by Figures 7 and 8. Figure 7 corresponds to the case of the varying parameter and fixed parameters , , , , , , , , , , , . Figure 8 corresponds to the case of the varying parameter and fixed parameters , , , , , , .


From Figures 7 and 8, we observe that the cost function and the mean number of replenishment exhibit convexity under the fixed parameters.
5. Conclusions
In this paper, we considered a QIS model with server’s multiple vacations and impatient customers, in which the new customer arrived waits for service with a probability when the server is off due to a vacation. We obtained the matrix geometric solution of the steady-state probability using the method of truncated approximation, and then derived some performance measures. We investigated the effect of the probability and the impatient rate on some performance measures using numerical analysis. The optimal policy and the optimal cost were computed using the genetic algorithm. The optimal service rate is derived. Our model can be further extended in the following directions: (i) studying non-Markovian QIS and (ii) investigating QIS with multiservers. Those directions become more challenging than the current one due to their analytical complexity.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China, grant no. 71971189, and the Key Project of Scientific Research in Higher Education of Hebei Province, China, grant no. ZD2018042.