Next Article in Journal
A Lindley-Type Distribution for Modeling High-Kurtosis Data
Next Article in Special Issue
Game Theory and an Improved Maximum Entropy-Attribute Measure Interval Model for Predicting Rockburst Intensity
Previous Article in Journal
Operators Obtained by Using Certain Generating Function for Approximation
Previous Article in Special Issue
A Novel Method for Remaining Useful Life Prediction of Bearing Based on Spectrum Image Similarity Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of a Stochastic Inventory Model on Random Environment with Two Classes of Suppliers and Impulse Customers

1
Department of Mathematics, Alagappa University, Karaikudi 630003, India
2
Ramanujan Centre for Higher Mathematics, Alagappa University, Karaikudi 630004, India
3
Ramanujan Institute for Advanced Study in Mathematics, University of Madras, Chennai 600005, India
4
Department of Electronic Engineering, Kwangwoon University, Seoul 01897, Korea
5
Department of Information and Communication Engineering and Convergence Engineering for Intelligent Drone, Sejong University, Seoul 05006, Korea
6
Department of Computer Science and Engineering, Sejong University, Seoul 05006, Korea
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(13), 2235; https://doi.org/10.3390/math10132235
Submission received: 29 April 2022 / Revised: 14 June 2022 / Accepted: 22 June 2022 / Published: 26 June 2022

Abstract

:
This paper explores the random environment with two classes of suppliers and impulse customers. The system’s greatest inventory size is S, and it has an infinitely large orbit. In this case, there are two categories of suppliers: temporary suppliers and regular suppliers. Whenever the inventory approaches r, we place on order Q 1 (= S r ) unit items to a temporary supplier. Similarly, when the inventory level drops to s (< Q 1 < r ), we place an order for Q 2 (= S s > s + 1 ) units of items to our regular supplier. Two types of suppliers’ lead times are considered to be exponentially distributed. Here, the customers who arrive from different states of the random environment (RE) are followed by the Markovian arrival process. If there is no inventory in the system when the customer arrives, they are automatically assigned to an orbit. The model was examined in steady state by using the matrix-analytic approach. Finally, the numerical examples for our structural model are discussed.

1. Introduction

Inventory management faces many replenishment difficulties. If the supplier’s activities are affected during replenishment, customers may be unable to get vital supplies. As a result, we use two suppliers, referred to in this paper as “regular suppliers” and “temporary suppliers”. For instance, let’s consider textile shops. The system always gives the order to the regular supplier, but the temporary supplier supplies trending dresses, so the system decides to make a trial for a temporary supplier at the same time to ensure the system doesn’t lose the regular supplier. Good connections with servers and customers are among these objectives, and customers should avoid being sent back due to inventory shortage and helping to minimise the total cost. Following that, all arriving customers have the option of purchasing or not purchasing their items after the detailed inquiries in N different states of random environment (RE) in this paper.
Soujanyaa and Vijaya Laxmi [1] investigated perishable inventory systems with dual supply chains and negative customers. They considered the lead time and perishable rate to be exponentially distributed, the arrival process followed by Poisson, and found a limiting time distribution. Finally, they demonstrate numerical examples for cost rate function and system performance parameters. Yang and Tseng [2] considered a three-echelon inventory model with a backorder and permissible delay in payments under a controllable lead time. The queueing-inventory system literature has more ordering policies such as [3,4,5,6].
The system’s behavior is determined by the current state of the Markovian random process with a finite state space. The technique of this operation in a fixed state is called a “random environment”. The system makes an effort to avoid crowds while also being mindful of not losing customers. The system also has several modalities for receiving the demand. This effort helps to satisfy all the customers without any customer accumulation. Kim and Dudin [7] considered a multi-server and impatient customer in the queueing system with a random environment. The service time is a phase-type distribution. They derived ergodicity conditions for this system with an infinite buffer. Some numerical results are presented. Kim et al. [8] investigated a multi-server queueing system and two types of customers with infinite buffers. The system is operating in a random environment, and the arrival stream is marked as a Markovian arrival process. Non-preemptive priority is given to type 1 customers over type 2 customers. The service time for two types of customers has a phase-type distribution with different parameters. The correlation in the arrival process is depicted numerically in the section.
In this article, we consider the arrival process as MAP under the fixed state of a random environment. Whenever the states of a random environment are changed, the parameters of the arrival flow and retrying flow are changed as well. Here, we give some reviews about the random environment (see [9,10,11,12,13,14]). In the 1970s, Neuts introduced the Versatile Markovian Point Process (VMPP) as the Markovian arrival process. The VMPP had extensive notations when it was first developed; however, it was considerably simplified by Lucantoni [15], and this process has been known as the Markovian arrival process ever since. The Markovian arrival process was greatly vetted in the literature [16,17,18,19,20,21]. See Krishnamurthy et al. [22] for details on the batch Markovian arrival process in the queueing-inventory system.
In this paper, all the arriving customers go into orbit at inventory level zero, and the arriving customers from orbit may or may not buy the product. For example, customers wait for the trending dresses at a nearby textile shop. At the textile shop, after coming across the trending dresses, the customer went to the shop, but the dress colour or size did not fit the customer. At that point, the customer didn’t buy the product and left the system. The responses to the last poll inspired our investigation, and there is no research into random environments with two supply chains to the best of our knowledge.

1.1. Research Gap

