Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 2005-0461(Print)
ISSN : 2287-7975(Online)
Journal of Society of Korea Industrial and Systems Engineering Vol.46 No.4 pp.93-100
DOI : https://doi.org/10.11627/jksie.2023.46.4.093

The Tradeoff of Bullwhip Effect with Inventory Costs in a Supply Chain

Heung-Kyu Kim†
School of Business Administration, Dankook University
Corresponding Author : heungkyu@dankook.ac.kr
30/11/2023 06/12/2023 06/12/2023

Abstract


In this paper, an alternative inventory policy that trades off the bullwhip effect at an upstream facility with cost minimization at a current facility, with the goal of reducing system wide total expected inventory costs, when external demand distributjon is autocorrelated, is considered. The alternative inventory policy has a form that is somewhere between one that completely neglects the autocorrleation and one that actively utilizes the autocorrelation. For this purpose, a mathematical model that allows us to evaluate system wide total expected inventory costs for a periodic review system is developed. This model enables us to identify an optimal inventory policy at a current facility that minimizes system wide total expected inventory costs by the best tradeoff of the bullwhip effect at an upstream facility with cost minimization at a current facility. From numerical experiments, it has been found that (i) when the autocorrelation is negative, the optimal policy is one that actively utilizes the autocorrelation, (ii) when the autocorrelation is small and positive, the optimal policy is one that neglects the autocorrelation, and (iii) when the autocorrelation is large and positive, the optimal policy is somewhere between one that actively utilizes the autocorrelation and one that neglect the autocorrelation.



공급사슬에서 채찍효과와 재고비용 사이의 상충

김흥규†
단국대학교 경영학부

