Next Article in Journal
Efficacy of Ultrasonic-Assisted Curing Is Dependent on Muscle Size and Ultrasonication System
Next Article in Special Issue
Scheduling Two Identical Parallel Machines Subjected to Release Times, Delivery Times and Unavailability Constraints
Previous Article in Journal
Accumulation of Airborne Toxic Elements and Photosynthetic Performance of Lolium multiflorum L. Leaves
Previous Article in Special Issue
Kinetic Parameter Determination for Depolymerization of Biomass by Inverse Modeling and Metaheuristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation-Based Optimization of a Two-Echelon Continuous Review Inventory Model with Lot Size-Dependent Lead Time

1
Industrial Engineering Department, College of Engineering, King Saud University, Riyadh 11421, Saudi Arabia
2
Industrial & Manufacturing Systems Engineering, College of Engineering& IT, Taiz University, Taiz 6803, Yemen
*
Author to whom correspondence should be addressed.
Submission received: 26 July 2020 / Revised: 11 August 2020 / Accepted: 17 August 2020 / Published: 19 August 2020
(This article belongs to the Special Issue Optimization Algorithms Applied to Sustainable Production Processes)

Abstract

:
This study analyzes a stochastic continuous review inventory system (Q,r) using a simulation-based optimization model. The lead time depends on lot size, unit production time, setup time, and a shop floor factor that represents moving, waiting, and lot size inspection times. A simulation-based model is proposed for optimizing order quantity (Q) and reorder point (r) that minimize the total inventory costs (holding, backlogging, and ordering costs) in a two-echelon supply chain, which consists of two identical retailers, a distributor, and a supplier. The simulation model is created with Arena software and validated using an analytical model. The model is interfaced with the OptQuest optimization tool, which is embedded in the Arena software, to search for the least cost lot sizes and reorder points. The proposed model is designed for general demand distributions that are too complex to be solved analytically. Hence, for the first time, the present study considers the stochastic inventory continuous review policy (Q,r) in a two-echelon supply chain system with lot size-dependent lead time L(Q). An experimental study is conducted, and results are provided to assess the developed model. Results show that the optimized Q and r for different distributions of daily demand are not the same even if the associated total inventory costs are close to each other.

1. Introduction

Inventory management plays a vital role in supply chains. It is one of the activities associated with managing customer and vendor relationships and plays a major role in facilitating the balancing of demand and supply. Effective inventory management in supply chains can lead to an increase in profits by cutting inventory-holding costs, and this indicates the need for a general model to manage inventories within a supply chain.
In classical continuous inventory review models ( Q , r ), demand and lead time are dealt with as constant input parameters or as random variables that follow a specific probability distribution. This could lead to an assumption that lead time is not subject to control and given as an input to the model. In practice, this assumption is not realistic since lead time can be controlled through some of its components, such as lot size, unit production cost, production time, and preparation time (setup, moving, and waiting times), as suggested by Kim and Benton [1]. In addition, considering lead time as a controllable factor allows for lead time reductions; this can lead to better customer satisfaction by lowering the risk of shortages and reducing safety stock [2]. Therefore, in many real-world systems, inventory models are stochastic and lead times are dependent and controllable. Several solution methodologies have been proposed to solve stochastic continuous inventory review models ( Q , r ) with a dependent lead time. These methodologies are classified into two main types: analytical models and simulation-based optimization models.
One of the first approximation models developed to study the impact of order quantities in a continuous review inventory policy ( Q , r ) on lead times was proposed by Kim and Benton [1]. The authors assumed that demand is normally distributed, order quantity Q is fixed, and the production lead time varies linearly with the order quantity. In addition, queue-waiting times were assumed to be a fixed fraction of the lead time. Based on an adjusted economic order quantity (EOQ), they developed an iterative algorithm to find the near-optimal order quantity ( Q ) and safety stock ( s ). Results showed that significant savings could be achieved if the interrelationship between order quantity and lead time is taken into account. Hariga [3] suggested a modification to the model by adjusting an expression related to annual backorder cost, in addition to proposing another relationship for the revised lot size. Both models assumed a constant processing time per unit. The model in [4] assumed random processing times and the processing time of a lot was correlated with its size. An exact model was developed, and two cases of lead time dependency were discussed. In the first, the lead time depends on lot size, and in the second, on the reserved capacity, whereby a customer can secure shorter processing times by reserving capacity at the supplier’s manufacturing facility. Ben-Daya and Hariga [5] proposed a model for an integrated production inventory problem for a vendor–buyer inventory system. A continuous review inventory policy was considered. Demand was assumed to be normally distributed and the lead time varied linearly with lot size. The lead time components were lot size-dependent run time and constant delay times, such as moving, waiting, and setup times. Moon and Cha [6] studied lead time reduction by controlling the manufacturer’s regular production rate and taking the risk of paying an additional cost into account. The relationship between the lot size and the production rate of the manufacturer expressed the lead time. An analytical model was proposed to optimize a continuous review inventory policy, based on the assumption that demand is normally distributed. Song et al. [7] proposed an integrated production inventory model for a supplier–buyer inventory system in which the buyer adopted a ( Q , r ) continuous review inventory policy. Demand was assumed to be normally distributed, while the supplier and buyer jointly controlled the lead time. The supplier determines the production rate while the buyer controls the lot size (order quantity). Song et al. [8] presented a Stackelberg game framework to model the interactions between a manufacturer and a retailer. Lead time and demand were distribution-free with only two known parameters: mean and variance. An iterative algorithm was developed to solve the proposed model. Sarkar et al. [9] proposed an integrated vendor–buyer model with shortages. Demand was assumed to be normally distributed and lead time was assumed to be a stochastic variable that depends on the order size and production rate. Sarkar and Giri [10] proposed an algorithm for finding the optimal solution of an integrated inventory model for a single-vendor single-buyer. The model assumed that demand follows a normal distribution and lead time is a linear function of batch size, setup time, and transport time.
Most of the proposed models were approximation models with the assumption that daily demand is normally distributed. However, it is extremely difficult to analytically solve the constrained optimization model that was developed by Hadley and Whitin [11]. Due to these limitations and difficulties, simulation-based optimization models are utilized to model and analyze complex inventory systems in order to obtain more practical results without the need to make restrictive assumptions. Al-harkan and Hariga [2] proposed a simulation-based optimization model for a single-supplier and single-retailer continuous review inventory system assuming a lead time dependent on lot size. The model was implemented for different daily demand distributions (normal, gamma, and log-normal). In addition, Tsai and Chen [12] proposed a simulation-based optimization model for a multi-objective single-supplier and single-retailer inventory system. The proposed model aims to find the values of reorder point and order quantity that simultaneously minimize the total inventory cost, average inventory level, and frequency of inventory shortages. The demand and lead time were assumed to be random and follow a uniform distribution. Avci and Selim [13] proposed a multi-objective simulation-based optimization approach for a three-echelon convergent supply chain system consisting of a customer, a manufacturer, and a number of suppliers. Kleijnen and Wan [14] studied an ( s , S ) inventory system to optimize reorder level s and maximum inventory level S . A simulation-based optimization model was developed using Arena and O p t Q u e s t . The metric to be minimized was the expected inventory costs, which include holding, ordering, and shortage costs while satisfying service (fill rate) constraints. Simulation-based optimization methods in stochastic inventory systems have been used to minimize total inventory costs while satisfying service level constraints as in [15,16,17,18]. A comprehensive survey of simulation-based optimization approaches to solve inventory replenishment problems is presented in [19].
Based on the literature review, we found that previous research solved the single-supplier and single-retailer continuous review inventory system. Most of the proposed models assumed that daily demand is normally distributed. Only studies that utilized simulation-based optimization models considered other distributions for daily demand. The lead time was considered as a variable and depended on several controllable factors in most previous studies. All the proposed models (analytical or simulation) aimed to minimize inventory costs.
In this paper, we solve a two-echelon three-stage supply chain stochastic continuous review inventory system with lot size-dependent lead time. The supply chain system has a supplier in the higher stage, a distributor in the middle stage, and two retailers in the lower stage. Different distributions for daily demand were considered. A simulation-based optimization model was developed to find near-optimal values of reorder point and order quantity that minimize the total inventory cost (holding, backlogging, and ordering costs). Hence, an efficient simulation-based optimization model is proposed for the first time for the problem considered herein.
The following section addresses the methodology. Section 3 presents the results and discussions. Finally, the conclusions are stated in Section 4.

