Research Article  Open Access
Zaiming Liu, Xuxiang Luo, Jinbiao Wu, "Analysis of an M/PH/1 Retrial QueueingInventory System with Level Dependent Retrial Rate", Mathematical Problems in Engineering, vol. 2020, Article ID 4125958, 10 pages, 2020. https://doi.org/10.1155/2020/4125958
Analysis of an M/PH/1 Retrial QueueingInventory System with Level Dependent Retrial Rate
Abstract
We analyze a queueinginventory system which can model airline and railway reservation systems. An arriving customer to an idle server joins for service immediately with exactly one item from inventory at the moment of service completion if there are some onhand inventory, or else he accesses to a buffer of varying size (the buffer capacity varies and equals to the number of the items in the inventory with maximum size ). When the buffer overflows, the customer joins an orbit of infinite capacity with probability or is lost forever with probability . Arrivals form a Poisson process, and service time has phase type distribution. The time between any two successive retrials of the orbiting customer is exponentially distributed with parameter depending on the number of customers in the orbit. In addition, the items have a common life time with exponentially distributed. Cancellation of orders is possible before their expiry and intercancellation times are assumed to be exponentially distributed. The stability condition and steadystate probability vector have been studied by Neuts–Rao truncation method using the theory of Level Dependent QuasiBirthDeath (LDQBD) processes. Several stationary performance measures are also computed. Furthermore, we provide numerical illustration of the system performance with variation in values of underlying parameters and analyze an optimization problem.
1. Introduction
In real life, people often purchase products for future use or cancel the purchased products due to some reasons. A typical case is that customers reserve flight/train/bus seats, in which the seats are regarded as inventory. When the fight/train/bus departs, the inventory will be invalidated immediately, namely, the inventory products have common life time. In addition, some customers choose to cancel their order because of schedule changes, resulting in the increase of the inventory. Therefore, managers need to consider the above factors and formulate corresponding management policies to achieve the maximum economic benefits.
There is a growing research interest in various queueing systems with an attached inventory under different inventory management policies and lost sales. Here, we take a brief look at works carried out so far in some related context. A significant contribution for the development of queueinginventory system is [1], where stationary distribution of joint queue length and inventory process was derived in explicit product form for various M/M/1system with inventory under different inventory management policies. Schwarz and Daduna [2] investigated an system with inventory management, continuous review, exponentially distributed lead times, and backordering. They computed performance measures and derived optimality conditions under different order policies. Zhao and Lian [3] considered a queueinginventory system with two classes of customers, and found a priority service rule to minimize the longrun expected waiting cost. Saffari et al. [4] researched an M/M/1 queueing system with inventory under policy and lost sales, in which lead times are random variables. Krishnamoorthy et al. [5] considered a single server supply chain model with a production inventory system. Rejitha and Jose [6] considered a stochastic inventory model with positive service time and retrial of customers. A detailed review on queueinginventory models with various cases can be seen in [7].
Retrial queueing system is widely used in many different fields such as modern communication system and supply chain management system. Artalejo [8] introduced the retry strategy into the stochastic inventory system for the first time and studied the numerical solution and optimal decision problem. On this basis, Ushakumari [9] studied a queueinginventory system with retrial and random lead time and obtained the optimal inventory ordering strategy . However, the service time is often neglected in some reported queueinginventory systems. In order to describe the actual situation, Berman et al. [10] considered positive service time firstly and formulated a model where both the demand and service rates are assumed to be constant and deterministic. Thereafter, Berman and Kim [11] extended the above model to the case where service time was a random variable. A discrete product inventory control system with positive service time was considered by Shajin et al. [12]. Krishnamoorthy et al. [13] made a detailed review of inventory models with positive service time and suggested possible directions for future research. One may refer this for the development taken place in queueinginventory with retrial customers which contains some of the finest contributions such as [14–18].
Concepts of reservation, cancellation, and common life time was first introduced in the queueinginventory system by Krishnamoorthy et al. [19] with a constant retrial rate. In their model, customers arrive according to a Poisson process and service time follows an exponentially distributed. Subsequently, they extended the above model to a type queueinginventory systems with postponed work [20]. A queueinginventory system with impatient customers and overbooking was analyzed by Shajin and Krishnamoorthy [21], in which customers arrive according to the Markovian arrival process and service time has phase type distribution. Furthermore, Shajin et al. [22] considered a single serve queueinginventory system with facility for reservation in advance for the next time frames ahead. For more works on queueinginventory models with reservation, cancellation, and common life time, we refer the reader to [23].
The aim of our research is to present stationary distribution and optimal inventory policies for service system with an attached inventory under reservation, cancellation of orders, common life time, and retrial. We analyze a single server queueing system of type with an attached inventory. Customers arrive according to a Poisson Process, and service time of each customer follows a PH distribution. The retrial times for any repeated customer are exponentially distributed. By rationing and canceling purchases, the inventory is replenished. In addition, we assume the items in the inventory have a common life time. Based on the above integrated model, we construct a multidimensional continuous time Markov chain, which is a Level Dependent QuasiBirthDeath process (LDQBD). Note that the replenishment policies include rationing and cancellation of purchases, and common life time of all items is considered in our model, which is quite different from [6]. We provide the stationary probabilities and stationary condition using the theory of LDQBD. Several performance measures are calculated and their sensitivity to system parameters is analyzed. According to the calculation results, an optimization problem is discussed. By comprehensively analyzing the influence of various factors on the profit, it is illustrated that the change trend of the profit on the system parameters is in line with the actual situation, which provides an effective and scientific basis for the inventory manager.
2. Model Description
We consider an M/PH/1 retrial queueinginventory system with reservation, cancellation of orders, common life time, and retrial. The customers arrive according to a Poisson process with rate . The service time for each customer follows a PH distribution with representation of order . It has the distribution function , and the mean rate of service , where is the column vector of 1’s with appropriate order. There are items in the inventory at the beginning of a cycle. They have a common life time, which is exponentially distributed with parameter . If an arriving customer finds the server idle and there is at least one item in the inventory, the customer starts service immediately. At the moment of service completion, the customer takes exactly one item from the inventory away, and hence the onhand inventory decreases by one. Otherwise, if an arriving customer finds the server busy, he joins a buffer of varying size depending on the number of items in the inventory (the buffer capacity varies and equals to the number of items in the inventory). If the buffer has reached the maximum capacity(that is, the number of customers in the buffer equals to the number of items in inventory), then an arriving customer either joins an orbit of infinite capacity with probability to seek the service again and again or leaves the system forever with complementary probability . Interretrial times are exponentially distributed with rate , when there are customers in the orbit. In addition, all customers in the buffer are assumed to leave the system when the common life time of the items is realized, and the customers in orbit do not leave the system on the contrary (that is, they are willing to wait for service in the next cycle). Thus, we define a cycle as the time duration from the epoch at which we start with S items at a replenishment epoch, to the moment when the common life time is realized.
Two types of inventory replenishment policies are considered. One is the instantaneous replenishment (zero lead time) that the inventory reaches its maximum level on realization of common life time for the next cycle. The other is cancellation of orders (that is, returning of a sold item) before their expiry. It occurs according to an exponentially distributed interoccurrence time with parameter when items in the inventory. By cancellation of purchases, inventory is increased until their expiry time. Note that inventory level will not go above by cancellation since the maximum possible number “in sold list” is when items are held in the inventory. Cancellation time is assumed to be negligible. Customers are served in FCFS discipline. The above processes are independent of each other. We start with the case of customers being flushed out from the finite buffer on realization of common life time. The graphical description of the proposed model is shown in Figure 1.
Throughout this paper, we use the following notation: = number of customers in the orbit at time = number of customers in the buffer at time = number of items in the inventory at time = identity matrix of appropriate size = zero matrix of appropriate size ,
The process is a continuous time Markov chain (CTMC), which is a LDQBD, with state space:
The subset of the state space consists of such states that the common life time is realized and the inventory level reaches its maximum in way of fresh replenishment with no customer in the buffer.
We enumerate the states of the above CTMC in lexicographic order. Thus, the infinitesimal generator matrix of this CTMC has the following structure:
The matrices , , , and are square matrices of order :where and of dimension with and , respectively, at the position and . For , denote the ordered sets:
Then, the matrices and are square matrices of order which are defined bywhere
The matrices and of size are represented bywhere
3. SteadyState Analysis
In this section, we study the stability condition of the queueinginventory system and calculate the steadystate probability vector using the theory of LDQBD processes.
Denote the steadystate probabilities of the continuous time Markov chain as
Let denote the steadystate probability vector of the generator . When the system is stable, we have
Partitioning as , where each of the subvectors is denoted by
Due to the multiple phases, it is hard to express in the closed form (e.g., generating function) even though is highly structured. In the following, we will adopt the Neuts–Rao truncation method using the theory of LDQBD processes to solve the system of equations given in (1). For more details about this theory, interested readers are referred to the book [24].
3.1. Stability Condition
Using the Neuts–Rao truncation method, we identify a positive integer suitably, such that defining , , and for , we have the following generator:
Let be the steady probability vector of :where , , , and .
Then,
The components of are obtained aswhere
From the equation,and the normalizing condition , we have
The stability condition of the queueinginventory system under study is given by the following theorem.
Theorem 1. The queueinginventory system under study is stable if and only if
Proof. The queueinginventory system in the model with the type generator given in (2) is stable if and only if (see [25])From the elements of and , we haveUsing (3) and (5), we acquire the stability condition (4).
3.2. SteadyState Probability Vector
Assuming that the stability condition (4) is satisfied, we see that the steadystate probability vectors , are given by [25]where is the minimal nonnegative solution of the matrix quadratic equation:and the vectors , are obtained by solving the following equations:
The normalizing condition iswherewith
4. Performance Measures
In this section, we derive some stationary performance measures and specific probabilities descriptions for the inventory system.(1)Expected number of customers in the orbit:(2)Expected number of customers in the buffer:(3)Expected number of items in the inventory before realization of common life time:(4)Expected cancellation rate before realization of common life time:(5)Expected purchase rate before realization of common life time: where represents the transfer rate from state to absorbing state .(6)Expected number of items in the inventory at the moment on realization of common life time (that is the lost items in inventory when the next cycle start–a typical case is flights/trains/buses departing with one or more vacant seats): where Here, , for instance, represents the probability of realization of common life time before an arriving or a cancellation occurs when the system states are .(7)Probability that the inventory is full at the time of realization of common life time (all seats remain vacant at the time flight takes off):(8)Probability that the inventory is depleted at the time of realization of common life time (all seats are sold out at departure time):(9)Expected number of purchases before realization of common life time:(10)Expected number of cancellations before realization of common life time:
5. Numerical Examples
In this section, we present some numerical results that show the system performance measures with variations in values of underlying parameters.
We assume that the parameters satisfied the stability condition (4) are set as the following: . The PH type distribution of service time with representation is described by the matrices:
Thus, the mean rate of service is .
Table 1 provides the joint stationary probabilities of the queueinginventory system. It can be seen that the largest value of the marginal probabilities (lateral and longitudinal) are equal to 0.1790 and 0.5488, respectively, indicating that there are no items in inventory and no repeated customers in orbit.

