1. Introduction
For the first time, the heading “Queuing-Inventory System” (QIS) was introduced in the papers [
1,
2]. These are the systems in which in order to service incoming consumer customers (
c-customers), along with idle server, some resources are required. These resources are called inventories.
The QIS theory has been intensively developed in the past three decades. The first publications in this direction are [
3,
4]. In [
4] an M/G/1 type QIS with exponentially distributed lead time and under light traffic is considered. They obtained product form solution for steady-state analysis. In [
3] a Markov Decision Process (MDP) approach to optimize the finite QIS is used where order size is a controllable parameter. They developed an exact and approximate (for large scale QIS) methods to solution of the constructed MDP.
The current state of QIS theory and its applications are described in detail in a recent review paper [
5]. Therefore, we indicate here only related works published after this review work, as well as those works that were not included in its bibliography.
This work is related to (i) QIS with destructible stocks, (ii) multi-parametric replenishment policies (RP) in double sources QIS and (iii) multiple type c-customers. In what follows, we shortly review the relevant literature.
Let us note that here we distinguish between the QISs with destructible and perishable stocks. In turn, in a class of perishable QISs two sub-classes of systems are differentiated: (1) QIS with individual life time (ILF) and (2) QIS with common life time (CLT). In QIS with ILT each item can perish independently of the others, while in QIS with CLT all items perish together, e.g., foods with the same expiry date, medicines manufactured with the same expiry date and so on. In considering perishable QIS with ILT, we draw from a vast stream of literature, see, e.g., [
6,
7,
8,
9,
10,
11,
12]. Models of QIS with CLT has been little studied. For the first time, such a model was examined in [
13], next in [
14] and recently, in [
15]. Some review of works in this direction can be found in the last paper. Note that QIS with CLT can be treated as QIS with catastrophes in inventory part (but not in a service facility).
Almost all studies on perishable QIS of both types are devoted to models in which items can perish within a certain time interval. However, in practice, there are QIS in which items can be destroyed instantly, e.g., (1) a shop selling breakable items in which destruction of the items may be the result of a worker’s careless actions; (2) a company selling electrical/electronic goods in which destruction of the item can occur even during the sale as a result of a power surge in the electrical network. Models of QIS with destructible inventory are similar to models of QIS with ILT, but their main difference is as follows: in QIS with ILT a perish rate of inventory is a linear function of inventory level (i.e., if inventory level is
m, then perish rate is
mγ, where
is a life time of individual item), but in QIS with destructible inventory, destruction rate is independent on inventory level, i.e., it is a constant quantity. In other words, QIS with destructible inventory can be treated as QIS with ILT in which perish rate is a constant quantity (as in models of retrial queues with linear and constant retrial rates from orbit). Models of QIS with destructible inventory have hardly been studied, see [
16]. In this paper, we consider models of QIS with destructible inventory.
In the vast majority of works, models of QIS with a single source are studied. One problem of optimal choosing of the source from a set of sources with different supply rates and costs was solved in [
17], where the optimality criterion was the total cost. In this work, it was assumed that in the entire period of the system operation, an order is made to only one source. However, the models of QIS with multiple sources in which an order can be sent to different sources depending on the current inventory level represents important scientific and practical interest. In other words, in such QISs the distribution of the order between a fast and an expensive supplier and a slow but inexpensive supplier are very important problems. Solutions to these problems require determining appropriate replenishment policies (RPs).
Recently, only a few papers have considered models of double-source QISs. Original RP in double-source QIS was recently proposed in [
18,
19]. Proposed in these papers, RP is defined as follows: whenever the inventory level falls to
r (>
S/2), an order with the temporary (regular) supplier for
items is placed,
where
is maximal capacity of store: whenever the inventory level reaches
s, an order of
items to the regular supplier is placed and immediately cancels the order for
items from the temporary supplier (to be fair, in [
18] do not explicitly state the cancelation procedure but in [
19] this fact is shown clearly). In both papers, the QIS models with infinite queue for
c-customers are considered and Neuts’ matrix-analytic method (MAM) [
20] is used.
Other kinds of RP in double-sources QISs were proposed in [
21,
22]. The main feature of the proposed policies is that they are based on classical RPs, i.e.,
and
policies are used for construction of the new RPs. In other words, in both RPs no order is placed until the inventory level is above the reorder point
.
In [
21] two models of double-source QIS with instantly damaging inventory were proposed: one model uses the (
s,
S)-policy, and the other model uses the (
s,
Q)-policy during the entire period of the system’s operation (for definiteness, we note that (
s,
S)-policy (sometimes it is called “Up to
S” policy) is a RP in which the order volume is as much as it is able to bring the level back to S at the replenishment epoch and in (
s,
Q)-policy the order volume is fixed and is equal to
Q =
S-
s). In both RPs, choice of the source is made depending on the current inventory level as follows: a regular order to SS is placed when the inventory level drops to the value
s, and an emergency order to FS is placed when the inventory level drops to some critical value
r, where 0 ≤
r <
s. Similar but hybrid RP in double-source QIS was proposed in [
22]. In hybrid RP, it is assumed that a regular order to the Slow Source (SS) is placed in accordance with (
s,
Q)-policy while an emergency order to the Fast Source (FS) is placed in accordance with (
s,
S)-policy. In all cases, at the time of placing an emergency order, the regular order is immediately canceled.
After the publications of papers [
21,
22], the authors became aware of works [
23,
24]. It appears that the first paper dealing with QIS with multiple sources was [
23]. In that paper, the service facility of examined QIS is a M/M/1/∞ queue, and the selection of source among two sources with different lead times is determined according to Bernoulli trials. Arriving customers are lost when inventory is out of stock (lost sales scheme) and fixed size order RP is used. The authors obtained explicit formulas for calculating the steady-state distribution of constructed 3D Markov chain. Note that the implementation of randomized source selection schemes in practice encounters certain difficulties, and sometimes becomes impossible. The following single-server Markovian QIS model with lost sales scheme and a deterministic rule for source choosing is explored in [
22]: a fixed quantity
is ordered to supplier
when on-hand inventory level reaches
where
. It is assumed that when a supplier is about to deliver an order, the existing outstanding order to the other supplier, if any, is automatically cancelled. By using the results of [
25] authors proof that the stationary distribution of queue length is independent of the inventory level, and is identical to the stationary distribution of the queue length in a classical M/M/1/distribution of, i.e., the joint stationary distribution of queue length and inventory level has multiplicative form. Note that in order to prove this fact, the key assumption is the use of the lost sales scheme. It is important to note that in the works [
21,
22], as well as in this paper, it is assumed that an arriving
c-customer can queue up even if the warehouse is empty, i.e., no lost sale scheme is used.
Classification of
c-customers in QIS can be done by using their various property. So, models of QIS with high and low priorities of
c-customers and with different RPs were considered in [
26,
27,
28,
29,
30,
31,
32,
33,
34].
In the indicated above papers the following admission scheme for different types of
c-customers is applied: if inventory level is exceeding some critical level then
c-customers of both types are serviced, but if it is lower than the critical one the service is given only to high priority
c-customers. Another admission scheme was proposed in [
35]: high-priority
c-customers are accepted if there is at least one free space in the buffer, and low priority
c-customers are only accepted when the total queue length is less than the given threshold value. Models of QISs in which
c-customers are required various size of inventory was examined in [
36].
Note that the classification of
c-customers can be carried out not at the instant of their arrival, but after the completion of the service process. In other words, in some QIS part of
c-customers may refuse to purchase the item after being served due to some reasons, i.e., they require only service and do not require items. Such a multi-class QIS model was first studied in [
37] and since then a series of papers have taken this phenomenon into account; see, e.g., the recent paper [
38] and its reference list. In this paper, we also take into account this phenomenon.
We refer the readers to [
14,
39,
40,
41,
42,
43,
44] for a discussion of more complex QIS models with Markovian Arrival Process (MAP) and phase-type distribution of service and/or lead times.
In all available papers devoted to double-source QIS, the models of systems with infinite queue for c-customers are considered and lost sale scheme is used. In this paper, we propose novel RPs in double-source QIS with finite room for c-customers and destructible stocks items to improve its performance measures by switching the order from slow source to fast one depending on the current level of inventory. Moreover, here unlike the known papers, no lost sale scheme for c-customers is used.
The main contributions of this paper can be summarized as follows:
A novel double-source QIS with finite room for waiting an impatient c-customers and destructible stocks under no lost sale scheme is formulated as a two-dimensional Markov process;
The novel replenishment policies in double-source QIS based on classical (s, S) and (s, Q) policies are proposed. An additional level-dependent parameter in these policies is defined that determine switching of order from slow to fast source;
An approximate method to derive closed-form expressions for the system’s steady-state probabilities and for its performance measures is developed;
By numerical experiments it is shown that under certain realistic assumptions related to the ratio of some system’s load parameters, the developed approximate method has high accuracy;
Assuming linear costs for each waiting and losing customer and each stored item, results of a cost analysis are provided, demonstrating how the optimal values of the replenishment policies is affected by the load parameters of the system.
The paper has the following structure. In
Section 2 we describe the models by spelling out the assumptions of the underlying variables and proposed replenishment policies. Exact and approximate approaches to steady-state analysis of the models including the list of main performance measures are given in
Section 3. Results of illustrative numerical examples along with optimization problems are demonstrated in
Section 4. Concluding remarks are given in
Section 5.
3. The Steady State Analysis of the Models
First consider the model QIS under policy. Mathematical model of this QIS is two dimensional Markov chain (2D MC) with states , where is the total number of c-customers in the system (in buffer or being serviced), and is on-hand inventory level, . State space of this 2D MC is given as .
Transition rates from state
to state
are denoted by
. So, the investigated 2D MC has an infinitesimal generator
with the following transition rates for
:
Hereinafter, is the indicator function of the event A, which is 1 if A is true and 0 otherwise.
Based on relations (1)–(8) we conclude that the constructed 2D MC is a reducible chain, i.e., for any positive values of the initial parameters, stationary mode is existing. Let
represent the probability of the state
. The system of global balance equations (SGBE) is expressed as follows:
In Equation (9) the following notations are used: is the subset of states of , that can be reached from state in one step; is the subset of states of , from which it is possible to go to the state in one step.
Normalizing condition should be added to SGBE (9):
Note that SGBE (9), (10) represent the system of linear algebraic equations with dimension . For moderate values of the parameters N and S, constructed SGBE (9), (10) can be solved by using existing software. However, certain computational difficulties arise in cases where the dimension of the SGBE is large. To overcome these difficulties, an approximate method for calculating the steady-state probabilities is proposed below.
The method developed below can be correctly applied to models of QIS in which the following assumption is fulfilled: the rate of c-customers exceeds both the rate of d-customers and the replenishment rate, . Note that this assumption corresponds to the working mode of real QISs. In addition, as a rule, in real QISs, the average service time of c-customers that receive inventory is much longer than the average service time of c-customers that receive only service, i.e., .
Under these assumptions, consider the following splitting of the initial state space
of the studied 2D MC:
where
Note 1. From relations (1)–(8) we conclude that when indicated above assumptions are fulfilled, the transitions intensities between the states within classesis greater than the transitions intensities between the states from different classes. Therefore, the application of the space merging method for the examined 2D MC is correct, see Appendix of the book [45].
Next, each subset of states is combined into one merged state , and the following merging function is determined in the initial state space . The merged states form the set
Below we show that when using the
-policy, the approximate values of the steady-state probabilities are calculated as
where
Note 2. Since the values of state probabilitiesfor casesdo not depend on parameter(see first row of the Formula (13)), in the rest of the paper in their notations subscriptfor casesis omitted for simplification. Additionally, note that in casesand/or, all state probabilitiesfor anyand.
Indeed, when using splitting (11) from relations (1)–(8), we conclude that the stationary distribution of all split models with state spaces
coincide with the stationary distribution of a single-server Markov queuing system M/M/1/N with a load
Stationary probability of state
within these split models are denoted by
. These quantities do not depend on
and they are calculated as
In addition, from relations (1)–(8) we conclude that the stationary distribution of a split model with state space
coincide with the stationary distribution of a single-server Markov queuing system M/M/1/N with a load
In other words, stationary probabilities of states
within this split model are calculated as
Combining Formulas (15) and (16), we obtain Formula (13).
Merged model represents 1D MC with state space
. Transition rate from merged state
to other merged state
is denoted by
. These quantities are calculated as
By taking into account relations (1)–(8), (15)–(17), after some algebra we obtain:
From (18) after some algebra we obtain that the stationary probabilities of merged states, , are calculated by (14).
Finally, according to the space merging algorithm, we obtain that the steady-state probabilities of the initial 2D MC are calculated from Formula (12).
Now consider the model under
-policy. For this model state space is also defined by the set
E. Here, the infinitesimal generator
of the appropriate 2D MC are determined via (1)–(6), but here Formulas (7) and (8) are replaced by the following formulas:
Similarly, we can prove that the when using the
-policy, the approximate values of the steady-state probabilities are calculated as
where
4. Performance Measures
Desired performance measures are calculated via steady-state probabilities as follows.
Average inventory level
is
Average supply from source-
i when using
-policy
are
Average supply from source-
i when using
-policy are
Average number of
c-customers in system
is
Average destruction rate of the stock
is
Average reorder rate of regular supply
is
Average reorder rate of emergency supply
is
Note 3. The average rate of emergency orders is equal to the average rate of cancellation (RC) of regular orders.
Loss probability of
c-customers
is
In (28), the first and second terms of the sum estimate the loss probabilities of c-customers when they enter the system, provided that at this moment the inventory level is zero and there is no free space in the buffer, respectively, the third and fourth terms determine the loss probabilities c-customers due to their impatience in the queue, provided that before the start of their service, the inventory level drops to zero.
Taking into account relations (12)–(14), after standard transformations, it can be shown that when using the
-policy, the approximate values of the performance measures are calculated as
In a similar way, taking into account relations (19) and (20), we can show that when using the
-policy, the approximate values of the performance measures are calculated via (29), (31)–(35), where the quantities
are replaced by the
, and an average supply from various sources are calculated as
Below, we consider the results of numerical experiments carried out using the obtained formulas.
5. Numerical Results
This section is threefold. Firstly, we study the accuracy of the developed approximate formulas to calculate the steady-state probabilities, secondly, we investigate the behavior of performance measures versus initial parameters of the model and thirdly, we consider the optimization problem.
Consider the first goal of the numerical experiments. Analytical investigation of the accuracy of the developed approximate formulas is impossible. Therefore, we consider solution of the problem via numerical experiments. Exact values of the steady-state probabilities are calculated via SGBE (9), (10). Some results of numerical experiments are demonstrated in
Table 1. From this table, we conclude that under both RPs the accuracy of the developed formulas to calculation of the steady-state probabilities are sufficiently high for practice. Results of numerical experiments show that the accuracy of the developed approximate formulas systematically increases with the increases in the dimension of the system, i.e., when
N and
S take on large values. The last fact is very important, since the approximate approach was developed specifically for large scale models.
Now, consider the second goal of the numerical experiments. The proposed approach makes it possible to carry out numerical experiments to study the behavior of performance measures (21)–(28) versus any initial system’s parameter. Below, in
Table 2,
Table 3,
Table 4,
Table 5,
Table 6,
Table 7,
Table 8,
Table 9 and
Table 10, we present the results of appropriate calculations where in each column the top row corresponds to the (
s,
S)-policy, and the lower one, to the (
s,
Q)-policy. Due to the limited size of the paper, conclusions regarding the behavior of performance measures are left to the reader.
Finally, consider the third aim of the numerical experiments, i.e., minimization of the expected total cost (TC), which is calculated as follows:
where
are the fixed costs of one order from Source-
i;
are the procurement costs per unit inventory from Source-
i, is the cost due to order cancellation from Source-1;
is the holding cost of an inventory unit per unit time;
is the cost due to damage per unit inventory;
is the cost of losing a
c-customer;
is the waiting cost of a
c-customer per unit time.
Let us assume that controllable are only parameters
s and
r. Then, the optimization problem is formulated as follows: it is required to find such pairs of optimal values
in order to minimize the functional (37). In other words, we need to solve the following task:
Some results of solving this problem for a model when using the (
s,
S) and (
s,
Q)-policies are shown in
Table 11. The optimal solution of the problem for both policies is