The authors above work with RE in the queuing system. This article examines RE in the inventory. The authors above handle the two-echelon chain, and both orders are obtained. In this article, we place an order with a temporary supplier, in case we have not received the order from the temporary supplier when the inventory level has reached some reordering point, we cancel that order and place the order with our regular supplier.

1.2. The Perspective of This Work

The model representation of the stochastic inventory model with two classes of suppliers and impulse customers on RE is portrayed in Section 2. In Section 3 the model analysis is portayed. The joint probability distribution of the number of customers in the orbit of infinite size, and the inventory level is evaluated in Section 4. We measure some important system peculiarities and construct the cost function in Section 5. Some numerical examples are provided in Section 6 and, finally, we describe a conclusion in Section 7.

2. Mathematical Framework of the Model

In this section, we discuss the random environment in the stochastic inventory system with two classes of suppliers named “temporary suppliers” and “regular suppliers”. In this system, an infinite-sized orbit is added, and here there is no waiting space attached to the system. The system has a maximum S units of items. Whenever the inventory level approaches r, we place an order with the temporary supplier for Q 1 items, where Q 1 = S r . Similarly, whenever the inventory level reaches s ( s < Q 1 < r ) , we will order Q 2 = S s ( > s + 1 ) items to the regular supplier and immediately cancel the order for Q 1 items from the temporary supplier. Both lead times of temporary and regular suppliers are exponentially distributed with β 1 and β 2 respectively.
The customers arrive according to a Markovian arrival process (MAP) from N different states of RE. Customers who require a single item have the option to choose any arrival mode (state) of the RE. The behaviour of the model depends on the state of the RE. The RE is used on the stochastic process J 2 ( t ) , which is an irreducible continuous time Markov chain with the state space { 1 , 2 , , N } and the infinitesimal generator H. The stationary row vector η 1 of the RE is to be obtained by using η 1 H = 0 , η 1 e = 1 . The underlying Markov chain J 4 ( t ) of the MAP has a generator D ( j 2 ) , which is a square matrix of dimension m with D ( j 2 ) = D 0 ( j 2 ) + D 1 ( j 2 ) in a fixed state j 2 ( = 1 , 2 , , N ) . In this case, D 0 ( j 2 ) denotes no arrival matrix of size m and D 1 ( j 2 ) denotes an arrival matrix of size m under fixed state j 2 . The arrival process of customer representation is ( D 0 ( j 2 ) , D 1 ( j 2 ) ) under fixed state j 2 . The average arrival rate λ ( j 2 ) of a customer is defined by λ ( j 2 ) = η 2 ( j 2 ) D 1 ( j 2 ) e , where stationary row vector η 2 ( j 2 ) of size 1 × m is to be obtained by using η 2 ( j 2 ) D ( j 2 ) = 0 and η 2 ( j 2 ) e = 1 under fixed state j 2 .
Not all incoming customers buy the product, those customers are referred to as impulse customers. Assume that customers buy the product with probability p ( j 2 ) and the complementary probability q ( j 2 ) under fixed state j 2 . The service process of the system is assumed to be instantaneous when the inventory level is positive. When the arriving customers find the stock stage is empty, they should enter into orbit. Even though the mode of approach to entering the system is different for each customer, once the customer enters the orbit, they will be considered retrial customers. Retrial customers have the option of choosing another arriving mode of RE during the retry without considering the previous choice of different states of RE. The retrial customers will get into service only if the stock stage is positive. The time between two successive retrials is exponentially distributed with a transition rate θ ( j 2 ) .

3. Model Analysis

In this sector, we construct the transition rate matrix on the stochastic inventory system. The Markov process of the form { ( J 1 ( t ) , J 2 ( t ) , J 3 ( t ) , J 4 ( t ) ) , t 0 } with state space
E = { ( j 1 , j 2 , j 3 , j 4 ) : j 1 0 ; 1 j 2 N ; 0 j 3 S ; 1 j 4 m } ,
where
  • J 1 ( t ) : The number of customers in the orbit of infinite size waiting place at time t.
  • J 2 ( t ) : The state of random environment at time t.
  • J 3 ( t ) : The number of items in the inventory at time t.
  • J 4 ( t ) : Phase of the arrival process at time t.
The process’s infinitesimal generator P is generated by
Mathematics 10 02235 i001
Mathematics 10 02235 i002
Here, j 2 = 1 , 2 , , N :
Mathematics 10 02235 i003
Mathematics 10 02235 i004
Here, j 2 = 1 , 2 , , N :
Mathematics 10 02235 i005
Mathematics 10 02235 i006
Mathematics 10 02235 i007
For k = 1 and 2 and j 2 = 1 , 2 , , N :
Mathematics 10 02235 i008
where b 1 = p ( j 2 ) D 1 ( j 2 ) and d i k = D 0 ( j 2 ) + δ ¯ i 0 q ( j 2 ) D 1 ( j 2 ) ( δ i 0 β 2 + δ i 1 β 2 + δ i 2 β 1 + δ k 2 δ ¯ i 0 θ ( j 2 ) + δ k 2 δ i 0 q ( j 2 ) θ ( j 2 ) ) I m , i = 0 , 1 , 2 and 3 .
It may be noted that the A 01 , A 10 , A 00 , and A 11 are all square matrices of dimension k 1 and the matrices A ( j 2 ) , B ( j 2 ) , C 1 ( j 2 ) and , C 2 ( j 2 ) are dimension k 2 , where j 2 = 1 , 2 , , N .

4. Joint Probability Distribution under Steady State