5.1. Effect of the Maximum Inventory Level on Various Performance Measures
From Table 2, we observe that with the increasing of the maximum inventory level , the measures, such as expected number of customers in the buffer, expected number of items in the inventory, expected cancellation rate before realization of common life time, and expected number of items in the inventory at the moment on realization, are increased as expected. However, the expected number of customers in the orbit experiences a decline owing to the increased buffer size.

5.2. Effect of the Arrival Rate on Various Performance Measures
Table 3 reveals that the expected number of customers in the orbit and in the buffer and expected cancellation rate before realization of common life time are increased when the arrival rate increases, whereas the expected number of items in the inventory before realization of common life time increases first with and then keeps going down, which may be caused by cancellation of orders. However, the expected number of items in the inventory immediately on realization of common life time shows a decreasing trend, resulting from an increasing demand arrival.

5.3. Effect of the Retrial Rate on Various Performance Measures
As shown in Table 4, we observe that an increase in retrial rate makes a sharp decrease in expected number of customers in orbit and a mild decrease in expected cancellation rate before realization of common life time. The former is attributed to large number of customers in the orbit getting into the buffer as increases. However, there is a slight increase in expected number of customers in the buffer and expected number of items in the inventory before realization of common life time as increases. To sum up, the retrial rate rarely affects the performance measures except the expected number of customers in orbit; this is true obviously.

