Next Article in Journal
Approach to the Numerical Study of Wave Processes in a Layered and Fractured Porous Media in a Two-Dimensional Formulation
Next Article in Special Issue
On the Control over the Distribution of Ticks Based on the Extensions of the KISS Model
Previous Article in Journal
Queueing Inventory System in Transport Problem
Previous Article in Special Issue
Random Motions at Finite Velocity on Non-Euclidean Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Double-Sources Queuing-Inventory Systems with Finite Waiting Room and Destructible Stocks

1
Institute of Control Systems, National Academy of Science, Baku AZ 1141, Azerbaijan
2
Faculty of Engineering, National Aviation Academy, Baku AZ 1045, Azerbaijan
3
Department of Informatics and Networks, Faculty of Informatics, University of Debrecen, 4032 Debrecen, Hungary
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(1), 226; https://doi.org/10.3390/math11010226
Submission received: 12 December 2022 / Revised: 29 December 2022 / Accepted: 30 December 2022 / Published: 2 January 2023

Abstract

:
Models of double-source queuing-inventory systems are studied in the presence of a finite buffer for waiting in the queue of consumer customers, where instant destruction of inventory is possible. It is assumed that the lead times of orders, as well as the cost of delivery from various sources, differ from each other. Replenishment of stocks from various sources is carried out according to the following scheme: if the inventory level drops to the reorder point s, then a regular order for the supply of inventory to a slow source is generated; if the inventory level falls below a certain threshold value r, where r < s, then the system instantly cancels the regular order and generates an emergency order to the fast source. Models of systems that use (s, S) or (s, Q) replenishment policies are studied. Exact and approximate methods for finding the performance measures of the models under study are proposed. The problems of minimizing the total cost are solved by choosing the appropriate values of the parameters s and r when using different replenishment policies. Numerical examples demonstrated the high accuracy of an approximate method as well as compared performance measures of the system under various replenishment policies.
MSC:
60K25; 90B22; 68M20

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 , where γ 1 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 Q 1 items is placed, Q 1 = S r , where S is maximal capacity of store: whenever the inventory level reaches s < Q 1 < r , s, an order of Q 2 = S s ( > s + 1 ) items to the regular supplier is placed and immediately cancels the order for Q 1 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., ( s , S ) and ( s , Q ) 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 s ,   s < S / 2 .
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 Q i is ordered to supplier i when on-hand inventory level reaches r i where max i ( r i ) min i ( Q i ) . 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.

2. The Models

Below we accept the following assumptions:
  • Single-server QIS contains the store of maximal capacity S ,   S < , and c-customers forms the Poisson flow with rate λ;
  • An incoming c-customer is accepted for servicing immediately if the server is idle at that moment and the inventory level is positive; if the inventory level is positive and the server is busy, then this c-customer is placed in a queue of finite length. The maximum buffer size for waiting of c-customers is N 1 ,   N < ; if upon arrival of the c-customer buffer is full, then this c-customer is lost with probability (w.p.) one;
  • An incoming c-customer can join the queue even when the inventory level is zero, i.e., if upon arrival of the c-customer an inventory level is zero, but there is free space in the buffer, then it either joins the queue w.p. φ 1 , or leaves the system w.p. φ 2 , where φ 1 + φ 2 = 1 ;
  • After the service is completed, each c-customer either does not receive items w.p. σ 1 , or receives items w.p. σ 2 , where σ 1 + σ 2 = 1 . In both cases, the service time of c-customers has exponential distribution function (d.f.) with different mean values, i.e., if the c-customer does not receive the items, then the average time of its service is equal to μ 1 1 ; otherwise, this time is equal to μ 2 1 ;
  • The c-customer at the head of the queue becomes impatient if the inventory level drops to zero, i.e., in such cases, the c-customer at the head of the queue waits for a random time, which has an exponential d.f. with an average τ 1 , and after this time it leaves the system without items;
  • Along with c-customers, the system receives a Poisson flow of destructive customers (d-customers) with rate κ , i.e., upon arrival of d-customers, the inventory level is instantly reduced by one. It is assumed that the d-customer may even destroy the item that is in the release status to the c-customer; if the inventory level is equal to zero, then the arrived d-customer does not affect the operation of the system;
  • Replenishment can be done either from the SS or from the FS. The lead time from each source has an exponential d.f., i.e., if an order is made to SS, then the average waiting time for the supply of stock is ν 1 1 , and if an order is made to FS, then the average waiting time for the supply of stock is ν 2 1 , where ν 2 > ν 1 ;
  • The system can use either (s, S) or (s, Q) replenishment policies. In both policies, it is assumed that if the stock level drops to the order point s ,   0 < s < ( S / 2 ) , then an order is made to SS, and when the stock level drops to the threshold value r ,   0 r < s , then the order is instantly canceled from SS and the order is sent to FS. The time required for a cancellation issue is negligible, i.e., this time is zero. It is also assumed that canceling an order from an SS results in certain penalties.
The problem is to find the joint distribution of the number of c-customers in the system and the inventory level of the system and to determine the main performance measures: average inventory level, average supply from each source, average reorder rate to each source, average destruction rate of inventory due to the receipt of d-customers, average length of queue of c-customers and loss probability of c-customers. In addition, the problems of optimizing the system by choosing the optimal values of the reorder point s and the threshold value r are of interest.

3. The Steady State Analysis of the Models