2. Methodology

2.1. Problem Formulation

This study considers a three-stage two-echelon inventory system consisting of a supplier at the higher stage, a distributor at the middle stage, and two retailers at the lower stage, as shown in Figure 1. The main aim of this study is to find the optimal or near-optimal order quantities, as well as reorder points for the distributor and the two retailers that minimize the total inventory costs. The inventory costs include ordering costs, holding costs, and shortage costs in the three-stage supply chain.
A stochastic continuous review inventory control policy ( Q , r ) is employed, which means that retailer 1, retailer 2, or the distributor will place an order whenever the inventory position drops below the reorder point r ; the order size is Q . An Arena model is developed to simulate the policy employed in the three-stage supply chain for ten years (2500 days). The optimization study was conducted using the developed Arena model, as well as the Arena O p t Q u e s t optimization tool. Assuming that the demand follows different probability distributions, several scenarios are considered in order to study the performance of the simulation model. Simulated results were validated with results obtained using an analytical method developed in [1]. The following assumptions are considered:
  • A customer’s demand is satisfied immediately if the on-hand inventory of the retailer is sufficient. Otherwise, it is considered as a backorder and will be filled by future deliveries.
  • A retailer’s order is satisfied immediately if the on-hand inventory of the distributor is sufficient. Otherwise, it is also considered as a backorder and will be filled by future deliveries.
  • There is always a sufficient number of items in the supplier’s storage.
  • The lead time depends on the order size Q ; this is explained below.
The replenishment order will be received after a period of time, referred to as the lead time, which depends on the order size ( Q ). Based on the research conducted by Kim and Benton [1], a linear relationship between lead time and lot size can be assumed, and can be expressed mathematically as follows:
L ( Q ) = ( θ + P ·   Q ) δ
where θ is the setup time (in days); P is the unit production time in days, and δ is the shop floor queuing factor ( δ > 1 ) . We find that the lead time also depends on production time (represented by P ) and other times related to the order preparation, such as setup (represented by θ ), moving, and waiting times (represented by the shop floor factor δ ).
In the system described above, the problem is determining the optimal or near-optimal order sizes and reorder points that minimize the total inventory cost in the system shown in Figure 1 under the given assumptions. Our objective function (the total inventory cost) can be represented as the sum of the holding, backorders, and ordering costs. Generally, it is difficult to create an exact formulation of the total inventory cost for a stochastic continuous review inventory system, and several approximate formulations have been proposed in the literature. One of these approximation treatments is the constrained model developed in [11]; the reader can refer to the original publication for more information. The authors reported that it is extremely difficult to analytically solve the constrained optimization model, even with a normally distributed daily demand.
Therefore, we propose a simulation-based-optimization methodology that can be utilized for daily demand distributions of any type. A time-persistent model is developed to calculate the simulated total inventory cost for a period of n days; it is formulated with the following variables:
  • n Number of days;
  • n Y e a r Number of days per year;
  • t Point in time;
  • R i Retailer i , where i = 1 , 2 ;
  • D C Distributor;
  • I ( t ) Inventory level at time t ;
  • I + ( t ) Physical inventory level at time t , I + ( t ) = max { 0 ,   I ( t ) } ;
  • I ( t ) Backlog amount at time t , I ( t ) = max { 0 ,   I ( t ) } ;
  • N ( n ) The number of replenishment orders triggered in ( 0 ,   n ] ;
  • H The holding cost per item per unit time;
  • B The backorder cost per item per unit time;
  • O The ordering cost for each replenishment order;
  • T C ( n ) Total inventory cost for a given number of days ( n );
  • T C ( n ) R i Total inventory cost for a given number of days ( n ) at retailer i ;
  • T C ( n ) D C Total inventory cost for a given number of days ( n ) at the distributor ( D C ).