We give the condition for stability through a theorem on the stochastic inventory system. Let us consider G = A 10 + A 11 + A 01 . It is easy to visualize that G is given by
Mathematics 10 02235 i009
Here j 2 = 1 , 2 , , N :
Mathematics 10 02235 i010
where b 2 = b 1 + p ( j 2 ) θ ( j 2 ) I m and f i = d i 2 + δ i 0 D 1 ( j 2 ) + q ( j 2 ) θ ( j 2 ) I m , i = 0 , 1 , 2 and 3 .
Here, Y denotes the steady-state probability vector of G.
where
Y = ( Y 1 , Y 2 , , Y N ) ,
Y j 2 = ( Y ( j 2 , 0 ) , Y ( j 2 , 1 ) , , Y ( j 2 , S ) ) , 1 j 2 N .
The below system of equations are derived from Y G = 0 :
j 2 = 1 N Y j 2 ( δ j 2 1 F ( 1 ) + H j 2 1 ) = 0 j 2 = 1 N Y j 2 ( δ j 2 2 F ( 2 ) + H j 2 2 ) = 0 j 2 = 1 N Y j 2 ( δ j 2 N F ( N ) + H j 2 N ) = 0 .
The steady-state probability vector Y i is obtained from the above system of equations and the normalizing condition ( Y e = 1 ).

4.1. Stability Condition

We give the condition for stability through theorem on the stochastic inventory system.
Theorem 1.
The stochastic inventory system under study is stable if and only if
j 2 = 1 N Y j 2 A ( j 2 ) e < j 2 = 1 N Y j 2 B ( j 2 ) e .
Proof. 
From the standard results of Neuts [23] on the positive recurrence of P we have
Y A 01 e < Y A 10 e ,
and by applying the structure of the matrices A 10 and A 01 and Y the declared results follow. □
It can be seen from the structure of the rate matrix P and from the Theorem 1, that the Markov process { ( J 1 ( t ) , J 2 ( t ) , J 3 ( t ) , J 4 ( t ) ) , t 1 } with the state space E is regular.

4.2. Steady-State Probability Vector

In this section, we calculate steady-state probability. Let Θ be the steady-state probability vector of the matrix P and it is obtained from Θ P = 0 , Θ e = 1 , where
Θ = ( Θ 0 , Θ 1 , Θ 2 , ) ,
where
Θ j 1 = ( Θ ( j 1 , 1 ) , Θ ( j 1 , 2 ) , Θ ( j 1 , 3 ) , , Θ ( j 1 , N ) ) , j 1 = 0 , 1 , 2 , , Θ ( j 1 , j 2 ) = ( Θ ( j 1 , j 2 , 0 ) , Θ ( j 1 , j 2 , 1 ) , Θ ( j 1 , j 2 , 2 ) , , Θ ( j 1 , j 2 , S ) ) , j 1 = 0 , 1 , 2 , , j 2 = 1 , 2 , , N , Θ ( j 1 , j 2 , j 3 ) = ( Θ ( j 1 , j 2 , j 3 , 1 ) , Θ ( j 1 , j 2 , j 3 , 2 ) , , Θ ( j 1 , j 2 , j 3 , m ) ) , j 1 = 0 , 1 , 2 , , j 2 = 1 , 2 , , N , j 3 = 0 , 1 , 2 , , S .
The relation Θ P = 0 leads to the following system of equations:
Θ 0 A 00 + Θ 1 A 10 = 0
Θ j 1 1 A 01 + Θ j 1 A 11 + Θ j 1 + 1 A 10 = 0 , ( j 1 1 ) .
The Markov process { ( J 1 ( t ) , J 2 ( t ) , J 3 ( t ) , J 4 ( t ) ) , t 0 } on the state space E and the limiting distribution
Θ ( j 1 , j 2 , j 3 , j 4 ) = lim t P r [ J 1 ( t ) = j 1 , J 2 ( t ) = j 2 , J 3 ( t ) = j 3 , J 4 ( t ) = j 4 | J 1 ( 0 ) , J 2 ( 0 ) , J 3 ( 0 ) , J 4 ( 0 ) ] ,
where Θ ( j 1 , j 2 , j 3 , j 4 ) is the steady-state probability for the state ( j 1 , j 2 , j 3 , j 4 ) , exists and is independent of the initial state.
Theorem 2.
When the stability condition { ( J 1 ( t ) , J 2 ( t ) , J 3 ( t ) , J 4 ( t ) ) , t 0 } holds well, the steady-state probability vector Θ is given by
Θ i = Θ 0 R i , i 1 ,
where the rate matrix R , it is defined by
Mathematics 10 02235 i011
and satisfies the matrix quadratic equation
R 2 A 10 + R A 11 + A 01 = 0 .
Proof. 
The proof follows from the well-known result on matrix geometric property given by Neuts [23]. □
Theorem 3.
The stationary probability vector Θ 0 is the unique solution of the system
Θ 0 ( A 00 + R A 10 ) = 0 ,
and subject to the normalising condition
Θ 0 ( I R ) 1 e = 1 .
Proof. 
The proof is based on Neuts’ [23] well-known matrix geometric property. And see [24] more details about the property. □
Once R is found, the boundary probability vector Θ 0 is computed by using the Equations (7) and (8), and the probability vectors Θ i , i 1 can be obtained from Equation (5).

5. A Few Significant System Peculiarities