First consider the model QIS under ( s , S ) policy. Mathematical model of this QIS is two dimensional Markov chain (2D MC) with states ( n , m ) , where n is the total number of c-customers in the system (in buffer or being serviced), n = 0 ,   1 , , N , and m is on-hand inventory level, m = 0 , 1 , , S . State space of this 2D MC is given as E = { 0 , 1 , , N } × { 0 , 1 , , S } .
Transition rates from state ( n 1 , m 1 ) E to state ( n 2 , m 2 ) E are denoted by q ( ( n 1 , m 1 ) ; ( n 2 , m 2 ) ) . So, the investigated 2D MC has an infinitesimal generator Q 1 = ( q 1 ( n 1 , m 1 ) ; ( n 2 , m 2 ) ) :   ( n 1 , m 1 ) , ( n 2 , m 2 ) E with the following transition rates for ( n 1 , m 1 ) E :
q 1 ( ( n 1 , m 1 ) ; ( n 1 + 1 , m 1 ) ) = λ φ 1 I ( n 1 < N ) · I ( m 1 = 0 ) ;
q 1 ( ( n 1 , m 1 ) ; ( n 1 + 1 , m 1 ) ) = λ I ( n 1 < N ) · I ( m 1 > 0 ) ;
q 1 ( ( n 1 , m 1 ) ; ( n 1 1 , m 1 ) ) = μ 1 σ 1 · I ( m 1 > 0 )
q 1 ( ( n 1 , m 1 ) ; ( n 1 1 , m 1 1 ) ) = μ 2 σ 2 · I ( m 1 > 0 )
q 1 ( ( n 1 , m 1 ) ; ( n 1 , m 1 1 ) ) = κ · I ( m 1 > 0 )
q 1 ( ( n 1 , m 1 ) ; ( n 1 1 , m 1 ) ) = τ · I ( n 1 > 0 ) · I ( m 1 = 0 )
q 1 ( ( n 1 , m 1 ) ; ( n 1 , S ) ) = ν 1 · I ( r < m 1 s )
q 1 ( ( n 1 , m 1 ) ; ( n 1 , S ) ) = ν 2 · I ( 0 m 1 r ) .
Hereinafter, I ( A ) 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 p ( n , m ) represent the probability of the state ( n , m ) E . The system of global balance equations (SGBE) is expressed as follows:
p ( n 1 , m 1 ) ( n 2 , m 2 ) E ( n 1 , m 1 ) + q ( ( n 1 , m 1 ) , ( n 2 , m 2 ) ) = ( n 2 , m 2 ) E ( n 1 . m 1 ) q ( ( n 2 , m 2 ) , ( n 1 , m 1 ) ) p ( n 2 , m 2 ) ,   ( n 1 , m 1 ) E
In Equation (9) the following notations are used: E ( n 1 , m 1 ) + is the subset of states of E , that can be reached from state ( n 1 , m 1 ) in one step; E ( n 1 , m 1 ) is the subset of states of E , from which it is possible to go to the state ( n 1 , m 1 ) in one step.
Normalizing condition should be added to SGBE (9):
( n , m ) E p ( n , m ) = 1 .
Note that SGBE (9), (10) represent the system of linear algebraic equations with dimension ( N + 1 ) ( S + 1 ) . 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, λ > m a x ( κ , ν ) . 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., μ 1 > μ 2 .
Under these assumptions, consider the following splitting of the initial state space E of the studied 2D MC:
E = m = 0 S E m , E m 1 E m 2 = ,   m 1 m 2 ,
where E m = { ( n , m ) E : n = 0 , 1 , , N }   ,   m = 0 , 1 , , S .
Note 1.
From relations (1)–(8) we conclude that when indicated above assumptions are fulfilled, the transitions intensities between the states within classes E m ,   m = 0 , 1 , ,   S , is 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 E m is combined into one merged state m , and the following merging function is determined in the initial state space E :   U ( n , m ) = m   ,   ( n , m ) E . The merged states form the set E ^ = { m : m = o , 1 , , S } .
Below we show that when using the ( s , S ) -policy, the approximate values of the steady-state probabilities are calculated as
p ( n , m ) ρ m ( n ) π 1 ( m ) ,
where
ρ m ( n ) = { θ n 1 θ 1 θ N + 1   , m > 0 ,   n = 0 , 1 , , N , ψ n 1 ψ 1 ψ N + 1   , m = 0 ,   n = 0 , 1 , , N ,
θ = λ / μ 1 σ 1 , ψ = φ 1 / τ ;
π 1 ( m ) = { α m π 1 ( 0 ) , 1 m r , β m π 1 ( 0 ) , r + 1 m s , β s π 1 ( 0 ) , s + 1 m S ,
π 1 ( 0 ) = ( 1 + m = 1 r α m + m = r + 1 s β m + ( S s ) β s ) 1 ;
α m = a 2 ( 1 + a 2 ) m 1   ,   β m = a 2 ( 1 + a 2 ) r ( 1 + a 1 ) m r 1   ; a k = ν k μ 2 σ 2 ( 1 ρ ( 0 ) ) + κ ,   k = 1 ,   2 ;
Note 2.
Since the values of state probabilities ρ m ( n ) for cases m > 0 do not depend on parameter m (see first row of the Formula (13)), in the rest of the paper in their notations subscript m for cases m > 0 is omitted for simplification. Additionally, note that in cases θ = 1 and/or ψ = 1 , all state probabilities ρ m ( n ) = 1 / ( N + 1 ) for any n ,   n = 0 , 1 , , N and m ,   m = 0 , 1 , , S .
Indeed, when using splitting (11) from relations (1)–(8), we conclude that the stationary distribution of all split models with state spaces E m ,   m = 1 ,   ,   S , coincide with the stationary distribution of a single-server Markov queuing system M/M/1/N with a load θ = λ / μ 1 σ 1 . Stationary probability of state ( n , m ) E m within these split models are denoted by ρ m ( n ) . These quantities do not depend on m = 1 ,   ,   S , and they are calculated as
ρ ( n ) = θ n 1 θ 1 θ N + 1   ,   n = 0 ,   1 ,   ,   N .
In addition, from relations (1)–(8) we conclude that the stationary distribution of a split model with state space E 0 coincide with the stationary distribution of a single-server Markov queuing system M/M/1/N with a load ψ = λ φ 1 / τ . In other words, stationary probabilities of states ( n , 0 ) E 0 within this split model are calculated as
ρ 0 ( n ) = ψ n 1 ψ 1 ψ N + 1   ,   n = 0 ,   1 ,   ,   N .
Combining Formulas (15) and (16), we obtain Formula (13).
Merged model represents 1D MC with state space E ^ . Transition rate from merged state m 1 E ^ to other merged state m 2 E ^ is denoted by q ( m 1 ,   m 2 ) . These quantities are calculated as
q ( m 1 ,   m 2 ) = ( n 2 , m 2 ) E m 2 q ( ( n 1 , m 1 ) ,   ( n 2 , m 2 ) ) ρ m 1 ( n 1 ) .
By taking into account relations (1)–(8), (15)–(17), after some algebra we obtain:
q ( m 1 ,   m 2 ) = { ν 2 , 0 m 1 r , m 2 = S ,   ν 1 , r + 1 m 1 s , m 2 = S , μ 2 σ 2 ( 1 ρ ( 0 ) ) + κ , m 1 > 0 ,   m 2 = m 1 1 .
From (18) after some algebra we obtain that the stationary probabilities of merged states, π 1 ( m ) ,   m E ^ , 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 ( s , Q ) -policy. For this model state space is also defined by the set E. Here, the infinitesimal generator
Q 2 = ( q 2 ( n 1 , m 1 ) ; ( n 2 , m 2 ) ) :   ( n 1 , m 1 ) , ( n 2 , m 2 ) E
of the appropriate 2D MC are determined via (1)–(6), but here Formulas (7) and (8) are replaced by the following formulas:
q 1 ( ( n 1 , m 1 ) ; ( n 1 , m 1 + S s ) ) = ν 1 · I ( r < m 1 s )
q 1 ( ( n 1 , m 1 ) ; ( n 1 , m 1 + S s ) ) = ν 2 · I ( 0 m 1 r )
Similarly, we can prove that the when using the ( s , Q ) -policy, the approximate values of the steady-state probabilities are calculated as
p ( n , m ) ρ m ( n ) π 2 ( m ) ,
where
π 2 ( m ) = { α m π 2 ( 0 ) , 1 m r , β m π 2 ( 0 ) , r + 1 m s , β s π 2 ( 0 ) , s + 1 m S s , η m π 2 ( 0 ) , S s + 1 m S ,
π 2 ( 0 ) = ( 1 + m = 1 r α m + m = r + 1 s β m + ( S 2 s ) β s + m = S s + 1 S η m ) 1 ,
η m = 1 μ 2 σ 2 ( 1 ρ ( 0 ) ) + κ i = m S + s s χ i x i ,   S s + 1 m S ;
χ i = { α i , 1 i r , β i , r + 1 i s ; x i = { ν 2 , 1 i r , ν 1 , r + 1 i s .

4. Performance Measures

Desired performance measures are calculated via steady-state probabilities as follows.
Average inventory level ( S a v ) is
S a v = m = 1 S m n = 0 N p ( n , m ) .
Average supply from source-i when using ( s , S ) -policy ( V a v ( i ) ) are
V a v ( 1 ) = m = r + 1 S ( S m ) n = 0 N p ( n , m ) ; V a v ( 2 ) = m = 0 r ( S m ) n = 0 N p ( n , m ) ;
Average supply from source-i when using ( s , Q ) -policy are
V a v ( 1 ) = ( S s ) m = r + 1 s n = 0 N p ( n , m ) ; V a v ( 2 ) = ( S s ) m = 0 r n = 0 N p ( n , m ) ;
Average number of c-customers in system ( L a v ) is
L a v = n = 1 N n m = 0 S p ( n , m ) ;
Average destruction rate of the stock ( D R S ) is
D R S = κ ( 1 n = 0 N p ( n , 0 ) ) ;
Average reorder rate of regular supply ( R R 1 ) is
R R 1 = κ p ( 0 , s + 1 ) + ( μ 2 σ 2 + κ ) n = 1 N p ( n , s + 1 ) ;
Average reorder rate of emergency supply ( R R 2 ) is
R R 2 = κ p ( 0 , r + 1 ) + ( μ 2 σ 2 + κ ) n = 1 N p ( n , r + 1 ) .
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 ( P L ) is
P L = φ 2 n = 0 N p ( n , 0 ) + m = 0 S p ( N , m ) + τ τ + λ φ 1 + ν 2 n = 1 N 1 p ( n , 0 ) + τ τ + ν 2 p ( N , 0 ) .
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 ( s , S ) -policy, the approximate values of the performance measures are calculated as
S a v m = 1 S m π 1 ( m ) ;
V a v ( 1 ) m = r + 1 s ( S m ) π 1 ( m ) ; V a v ( 2 ) m = 0 r ( S m ) π 1 ( m ) ;
L a v n = 1 N n ( ( 1 π 1 ( 0 ) ) ρ ( n ) + π 1 ( 0 ) ρ 0 ( n ) ) ;
D R S κ ( 1 π 1 ( 0 ) ) ;
R R 1 π 1 ( s + 1 ) ( κ ρ ( 0 ) + ( μ 2 σ 2 + κ ) ( 1 ρ ( 0 ) ) ) ;
R R 2 π 1 ( r + 1 ) ( κ ρ ( 0 ) + ( μ 2 σ 2 + κ ) ( 1 ρ ( 0 ) ) ) ;
P L φ 2 π ( 0 ) ( 1 ρ 0 ( N ) ) + τ τ + λ φ 1 + ν 2 π ( 0 ) ( 1 ρ 0 ( 0 ) ρ 0 ( N ) ) + τ τ + ν 2 ρ 0 ( N ) π ( 0 ) ( 1 π ( 0 ) ) ρ ( N ) + π ( 0 ) ρ 0 ( N ) .
In a similar way, taking into account relations (19) and (20), we can show that when using the ( s , Q ) -policy, the approximate values of the performance measures are calculated via (29), (31)–(35), where the quantities π 1 ( m ) are replaced by the π 2 ( m ) , and an average supply from various sources are calculated as
V a v ( 1 ) = ( S s ) m = r + 1 s π 2 ( m ) ;           V a v ( 2 ) = ( S s ) m = 0 r π 2 ( m ) ;
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.
In all numerical experiments (in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9), it is assumed that σ 1 = 0.4 ,   φ 1 = 0.6 and S = 30 ,   N = 100 . The values of other parameters are indicated in the title of the tables.
Finally, consider the third aim of the numerical experiments, i.e., minimization of the expected total cost (TC), which is calculated as follows:
TC ( s , r ) = i = 1 2 ( K i + c r ( i ) V a v ( i ) ) R R i + c c R R 2 + c h S a v + c d D R S + c l λ P L + c w L a v ,
where K i are the fixed costs of one order from Source-i; c r ( i ) are the procurement costs per unit inventory from Source-i,  i = 1 , 2 ;   c c is the cost due to order cancellation from Source-1; c h is the holding cost of an inventory unit per unit time; c d is the cost due to damage per unit inventory; c l is the cost of losing a c-customer; c w 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 ( s * , r * ) in order to minimize the functional (37). In other words, we need to solve the following task:
( s * , r * ) = arg m i n ( s , r ) X T C ( s , r ) .
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 ( s * , r * ) = ( 6 ,   1 ) .

6. Conclusions

Models of double-source QIS are studied, where the delivery rates from various sources differ from each other. They use two replenishment policies, i.e., with a fixed supply volume and with a variable supply volume. In each policy, along with the order point, a threshold parameter is introduced that determines the moment when an order is canceled from a slow source and at the same time a new order is generated to a fast source. Models with a finite buffer for c-customers are studied in which c-customers are queued even if inventory level is zero and there is free space in the waiting room. However, c-customers in the queue can be impatient when the inventory level drops to zero before their service is started. After the service is completed, part of the c-customers, according to the Bernoulli trials, either leaves the system without receiving items, or receives items. Along with c-customers, d-customers also enter the system, the appearance of which leads to an instant decrease in the system’s inventory level by one.
Exact and approximate methods for calculating the performance measures of the considered QIS are proposed. The exact method is based on the balance equations of the constructed 2D MC, while the approximate method uses the space merging approach for the 2D MC. When using the approximate method, explicit formulas for calculating both the steady-state probabilities and the performance measures are obtained. The conducted numerical experiments demonstrate the high accuracy of the developed approximate formulas. By the numerical experiments, the behavior of the performance measures versus initial parameters of the models was also studied. In addition, the problems of finding the optimal values of the parameters of the used RPs to minimize the total cost were solved.
Suggestions for future directions of research are:
  • Development and investigation of replenishment policies in double source QIS in which switching the order from slow source to fast depends both on inventory level and number of c-customers in the system;
  • Investigation of similar QIS models with more complex d.f. for c- and d-customers, service time of c-customers, lead times etc.;
  • Investigation of similar QIS models with common life times of the items.

Author Contributions

Conceptualization, A.M., R.M. and J.S.; methodology, A.M., R.M. and J.S.; software, R.M.; validation, R.M.; formal analysis, A.M., R.M. and J.S.; investigation, A.M., R.M. and J.S.; writing-original draft preparation, A.M., R.M. and J.S.; writing-review and editing, A.M., R.M. and J.S.; supervision, A.M.; project administration, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Study did not require ethical approval.

Informed Consent Statement

Not applicable.

Data Availability Statement

Did not report any data.

Acknowledgments

The authors are very grateful to the reviewers for their valuable comments and suggestions, which improved the quality and the presentation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schwarz, M.; Daduna, H. Queuing Systems with Inventory Management with Random Lead Times and with Backordering. Math. Meth. Oper. Res. 2006, 64, 383–414. [Google Scholar] [CrossRef]
  2. Schwarz, M.; Sauer, C.; Daduna, H.; Kulik, R.; Szekli, R. M/M/1 Queuing Systems with Inventory. Queuing Syst. 2006, 54, 55–78. [Google Scholar] [CrossRef]
  3. Melikov, A.Z.; Molchanov, A.A. Stock optimization in transportation/storage systems. Cybern. Syst. Anal. 1992, 28, 484–487. [Google Scholar] [CrossRef]
  4. Sigman, K.; Simchi-Levi, D. Light traffic heuristic for anM/G/1 queue with limited inventory. Ann. Oper. Res. 1992, 40, 371–380. [Google Scholar] [CrossRef]
  5. Krishnamoorthy, A.; Shajin, D.; Narayanan, W. Inventory with Positive Service Time: A Survey. In Advanced Trends in Queueing Theory; Series of Books “Mathematics and Statistics” Sciences; Anisimov, V., Limnios, N., Eds.; ISTE & Wiley: London, UK, 2021; pp. 201–238. [Google Scholar]
  6. Hanukov, G.; Avinadav, T.; Chernonog, T.; Spiegel, U.; Yechiali, U. A queueing system with decomposed service and inventoried preliminary services. Appl. Math. Model. 2017, 47, 276–293. [Google Scholar] [CrossRef]
  7. Hanukov, G.; Avinadav, T.; Chernonog, T.; Spiegel, U.; Yechiali, U. Improving efficiency in service systems by performing and storing “preliminary services”. Int. J. Prod. Econ. 2018, 197, 174–185. [Google Scholar] [CrossRef]
  8. Hanukov, G.; Avinadav, T.; Chernonog, T.; Yechiali, U. Performance improvement of a service system via stocking perishable preliminary services. Eur. J. Oper. Res. 2018, 274, 1000–1011. [Google Scholar] [CrossRef]
  9. Hanukov, G.; Avinadav, T.; Chernonog, T.; Yechiali, U. A multi-server QIS with stock-dependent demand. IFAC PapersOnLine 2019, 52, 671–676. [Google Scholar] [CrossRef]
  10. Hanukov, G.; Avinadav, T.; Chernonog, T.; Yechiali, U. A service system with perishable products where customers are either fastidious or strategic. Int. J. Prod. Econ. 2020, 228, 107696. [Google Scholar] [CrossRef]
  11. Ko, S.S. A Nonhomogeneous Quas-Birth Process Approach for an (s, S) Policy for a Perishable Inventory System with Retrial Demands. J. Ind. Manag. Opt. 2020, 16, 1415–1433. [Google Scholar] [CrossRef]
  12. Melikov, A.; Shahmaliyev, M.; Nair, S.S. Matrix-geometric Method to Study Queuing System with Perishable Inventory. Autom. Remote Control. 2021, 82, 2168–2181. [Google Scholar] [CrossRef]
  13. Lian, Z.; Liu, L.; Neuts, M.F. A Discrete-Time Model for Common Lifetime Inventory Systems. Math. Oper. Res. 2005, 30, 718–732. [Google Scholar] [CrossRef]
  14. Chakravarthy, S.R. An inventory system with Markovian demands, phase type distributions for perishability and replenishment. Opsearch 2010, 47, 266–283. [Google Scholar] [CrossRef]
  15. Shajin, D.; Krishnamoorthy, A.; Manikandan, R. On a queueing-inventory system with common life time and Markovian lead time process. Oper. Res. 2020, 1–34. [Google Scholar] [CrossRef]
  16. Melikov, A.; Aliyeva, S.; Nair, S.S.; Kumar, B.K. Retrial Queuing-Inventory Systems with Delayed Feedback and Instantaneous Damaging of Items. Axioms 2022, 11, 241. [Google Scholar] [CrossRef]
  17. Melikov, A.; Krishnamoorthy, A.; Shahmaliyev, M. Numerical Analysis and Long Run Total Cost Optimization of Perishable Queuing Inventory Systems with Delayed Feedback. Queuing Model. Serv. Manag. 2019, 2, 83–111. [Google Scholar]
  18. Soujanya, M.L.; Laxmi, P.V. Analysis on Dual Supply Inventory Model Having Negative Arrivals and Finite Lifetime Inventory. Reliab. Theory Applications 2021, 16, 295–301. [Google Scholar]
  19. Vinitha, V.; Anbazhagan, N.; Amutha, S.; Jeganathan, K.; Shrestha, B.; Song, H.-K.; Joshi, G.P.; Moon, H. Analysis of a Stochastic Inventory Model on Random Environment with Two Classes of Suppliers and Impulse Customers. Mathematics 2022, 10, 2235. [Google Scholar] [CrossRef]
  20. Neuts, M.F. Matrix-Geometric Solutions in Stochastic Models. An Algorithmic Approach; Johns Hopkins University Press: Baltimore, MD, USA, 1981. [Google Scholar]
  21. Melikov, A.Z.; Mirzayev, R.R.; Nair, S.S. Numerical Study of a Queuing-Inventory System with Two Supply Sources and Destructive Customers. J. Comput. Syst. Sci. Int. 2022, 61, 581–598. [Google Scholar] [CrossRef]
  22. Melikov, A.; Mirzayev, R.; Nair, S.S. Double Sources Queuing-Inventory System with Hybrid Replenishment Policy. Mathematics 2022, 10, 2423. [Google Scholar] [CrossRef]
  23. Saffari, M.; Haji, R.; Hassanzadeh, F. A queueing system with inventory and mixed exponentially distributed lead times. Int. J. Adv. Manuf. Technol. 2010, 53, 1231–1237. [Google Scholar] [CrossRef]
  24. Saffari, M.; Sajadieh, M.S.; Hassanzadeh, F. A queuing system with inventory and competing suppliers. Eur. J. Ind. Eng. 2019, 13, 420–433. [Google Scholar] [CrossRef]
  25. Saffari, M.; Asmussen, S.; Haji, R. The M/M/1 queue with inventory, lost sale, and general lead times. Queueing Syst. 2013, 75, 65–77. [Google Scholar] [CrossRef] [Green Version]
  26. Dekker, R.; Hill, R.; Kleijn, M.; Teunter, R. On the (S-1, S) lost sales inventory model with priority demand classes. Nav. Res. Logist. (NRL) 2002, 49, 593–610. [Google Scholar] [CrossRef] [Green Version]
  27. Ha, A.Y. Stock-rationing policy for a make-to-stock production system with two priority classes and backordering. Nav. Res. Logist. (NRL) 1997, 44, 457–472. [Google Scholar] [CrossRef]
  28. Ha, A.Y. Inventory rationing in a make-to-stock production system with two priority classes and backordering. Manag. Sci. 1997, 43, 1093–1103. [Google Scholar] [CrossRef]
  29. Ha, A.Y. Stock rationing in an M/Ek/1 make-to-stock queue. Manag. Sci. 2000, 46, 77–87. [Google Scholar] [CrossRef]
  30. Isotupa, K.P.S. An (S, Q) Markovian inventory system with lost sales and two demand classes. Math. Comp. Model. 2006, 43, 687–694. [Google Scholar] [CrossRef]
  31. Isotupa, K.P.S. An (S, Q) inventory system with two demand classes of customers. Int. J. Oper. Res. 2011, 12, 12–19. [Google Scholar] [CrossRef]
  32. Karthick, T.; Sivakumar, B.; Arivarignan, G. An inventory system with two types of customers and retrial demands. Int. J. Syst. Sci. Oper. Logist. 2015, 2, 90–112. [Google Scholar] [CrossRef]
  33. Kranenburg, A.; van Houtum, G. Cost optimization in the (S−1, S) lost sales inventory model with multiple demand classes. Oper. Res. Lett. 2007, 35, 493–502. [Google Scholar] [CrossRef]
  34. Melikov, A.Z.O.; Ponomarenko, L.A.; Aliyev, I.A.O. Markov Models of Queuing-Inventory Systems with Different Types of Retrial Customers. J. Autom. Inf. Sci. 2019, 51, 1–15. [Google Scholar] [CrossRef]
  35. Melikov, A.Z.O.; Ponomarenko, L.A.; Aliyev, I.A.O. Analysis and Optimization of Models of Queuing-Inventory Systems with Two Types of Requests. J. Autom. Inf. Sci. 2018, 50, 34–50. [Google Scholar] [CrossRef]
  36. Melikov, A.; Fatalieva, M. Situational Inventory in Counter-Stream Serving Systems. Eng. Simul. 1998, 15, 839–848. [Google Scholar]
  37. Krishnamoorthy, A.; Manikandan, R.; Lakshmy, B. A revisit to queueing-inventory system with positive service time. Ann. Oper. Res. 2013, 233, 221–236. [Google Scholar] [CrossRef]
  38. Shajin, D.; Krishnamoorthy, A.; Melikov, A.Z.; Sztrik, J. Multi-Server Queuing Production Inventory System with Emergency Replenishment. Mathematics 2022, 10, 3839. [Google Scholar] [CrossRef]
  39. Chakravarthy, S.R. Queueing-Inventory Models with Batch Demands and Positive Service Times. Autom. Remote Control 2020, 81, 713–730. [Google Scholar] [CrossRef]
  40. Chakravarthy, S.R. Analysis of a Queueing Model with M AP Arrivals and Heterogeneous Phase-Type Group Services. Mathematics 2022, 10, 3575. [Google Scholar] [CrossRef]
  41. Chakravarthy, S.; Rao, B. Queuing-Inventory Models with MAP Demands and Random Replenishment Opportunities. Mathematics 2021, 9, 1092. [Google Scholar] [CrossRef]
  42. Chakravarthy, S.R.; Maity, A.; Gupta, U.C. Modeling and Analysis of Bulk Service Queues with an Inventory under (s, S) Policy. Ann. Oper. Res. 2017, 258, 263–283. [Google Scholar] [CrossRef]
  43. Chakravarthy, S.R.; Rumyantsev, A. Analytical and simulation studies of queueing-inventory models with MAP demands in batches and positive phase type services. Simul. Model. Pract. Theory 2020, 103, 102092. [Google Scholar] [CrossRef]
  44. Choudhary, A.; Chakravarthy, S.R.; Sharma, D.S. Analysis of M AP/PH/1 Queueing System with Degrading Service Rate and Phase Type Vacation. Mathematics 2021, 9, 2387. [Google Scholar] [CrossRef]
  45. Ponomarenko, L.; Kim, C.S.; Melikov, A. Performance Analysis and Optimization of Multi-Traffic on Communication Networks; Springer: London, UK, 2010. [Google Scholar]
Table 1. Dependence of the absolute error of the steady-state probabilities on parameter s ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , S = 22 ,   N = 100 .
Table 1. Dependence of the absolute error of the steady-state probabilities on parameter s ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , S = 22 ,   N = 100 .
sMax of Error
(s, S)-Policy(s, Q)-Policy
61.06 × 10−21.31 × 10−2
71.15 × 10−21.66 × 10−2
81.26 × 10−21.98 × 10−2
91.38 × 10−22.26 × 10−2
101.37 × 10−22.54 × 10−2
Table 2. Dependence of the performance measures on parameter λ ; μ 1 = 45 , μ 2 = 15 , κ = 2 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 2. Dependence of the performance measures on parameter λ ; μ 1 = 45 , μ 2 = 15 , κ = 2 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
λ Vav(1)Vav(2)SavLavDRSRR1RR2PL
810.9680.26218.7700.8152.0000.2588.1520 × 10−28.1520 × 10−2
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.2171 × 10−5
910.9930.30618.7011.0212.0000.2789.4923 × 10−29.4923 × 10−2
3.3990.0955.7831.0072.0000.0862.9353 × 10−26.4781 × 10−5
1011.0130.35218.6351.2801.9990.2971.0882 × 10−11.0882 × 10−1
3.5010.1125.9231.2592.0000.0953.4592 × 10−29.4994 × 10−5
1111.0300.39918.5721.6111.9990.3171.2315 × 10−11.2315 × 10−1
3.5960.1306.0561.5842.0000.1034.0156 × 10−21.3405 × 10−4
1211.0420.44718.5102.0531.9990.3371.3785 × 10−11.3785 × 10−1
3.6870.1496.1812.0182.0000.1124.6028 × 10−21.8316 × 10−4
Table 3. Dependence of the performance measures on parameter κ ; λ = 8 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 3. Dependence of the performance measures on parameter κ ; λ = 8 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
κ Vav(1)Vav(2)SavLavDRSRR1RR2PL
210.9680.26218.7700.8152.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.2171 × 10−5
311.0130.35218.6350.8302.9990.2971.0882 × 10−12.8717 × 10−4
3.5010.1125.9230.8093.0000.0953.4592 × 10−29.1282 × 10−5
411.0420.44718.5100.8533.9980.3371.3785 × 10−15.1361 × 10−4
3.6870.1496.1810.8183.9990.1124.6028 × 10−21.7150 × 10−4
511.0580.54818.3940.8864.9960.3751.6813 × 10−18.3369 × 10−4
3.8560.1916.4130.8304.9980.1315.8624 × 10−22.9069 × 10−4
611.0640.65218.2830.9305.9920.4131.9935 × 10−11.2579 × 10−3
4.0110.2376.6280.8475.9970.1507.2263 × 10−24.5597 × 10−4
Table 4. Dependence of the performance measures on parameter μ 1 ; λ = 8 , κ = 2 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 4. Dependence of the performance measures on parameter μ 1 ; λ = 8 , κ = 2 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
μ 1 Vav(1)Vav(2)SavLavDRSRR1RR2PL
4510.9680.26218.7700.8152.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.2171 × 10−5
5010.9440.22918.8270.6772.0000.2427.1212 × 10−21.0050 × 10−4
3.1990.0675.5040.6702.0000.0712.0819 × 10−22.9382 × 10−5
5510.9220.20218.8760.5792.0000.2286.3090 × 10−27.4425 × 10−5
3.1200.0585.3920.5742.0000.0651.8023 × 10−22.1261 × 10−5
6010.9010.18118.9180.5062.0000.2175.6559 × 10−25.6757 × 10−5
3.0510.0515.2940.5022.0000.0611.5828 × 10−21.5883 × 10−5
6510.8820.16418.9540.4492.0000.2085.1217 × 10−24.4383 × 10−5
2.9890.0455.2070.4462.0000.0571.4070 × 10−21.2192 × 10−5
Table 5. Dependence of the performance measures on parameter μ 2 ; λ = 8 , κ = 2 , μ 1 = 45 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 5. Dependence of the performance measures on parameter μ 2 ; λ = 8 , κ = 2 , μ 1 = 45 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
μ 2 Vav(1)Vav(2)SavLavDRSRR1RR2PL
1510.9680.26218.7700.8152.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.2171 × 10−5
2011.0250.38318.5920.8371.9990.3101.1833 × 10−13.5303 × 10−4
3.5650.1246.0120.8122.0000.1003.8266 × 10−21.1416 × 10−4
2511.0540.51418.4320.8741.9980.3621.5792 × 10−17.1588 × 10−4
3.8010.1776.3380.8251.9990.1255.4305 × 10−22.4617 × 10−4
3011.0640.65218.2830.9301.9970.4131.9935 × 10−11.2579 × 10−3
4.0110.2376.6280.8471.9990.1507.2263 × 10−24.5597 × 10−4
3511.0600.79618.1441.0071.9960.4642.4198 × 10−11.9972 × 10−3
4.2010.3026.8910.8791.9980.1769.1906 × 10−27.5854 × 10−4
Table 6. Dependence of the performance measures on parameter τ ;   λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 6. Dependence of the performance measures on parameter τ ;   λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
τ Vav(1)Vav(2)SavLavDRSRR1RR2PL
1.510.9680.26218.7700.8152.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.2171 × 10−5
310.9680.26218.7700.8142.0000.2588.1520 × 10−21.2463 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.1657 × 10−5
4.510.9680.26218.7700.8122.0000.2588.1520 × 10−21.0370 × 10−4
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.1100 × 10−5
610.9680.26218.7700.8002.0000.2588.1520 × 10−29.6548 × 10−5
3.2910.0795.6330.8042.0000.0772.4463 × 10−24.0504 × 10−5
7.510.9680.26218.7700.8002.0000.2588.1520 × 10−29.3768 × 10−5
3.2910.0795.6330.8042.0000.0772.4463 × 10−23.9870 × 10−5
Table 7. Dependence of the performance measures on parameter ν 1 ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 2 = 8 , r = 5 , s = 10 .
Table 7. Dependence of the performance measures on parameter ν 1 ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 2 = 8 , r = 5 , s = 10 .
ν 1 Vav(1)Vav(2)SavLavDRSRR1RR2PL
210.9680.26218.7700.8152.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0790.0790.8042.0000.0772.4463 × 10−24.2171 × 10−5
310.7820.16819.0500.8092.0000.2645.2089 × 10−28.9793 × 10−5
2.7740.0430.0430.8022.0000.0681.3401 × 10−22.3102 × 10−5
410.6500.11219.2380.8062.0000.2683.4710 × 10−25.9835 × 10−5
2.4190.0250.0250.8012.0000.0617.8834 × 10−31.3590 × 10−5
510.5530.07719.3700.8042.0000.2712.3967 × 10−24.1315 × 10−5
2.1560.0160.0160.8012.0000.0554.8959 × 10−38.4397 × 10−6
610.4790.05519.4660.8032.0000.2731.7058 × 10−22.9405 × 10−5
1.9500.0100.0100.8012.0000.0513.1749 × 10−35.4731 × 10−6
Table 8. Dependence of the performance measures on parameter ν 2 ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , r = 5 , s = 10 .
Table 8. Dependence of the performance measures on parameter ν 2 ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , r = 5 , s = 10 .
ν 2 Vav(1)Vav(2)SavLavDRSRR1RR2PL
810.9680.26218.7700.8152.0002.5764 × 10−18.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0007.7317 × 10−22.4463 × 10−24.2171 × 10−5
910.9800.23318.7870.8092.0002.5793 × 10−18.1612 × 10−28.7411 × 10−5
3.2850.0705.6200.8032.0007.7163 × 10−22.4415 × 10−22.6150 × 10−5
1010.9900.20918.8010.8062.0002.5817 × 10−18.1686 × 10−25.6393 × 10−5
3.2800.0625.6100.8022.0007.7041 × 10−22.4376 × 10−21.6828 × 10−5
1110.9980.19018.8120.8042.0002.5836 × 10−18.1747 × 10−23.7529 × 10−5
3.2750.0575.6020.8012.0007.6942 × 10−22.4345 × 10−21.1176 × 10−5
1211.0050.17418.8210.8032.0002.5852 × 10−18.1798 × 10−22.5653 × 10−5
3.2720.0525.5960.8012.0007.6860 × 10−22.4319 × 10−27.6268 × 10−6
Table 9. Dependence of the performance measures on parameter s ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 .
Table 9. Dependence of the performance measures on parameter s ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 .
sVav(1)Vav(2)SavLavDRSRR1RR2PL
67.6620.75013.5880.8601.9990.3380.3385.8271 × 10−4
7.7600.73913.5020.8591.9990.3330.3335.7427 × 10−4
77.6250.57613.7990.8461.9990.3470.2604.4808 × 10−4
7.8870.57513.5380.8461.9990.3460.2594.4733 × 10−4
87.4850.44714.0680.8361.9990.3590.2023.4767 × 10−4
7.8680.44913.6830.8361.9990.3600.2023.4895 × 10−4
97.2700.35014.3800.8281.9990.3740.1582.7185 × 10−4
7.7530.35113.8960.8281.9990.3750.1582.7264 × 10−4
107.0000.27514.7250.8222.0000.3920.1242.1406 × 10−4
7.5820.27414.1440.8222.0000.3910.1242.1338 × 10−4
Table 10. Dependence of the performance measures on parameter N ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
Table 10. Dependence of the performance measures on parameter N ; λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , r = 5 , s = 10 .
NVav(1)Vav(2)SavLavDRSRR1RR2PL
5010.9680.26218.7700.8072.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6320.7992.0002.0002.4457 × 10−22.0923 × 10−4
7010.9680.26218.7700.8102.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8022.0002.0002.4463 × 10−24.2171 × 10−5
9010.9680.26218.7700.8132.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8032.0002.0002.4463 × 10−24.2171 × 10−5
11010.9680.26218.7700.8162.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8042.0002.0002.4463 × 10−24.2171 × 10−5
13010.9680.26218.7700.8192.0000.2588.1520 × 10−21.4053 × 10−4
3.2910.0795.6330.8062.0002.0002.4463 × 10−24.2171 × 10−5
Table 11. Results of the task (38); λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , S = 22   , K 1 = 100 ,   K 2 = 200 ,   c r ( 1 ) = 50 ,   c r ( 2 ) = 100 ,   c c = 50 ,   c h = 35 ,   c d = 75 ,   c l = 200 ,   c w = 50 . The minimal value of TC is indicated in bold.
Table 11. Results of the task (38); λ = 8 , κ = 2 , μ 1 = 45 , μ 2 = 15 , τ = 1.5 , ν 1 = 2 , ν 2 = 8 , S = 22   , K 1 = 100 ,   K 2 = 200 ,   c r ( 1 ) = 50 ,   c r ( 2 ) = 100 ,   c c = 50 ,   c h = 35 ,   c d = 75 ,   c l = 200 ,   c w = 50 . The minimal value of TC is indicated in bold.
Total Cost
s1098765432
r
1896.90885.21875.95869.80867.73871.11881.95903.24939.62
870.62850.90835.46824.94820.45823.69837.22864.91912.74
2896.22884.54875.36869.44867.88872.30885.20910.49
872.35853.23838.68829.49827.00833.34851.80887.57
3899.59889.07881.50877.89879.71889.26910.23
879.47862.48850.85845.75849.10863.98895.31
4906.67898.44894.07895.01903.51923.23
892.02878.78872.33874.49888.29918.65
5917.85913.31914.16922.68942.6
910.98903.58905.28919.10950.08
6934.27935.39944.45965.31
938.43939.84954.14986.50
7958.00967.89990.1
977.88992.761026.87
8992.671016.55
1035.161071.21
91044.60
1120.16
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

Melikov, A.; Mirzayev, R.; Sztrik, J. Double-Sources Queuing-Inventory Systems with Finite Waiting Room and Destructible Stocks. Mathematics 2023, 11, 226. https://doi.org/10.3390/math11010226

AMA Style

Melikov A, Mirzayev R, Sztrik J. Double-Sources Queuing-Inventory Systems with Finite Waiting Room and Destructible Stocks. Mathematics. 2023; 11(1):226. https://doi.org/10.3390/math11010226

Chicago/Turabian Style

Melikov, Agassi, Ramil Mirzayev, and Janos Sztrik. 2023. "Double-Sources Queuing-Inventory Systems with Finite Waiting Room and Destructible Stocks" Mathematics 11, no. 1: 226. https://doi.org/10.3390/math11010226

APA Style

Melikov, A., Mirzayev, R., & Sztrik, J. (2023). Double-Sources Queuing-Inventory Systems with Finite Waiting Room and Destructible Stocks. Mathematics, 11(1), 226. https://doi.org/10.3390/math11010226

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