Then, the simulated total cost for n days, based on one simulation replication of the system, can be modeled as follows:
T C ( n ) = T C ( n ) R 1 + T C ( n ) R 2 + T C ( n ) D C
where
T C ( n ) R 1 = H R 1 ·   0 n I + R 1 ( t ) d t + B R 1 ·   0 n I R 1 ( t ) d t + O R 1 · N ( n ) R 1
T C ( n ) R 2 = H R 2 ·   0 n I + R 2 ( t ) d t + B R 2 ·   0 n I R 2 ( t ) d t + O R 2 . N ( n ) R 2
T C ( n ) D C = H D C ·   0 n I + D C ( t ) d t + B D C ·   0 n I D C ( t ) d t + O D C · N ( n ) D C
When conducting a simulation study, a certain number of replications are needed to estimate a performance measure. These replications are required to find the system’s steady-state and should be determined before the actual run of the simulation model. Several pilot runs, as well as some calculation steps, are conducted using the sequential procedure proposed in [20] to find the number of replications required. Based on the calculated number of replications, K , the following model is developed to calculate the expected annual total cost, which is formulated with the following variables:
  • β The source of system randomness, which represents the sequence of independent
  • U ( 0   ,   1 ) random numbers that activate a certain number of simulation replications;
  • β 1 ,   β 2 ,   ,   β K Independent realizations from the distribution of β ;
  • T C ( n ,   β k ) The observed total inventory cost ( TC ) in the interval ( 0 , n ] for the k t h simulation replication.
Then, the expected annual total inventory cost (ATC) can be modeled as follows:
E [ A T C ] = E [ T C ( n ) ] n * n Y e a r
where, n Y e a r is the number of days per year and
E [ T C ( n ) ] = lim K 1 K k = 1 K T C ( n ,   β k )
In order to optimize the ordering policy ( Q , r ), a mixed integer-programming model is formulated with the following variables:
  • r L R i The lower bound of the reorder point for retailer i , i = 1 ,   2 ;
  • r L D C The lower bound of the reorder point of the distributor;
  • r U R i The upper bound of the reorder point for retailer i ;
  • r U D C The upper bound of the reorder point for the distributor;
  • Q L R i The lower bound of the order quantity for retailer i ;
  • Q L D C The lower bound of the order quantity for the distributor;
  • Q U R i The upper bound of the order quantity for retailer i ;
  • Q U D C The upper bound of the order quantity for the distributor.
The objective function is given by
Min z = E [ A T C ]
subject to
r L R i r R i r U R i i = 1 , 2
r L D C r D C r U D C
Q L R i Q R i Q U R i i = 1 , 2
Q L D C Q D C Q U D C
Constraints (9)–(12) ensure that the reorder points and order quantities are within the corresponding lower and upper bounds.

2.2. Simulation-Optimization Model

As previously mentioned, in real-world problems, solving inventory systems analytically is challenging or extremely complex. This complexity is mainly due to the variability (or randomness) in demand and lead time. Based on the results reported in [11], most inventory models with a stochastic demand and lead time are a heuristic approximation of real-world systems. In stochastic continuous review systems ( Q , r ), analytical solutions can be obtained under some assumptions, such as normally distributed daily demand, independent lead time, and the out-of-stock time in a cycle being smaller than the cycle length [2]. However, using simulation-based optimization models, these assumptions can be easily met and included. Therefore, simulation is utilized to model and analyze complex inventory systems to obtain more practical results, which is the aim of this study.

2.2.1. Arena Model

As stated above, the system consists of a supplier and a distributor in the higher echelon, and two retailers in the lower echelon (see Figure 1). An Arena network model for the system is shown in Figure 2. The modeled network has five parts: each one represents a complete function of the system. In part 1, several variables related to retailers (1 and 2) are initialized, calculated, and updated. These variables are demand (customers’ orders) creations, inventory levels, inventory positions, and backorders (if any). It should be noted that the simulation model does not accept negative demand values that could be generated in the case of normal distribution. Instead, the negative demand values were regenerated again in a similar way to that of Al-harkan and Hariga [2]. In part 2, decisions related to placing orders (for both retailers) are made. In addition, inventory positions for the distributor and the associated retailer are updated. Then, in parts 3 and 4, replenishments are created for retailers 1 and 2, ordering costs are calculated for retailers 1 and 2, and backorders, if any, are updated. Finally, distributor replenishments are created, distributor ordering costs are calculated (if an order is made), and backorders (if assigned) are updated in part 5. Parameters were defined as variables in the “Variable module” as shown in Table 1. In addition, Table 2 presents the calculations of the costs reported in Equations (2)–(7) which were calculated using the “Statistic module”.
The model is checked to ensure that it is error-free and running as per requirement. Several pilot runs are also conducted using the required simulation statistics to determine the required number of replications K , which is needed to calculate the expected annual total inventory cost in Equations (6) and (7). Results from the pilot runs were also used in the model validation (refer to Section 2.3.2). Based on the selected performance measures of the system, we verify the accuracy of the simulation model in Section 2.3.2.