In the segment, we examine a few significant peculiarities measures.
  • Mean inventory level, η I is given by
    η I = j 1 = 0 j 2 = 1 N j 3 = 1 S j 3 Θ ( j 1 , j 2 , j 3 ) e .
  • Mean reorder for temporary supplier, η R 1 is given by
    η R 1 = j 1 = 0 j 2 = 1 N p ( j 2 ) D 1 ( j 2 ) [ Θ ( j 1 , j 2 , r + 1 ) ] e + j 1 = 1 j 2 = 1 N p ( j 2 ) θ ( j 2 ) [ Θ ( j 1 , j 2 , r + 1 ) ] e .
  • Mean reorder for regular supplier, η R 2 is given by
    η R 2 = j 1 = 0 j 2 = 1 N p ( j 2 ) D 1 ( j 2 ) [ Θ ( j 1 , j 2 , s + 1 ) ] e + j 1 = 1 j 2 = 1 N p ( j 2 ) θ ( j 2 ) [ Θ ( j 1 , j 2 , s + 1 ) ] e .
  • Mean number of customers in the orbit
    η O C = Θ 0 R ( I R ) 2 e .
  • Mean loss rate of arrival of impulse customers
    η L C = j 1 = 0 j 2 = 1 N j 3 = 1 S q ( j 2 ) D 1 ( j 2 ) Θ ( j 1 , j 2 , j 3 ) e + j 1 = 1 j 2 = 1 N j 3 = 1 S q ( j 2 ) θ ( j 2 ) Θ ( j 1 , j 2 , j 3 ) e .
  • Expected number of time the replenishment is to be done from temporary supplier
    η T S = j 1 = 0 j 2 = 1 N j 3 = s + 1 r β 1 Θ ( j 1 , j 2 , j 3 ) e .
  • Expected number of time the replenishment is to be done from regular supplier
    η R S = j 1 = 0 j 2 = 1 N j 3 = 0 s β 2 Θ ( j 1 , j 2 , j 3 ) e .
  • Overall retrial rate
    η O R = j 1 = 1 j 2 = 1 N θ ( j 2 ) Θ ( j 1 , j 2 ) e .
  • Success retrial rate
    η S R = j 1 = 1 j 2 = 1 N j 3 = 1 S θ ( j 2 ) Θ ( j 1 , j 2 , j 3 ) e .
  • Fractional success retrial rate
    η F S R = η S R η O R .

Construction of the Cost Feature

To construct the expected total cost function per unit time is provided by
C ( S , s ) = c h η I + c w η O C + c T S η R 1 + c r s η R 2 + c i l η L C , where c h : The inventory carrying cost per unit time.
c w : Waiting cost of a customer in the orbit per unit time.
c t s : Setup cost per order for temporary supplier.
c r s : Setup cost per order for regular supplier.
c i l : Lost cost of a impulse customer per unit time.

6. Numerical Analysis

We give a few descriptive numerical examples that expose the convexity of the expected cost rate and consider that orders are received with 2 different states of RE; these are confined by the infinitesimal generator
H = 0.10 0.10 0.01 0.01 .
The MAP for appearance of customers are
  • Hyper-exponential(HEX):
    D 0 ( 1 ) = 1.90 0 0 0.19 , D 1 ( 1 ) = 1.710 0.190 0.171 0.019 ;
    D 0 ( 2 ) = 2.85 0 0 0.285 , D 1 ( 2 ) = 2.565 0.285 0.2565 0.0285 .
  • Negative Correlation (NC):
    D 0 ( 1 ) = 1.00222 1.00222 0 0 1.00222 0 0 0 225.75 , D 1 ( 1 ) = 0 0 0 0.01002 0 0.99220 223.49250 0 2.25750 ;
    D 0 ( 2 ) = 2.35 2.35 0 0 2.35 0 0 0 3.5 , D 1 ( 2 ) = 0 0 0 0.0235 0 2.3265 3.465 0 0.035 .
  • Positive Correlation (PC):
    D 0 ( 1 ) = 1.00222 1.00222 0 0 1.00222 0 0 0 225.75 , D 1 ( 1 ) = 0 0 0 0.99220 0 0.01002 2.25750 0 223.49250 ; D 0 ( 2 ) = 2.35 2.35 0 0 2.35 0 0 0 3.5 , D 1 ( 2 ) = 0 0 0 2.3265 0 0.0235 0.035 0 3.465 .
The demand process has negative (positive) correlated arrival with coefficient of variance c v a r ( 1 ) = 2 λ ( 1 ) η 2 ( 1 ) ( D 0 ( 1 ) ) 1 e 1 = 1.986 ( 1.986 ) and coefficient of correlation c c o r ( 1 ) = ( λ ( 1 ) η 2 ( 1 ) ) ( D 0 ( 1 ) ) 1 D 1 ( 1 ) ( D 0 ( 1 ) ) 1 ) e 1 ) / c v a r ( 1 ) = 0.4889 ( 0.4889 ) with an arrival rate λ ( 1 ) = 1.0 .
The demand process has negative (positive) correlated arrival with coefficient of variance c v a r ( 2 ) = 2 λ ( 2 ) η 2 ( 2 ) ( D 0 ( 2 ) ) 1 e 1 = 0.9342 ( 0.9342 ) and coefficient of correlation c c o r ( 2 ) = ( λ ( 2 ) η 2 ( 2 ) ) ( D 0 ( 2 ) ) 1 D 1 ( 2 ) ( D 0 ( 2 ) ) 1 ) e 1 ) / c v a r ( 2 ) = 0.2595 ( 0.2595 ) with an arrival rate λ ( 2 ) = 1.7594 .
Some discussion about numerical examples of our perspective model and its obtainment with the parameters p ( 1 ) = 0.8 , q ( 1 ) = 1 p ( 1 ) , β 1 = 0.2 , λ 1 ( 1 ) = 1 , β 2 = 0.9 , θ ( 1 ) = 5 , p ( 2 ) = 0.6 , q ( 2 ) = 1 p ( 2 ) , θ ( 2 ) = 4.8 , λ 1 ( 2 ) = 1.5 , R = 2 and cost values are c h = 0.1 , c w = 10 , c t s = 28 , c r s = 15 , c i l = 10 are listed in the Results and Discussion section.