5.4. Effect of Variation in on
We assign the following values to the parameters: . Keeping fixed, the effect of on measures are depicted through Figure 2. As shown in Figure 2(a), we observe that the probability for the inventory being full at the time of realization of common life time (at ) increases sharply as increases. Similarly, the probability for the inventory being depleted (at 0) decreases slowly with increasing value of . It can be seen that common life time has a great influence on the probability but almost does not affect the probability . The results accord with the expectations. Figure 2(b) shows that both expected number of purchases and expected number of cancellations before realization of common life time decrease with increasing value of for fixed . This is attributed to faster realization of common life time.
(a)
(b)
5.5. Effect of Variation in on
We assign the same values as above. Keeping fixed, the effect of on measures are displayed in Figure 3. Probability for the inventory being full at the time of realization of common life time (at ) is a linearly increasing function of (see Figure 3(a)). Conversely, the increase in results in a decrease in probability for the inventory being zero (at 0) at the time of realization of common life time. The former phenomenon can be explained that the set of states with inventory level at the maximum could be revisited several times due to cancellations in a cycle. The latter is due to the adding inventory as the value increases. Figure 3(b) indicates that an increase in makes expected number of purchases before realization of common life time increase. Meanwhile, it is strange that expected number of cancellations before realization of common life falls off. The reason accounting for this phenomenon can be expressed by the provided in Section4. The value of depends on the multiplication factor and the probabilities. The former decreases as increases. On the contrary, the probabilities are risen rapidly with increasing value of . Thus, their combined effect leads to a decline in the value.
(a)
(b)
6. Optimization
Based on the above performance measures, we construct a cost function for checking the optimality of the inventory capacity . It may be noted that the system profit will increase when cancellation takes place before realization of common life time since there is a cancellation penalty imposed on the customer. Hence, we define profit function aswhere represents the revenue to the system per unit cancellation of inventory purchased; represents the revenue to the system per customer per unit time in the buffer (this is an income provided by being sold items in account of inventory for sure); and represents the holding cost per unit time per inventory.
The objective of this section is to find the optimal inventory level that maximizes expected profit function and analyze the effect of system parameters , respectively, on expected profit. The cost coefficients are assumed that per cancellation, per customer per time, and per item per time. Table 5 shows the optimal inventory level and the maximum profit when .