2.2.2. OptQuest

We use O p t Q u e s t , which is a useful tool embedded in the Arena software, to find the best (near-optimal) order sizes and reorder points for the system. O p t Q u e s t can be used to optimize complex systems that are difficult to formulate or solve analytically using classical optimization tools. It combines several metaheuristic approaches into a single search heuristic. These approaches are based on the tabu search algorithm, neural networks, and scatter search algorithm [14]. This advanced search framework is the key strength of O p t Q u e s t [21,22]. O p t Q u e s t performs a special (non-monotonic) search, where different evaluations are produced by the progressively generated inputs; these may not always be an improvement, but they provide an efficient path to the best solutions over time. In other words, if a candidate solution is not feasible, O p t Q u e s t eliminates it and investigates alternatives that are potentially better [20].
The optimization approach used by O p t Q u e s t requires a set of inputs to start the search process for better solutions in the model under investigation. These inputs are related to the model’s decision variables, objective function, and constraints. In our case, decision variables are the order size Q and reorder points r for the two retailers and the distributor (refer to Equations (9)–(12)). O p t Q u e s t needs to specify the initial, lower, and upper values of variables to be optimized. The suggested initial values state the search starting point of each optimized variable and have a significant effect on the efficiency and effectiveness of the search process. The related constraints are the lower and upper bound values for each decision variable to be optimized. These lower and upper bound values are based on the analytical solutions obtained in Section 2.3.2 and preliminary experiments. The objective function to be minimized is the annual total cost ( A T C ) (see Equation (8)). For practical cases, O p t Q u e s t allows simulation inputs to be defined as integers or continuous inputs, based on their nature. Moreover, a stopping criterion needs to be specified; in our case, the number of simulated candidate solutions is used as the stopping criterion.

2.3. Empirical Evaluation

The performance of the simulation-based optimization model is assessed in this section. The model configurations are first illustrated using an analytically solvable example (the case of normally distributed daily demands) that is used to validate the developed model. An experimental study is then conducted to solve the considered system. The effects of the variability of daily demand and lead time, as well as of the types of fitted distributions for the daily demand, on the inventory system’s performance are investigated.

2.3.1. Configurations

In order to clarify the developed simulation-based optimization model described in the previous section, we fitted the solved distributions into the system. In the existing literature, a similar case (two-stage supply chain with a supplier in the higher stage and a retailer in the lower stage) was solved analytically by [1,3]. Moreover, the same case was solved in [2] using a simulation optimization model for three cases of daily demand distributions: normal, gamma, and log-normal. The authors reported that the approximate problem is analytically solvable only for the normal distribution case, and it is hard to formulate and compute for other distributions. Therefore, we used the approximate analytical model to solve the system under the assumption that daily demands are normally distributed. We also assumed that the two retailers are identical and have the same average demand as the retailer in [1,3]. As a result, the average demand at the distributor was the sum of the demand at the two retailers. The rest of the required data were considered in a way similar to that used in [1,3]. A summary of the data is presented in Table 3. In order to verify the developed simulation-based optimization model, the results of the developed model were compared with the results of the analytical model.
The simulation-based optimization model has several configurations that need to be set before solving the special case under consideration. The lower bounds of the decision variables Q and r for the two retailers ( R 1 and R 2 ) and the distributor ( D C ) are (30, 30, 30) and (300, 300, 700), respectively. The upper bounds for the retailers and distributor are (150, 150, 300) and (800, 800, 1700), respectively. These lower and upper bound values of the Q and r are based on the analytical solutions obtained in Section 2.3.2 and preliminary experiments. The developed model is simulated and optimized for a period of 10 years (2500 days). The daily demand is normally distributed with a mean μ = 20 and standard deviation (SD)   σ = 5 ,   8 ,   10 ,   15 ,   and   22 . We used different values for the SD to study the effect of the variability in demand and lead time on the simulated annual total inventory cost for the specified period (2500 days). The number of replications K is estimated (based on the simulated annual total cost as the performance measure) using the sequential procedure proposed by Law and Kelton [20]. The relative error considered is γ = 0.15 with a confidence level of 95%. We found that 10 replications are enough to obtain the relative error. However, the actual simulation runs are conducted with 20 replications to obtain a better estimate of the annual total cost.

2.3.2. Model Validation

The results of the proposed simulation-based optimization model and the analytical approximation model are reported in Table 4. The terms used in Table 4 are order size (Q), reorder point (R), annual ordering cost (AOC), annual holding cost (AHC), annual backorders cost (ABC), the annual total cost for retailer 1 (ATC_R1), the annual total cost for retailer 2 (ATC_R2), the annual total cost for the distributor (ATC_D), and the total annual cost for the system (ATC).
A comparison between the two models was conducted using the results reported in Table 4. The comparison demonstrated that the proposed model provided accurate results. In fact, the simulation-based optimization model showed (for different values of SD) a slight deviation compared to the analytical model results. The maximum absolute ATC deviation is 9.74 % when σ = 15 . In almost all scenarios, which are different values of SD, the optimization-based simulation model provided smaller ATC values than the approximate analytical model. Results showed that the optimization-based simulation model outperformed the approximate analytical model, except for σ = 5 , when the absolute ATC deviation was 0.96 % , which is a very small value.
Figure 3 also represents the results of the two models (analytical and simulation). The terms used in Figure 3 are the same as those in Table 4. ATC increases with the increase in the demand variability ( σ ), which demonstrates the significance of the effect of demand variability on the ATC. Moreover, a paired t-confidence interval was used to analyze the differences in the means of the ATC generated by the analytical and simulation models. At a confidence level of α = 95 % , the reported paired t-confidence interval was ( 13994 ,   1067 ) and the p-value was 0.079 . Based on the reported interval and p -value, we can perceive that zero is included in the interval and the p-value is larger than 0.05; this means that there is no significant difference between the results generated by the two models. Thus, we can report that our model is accurate and valid.