Results and Discussion

  • We discuss the behavior of the cost function of two variables, C ( S , s ) , under hyper-exponential distribution. The values are divulged in bold in each column to indicate the minimum cost rate, whereas the least cost rate is specified in each row by underlining the values. As a result, a value (bold and underlined) represents the local minimum of the function C ( S , s ) . At S * = 28 and s * = 5, the optimal cost value C * ( S , s ) = 5.7025. The function C ( S , s ) is convex, as shown in Table 1 and Figure 1. Figure 2 depicts a contour plot of the total cost function, which also demonstrates that the function C ( S , s ) is convex.
  • Figure 3 and Figure 4 show that the mean reorder for temporary supplier ( η R 1 ) and regular supplier ( η R 2 ) compare with s and r respectively. Figure 3 demonstrates that η R 1 decreases and η R 2 increases, whenever s increases. Figure 4 shows that η R 1 is increased and regular suppliers η R 2 is decreased whenever r is increased.
  • Figure 5 and Figure 6 show that the mean number of times the replenishment is to be done from a temporary supplier ( η T S ) and regular supplier ( η R S ) . Here, Figure 5 shows that η T S decreases and η R S increases whenever s increases. Figure 6 shows that when r increases, η T S increases and η R S decreases.
  • Figure 7 compares the lead time rate of temporary ( β 1 ) and regular ( β 2 ) suppliers with their total expected cost value. Here, the total expected cost value decreases whenever β 1 and β 2 rates are increased.
  • Table 2 shows that when increasing the lead time rate of the temporary supplier ( β 1 ) increases the optimal cost value C * ( S , s ) , whereas when increasing the lead time rate of the regular supplier ( β 2 ) decreases the C * ( S , s ) . The C * ( S , s ) for the lead time rates of two suppliers decreases as the maximum inventory level (S) increases.
  • When the setup cost for temporary supplier ( c t s ) and regular supplier ( c r s ) increases, the value of C * ( S , s ) increases in Table 3. When S increases, the value of C * ( S , s ) decreases. When the values of c h , c w , c i l increase, so does the value of C * ( S , s ) as shown in Table 4.
  • As shown in Table 5, when s and s values rise, the mean inventory level ( η I ) rises, but the mean number of customers in the orbit ( η O C ) and the mean loss rate of arrival of impulse customers ( η L C ) decreases.

7. Conclusions

In this paper, we discussed two classes of suppliers in the stochastic inventory system on the RE. These two classes of suppliers have a huge impact on the management of the inventory system and also reduce replenishment difficulties. We showed the minimised total expected cost rate (refer to Table 1 and Figure 1) and showed it via a contour plot (refer Figure 2). When the lead time rate of both suppliers is smaller, the total expected cost value of the regular supplier is greater than the total expected cost value of the temporary supplier. Moreover, when the lead time rate of both suppliers goes up, the total expected cost value of the regular supplier is smaller than the total expected cost value of the temporary supplier (refer Table 2). Similarly, when the setup cost value is small, the total expected cost value of the temporary supplier is smaller than the total expected cost value of the regular supplier. Furthermore, when the setup cost value goes up, the total expected cost value of the temporary supplier is greater than the total expected cost value of the regular supplier (refer Table 3). Finally, we showed that both the value of the setup cost and the lead time rate are small, the best supplier is a temporary supplier, and both have higher value, making the best supplier a regular supplier.

7.1. Limitations

  • This study deals with impulse customers and the probability of customers who may buy an item, p j 2 and the complementary probability q j 2 . Assuming the probability q j 2 value is zero, in this case, not all incoming customers will purchase the item, which does not always happen in real-life situations.
  • The sum of the fixed probability distribution values must be one.

7.2. Future Directions

  • We will discuss multi-server with phase-type distribution for service.
  • We plan to investigate multi-suppliers.
  • We plan to study RE use with payment mode.

Author Contributions

Conceptualization, N.A. and V.V.; data curation, N.A., V.V., and K.J.; formal analysis, N.A. and V.V.; funding acquisition, H.-K.S., G.P.J., and H.M.; methodology, S.A. and K.J.; project administration, H.-K.S. and H.M.; resources, B.S. and H.M.; software, S.A.; visualization, B.S. and G.P.J.; writing—original draft, N.A.; writing—review & editing, G.P.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1A6A1A03038540) and National Research Foundation of Korea (NRF) grant funded by the Korean government, Ministry of Science and ICT (MSIT) (2021R1F1A1046339) and by a grant (20212020900150) from “Development and Demonstration of Technology for Customers Bigdata-based Energy Management in the Field of Heat Supply Chain” funded by Ministry of Trade, Industry and Energy of Korean government.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Anbazhagan and Amutha would like to thank RUSA Phase 2.0(F 24-51/2014-U), DST-FIST( SR/FIST/MS-I/2018/17), DST-PURSE 2nd Phase programme (SR/PURSE Phase 2/38), Government of India.