초록


    1. Introduction

    In industry, actual demand distribution of many consumer goods have often been found to be autocorrelated [1,2,4]. However, industry practitioners are not aware of accurate inventory control models since most inventory control models are based on the assumption that the demand is independent from one period to the next. Therefore, they do not make use of accurate inventory control models [6].

    Urban [6] determined accurate reorder levels for a continuous review system and Zinn et al. [7] conducted simulation studies to evaluate the impact of autocorrelation on safety stocks. Recently, Kim [2] developed a mathematical model that, for a periodic review system, compared accurate inventory control models with traditional inventory control models under autocorrelated demand in a supply chain in which there are two participants, a retailer and a manufacturer. From which, Kim [2] found that knowledge of the actual demand process is always beneficial to the retailer in terms of reducing long run average total inventory costs per period, given that the retailer knows how to use it, and that the retailer’s knowledge of the actual demand process may induce bullwhip effect, thus increasing the long run average total inventory costs per period at the manufacturer. Therefore, it is possible that sometimes a supply chain in which the retailer uses traditional inventory control models and the manufacturer uses accurate inventory control models has a lower system wide long run average total costs than a supply chain in which both the retailer and the manufacturer use accurate inventory control models.

    The above observation naturally leads us to consider a policy for the retailer that trades off reduction of the bullwhip effect at the manufacturer with cost minimization at the retailer, which is the goal of this paper. Therefore, this paper can be regarded as a sequel to Kim [2].

    The remainder of this paper is organized as follows: Section 2 reviews the supply chain model of Kim [2] and Section 3 develops a policy for the retailer that trades off reduction of the bullwhip effect at the manufacturer with cost minimization at the retailer. Then, in Section 4, through numerical studies, managerial insights gleaned are explained. Final remarks are addressed in Section 5.

    2. Model Review

    As mentioned in Section 1, the model developed by Kim [2], is used in this paper. For the clear description of this paper, a special case of the model, where the review period (lead time) of the retailer’s is 1 (0) and the review period (lead time) of the manufacturer’s is 1 (1), is reviewed in this section. For general cases, see Kim [2].

    2.1 Model Structure

    A supply chain consisting of a single retailer and a single manufacture is considered where external demand distribution is an autocorrelated one. If we let dt be the demand faced by the retailer during period t, t∈{1, 2, 3,…}, then we can write

    d t = μ + ρ d t 1 + t
    (1)

    where μ > 0 , - 1 < ρ < 1, and ϵtN(0, σ2).

    The retailer, at the start of every review period t, observes the inventory level and the previous demands, and calculates the order-up-to level yj, t , j = s, n, from which the retailer determines the order quantity qj, t , j = s, n, to place to the manufacturer, and the shipment of which the retailer receives right away. The order the retailer receives, in turn, is used to meet the previous backlogged demands and the current period’s demand. Here the subscript ‘s’, short for ‘smart’, refers to the retailer who is aware that the demand distribution is an autocorrelated one and thus takes advantage of this knowledge to determine the order-up-to level and the subscript ‘n’, short for ‘naïve’, refers to the retailer who is not aware that the demand distribution is an autocorrelated one and thus resort to an inventory policy based on the assumption of i.i.d. demand to determine the order-up-to level.

    The manufacturer, at the start of period t, receives and ships the order quantity qj, t, j = s, n, to the retailer. If the manufacturer does not have enough stock on hand to fill the order quantity, the manufacturer can always find an alternative source to borrow from and that the borrowed items are returned to the source when the next replenishment arrives. The manufacturer places an order at the start of period t, right before the retailer orders. The order arrives at the start of period t + 1. The supplier from which the manufacturer orders is assumed to have infinite capacity so that the manufacturer’s order is always satisfied. As for the manufacturer, whether the retailer acts smart or naïve (j = s, n) determines the demand stream faced by the manufacturer. In addition, it is assumed that the manufacturer acts smart and knows the demand stream faced by the retailer, which is made possible by information sharing with the retailer.

    2.2 The Retailer

    Smart Retailer From Equation (1), the actual effective lead time (=1) demand follows a normal distribution with mean

    E [ d t | d t 1 ] = μ + ρ d t 1

    and variance

    V ( d t | d t 1 ) = σ 2 .

    Next, the order-up-to level for period t, ys, t will be calculated as

    y s , t = μ + ρ d t 1 + z σ
    (2)

    where z is a constant, whose value will be described later. Next, to calculate the expected holding cost per period, an approximation for the average inventory level, which is y s , t E [ d t | d t 1 ] + E [ d t | d t 1 ] 2 , given by Silver and Peterson [5], is used. Then the average inventory level over the period t is approximated as follows:

    i n υ s = z σ + μ + ρ d t 1 2 .

    Therefore, the expected average inventory level per period can be represented as follows:

    i n υ s = z σ + 1 2 ( μ 1 ρ )
    (3)

    Next, notice that the average number of stockouts for period t for an order-up-to level of ys, t can be written as

    y s , t ( d t | d t 1 y s , t ) d F ( d t | d t 1 ) ,

    where F (dtdt - 1) is the cumulative distribution function of demand for the period t, given dt - 1. Since dtdt - 1 follows a normal distribution, the above formula can be simplified to

    σ z s ( x z s ) ϕ ( x ) d x = σ [ ϕ ( z s ) z s ( 1 Φ ( z s ) ) ]
    (4)

    where x is a standard normal random variable, z s = y s , t ( μ + ρ d t 1 ) σ ( = z ) is the standardized value of the order- up-to level, ϕ(∙) is the probability distribution for the standard normal random variable, and Φ(∙) is the cumulative distribution function for the standard normal random variable.

    Therefore, if we denote the holding cost per unit per unit time by h and the penalty cost per unit associated with backlogged demand by p, the long run average inventory cost per period, gs , can be written as

    g s = h i n υ s + p σ [ ϕ ( z s ) z s ( 1 Φ ( z s ) ) ] ,
    (5)

    where z = Φ 1 ( p h p ) .

    Finally, the order quantity for period t, which becomes the manufacturer’s demand for that period, can be written as

    q s , t = y s , t y s , t 1 + d t 1 = ( 1 + ρ ) d t 1 ρ d t 2
    (6)

    Naive Retailer A naïve retailer believes that the effective lead time demand is i.i.d. from a normal distribution with a mean E [ d t ] = μ 1 ρ and a variance V [ d t ] = σ 2 1 ρ 2 . which becomes the perceived effective lead time demand,

    Next, the order-up-to level for period t, yn, t, which is constant from period to period, will be calculated as

    y n , t = μ 1 ρ + z σ 1 ρ 2
    (7)

    Notice that the expected average inventory level per period can be represented as follows:

    i n υ n = z σ 1 ρ 2 + 1 2 ( μ 1 ρ )
    (8)

    Next, since z n = y n , t ( μ + ρ d t 1 ) σ is a random variable due to the dependence of dt - 1, the long run average number of stockouts per period can be written as

    E [ σ [ ϕ ( z n ) z n ( 1 Φ ( z n ) ) ] ] ,

    where the expectation is taken over zn.

    Since zn follows a normal distribution with mean z 1 ρ 2 and variance ρ 2 1 ρ 2 , the above expression for the long run average number of stockouts per period can be simplified as

    σ [ h ( E [ z n ] ) ( 1 + V ( z n ) ) E [ z n ] ( 1 H ( E [ z n ] ) ) ] ,
    (9)

    where h(∙) and H(∙) are the pdf and cdf for a normal distribution with mean 0 and variance 1 +V(zn). Fore derivation, see Kim and Ryan [3].

    Therefore, the long run average inventory cost per period, gs , can be written as

    g n = h i n υ n + p [ h ( E [ z n ] ) ( 1 + V ( z n ) ) E [ z n ] ( 1 H ( E [ z n ] ) ) . ]
    (10)

    Finally, the order quantity for period t, which becomes the manufacturer’s demand for that period, can be written as

    q n , t = y n , t y n , t 1 + d t 1 = d t 1
    (11)

    The long run average total inventory costs per period for a naive retailer, gs , will be always larger than or equal to the long run average total inventory costs per period for a smart retailer, gs . For details, see Kim [2].

    2.3 The Manufacturer

    With Smart Retailer First, notice that the actual effective lead time (=2) demand faced by the manufacturer with smart retailer can be written as

    D s , t = u = 0 1 q s , t + u | d t 1 , d t 2 , = ( 1 + ρ ) μ + ( 1 + ρ + ρ 2 ) d t 1 ρ d t 2 + ( 1 + ρ ) t

    Therefore, the actual effective lead time demand faced by the manufacturer with the smart retailer follows a normal distribution with mean

    E [ D s , t ] = ( 1 + ρ ) μ + ( 1 + ρ + ρ 2 ) d t 1 ρ d t 2 ,

    and variance

    V [ D s , t ] = ( 1 + ρ 2 ) σ 2 .

    Next, the order-up-to level at the start of period t, Ys, t will be calculated as

    Y s , t = ( 1 + ρ ) μ + ( 1 + ρ + ρ 2 ) d t 1 ρ d t 2 + Z ( 1 + ρ ) σ
    (12)

    where Z is a constant, whose value will be described later.

    Notice that the expected average inventory level per period can be represented as follows:

    I N V s = Z ( 1 + ρ ) σ + 1 2 ( μ 1 ρ ) .
    (13)

    Next, notice that, since Z s = Y s , t { ( 1 + ρ ) μ + ( 1 + ρ + ρ 2 ) d t 1 ρ d t 2 } ( 1 + ρ ) σ ( = Z ) is the standardized value of the order-up-to level, the average number of stockouts for period t for an order-up-to level of Ys, t can be written as

    ( 1 + ρ ) σ [ ϕ ( Z s ) Z s ( 1 Φ ( Z s ) ) ]
    (14)

    Therefore, if we denote the holding cost per unit per unit time, by H , and the penalty cost per unit associated with backlogged demand by P , the long run average inventory cost per period, Gs , can be written as

    G s = H I N V s + P ( 1 + ρ ) σ [ ϕ ( Z ) Z ( 1 Φ ( Z ) ) ] ,
    (15)

    where Z = Φ 1 ( P H P ) .

    With Naive Retailer First, notice that the actual effective lead time demand faced by the manufacturer with naive retailer can be shown as

    D n , t = u = 0 1 q n , t + u | d t 1 , d t 2 , = μ + ( 1 + ρ ) d t 1 + t

    Therefore, the actual effective lead time demand faced by the manufacturer with the naive retailer follows a normal distribution with mean

    E [ D n , t ] = μ + ( 1 + ρ ) d t 1 ,

    and variance

    V ( D n , t ) = σ 2 .

    Next, the order-up-to level at the start of period t, Yn, t will be calculated as

    Y n , t = μ + ( 1 + ρ ) d t 1 + Z σ .
    (16)

    Notice that the expected average inventory level per period can be represented as follows:

    I N V n = Z σ + 1 2 ( μ 1 ρ ) .
    (17)

    Next, notice that, since Z n = Y n , t { μ + ( 1 + ρ ) d t 1 } σ ( = Z ) is the standardized value of the order-up-to level, the average number of stockouts for period t for a fixed value of Yn, t can be written as

    σ [ ϕ ( Z n ) Z n ( 1 Φ ( Z n ) ) ] .
    (18)

    Therefore, the long run average inventory cost per period, Gn , can be written as

    G n = H I N V n + P σ [ ϕ ( Z ) Z ( 1 Φ ( Z n ) ) ] ,
    (19)

    where Z = Φ 1 ( P H P ) .

    Notice that when, ρ > 0 and the retailer is navie, the manufacturer faces less severe bullwhip effect, i.e., the variance of the orders placed by the naive retailer, V(Dn, t), is less than that of the orders placed by the smart retailer, V(Ds, t). Therefore, when ρ > 0, the manufacturer is better off with a naive retailer than with a smart retailer. That is, when ρ > 0, Gs > Gn .

    Now that we have briefly reviewed the model in Kim [2] for calculating the long run average inventory cost per period incurred at each participant in the supply chain, we are now ready to consider an alternative policy for the retailer that trades off reduction of the bullwhip effect at the manufacturer with cost minimization at the retailer.

    3. An Alternative Policy at the Retailer

    3.1 Relationship between Retailer and Manufacturer

    Since we are interested in system wide expected total inventory costs, we define the following function:

    T C j = g j + G j , j = s , n ,
    (20)

    where j, j = s, n, denotes a form of the retailer’s order-up-to level.

    As demonstrated in Section 2, when ρ > 0, we have V(Ds, t) > V(Dn, t) , which causes Gs > Gn, which, in turn, sometimes leads us to observe TC s > TC n . In other words, the smart retailer causes more bullwhip effect at the manufacturer, which leads to higher costs at the manufacturer. Thus, sometimes, depending on the system parameter values, we have a system with naive retailer perform better than a system with smart retailer.

    This observation leads us to the following question: “What is the form for an order-up-to level at the retailer that trades off reduction of the bullwhip effect at the manufacturer with cost minimization at the retailer, with the goal of reducing system wide expected total inventory costs?” In other words, we want to find an order-up-to policy for the retailer that provides lower system wide expected total inventory costs than TC s or TC n , i.e., one that provides lower system wide expected total inventory costs than a system with a smart or a system with a niave retailer.

    As noted above, the policy we seek will trade off reduction of the bullwhip effect at the manufacturer with cost minimization at the retailer. Notice that order-up-to level for period t, which we denote by yα, t, of form

    y α , t = α y n , t ( 1 α ) y s , t ,
    (21)

    where 0 ≤ α ≤ 1, allows us to make this trade off.

    From Equation (21), if α is close to 0 and ρ is positive, then yα, tys, t. Therefore, minimizing costs at the retailer becomes a priority. If α is close to 1, then yα, tyn, t . Therefore, reducing the bullwhip effect becomes a priority. Notice that α is also a measure of the amount of neglect to the most recent demand since ys, t (when α = 0) refers to an order-up-to policy that reacts to the most recent demand information and yn, t (when α = 1) refers to an order-up-to policy that completely neglects the most recent demand information. The smaller value of α indicates more reaction to the most recent demand information.

    We define TC α, i.e., the system wide expected total inventory costs, as follows:

    T C α = g α + G α .

    Next, we have a following proposition:

    Proposition 1minαTCαmin {TC n, TC s}, where 0 ≤ α ≤ 1.

    Proof. Proposition 1 is self-explanatory.

    3.2 Optimal Value of α

    Given a set of order-up-to levels of this form, we next consider the optimal value of α , i.e., the value of α that minimizes TC α, as a function of the autocorrelation.

    Retailer In this case, the order-up-to level, at the start of period t, for the retailer, becomes:

    y α , t = α y n , t + ( 1 α ) y s , t + z ( α σ 1 ρ 2 + ( 1 α ) σ ) , = α ( μ 1 ρ ) + ( 1 α ) ( μ + ρ d t 1 )
    (22)

    where z = Φ 1 ( p h p ) . In addition, the average inventory level per period, inυα, becomes

    i n υ α = z ( α ( σ 1 ρ 2 ) + ( 1 α ) σ ) + 1 2 ( μ 1 ρ ) .
    (23)

    Therefore, the expected total inventory costs per period for the retailer, gα, can be written as:

    g α = h i n υ α + p σ [ h ( E [ z α ] ) ( 1 + V ( z α ) ) E [ z α ] ( 1 H ( E [ z α ] ) ) ] ,
    (24)

    where E [ z α ] = z ( α 1 1 ρ 2 + ( 1 α ) ) , V ( z α ) = α 2 ρ 2 1 ρ 2 , and h(∙) and H(∙) are the pdf and cdf for a normal distribution with mean 0 and variance 1 +V(zα). Finally, the order quantity placed by this retailer, at the start of period t, can be written as:

    q α , t = y α , t y α , t 1 + d t 1 = ( 1 α ) ρ ( d t 1 d t 2 ) + d t 1 .
    (25)

    Notice that qα, t is a function of the decision variable, α , in addition to the previous demands at the retailer.

    Manufacturer Next, the actual effective lead time demand faced by the manufacturer at the start of period t, becomes:

    D α , t = u = 0 1 q α , t + u | d t 1 , d t 2 , = q α , t | d t 1 , d t 2 + q α , t + 1 | d t 1 , d t 2 = ( 1 + ( 1 α ) ρ ) μ + ( 1 + ρ + ( 1 α ) ρ 2 ) d t 1 ( 1 α ) ρ d t 2 + ( 1 + ( 1 α ) ) t
    (26)

    From Equation (26), we have

    E [ D α , t ] = ( 1 + ( 1 α ) ρ ) μ + ( 1 + ρ + ( 1 α ) ρ 2 ) d t 1 ( 1 α ) ρ d t 2

    and

    V ( D α , t ) = ( 1 + ( 1 α ) ) 2 σ 2 .

    Notice that, if α = 0 (α = 1), these expressions are identical to those presented in Section 2 for the smart (naive) retailer.

    Therefore, we have

    Y α , t = E [ D α , t ] + Z V ( D α , t ) = ( 1 + ( 1 α ) ρ ) μ + ( 1 + ρ + ( 1 α ) ρ 2 ) d t 1 ( 1 α ) ρ d t 1 + Z ( 1 + ( 1 α ) ρ ) σ ,

    where Z = Φ 1 ( P H P ) .

    Now that we have written the order-up-to level for the manufacturer as a function of α, we have:

    I N V α = Z ( 1 + ( 1 α ) ρ ) σ + 1 2 ( μ 1 ρ ) .
    (27)

    Therefore, the expected total inventory costs per period for the manufacturer, given α , can be written as follows:

    G α = H ( Z ( 1 + ( 1 α ) ρ ) σ + 1 2 ( μ 1 ρ ) ) + P ( 1 + ( 1 α ) ρ ) σ [ ϕ ( Z ) Z ( 1 Φ ( Z ) ) ] ,
    (28)

    where Z = Φ 1 ( P H P ) .

    Thus we can write the system wide expected total inventory costs as TC α = gα +Gα. Notice that this system wide expected total inventory costs is a function of α only, the decision variable, once the system parameters such as μ, ρ, σ, etc. are given.

    4. Numerical Experiment

    Now that we have developed expressions for the expected total inventory costs per period for the retailer and the manufacturer as a function of α , we are ready to determine the optimal value of α , the value of α that leads to the minimum system wide expected total inventory costs, either numerically or analytically.

    We now present a numerical example to examine how the optimal value of α varies as a function of ρ. In this example, the demand process is specified by μ = 600, σ = 15. The cost parameters for the retailer are h = 2, p = 50, and those for the manufacturer are H = 1, P = 25. In Figure 1 (a), (b), (c) and (d), the y-axis refers to the system wide expected total inventory costs per period whereas the x-axis refers to the values of α which vary between 0 and 1 when ρ =- 0.7, - 0.3, 0.3, 0.7, respectively.

    From <Figure 1>, some interesting results are found. They are:

    1. When ρ is negative, the optimal value of α is always 0. This implies that ys, t is optimal, i.e., a system with a smart retailer is optimal.

    2. When ρ is small and positive, the optimal value of α is close to 1. This implies that more weight should be placed on the form of order-up-to level that reduces bullwhip effect, which is yn, t .

    3. When ρ is large and positive, the optimal value of α is somewhere between 0 and 1. This implies that an appropriate balance is needed between reduction of bullwhip effect at the manufacturer and cost minimization at the retailer.

    We can interpret these results as follows. Result 1 is due to the fact that, when ρ is negative, ys, t not only minimizes the retailer's costs but also causes less bullwhip effect. Notice that, when ρ is negative,

    V ( q s , t ) = V ( ( 1 + ρ ) d t 1 ρ d t 2 ) = ( 1 + 2 ( 1 ρ 2 ) ) σ 2 1 ρ 2 < σ 2 1 ρ 2 = V ( d t ) = V ( q n , t ) .

    Result 2 is due to the fact that, when ρ is small and positive, the previous demand, dt - 1, is not a good indicator for the current demand, dt . Therefore, by adapting to the most recent demand information, the retailer causes more bullwhip effect without significantly improving its forecasts and its own costs. So, in this case, reducing the bullwhip effect is more critical than reducing costs at the retailer.

    Finally, result 3 is due to the fact that, when ρ is large, the previous demand becomes a good predictor for the current demand. Therefore, by adapting the order-up-to level to the most recent demand information, the reduction in the expected total inventory cost at the retailer is great, but causes an increase in the expected total inventory cost at the manufacturer due to the bullwhip effect. In this case, reducing costs at the retailer is critical as well as reducing the bullwhip effect at the manufacturer.

    <Figure 2> shows the optimal value of α as a function of ρ. In this figure, the system parameters are the same as in Figure 1. The x-axis refers to the value of ρ and the y-axis refers to the optimal value of α . Notice that when ρ = 0, we have ys, t = yn, t . Therefore, when ρ = 0, the optimal value of α becomes any number between 0 and 1.

    5. Final Remarks

    The alternative order-up-to level developed in Section 3 can be used to demonstrate the fundamental tradeoff of the bullwhip effect at an upstream facility with cost minimization at a current facility, with the goal of reducing system wide total expected inventory costs. Notice that this policy allows us to trade off the costs/benefits of not reacting to the most recent demand information vs. reacting to the most recent demand information. That is, the optimal value of α measures the optimal degree of neglect. For instance, when the optimal value of α is small (large), more (less) reaction to the most recent demand information is better in terms of minimizing system wide expected total inventory costs.

    Acknowledgement

    This study was supported by Research Fund of Dankook University in 2022.

    Figure

    JKSIE-46-4-93_F1.gif

    Costs as a Function of α

    JKSIE-46-4-93_F2.gif

    Optimal α as a Function of ρ

    Table

    Reference

    1. Erkip, N., Hausman, W.H., and Nahmias, S., Optimal Centralized Ordering Policies in Multi-Echelon Inventory Systems with Correlated Demands, Management Science, 1990, Vol. 36, No. 3, pp. 381-392.
    2. Kim, H., The Cost Impact of Incorrect Assumption in a Supply Chain, International Journal of Management Science, 2004, Vol. 10, No. 2, pp. 29-52.
    3. Kim, H. and Ryan J.K., The Cost Impact of Using Simple Forecasting Techniques in a Supply Chain, Naval Research Logistics, 2003, Vol. 50, pp. 388-411.
    4. Lee, H.L., So, K.C., and Tang, C.S., The Value of Information Sharing in a Two-Level Supply Chain, Management Science, 2000, Vol. 46, No. 5, pp. 626-643.
    5. Silver, E. and Peterson, S., Decision Systems for Inventory Management and Production Planning, 2nd Edition, John Wiley and Sons, New York, 1985.
    6. Urban, T.L., Reorder Level Determination with Serially-Correlated Demand, Journal of the Operational Research Society, 2000, Vol. 51, No. 6, pp. 762-768.
    7. Zinn, W., Marmorstein, H., and Charnes, J., The Effect of Autocorrelated Demand on Customer Service, Journal of Business Logistics, 1992, Vol. 13, No. 1, pp. 173-192.