3. Results and Discussion

As discussed in the previous section, a special case of the problem under consideration was solved analytically to validate the developed model. The solved case was a two-stage supply chain (one supplier and one retailer) with a lot size-dependent lead time. Moreover, it was assumed that the daily demand is normally distributed. However, it has been reported that the analytical solution to this special case is extremely difficult to arrive at [2,6]. Therefore, a simulation-based optimization model was developed to solve more complex real-world problems. In this study, we considered the stochastic inventory continuous review policy ( Q , r ) in a three-stage two-echelon supply chain system with lot size-dependent lead time,   L ( Q ) . The system has a supplier, a distributor, and two identical retailers (see Figure 1).
In this section, an experimental study was conducted to solve the system. The developed model was run for different daily demand distributions to study the effects of the variability of daily demand and lead time on the simulated ATC of the inventory system. We used three distributions for daily demand: normal, gamma, and log-normal. The simulation runs are conducted using different values of the coefficient of variance c v ( c v = 0.25 ,   0.4 ,   0.5 ,   0.75 ,   and   1.1 ) for each of the three distributions of daily demand. Based on the relative conversion formulas reported in [23] and the reported parameter values in [2], the used parameters of the gamma and log-normal distributions were chosen so that they provide the same values for the normal distribution parameters ( μ and σ ). Table 5 provides a summary of the parameters used for the three distributions. The results of the experimental study are summarized in Table 6 and Figure 4.
As shown in Figure 4, the ATC increases with the increase in the coefficient of variance ( c v ) for all the distributions and increases the most when the daily demand follows log-normal. This indicates a significant effect of the c v and the distribution type on the ATC. The results in Table 6 demonstrate that the solutions generated (order sizes and reorder points) for different distributions of daily demand are not the same even if the associated total costs are close to each other. This explains the difficulty in solving inventory systems with a stochastic daily demand.
Moreover, results have been tested statistically using the analysis of variance (ANOVA) with 0.95 confidence to examine whether the c v or the type of the daily demand distribution has a significant effect on the simulated ATC. p-values less than 0.05 indicate model terms are significant.
From Table 7, ANOVA test results revealed a significant effect of c v on ATC since its p-value is 0.00, less than 0.05. Overall, the results showed that c v has the most significant effect on ATC. Figure 5a,b illustrate the interaction between the distribution types and c v . From these interaction plots, it is evident that the ATC results difference between the three distributions is high with higher values of c v (Figure 5a). In addition, it is clear that the ATC increases with an increase in c v (Figure 5b).

4. Conclusions

In most inventory models, inventory costs are optimized with the assumption that the lead time is constant. Unfortunately, this is not the case in many real-world inventory systems. In practice, inventory systems are complex and difficult to formulate because of the stochastic nature of demand and lead times. In some cases, the lead time could have a few controllable components, such as lot size, production time, setup time, waiting and moving times, and lot size inspection time, which are related to order preparation and production. This study aimed to solve a complex inventory system that consists of two identical retailers, a distributor, and a supplier. The system is a stochastic inventory continuous review system ( Q , r ) with the lead time depending on lot size, unit production time, setup time, and a shop floor factor that represents moving, waiting, and lot size inspection times. A simulation-based optimization model was developed using Arena and O p t Q u e s t . A special case of the system was solved analytically to validate the model. Results indicate that the developed model was more accurate. Moreover, the model was tested for a more complex inventory system with different distributions of the system’s daily demand.
The results demonstrated that the optimized solutions (order sizes and reorder points) for different distributions of daily demand are not the same even if the associated total costs are close to each other. This explains the difficulty in solving inventory systems with a stochastic daily demand. Therefore, simulation-based optimization models can be considered as a powerful tool to optimize complex inventory systems without the need to make restrictive assumptions. Extensions of this work can be made by using different inventory control policies, taking into account supplier capacity, multi-products, and finally, cooperation between retailers.

Author Contributions

Conceptualization, I.A. and M.S.; methodology, M.S., A.F. and A.B.; formal analysis, I.A., M.S., M.G., A.F. and A.B. investigation, I.A., M.S., A.F. and A.B.; validation, M.S., A.F. and A.B.; visualization, M.S. and A.F.; data curation, M.S. and A.F; software, M.S. and A.F.; writing—original draft M.S., A.F., M.G. and A.B.; writing—review and editing, I.A., M.S. and M.G.; supervision, I.A.; resources, I.A.; funding acquisition, I.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Vice Deanship of Scientific Research Chairs (DSRVCH) and the APC was funded by Vice Deanship of Scientific Research Chairs (DSRVCH).

Acknowledgments