Conflicts of Interest

The authors declare no conflict of interest.

Notations and Abbreviations

The following Notations and Abbreviations are used in this manuscript:
[ A ] i j The element of submatrix at (i,j) the position of A.
0 Zero matrix.
I Identity matrix of appropriate dimension.
I m Identity matrix of dimension m.
e A column vector of 1’s appropriate dimension.
A B Kronecker product of matrices A and B.
A B Kronecker sum of matrices A and B.
RERandom Envirnonment.
MAPMarkovian Arrival Process.
δ ¯ i j 1 δ i j .
k 1 N ( S + 1 ) m
k 2 ( S + 1 ) m

References

  1. Soujanyaa, M.L.; Laxmi, P.V. Analysis on Dual Supply Inventory Model having Negative Arrivals and Finite Life Time Inventory. Reliab. Theory Appl. 2021, 16, 295–301. Available online: http://www.gnedenko.net/RTA/index.php/rta/article/view/778 (accessed on 28 April 2022).
  2. Yang, M.F.; Tseng, W.C. Three-echelon inventory model with permissible delay in payments under controllable lead time and backorder consideration. Math. Probl. Eng. 2014, 2014, 809149. [Google Scholar] [CrossRef]
  3. Anbazhagan, N.; Arivarignan, G. Two-commodity inventory system with individual and joint ordering policies. Int. J. Manag. Syst. 2003, 19, 129–144. [Google Scholar]
  4. Sivakumar, B.; Anbazhagan, N.; Arivarignan, G. Two-commodity inventory system with individual and joint ordering policies and renewal demands. Stoch. Anal. Appl. 2007, 25, 1217–1241. [Google Scholar] [CrossRef]
  5. Anbazhagan, N.; Kumaresan, V. Coordinated and joint ordering policies for two commodity inventory system. Iran. J. Oper. Res. 2012, 3, 94–105. [Google Scholar]
  6. Anbazhagan, N.; Arivarignan, G. Two commodity continuous review inventory system with coordinated reorder policy. Int. J. Inf. Manag. Sci. 2000, 11, 19–30. [Google Scholar]
  7. Kim, C.; Dudin, S. Analysis of BMAP(r)/M(r)/N(r) Type Queueing System Operating in Random Environment. J. Korean Inst. Ind. Eng. 2016, 42, 30–37. [Google Scholar]
  8. Kim, C.; Dudin, A.; Dudin, S.; Dudina, O. Analysis of an MMAP/PH1,PH2/N/ queueing system operating in a random environment. Int. J. Appl. Math. Comput. Sci. 2014, 24, 485–501. [Google Scholar] [CrossRef] [Green Version]
  9. Dudin, A.N.; Nazarov, A.A. The MMAP/M/R/0 queueing system with reservation of servers operating in a random environment. Probl. Inf. Transm. 2015, 51, 289–298. [Google Scholar] [CrossRef]
  10. Kim, C.; Dudin, A.; Dudina, O.; Kim, J. Queueing system operating in random environment as a model of a cell operation. Ind. Eng. Manag. Syst. 2016, 15, 131–142. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, C.S.; Dudin, A.; Klimenok, V.; Khramova, V. Performance analysis of multi-server queueing system operating under control of a random environment. In Trends in Telecommunications Technologies; Christos, J.B., Ed.; IntechOpen: London, UK, 2010; ISBN 978-953-307-072-8. [Google Scholar]
  12. Kim, C.S.; Dudin, A.; Klimenok, V.; Khramova, V. Erlang loss queueing system with batch arrivals operating in a random environment. Comput. Oper. Res. 2009, 36, 674–697. [Google Scholar] [CrossRef]
  13. Kim, C.; Klimenok, V.; Lee, S.; Dudin, A. The BMAP/PH/1 retrial queueing system operating in random environment. J. Stat. Plan. Inference 2007, 137, 3904–3916. [Google Scholar] [CrossRef]
  14. Dudin, A.; Kim, C.S.; Dudin, S.; Dudina, O. Priority Retrial Queueing Model Operating in Random Environment with Varying Number and Reservation of Servers. Appl. Math. Comput. 2015, 269, 674–690. [Google Scholar] [CrossRef]
  15. Lucantoni, D. New results on the single server queue with a batch Markovian arrival process. Commun. Stat. Stoch. Model. 1991, 7, 1–46. [Google Scholar] [CrossRef]
  16. Nair, S.S.; Jose, K.P.A. PH Distributed Production Inventory Model with Different Modes of Service and MAP Arrivals. In Applied Probability and Stochastic Processes; Infosys Science Foundation Series; Springer: Singapore, 2020. [Google Scholar]
  17. Punalal, J.; Babu, S. The MAP/(PH, PH, PH)/1 Model with Self-Generation of Priorities, Customer Induced Interruption and Retrial of Customers. In Applied Probability and Stochastic Processes; Infosys Science Foundation Series; Springer: Singapore, 2020. [Google Scholar]
  18. Lee, S.; Dudin, S.; Dudina, O.; Kim, C.; Klimenok, V. A Priority Queue with Many Customer Types, Correlated Arrivals and Changing Priorities. Mathematics 2020, 8, 1292. [Google Scholar] [CrossRef]
  19. Klimenok, V.; Dudin, A.; Dudina, O.; Kochetkova, I. Queuing System with Two Types of Customers and Dynamic Change of a Priority. Mathematics 2020, 8, 824. [Google Scholar] [CrossRef]
  20. Kim, C.; Dudin, S.; Taramin, O.; Baek, J. Queueing system MAP|PH|N|N+R with impatient heterogeneous customers as a model of call center. Appl. Math. Model. 2013, 37, 958–976. [Google Scholar] [CrossRef]
  21. He, Q.-M. The Versatility of MMAP[K] and the MMAP[K]/G[K]/1 Queue. Queueing Syst. 2001, 38, 397–418. [Google Scholar] [CrossRef]
  22. Krishnamoorthy, A.; Joshua, A.N.; Kozyrev, D. Analysis of a Batch Arrival, Batch Service Queuing-Inventory System with Processing of Inventory While on Vacation. Mathematics 2021, 9, 419. [Google Scholar] [CrossRef]
  23. Neuts, M.F. Matrix-Geometric Solutions in Stochastic Models: An Algorithmic Approach; Dover Publication: New York, NY, USA, 1994. [Google Scholar]
  24. Latouche, G.; Ramaswami, V. Introduction to Matrix Analytic Methods in Stochastic Modeling; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 1999. [Google Scholar] [CrossRef]