In Table 5, we note that the optimal inventory size , and the maximum expected profit is $78.8495. Figure 4 presents expected profit with the changing of system parameters , respectively, when other parameters are set the same as Table 5. In this figure, the expected profit is nondecreasing in and nonincreasing in as expected. Those phenomenons can be easily explained that there are more arriving demands and more inventory items due to the larger arrival rate and higher cancellation rate, respectively. In addition, the shorter the common life time, the lesser the arriving demands and the number of cancellations, which results in the latter phenomenon. It is noteworthy that the retrial rate has little impact on the expected profit. The appearance of this phenomenon may be based on the fact that the mean number of repeated customers is very small compared to the arriving customers.
7. Conclusions
This paper investigates a queueinginventory system with reservation in advance, cancellation of orders, common life time, and retrial. Customers arrive according to a Poisson process and service time for each customer has phase type distribution. The time between any two successive retrials of the orbiting customer is distributed as an exponential with parameter depending on the number of customers in the orbit. We have constructed a multidimensional continuoustime Markov chain with infinitesimal generator which is level dependent. The stability condition and steadystate probability vector have been studied by the Neuts–Rao truncation method using . Several stationary performance measures have been obtained. In addition, we have established an optimization model for the inventory level and analyzed the effect of system parameters on various performance measures and expected profit.
This research can be generalized to the situation that each customer demands more than one unit items (one may order more than one ticket). Vacation policy and priority customers can also be introduced in this model.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have on conflicts of interest.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (Grant no.11671404).
References
 M. Schwarz, C. Sauer, H. Daduna, R. Kulik, and R. Szekli, “M/M/1 queueing systems with inventory,” Queueing Systems, vol. 54, no. 1, pp. 55–78, 2006. View at: Publisher Site  Google Scholar
 M. Schwarz and H. Daduna, “Queueing systems with inventory management with random lead times and with backordering,” Mathematical Methods of Operations Research, vol. 64, no. 3, pp. 383–414, 2006. View at: Publisher Site  Google Scholar
 N. Zhao and Z. Lian, “A queueinginventory system with two classes of customers,” International Journal of Production Economics, vol. 129, no. 1, pp. 225–231, 2011. View at: Publisher Site  Google Scholar
 M. Saffari, S. Asmussen, and R. Haji, “The M/M/1 queue with inventory, lost sale, and general lead times,” Queueing Systems, vol. 75, no. 1, pp. 65–77, 2013. View at: Publisher Site  Google Scholar
 A. Krishnamoorthy, D. Shajin, and B. Lakshmy, “Product form solution for some queueinginventory supply chain problem,” Opsearch, vol. 53, no. 1, pp. 85–102, 2016. View at: Publisher Site  Google Scholar
 K. R. Rejitha and K. P. Jose, “A stochastic inventory system with two modes of service and retrial of customers,” Opsearch, vol. 55, no. 1, pp. 134–149, 2018. View at: Publisher Site  Google Scholar
 K. Karthikeyan and R. Sudhesh, “Recent review article on queueing inventory systems,” Research Journal of Pharmacy and Technology, vol. 9, no. 11, pp. 2056–2066, 2016. View at: Publisher Site  Google Scholar
 J. R. Artalejo, A. Krishnamoorthy, and M. J. LopezHerrero, “Numerical analysis of (s, S) inventory systems with repeated attempts,” Annals of Operations Research, vol. 141, no. 1, pp. 67–83, 2006. View at: Publisher Site  Google Scholar
 P. V. Ushakumari, “On (s, S) inventory system with random lead time and repeated demands,” Journal of Applied Mathematics and Stochastic Analysis, vol. 2006, pp. 1–22, 2006. View at: Publisher Site  Google Scholar
 O. Berman, E. H. Kaplan, and D. G. Shevishak, “Deterministic approximations for inventory management at service facilities”,” IIE Transactions, vol. 25, no. 5, pp. 98–104, 1993. View at: Google Scholar
 O. Berman and E. Kim, “Stochastic models for inventory management at service facilities,” Communications in Statistics. Stochastic Models, vol. 15, no. 4, pp. 695–718, 1999. View at: Publisher Site  Google Scholar
 D. Shajin, B. Benny, R. V. Razumchik, and A. Krishnamoorthy, “Discrete product inventory control system with positive service time and two operation modes,” Automation and Remote Control, vol. 79, no. 9, pp. 1593–1608, 2018. View at: Publisher Site  Google Scholar
 A. Krishnamoorthy, B. Lakshmy, and R. Manikandan, “A survey on inventory models with positive service time,” Opsearch, vol. 48, no. 2, pp. 153–169, 2011. View at: Publisher Site  Google Scholar
 T. Yang and J. G. C. Templeton, “A survey on retrial queues,” Queueing Systems, vol. 2, no. 3, pp. 201–233, 1987. View at: Publisher Site  Google Scholar
 B. K. Kumar, A. Vijayakumar, and D. Arivudainambi, “An M/G/1 retrial queueing system with twophase service and preemptive resume,” Annals of Operations Research, vol. 113, no. 14, pp. 61–79, 2002. View at: Publisher Site  Google Scholar
 B. Sivakumar, “A perishable inventory system with retrial demands and a finite population,” Journal of Computational and Applied Mathematics, vol. 224, no. 1, pp. 29–38, 2009. View at: Publisher Site  Google Scholar
 R. P. Hannah and M. G. Subramanian, “Two server (s, S) inventory system with positive service time, positive lead time, retrial customers and negative arrivals,” International Journal of Computer Applications, vol. 62, no. 10, pp. 9–13, 2013. View at: Google Scholar
 T. Karthick, B. Sivakumar, and G. Arivarignan, “An inventory system with two types of customers and retrial demands,” International Journal of Systems Science: Operations & Logistics, vol. 2, no. 2, pp. 90–112, 2015. View at: Publisher Site  Google Scholar
 A. Krishnamoorthy, D. Shajin, and B. Lakshmy, “On a queueinginventory with reservation, cancellation, common life time and retrial,” Annals of Operations Research, vol. 247, no. 1, pp. 365–389, 2016. View at: Publisher Site  Google Scholar
 A. Krishnamoorthy, D. Shajin, and B. Lakshmy, “GI/M/1 type queueinginventory systems with postponed work, reservation, cancellation and common life time,” Indian Journal of Pure and Applied Mathematics, vol. 47, no. 2, pp. 357–388, 2016. View at: Publisher Site  Google Scholar
 D. Shajin and A. Krishnamoorthy, “On a queueinginventory system with impatient customers, advanced reservation, cancellation, overbooking and common life time,” Operational Research, vol. 47, pp. 1–25, 2019. View at: Google Scholar
 D. Shajin, A. Krishnamoorthy, A. N. Dudin, V. C. Joshua, and V. Jacob, “On a queueinginventory system with advanced reservation and cancellation for the next K time frames ahead: the case of overbooking,” Queueing Systems, vol. 94, no. 12, pp. 3–37, 2020. View at: Publisher Site  Google Scholar
 D. Shajin, On QueueingInventory Models  Product Form Solution; Reservation, Cancellation and Common Life Time, Cochin University of Science and Technology, Kochi, Kerala, 2015.
 G. Latouche and V. Ramaswami, Introduction to Matrix Analytic Methods in Stochastic Modeling, Siam, New Delhi, India, 1999.
 M. F. Neuts, Matrixgeometric Solutions in Stochastic Models : An Algorithmic Approach, The Johns Hopkins University Press, Baltimore, ML, USA, 1981.
Copyright
Copyright © 2020 Zaiming Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.