The authors are grateful to the Deanship of Scientific Research, King Saud University for funding this research project through Vice Deanship of Scientific Research Chairs (DSRVCH) and the authors thank RSSU at King Saud University for their technical support.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Kim, J.S.; Benton, W.C. Lot size dependent lead times in a Q,R inventory system. Int. J. Prod. Res. 1995, 33, 41–58. [Google Scholar] [CrossRef]
  2. Al-harkan, I.M.; Hariga, M.; Arabia, S. A simulation optimization solution to the inventory continuous review problem with lot size dependent lead time. Arab. J. Sci. Eng. 2007, 32, 327–338. [Google Scholar]
  3. Hariga, M.A. A stochastic inventory model with lead time and lot size interaction. Prod. Plan. Control 1999, 10, 434–438. [Google Scholar] [CrossRef]
  4. Çakanyildirim, M.; Bookbinder, J.H.; Gerchak, Y. Continuous review inventory models where random lead time depends on lot size and reserved capacity. Int. J. Prod. Econ. 2000, 68, 217–228. [Google Scholar] [CrossRef]
  5. Ben-Daya, M.; Hariga, M. Integrated single vendor single buyer model with stochastic demand and variable lead time. Int. J. Prod. Econ. 2004, 92, 75–80. [Google Scholar] [CrossRef]
  6. Moon, I.K.; Cha, B.C. A continuous review inventory model with the controllable production rate of the manufacturer. Int. Trans. Oper. Res. 2005, 12, 247–258. [Google Scholar] [CrossRef]
  7. Song, H.; Yang, H.; Luo, J. Integrated inventory model with lot size, production rate and lead time interactions. Int. J. Manag. Sci. Eng. Manag. 2010, 5, 141–148. [Google Scholar] [CrossRef]
  8. Song, H.-M.; Yang, H.; Bensoussan, A. Optimizing production and inventory decisions in a supply chain with lot size, production rate and lead time interactions. Appl. Math. Comput. 2013, 224, 150–165. [Google Scholar]
  9. Sarkar, S.; Giri, B.C.; Sarkar, A.K. A vendor–buyer inventory model with lot-size and production rate dependent lead time under time value of money. RAIRO Oper. Res. 2020, 54, 961–979. [Google Scholar] [CrossRef] [Green Version]
  10. Sarkar, S.; Giri, B.C. A vendor–buyer integrated inventory system with variable lead time and uncertain market demand. Oper. Res. 2020, 20, 491–515. [Google Scholar] [CrossRef]
  11. Hadley, G.F.; Whitin, T.M. Analysis of Inventory; Prentice-Hall: Upper Saddle River, NJ, USA, 1963. [Google Scholar]
  12. Tsai, S.C.; Chen, S.T. A simulation-based multi-objective optimization framework: A case study on inventory management. Omega 2017, 70, 148–159. [Google Scholar] [CrossRef]
  13. Avci, M.G.; Selim, H. A multi-objective simulation-based optimization approach for inventory replenishment problem with premium freights in convergent supply chains R. Omega 2017, 80, 153–165. [Google Scholar] [CrossRef]
  14. Kleijnen, J.P.C.; Wan, J. Optimization of simulated systems: OptQuest and alternatives. Simul. Model. Pract. Theory 2007, 15, 354–362. [Google Scholar] [CrossRef]
  15. Chu, Y.; You, F.; Wassick, J.M.; Agarwal, A. Simulation-based optimization framework for multi-echelon inventory systems under uncertainty. Comput. Chem. Eng. 2015, 73, 1–16. [Google Scholar] [CrossRef]
  16. Tsai, S.C.; Liu, C.H. A simulation-based decision support system for a multi-echelon inventory problem with service level constraints. Comput. Oper. Res. 2015, 53, 118–127. [Google Scholar] [CrossRef]
  17. Tsai, S.C.; Zheng, Y.-X. A simulation optimization approach for a two-echelon inventory system with service level constraints. Eur. J. Oper. Res. 2013, 229, 364–374. [Google Scholar] [CrossRef]
  18. Ye, W.; You, F. A computationally efficient simulation-based optimization method with region-wise surrogate modeling for stochastic inventory management of supply chains with general network structures. Comput. Chem. Eng. 2016, 87, 164–179. [Google Scholar] [CrossRef]
  19. Jalali, H.; Nieuwenhuyse, I.V. Simulation optimization in inventory replenishment: A classification. IIE Trans. (Inst. Ind. Eng). 2015, 47, 1217–1235. [Google Scholar] [CrossRef]
  20. Law, A.M.; Kelton, W.D.A. Simulation Modeling and Analysis; McGraw-Hill Education: New York, NY, USA, 2000. [Google Scholar]
  21. Hong, L.J.; Nelson, B.L. Discrete Optimization via Simulation Using COMPASS. Oper. Res. 2006, 54, 115–129. [Google Scholar] [CrossRef]
  22. Adenso-Díaz, B.; Laguna, M. Fine-Tuning of Algorithms Using Fractional Experimental Designs and Local Search. Oper. Res. 2006, 54, 99–114. [Google Scholar] [CrossRef]
  23. Leemis, L.M.; McQueston, J.T. Univariate distribution relationships. Am. Stat. 2008, 62, 45–53. [Google Scholar] [CrossRef]