Figure 1. Total expected cost rate as a function of s and S.
Figure 1. Total expected cost rate as a function of s and S.
Mathematics 10 02235 g001
Figure 2. Contour plot of total expected cost rate.
Figure 2. Contour plot of total expected cost rate.
Mathematics 10 02235 g002
Figure 3. Mean reorder rate for two suppliers vs. s.
Figure 3. Mean reorder rate for two suppliers vs. s.
Mathematics 10 02235 g003
Figure 4. Mean reorder rate for two suppliers vs. r.
Figure 4. Mean reorder rate for two suppliers vs. r.
Mathematics 10 02235 g004
Figure 5. η T S and η R S vs. s.
Figure 5. η T S and η R S vs. s.
Mathematics 10 02235 g005
Figure 6. η T S and η R S vs. r.
Figure 6. η T S and η R S vs. r.
Mathematics 10 02235 g006
Figure 7. The total expected cost value for different values of β 1 and β 2 .
Figure 7. The total expected cost value for different values of β 1 and β 2 .
Mathematics 10 02235 g007
Table 1. Total expected cost rate as a function of s and S.
Table 1. Total expected cost rate as a function of s and S.
S/s2345678
245.82515.78875.76915.76735.78485.82395.8884
265.78005.74625.72605.72005.72905.75445.7985
285.76095.73025.71065.70255.70635.72305.7540
305.76345.73575.71735.70835.70905.72005.7422
325.78385.75905.74215.73295.73175.73895.7551
345.81945.79735.78185.77295.77075.77535.7872
Table 2. The optimal cost value for different values of lead times.
Table 2. The optimal cost value for different values of lead times.
S = 14S = 16S = 18S = 20
β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s )
0.16.67300.76.81890.16.35260.76.48030.16.12780.76.23800.15.96900.76.0646
0.26.77300.86.79550.26.43170.86.45520.26.19050.86.21320.26.02020.86.0414
0.36.83340.96.77300.36.47400.96.43170.36.22070.96.19050.36.04340.96.0202
0.46.86891.06.75160.46.49521.06.40990.46.23361.06.16970.46.05271.06.0009
0.56.88921.16.73150.56.50481.16.38970.56.23811.16.15050.56.05581.15.9832
0.66.90021.26.71250.66.50811.26.37100.66.23851.26.13300.66.05621.25.9671
S = 22S = 24S = 26S = 28
β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s ) β 1 C * ( S , s ) β 2 C * ( S , s )
0.15.85830.75.94310.15.78390.75.86190.15.73770.75.81290.15.71420.75.7903
0.25.90240.85.92170.25.82510.85.84260.25.78000.85.79560.25.76090.85.7748
0.35.92300.95.90240.35.84640.95.82510.35.80440.95.78000.35.79020.95.7609
0.45.93201.05.88480.45.85731.05.80930.45.81851.05.76590.45.80831.05.7484
0.55.93601.15.86890.55.86351.15.79500.55.82771.15.75310.55.82071.15.7371
0.65.93791.25.85430.65.86781.25.78200.65.83451.25.74150.65.83001.25.7268
Table 3. The optimal cost value for different setup cost values of two suppliers.
Table 3. The optimal cost value for different setup cost values of two suppliers.
S = 14S = 16S = 18S = 20
c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s )
256.4712106.5893256.1603106.2936255.9435106.0855255.7932105.9397
266.5718116.6260266.2508116.3212266.0258116.1065265.8688115.9558
276.6724126.6628276.3412126.3488276.1082126.1275275.9445125.9719
286.7730136.6995286.4317136.3765286.1905136.1485286.0202135.9880
296.8736146.7363296.5221146.4041296.2729146.1695296.0958146.0041
306.9742156.7730306.6126156.4317306.3552156.1905306.1715156.0202
317.0749166.8098316.7030166.4593316.4376166.2115316.2472166.0363
327.1755176.8465326.7935176.4869326.5199176.2325326.3228176.0524
337.2761186.8832336.8839186.5145336.6023186.2535336.3985186.0685
347.3767196.9200346.9744196.5421346.6846196.2745346.4742196.0846
357.4773206.9567357.0648206.5697356.7670206.2955356.5498206.1007
S = 22S = 24S = 26S = 28
c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s ) c t s C * ( S , s ) c r s C * ( S , s )
255.6923105.8403255.6297105.7769255.5974105.7425255.5896105.7316
265.7624115.8527265.6948115.7866265.6582115.7500265.6467115.7375
275.8324125.8651275.7600125.7962275.7191125.7575275.7038125.7433
285.9024135.8775285.8251135.8058285.7800135.7650285.7609135.7492
295.9724145.8900295.8903145.8155295.8409145.7725295.8180145.7551
306.0424155.9024305.9554155.8251305.9018155.7800305.8751155.7609
316.1124165.9148316.0206165.8348315.9626165.7875315.9323165.7668
326.1825175.9273326.0857175.8444326.0235175.7950325.9894175.7727
336.2525185.9397336.1508185.8540336.0844185.8025336.0465185.7785
346.3225195.9521346.2160195.8637346.1453195.8100346.1036195.7844
356.3925205.9645356.2811205.8733356.2062205.8175356.1607205.7903
Table 4. Total expected cost values for different cost values.
Table 4. Total expected cost values for different cost values.
c h = 0.1 c h = 0.2
c w 8910111289101112
c i l 54.65604.65664.65724.65784.65846.36736.36796.36856.36916.3697
64.88064.88124.88184.88244.88306.59196.59256.59316.59366.5942
75.10515.10575.10635.10695.10756.81646.81706.81766.81826.8188
85.32975.33035.33095.33155.33217.04107.04167.04227.04287.0433
95.55435.55495.55545.55605.55667.26557.26617.26677.26737.2679
105.77885.77945.78005.78065.78127.49017.49077.49137.49197.4925
116.00346.00406.00466.00516.00577.71467.71527.71587.71647.7170
126.22796.22856.22916.22976.23037.93927.93987.94047.94107.9416
136.45256.45316.45376.45436.45488.16388.16448.16498.16558.1661
146.67706.67766.67826.67886.67948.38838.38898.38958.39018.3907
156.90166.90226.90286.90346.90408.61298.61358.61418.61468.6152
c h = 0.3 c h = 0.4
c w 8910111289101112
c i l 58.07868.07928.07988.08048.08109.78999.79059.79119.79169.7922
68.30328.30378.30438.30498.305510.014410.015010.015610.016210.0168
78.52778.52838.52898.52958.530110.239010.239610.240210.240810.2414
88.75238.75298.75348.75408.754610.463510.464110.464710.465310.4659
98.97688.97748.97808.97868.979210.688110.688710.689310.689910.6905
109.20149.20209.20269.20319.203710.912710.913210.913810.914410.9150
119.42599.42659.42719.42779.428311.137211.137811.138411.139011.1396
129.65059.65119.65179.65239.652811.361811.362411.362911.363511.3641
139.87509.87569.87629.87689.877411.586311.586911.587511.588111.5887
1410.099610.100210.100810.101410.102011.810911.811511.812111.812611.8132
1510.324110.324710.325310.325910.326512.035412.036012.036612.037212.0378
Table 5. Ramification of some measures.
Table 5. Ramification of some measures.
S = 20S = 22S = 24
s η I η O C η L C s η I η O C η L C s η I η O C η L C
213.01960.00130.2345214.36430.00100.2309215.72930.00080.2276
313.23350.00090.2285314.55810.00070.2256315.90290.00050.2229
413.46080.00060.2232414.76550.00040.2210416.09010.00030.2188
513.70390.00040.2187514.98960.00030.2170516.29420.00020.2153
613.96150.00030.2150615.23050.00020.2137616.51610.00020.2124
714.22860.00020.2120715.48550.00010.2109716.75480.00010.2100
814.49680.00010.2095815.74940.00010.2087817.00710.00010.2080
S = 26S = 28S = 30
s η I η O C η L C s η I η O C η L C s η I η O C η L C
217.11280.00060.2246218.51310.00050.2218219.92820.00040.2193
317.26710.00040.2204318.64900.00030.2182320.04720.00020.2161
417.43430.00030.2168418.79720.00020.2150420.17760.00020.2133
517.61810.00020.2137518.96120.00010.2123520.32280.00010.2109
617.82010.00010.2112619.14310.00010.2100620.48500.00010.2089
718.04010.00010.2090719.34340.00010.2081720.66540.00000.2073
818.27660.00010.2072819.56160.00000.2065820.86420.00000.2059
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/math10132235

AMA Style

Vinitha V, Anbazhagan N, Amutha S, Jeganathan K, Shrestha B, Song H-K, Joshi GP, Moon H. Analysis of a Stochastic Inventory Model on Random Environment with Two Classes of Suppliers and Impulse Customers. Mathematics. 2022; 10(13):2235. https://doi.org/10.3390/math10132235

Chicago/Turabian Style

Vinitha, V., N. Anbazhagan, S. Amutha, K. Jeganathan, Bhanu Shrestha, Hyoung-Kyu Song, Gyanendra Prasad Joshi, and Hyeonjoon Moon. 2022. "Analysis of a Stochastic Inventory Model on Random Environment with Two Classes of Suppliers and Impulse Customers" Mathematics 10, no. 13: 2235. https://doi.org/10.3390/math10132235

APA Style

Vinitha, V., Anbazhagan, N., Amutha, S., Jeganathan, K., Shrestha, B., Song, H. -K., Joshi, G. P., & Moon, H. (2022). Analysis of a Stochastic Inventory Model on Random Environment with Two Classes of Suppliers and Impulse Customers. Mathematics, 10(13), 2235. https://doi.org/10.3390/math10132235

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