Figure 1. Three-stage two-echelon supply chain.
Figure 1. Three-stage two-echelon supply chain.
Processes 08 01014 g001
Figure 2. Arena network model of the ( Q , r ) two-echelon stochastic inventory system.
Figure 2. Arena network model of the ( Q , r ) two-echelon stochastic inventory system.
Processes 08 01014 g002
Figure 3. Total annual cost for the system (ATC) of R1, R2, and DC.
Figure 3. Total annual cost for the system (ATC) of R1, R2, and DC.
Processes 08 01014 g003
Figure 4. ATC of R1, R2, and DC with normal, gamma, and log-normal distribution of daily demand.
Figure 4. ATC of R1, R2, and DC with normal, gamma, and log-normal distribution of daily demand.
Processes 08 01014 g004
Figure 5. Interaction of distribution types (pdf) and coefficient of variance ( c v ) : (a) cv * pdf and (b) pdf * cv.
Figure 5. Interaction of distribution types (pdf) and coefficient of variance ( c v ) : (a) cv * pdf and (b) pdf * cv.
Processes 08 01014 g005
Table 1. Parameter definitions in the “Variable module”.
Table 1. Parameter definitions in the “Variable module”.
No.VariableDefinitionNo.VariableDefinition
1Setup_CostFixed ordering cost15STOCK_DStock at DC
2Short_CostStock out cost16BACKO_R1Backorder R1
3Unit_Holding_Cost_dayInventory holding cost17BACKO_R2Backorder R2
4Prod_TimeUnit Preparation time18BACKO_DBackorder DC
5ThetaSetup time19Order_Size_R1Order size (R1)
6ZetaShop floor queuing factor20Order_Size_R2Order size (R2)
7n_runsNo. of days (n)21Order_Size_DOrder size (DC)
8n_yearNo. of days per year ( n Y e a r ) 22Reorder_Point_R1Reorder point (R1)
9No_YearsNo. of years ( n / n Y e a r ) 23Reorder_Point_R2Reorder point (R2)
10POS_R1Inventory position (R1)24Reorder_Point_DReorder point (DC)
11POS_R2Inventory position (R2)25Total_Ordering_Cost_R1Order cost of R1
12POS_DInventory position (DC)26Total_Ordering_Cost_R2Order cost of R2
13STOCK_R1Stock at R127Total_Ordering_Cost_DOrder cost of DC
14STOCK_R2Stock at R228Total_Ordering_CostTotal ordering cost
Table 2. Cost calculations using the “Statistic module”.
Table 2. Cost calculations using the “Statistic module”.
NameTypeExpression
Holding CostTime-Persistentn_runs*Unit_Holding_Cost_day*MX(STOCK_R1,0)+n_runs* Unit_Holding_Cost_day*MX(STOCK_R2,0)+n_runs*Unit_Holding_Cost_day*MX(STOCK_D,0)
Annual Holding CostOutputDAVG(Holding Cost)/No_Years
Backlog CostTime-Persistentn_runs*Short_Cost * BACKO_R1+n_runs*Short_Cost * BACKO_R2 + n_runs*Short_Cost * BACKO_D
Annual Backlog CostOutputDAVG(Backlog Cost)/No_Years
Annual Order CostOutput (Total_Ordering_Cost_R1)/No_Years+ (Total_Ordering_Cost_R2)/No_Years + (Total_Ordering_Cost_D)/No_Years
Annual Total CostOutputOvalue(Annual Order Cost) +Ovalue(Annual Holding Cost)+ Ovalue(Annual Backlog Cost)
Table 3. Data of Retailer 1, Retailer 2, and Distributor.
Table 3. Data of Retailer 1, Retailer 2, and Distributor.
Retailer 1 and Retailer 2
Average demand per year D R 1 , R 2 = 5000
Average demand per day (250 days in a year) μ R 1 , R 2 = 20
Inventory holding cost H R 1 , R 2 = $ 0.5 per unit per day
Stock out cost B R 1 , R 2 = 600 per unit (out of stock)
Setup time θ R 1 , R 2 = 0.125 days
Ordering cost O R 1 , R 2   = $ 125 per order
Unit preparation time p R 1 , R 2 = 0.025 day per unit
Queueing factor δ R 1 , R 2 = 10
Distributor DC
Average demand per year D D C = 10000
Average demand per day (250 days in a year) μ D C = 40
Inventory holding cost H D C = $ 0.5 per unit per day
Stock out cost B D C = 600 per unit (out of stock)
Setup time θ D C = 0.125 days
Ordering cost O D C   = $ 125 per order
Unit production time p D C = 0.025 day per unit
Queueing factor δ D C = 10
Table 4. Results of the simulation-optimization model and the analytical method.
Table 4. Results of the simulation-optimization model and the analytical method.
σ 58101522
Model AnalyticalSimulation AnalyticalSimulation Analytical Simulation Analytical Simulation Analytical Simulation
Retailer 1Q77.7583.0068.5367.0063.5668.0053.9655.0044.8945.00
R475.67494.00462.48460.00458.24480.00458.43483.00474.70500.00
AOC8038.517536.259119.439363.759833.489446.8811,582.1112,829.3813,923.2618,874.38
AHC12,599.2012,143.4016,134.1916,699.9418,402.7817,527.3023,824.8121,619.2830,962.9822,640.41
ABC863.48766.501288.83224.441545.001468.312120.88772.102818.768632.84
ATC_R121,501.1920,446.1526,542.4526,288.1429,781.2628,442.4837,527.8035,220.7547,705.0050,147.62
Retailer 2Q77.7579.0068.5371.0063.5663.0053.9652.0044.8945.00
R475.67474.00462.48465.00458.24446.00458.43467.00474.70499.00
AOC8038.517877.509119.438847.509833.4810,143.7511,582.1113,684.3813,923.2618,720.63
AHC12,599.2012,140.8716,134.1915,022.7718,402.7816,409.8623,824.8121,192.3730,962.9822,867.87
ABC863.48509.221288.831308.881545.001738.692120.881283.832818.768279.70
ATC_R221,501.1920,527.5926,542.4525,179.1529,781.2628,292.3037,527.8036,160.5847,705.0049,868.19
DistributorQ113.40111.00103.14107.0097.86107.0088.17103.0079.4979.00
R1292.681297.001249.321310.001234.601368.001229.101490.001265.001476.00
AOC11,022.9711,247.5012,120.0311,741.2512,773.0211,989.3814,177.2213,788.7515,723.5421,449.38
AHC20,673.4525,025.5327,442.4130,489.1131,863.1634,770.5842,685.8740,124.6257,469.7947,645.67
ABC2564.72762.004311.95681.005540.191038.008778.421698.0013,598.611521.00
ATC_D34,261.1437,035.0343,874.3942,911.3650,176.3747,797.9665,641.5155,611.3786,791.9470,616.04
ATC77,263.5278,008.7796,959.2994,378.65109,738.90104,532.74140,697.10126,992.70182,201.90170,631.85
Table 5. The parameters of the three distributions.
Table 5. The parameters of the three distributions.
c v Normal
( μ , σ )
Gamma
( α , β )
Log-Normal
( α , β )
0.25 ( 20 ,   5 ) (16.00, 1.25)(2.965, 0.246)
0.40 ( 20 ,   8 ) (6.25, 3.20)(2.9215, 0.385)
0.50 ( 20 ,   10 ) (4.00, 5.00)(2.884, 0.472)
0.75 ( 20 ,   15 ) (1.78, 11.25)(2.773, 0.668)
1.10 ( 20 ,   22 ) (0.83, 24.20)(2.599, 0.891)
Table 6. Results for different types of normal, gamma, and log-normal demand.
Table 6. Results for different types of normal, gamma, and log-normal demand.
Inventorycv = σ/μ0.250.40.50.751.1
pdfNormalGammaLog-NormalNormalGammaLog-NormalNormalGammaLog-NormalNormalGammaLog-NormalNormalGammaLog-Normal
Retailer 1Q834683674967684671555562454542
R494302491460363463480366517483493582500498750
AOC7536.313,516.37516.39363.812,696.99313.89446.913,493.88772.512,829.411,300.010,091.318,874.413,882.514,979.4
AHC12,143.49104.911,907.716,699.914,977.917,316.817,527.317,092.421,910.021,619.327,828.834,772.822,640.433,342.865,892.3
ABC766.5822.51590.7224.41086.21194.61468.3757.61750.9772.1987.06041.08632.816,519.930,001.2
ATC_R120,446.223,443.621,014.626,288.128,761.027,825.228,442.531,343.732,433.435,220.840,115.850,905.150,147.663,745.2110,872.9
Retailer 2Q795571715165635191525053454547
R474352434465358455446397629467474533499500750
AOC7877.511,306.38782.58847.512,192.59580.010,143.812,228.16886.913,684.412,496.311,750.618,720.613,912.513,273.1
AHC12,140.910,298.311,507.815,022.813,281.217,504.216,409.918,005.924,430.921,192.428,034.733,798.722,867.933,536.763,723.3
ABC509.2910.41140.41308.91521.01424.71738.71096.2655.01283.81045.02317.18279.715,886.4726.7
ATC_R220,527.622,515.021,430.725,179.226,994.828,508.928,292.331,330.231,972.836,160.641,575.947,866.449,868.263,335.577,723.2
DistributorQ11110311110783107107911101039795799375
R129711961304131010421313136811551397149013451343147615001311
AOC11,247.512,067.511,246.311,741.314,953.811,660.011,989.413,660.011,378.113,788.812,860.013,153.121,449.413,482.516,713.1
AHC25,025.522,417.925,986.930,489.126,783.331,858.134,770.631,345.938,374.640,124.647,547.049,432.847,645.771,014.868,653.6
ABC762.0435.0909.0681.01215.01092.01038.03270.04582.51698.04360.57125.01521.03660.04290.0
ATC_D37,035.034,920.438,142.142,911.442,952.144,610.147,798.048,275.954,335.255,611.464,767.569,711.070,616.088,157.389,656.7
ATC78,008.880,879.080,587.494,378.798,707.8100,944.2104,532.7110,949.9118,741.4126,992.7146,459.2168,482.4170,631.9215,238.1278,252.8
Table 7. ANOVA table for the ATC.
Table 7. ANOVA table for the ATC.
SourceDFSeq SSContributionAdj SSAdj MSF-Valuep-Value
CV43.76 × 101084.60%3.76 × 10109.39 × 10919.500.000
PDF22.98 × 1096.72%2.98 × 1091.49 × 1093.100.101
Error83.85 × 109 8.68%3.85 × 1094.82 × 108
Total144.44 × 1010100.00%

Share and Cite

MDPI and ACS Style

Alharkan, I.; Saleh, M.; Ghaleb, M.; Farhan, A.; Badwelan, A. Simulation-Based Optimization of a Two-Echelon Continuous Review Inventory Model with Lot Size-Dependent Lead Time. Processes 2020, 8, 1014. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8091014

AMA Style

Alharkan I, Saleh M, Ghaleb M, Farhan A, Badwelan A. Simulation-Based Optimization of a Two-Echelon Continuous Review Inventory Model with Lot Size-Dependent Lead Time. Processes. 2020; 8(9):1014. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8091014

Chicago/Turabian Style

Alharkan, Ibrahim, Mustafa Saleh, Mageed Ghaleb, Abdulsalam Farhan, and Ahmed Badwelan. 2020. "Simulation-Based Optimization of a Two-Echelon Continuous Review Inventory Model with Lot Size-Dependent Lead Time" Processes 8, no. 9: 1014. https://0-doi-org.brum.beds.ac.uk/10.3390/pr8091014

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