MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/4391970 4391970 Research Article Mathematical Model of (R,Q) Inventory Policy under Limited Storage Space for Continuous and Periodic Review Policies with Backlog and Lost Sales Singha Kanokwan 1 http://orcid.org/0000-0003-1028-0950 Buddhakulsomsiri Jirachai 1 Parthanadee Parthana 2 Hanne Thomas 1 School of Manufacturing Systems and Mechanical Engineering Sirindhorn International Institute of Technology Thammasat University Pathum Thani 12121 Thailand tu.ac.th 2 Department of Agro-Industrial Technology Faculty of Agro-Industry Kasetsart University Bangkok 10900 Thailand ku.ac.th 2017 10122017 2017 31 05 2017 25 10 2017 10122017 2017 Copyright © 2017 Kanokwan Singha et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper involves developing new mathematical expressions to find reorder point and order quantity for inventory management policies that explicitly consider storage space capacity. Both continuous and periodic reviews, as well as backlogged and lost demand during stockout, are considered. With storage space capacity, when on-hand inventory exceeds the capacity, the over-ordering cost of storage at an external warehouse is charged on a per-unit-period basis. The objective is to minimize the total cost, consisting of ordering, shortage, holding, and over-ordering costs. Demand and lead time are stochastic and discrete in nature. Demand during varying lead time is modeled using an empirical distribution so that the findings are not subject to assumptions of demand and lead time probability distributions. Due to the complexity of the developed mathematical expressions, the problems are solved using an iterative method. The method is tested with problem instances that use real data from industry. Optimal solutions of the problem instance are determined by performing exhaustive search. The proposed method can effectively find optimal solutions for continuous review policies and near optimal solutions for periodic review policies. Fundamental insights about the inventory policies are reported from a comparison between continuous review and periodic review solutions, as well as a comparison between backlog and lost sales cases.

1. Introduction

This paper focuses on the widely used (R,Q) inventory policy in distribution networks, where an order quantity Q is placed when the inventory position falls on or below the reorder point R at the time of a review. Inventory can be reviewed either continuously in real time or periodically at the end of some period. Continuous reviewing mostly requires a warehouse management system that supports it, whereas periodic reviewing can be performed without such a system. Mostly, the (R,Q) policy is designed for a continuous review, whereas an order-up-to level policy or basestock policy is for a periodic review. However, even without the technology to support a continuous review, some industrial users prefer the (R,Q) policy and use it with a periodic review. This is because some users would not like to change their order quantity every order. Also, their order quantity has to be a multiple of some pack sizes or container.

There are two types of behaviors of demands during stockout: lost sales, where the demands are lost; backlog, where the demands are fulfilled when inventories become available. This behavior depends on the characteristics of products, market, and relationships between suppliers and customers. Both types of demands during stockout are considered in this paper.

The inventory system under study is motivated by and therefore modeled after real distribution networks in various industries in Thailand. In such networks, internal storage space of an item in a warehouse or distribution center (DC) is preassigned and limited. When a replenishment order from a supplier arrives, there are instances when on-hand inventory exceeds the internal storage space. When this occurs, external space can be rented with an external storage cost charged on a per-unit per-period basis. In this system, the total cost of inventory management consists of four components: fixed ordering cost, inventory holding cost at internal storage, inventory holding cost at external storage (called over-ordering cost in this paper), and shortage cost. Item demands and replenishment lead time are stochastic and discrete in nature; therefore, the demands during varying lead time (i.e., the period of shortage risk) are of the same nature. Instead of assuming probability distribution(s) of demand, lead time, and demand during lead time, they are modeled using an empirical distribution. This is to overcome a problem found in real systems that many item demands do not follow widely used probability distributions.

The objective of the paper is to develop mathematical expressions for determining optimal or near optimal R and Q for the inventory system under study, so as to minimize the average total cost per period. The scope of the problem covers a wide range of real problems. The system considers four cost components (as mentioned above), discrete and varying demand and lead time that are not subject to assumptions of probability distribution, and four cases of review policies and demands during stockout. To the best of our knowledge, no mathematical expressions for a problem with this scope have been reported in the literature. The contributions of the paper are, therefore, the developed mathematical expressions, a method that implements them to solve a problem, and important fundamental insights gained from considering four cases of the problem.

There is a vast literature on inventory policy. Therefore, only relevant studies are included, that is, research studies that have involved storage space capacity. Those studies can be categorized into three groups: (1) storage space capacity is considered as a hard constraint, which means over-ordering storage is not allowed; (2) storage space capacity is a hard constraint, but with an additional cost of returning the over-ordered quantity to the supplier; and (3) storage space is a soft constraint, where an over-ordered amount is stored at an external, rented warehouse. These studies vary considerably in other aspects of the problem, such as inventory policies used, cost components considered, nature of demands during stockout, and the solution methodology.

In the first category, recent studies include Mandal et al. , Zhao et al. , Chou et al. , Zhao et al. , Zhong and Zhou , Pan et al. , Ghosh et al. , Beemsterboer et al. , and Zhang and Rajaram . Among them, studies that focused on (R,Q) policy are Zhao et al. , Zhao et al. , and Pan et al. . These studies examined a single product, continuous review (R,Q) policy, where demands are stochastic, shortages are backordered, and lead time is assumed to be constant. Among the three studies, only Zhao et al.  extended the problem to the case of multiple products. Heuristic algorithms were proposed to solve the problem: polynomial time algorithm in Zhao et al. , genetic algorithm in Pan et al. , and a solution approach consisting of an iterative part and local search in Zhao et al. . For this set of problems, since storage space is a hard constraint, this implies that the solution methodology proposed in these studies would attempt to set the inventory parameters to be within the available space, which is not the same as this paper.

For the second category where over-ordered quantity is returned to the supplier, Hariga  considered the following problem: an (R,Q) policy for a single item problem using continuous review: demands during stockout are backordered, and cycle length does not include the stockout period. The cost function was derived and solved using a simple economic order quantity (EOQ) based heuristic solution.

Studies that allowed external storage include Huang , Huang et al. , Hariga , Zhou et al. , Huang , Ouyang et al. , and Sana . The inventory holding cost at an owned warehouse is normally smaller than at a rented warehouse. Huang  and Huang et al.  developed a retailer’s inventory model when demand is known and constant and considered ordering cost and holding cost while shortages are not considered. The inventory policy was to find optimal time between order and order quantity, a (Q,T) policy. Huang  later extended their previous problem by adding the perishable characteristics of the product. Hariga  considered multiwarehouse systems to choose where to store excessive stock. The optimal solution consisted of the optimal order quantity and leased storage space. Some studies that focused on a problem where a supplier offers trade credit to a buyer to create an incentive to place larger orders are Zhou et al.  and Ouyang et al. . For this problem, the optimal decision for the buyer is the order quantity. The optimal decisions for the supplier are the optimal trade credit period to offer and the corresponding order quantity. Sana  presented an EOQ model for stochastic demand without considering lead time and with storage space capacity. Demands during stockout are lost. It can be seen that among recent studies where external rented storage is allowed, no studies have focused on determining (R,Q) under the same settings as this paper.

2. Problem Characteristics and Notation

The mathematical model developed in this paper adheres to the following problem characteristics or assumptions.

The inventory policy is an (R,Q) policy for a single continuously stocked item.

Two review policies are considered: continuous review and periodic review.

At the time of a review, if the inventory position (IP) falls on or below R, a replenishment order Q is placed to the supplier, and a fixed ordering cost is charged.

In the periodic review case, it is assumed that the IP is reviewed frequently enough so that a review period is shorter than the replenishment lead time.

The order arrives after a replenishment lead time that is stochastic and discrete.

Demands are stochastic and discrete in nature.

The item has a preassigned limited storage capacity. At the beginning of a cycle when a replenishment order arrives, if the on-hand inventory (OH) exceeds the storage space, over-ordered inventories will be kept at an external warehouse.

The over-ordered amount stored at an external warehouse is charged on a per-unit per-period basis.

It is assumed that the rate of external storage cost is no less than the rate of internal storage cost.

When demand arrives, if OH is less than the demand, shortage cost is charged on a per-unit basis.

Two types of shortages are considered. If there exist shortages, demands during stockout may be backlogged and may be lost.

An inventory cycle covers the elapsed time between replenishment orders, including the time the system has stock and the time the system is out of stock.

The objective is to determine the optimal or near optimal R and Q which minimizes the total inventory management cost that consists of fixed ordering cost, inventory holding cost for internal storage space, shortage cost, and storage cost at an external warehouse (i.e., the so-called over-ordering cost in this paper).

Given data include historical demand data, lead time data, all cost parameters, and internal storage space capacity. External storage space is always available.

Since both continuous review and periodic reviews as well as backlog and lost sales during stockout are considered, there are a total of four cases of the problem.

2.1. Notations

Q : replenishment order quantity, units

R : reorder point, units

C P : fixed ordering cost, THB/order

C H : inventory holding cost, THB/(unit-period)

C S : shortage cost, THB/unit

C O : over-ordering cost, THB/(unit-period)

X D : random demand, units/period

μ D : average demand, units/period

f ( x d ) : probability mass function of XD

X L : random lead time, periods

μ L : average lead time, periods

f ( x l ) : probability mass function of XL

X : random demand during varying lead time, units

μ : average demand during varying lead time, units, where μ=μDμL

f ( x ) : probability mass function of X

W : storage space capacity, units

E S C , ESP: expected number of shortages, for a given R, under continuous and periodic reviews, respectively, units

P S C , PSP: probability of a shortage, for a given R, under continuous and periodic reviews, respectively, in cases that demand during stockout is backlogged and is lost, respectively, %

O H e n d , OHordering, and OHbegin: random level of on-hand inventory at the end of a cycle, at the time of ordering, and at the beginning of a cycle, respectively, units

E O H e n d , EOHordering, and EOHbegin, EOH: the expected OH at the end of a cycle, at the time of ordering, at the beginning, and over the length of time the item is in stock in a cycle, respectively

IP: random level of inventory position at the time of ordering, units

E O C , L , EOP,L, EOC,B, and EOP,B: expected over-ordered amount at the beginning of an order cycle for different cases: continuous (C) and periodic (P) reviews, in combination with lost sales (L) and backlog (B), units

P O C , L , POP,L, POC,B, and POP,B: probability of over-ordering at the beginning of an order cycle for the four cases, %

T C C , L , TCP,L, TCC,B, and TCP,B: total inventory management cost per period for the four different cases, THB.

3. Methodology 3.1. Probability Distribution of Demand during Varying Lead Time

The stochastic and discrete behaviors of the demand and lead time are modeled using the distribution of demand during varying lead time. Given historical demand data and lead time data, one can list all possible values of random demand, XD=xd, and random lead time, XL=xl, and empirically derive their probability mass functions, f(xd) andf(xl). Then, the demand during a given lead time xl(X=xXL=xl) is the sum of the demands that occur during xl, where X=x=xd,1+xd,2++xd,xl=i=1xlxd,i, and xd,i are independently and identically distributed random variables with f(xd). In addition, its probability is the product of the probability of all demands P(X=xXL=xl)=f(xd,1)f(xd,2)f(xd,xl)=i=1lf(xd,i). Combining with f(xl) gives the probability of demand during varying lead time: f(x)=lP(X=xXL=xl)P(XL=xl).

3.2. Probability of Stockout and Probability of Over-Ordering

Given the values of Q, R, and W, in a replenishment cycle IP, OH, shortage amount, and over-ordered amount are shown in Figure 1, while the probability of stockout and the probability of over-ordering are shown in Figure 2.

Inventory level in a replenishment cycle.

Probabilities of stockout and over-ordering under distribution of demand during lead time.

In Figure 2, for illustration, suppose f(x) is normal. The probability of stockout would occur at the right tail of f(x) and towards the end of an order cycle, whereas the probability of over-ordering would occur at the left tail of f(x), at the beginning of the next order cycle, and depends on the value of X in the previous cycle.

3.3. Inventory Review and Demands during Stockout

The key difference between continuous and periodic reviews is that in a continuous review a replenishment order is placed when OH falls on or below R, whereas in a periodic review those events usually occur between review periods, so there is an average delay of half a period before the order is placed. That delay leads to an undershoot amount of μD/2 units, that is, the difference between R and OH at the point of ordering (i.e., the review time when an order is placed). The undershoot amount affects the probability of stockout and expected shortage as follows. With (1)–(4), it can be seen that both ES and PS are functions of R:(1)PSC=PX>EOHordering=PX>R=x=R+1xmaxfx,(2)ESC=x=R+1xmaxx-Rfx,(3)PSC=PX>EOHordering=PX>R-μD2=x=R-μD/2+1xmaxfx,(4)ESP=x=R-μD/2+1xmaxx-R+μD2fx.

The two cases of demands during stockout, lost sales and backlog, are considered. The difference between the two cases is the amount of OH that is replenished at the beginning of a cycle, which would be Q for lost sales and Q-ES for backlog.

3.4. Expected Over-Ordered Amount

To derive EO, first consider the expected OH at the end of a cycle, EOHend, for continuous review: EOHend=x=0R(R-x)f(x)=x=0xmax(R-x)f(x)-x=Rxmax(R-x)f(x)=R-μ+ESC. For a periodic review, EOHend is adjusted with the undershoot amount, EOHend=R-μD/2-μ+ESP. Hence, a random OH at different points in an order cycle can be derived by simply replacing μ with X. Expressions of the random OH for the four cases under different review policies and demands during stockout are shown in Table 1.

On-hand inventory in a cycle for four cases of (R,Q) inventory policies.

Point in a cycle Continuous review Periodic review
OH falls on or below R O H = R O H = R
At the time of ordering O H o r d e r i n g = R O H o r d e r i n g = R - μ D 2
At the end of a cycle O H e n d = R - X + E S C O H e n d = R - μ D 2 - X + E S P
At the beginning of the next cycle
Lost sales O H b e g i n = Q + R - X + E S C O H b e g i n = Q + R - μ D 2 - X + E S P
Backlog O H b e g i n = Q + R - X O H b e g i n = Q + R - μ D 2 - X

Now, EO can be derived from the difference between the random OH at the beginning of a cycle (OHbegin) and W, for the four cases of policies in (5). The expected inventory position at the time of ordering (IP) is used to generalize the formula for the four cases. Similarly, PO is generalized according to (6):(5)EO=0IP-W0x=0IP-WIP-W-xfxIP-W>0,(6)PO=0IP-W0x=0IP-WfxIP-W>0EO and PO for the four cases use the corresponding IP, where IPC,L=Q+R+ESC, IPP,L=Q+R-μD/2+ESP, IPC,B=Q+R, and IPP,B=Q+R-μD/2. With (5) and (6), it can be seen that both EO and PO are functions of Q and R.

3.5. Average Inventory Management Cost per Cycle

In a cycle, the inventory management cost has four components: ordering, shortage, holding, and over-ordering costs. The ordering cost and shortage cost are simply CP and CSES, respectively. The over-ordering cost is the storage cost at an external warehouse, which occurs at the beginning of an order cycle, and is subject to CO. Based on EO, that is, the expected over-ordered amount at the beginning of a cycle, the over-ordering cost is charged on average EO/2 for the length of EO/μD periods. Therefore, the expected over-ordering cost is CO(EO2/2μD). Finally, the holding cost at internal storage is charged, based on the average OH at the beginning and at the end of a cycle for the length of Q/μD, taking into account the expected over-ordered amount. The sum of the four cost components becomes the total inventory management cost. Dividing the total cost by the average cycle length (Q+ES)/μD, which takes into account the length of shortage period, gives the average total cost per period for the four cases in (7). To save space, EOH in the formula is used to denote the expected on-hand inventory during the period of Q/μD for the inventory holding cost component:(7)TC=μDCP+CSES+CHQ/μDEOH+CO-CHEO2/2μDQ+ES,where TC for each policy needs ES to be computed using (2) or (4), EO using (5), EOHC,L=Q/2+R-μ+ESC, EOHP,L=Q/2+R-μD/2-μ+ESP, EOHC,B=Q/2+R-μ, and EOHP,B=Q/2+R-μD/2-μ.

3.6. Mathematical Model to Find <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M132"><mml:mrow><mml:mi>R</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M133"><mml:mrow><mml:mi>Q</mml:mi></mml:mrow></mml:math></inline-formula>

The mathematical model to determine R and Q consists of two mathematical expressions that are derived by taking the first derivative of TC with respect to both R and Q, setting both equations to zero, and then solving them. The general form of dTC/dQ for the four cases is shown in (8). Solving the equations results in the general expression for Q in (9). Note that, from (5) and (6), since EOs are functions of IPs, which are functions of Q, we have dEO/dQ=PO:(8)dTCdQ=μDCH/μDQ/2+EOH+CO-CHEO/μDPO-TCQ+ES=0,(9)Q=2μDCHCP+CSES+CHESμDG+CO-CH1+CO-CHPO2CHEO22μD-ES-CO-CHPOCHEO.To save space, an expression G is used in (9) for the four cases: GC,L=μ-R-ESC/2, GP,L=μ+μD/2-R-ESP/2, GC,B=μ-R+ESC/2, and GP,B=μ+μD/2-R+ESP/2.

Similarly, the dTC/dR for the four cases is shown in (10). Solving them results in the expression for R in (11). Note that dTC/dR is based on the results of three derivatives: (1)dES/dR=-PS, which is based on (1)–(4) where PS and ES are both functions of R, (2)dEO/dR=PO(1-PS) for lost sales, which is based on (5) and (6) where EO (or IP) and ES are functions of R, and (3)dEO/dR=PO for backlog, where their IPs in (5) and (6) do not contain ES:(10)dTCdR=μDQ+ES-CSPS+CHH+CO-CHI+PSQ+ESTC=0.To save space, expressions H and I are used in (10), where HL=(Q/μD)(1-PS) for lost sales; HB=Q/μD for backlog; IL=PO1-PS(EO/μD) for lost sales; and IB=PO(EO/μD) for backlog:(11)R=a-b,where a=(μDQ+ES/PSCHQ)CSPS-CHJ-(CO-CH)K. Substitute J and K with JL=(Q/μD)(1-PS) and KL=PO(1-PS)(EO/μD) for lost sales, JB=Q/μD and KB=PO(EO/μD) for backlog, and b=(μD/CHQ)CP+CSES+(CHQ/μD)N+(CO-CH)(EO2/2μD). Substitute N with NC,L=Q/2-μ+ESC, NP,L=Q/2-μD/2-μ+ESP, NC,B=Q/2-μ, and NP,B=Q/2-μD/2-μ, accordingly.

3.7. Solving the Model for <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M175"><mml:mrow><mml:mi>R</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M176"><mml:mrow><mml:mi>Q</mml:mi></mml:mrow></mml:math></inline-formula>

Given all the input data, CP, CS, CH, CO, f(x), W, μD, μ, the two mathematical expressions for R and Q are recursive functions of each other. Therefore, they must be solved iteratively until their values converge. In each iteration, given an initial value Q, solving for a new R requires a search in the range of possible values of R because some input parameters in (11), including ES, PS, EO, and PO, are themselves functions of R. In other words, the search for Rnew begins with computing ES, PS, EO, and PO using R (a variable) and Q (fixed at an initial value for this search step). Then, compute Rnew using (11) and compute R-Rnew. It is important to note that the calculations of ES, PS, EO, and PO must be done in a sequence because EO uses ES as one of the inputs. While holding Q fixed and moving R in the direction that reduces R-Rnew, the best Rnew is the value that minimizes R-Rnew. Solving for Qnew using (10) is also performed in a similar manner. We now state the algorithm to solve for R and Q.

Step 1 (initialization).

Input data are CP, CS, CH, CO, W, d, l; initialize TC = M (a large number).

From xd and xl, compute μD=x-d, μL=x-l, and fit f(xd) and f(xl).

Compute all x for each xl and the probability: x=i=1xlxd,i, P(X=xXL=xl)=i=1lf(xd,i).

Compute μ=μDμL, xmax, and f(x)=lP(X=xXL=xl)P(XL=xl).

Compute initial Q=2μDCP/CH.

Step 2.

Compute new R

For R = 0 to xmax, the following steps are considered.

Compute PSC=x=R+1xmaxf(x) and ESC=x=R+1xmax(x-R)f(x), if it is a continuous review; or PSP=x=R-μD/2+1xmaxf(x) and ESP=x=R-μD/2+1xmax(x-R+μD/2)f(x), if it is a periodic review.

Compute EO=x=0IP-W(IP-W-x)f(x) if IP-W>0 and PO=x=0IP-Wf(x) if IP-W>0, where IP takes on one of the following: IPC,L=Q+R+ESC, IPP,L=Q+R-μD/2+ESP, IPC,B=Q+R, and IPP,B=Q+R-μD/2, depending on the case.

Compute Rnew(R)=a-b, where a=(μDQ+ES/PSCHQ)CSPS-CHJ-(CO-CH)K, JL=(Q/μD)(1-PS), and KL=PO(1-PS)(EO/μD) for lost sales, or JB=Q/μD and KB=PO(EO/μD) for backlog; b=(μD/CHQ)CP+CSES+(CHQ/μD)N+(CO-CH)(EO2/2μD), and N takes on one of the following: NC,L=Q/2-μ+ESC, NP,L=Q/2-μD/2-μ+ESP, NC,B=Q/2-μ, and NP,B=Q/2-μD/2-μ, depending on the case.

Compute R-Rnew(R).

Update R=Rnew(R) for minR-Rnew(R).

Step 3.

Compute new Q.

Compute ES with R from Step 2.

For Q = 0 to xmax, the following are considered.

Compute PO and EO, and compute Q according to (9).

Compute Q-Qnew(Q).

Update Q=Qnew(Q) with minQ-Qnew(Q), and update PO, EO.

Update TC=μDCP+CSES+(CHQ/μD)EOH+(CO-CH)(EO2/2μD)/(Q+ES) where EOH is one of the following: EOHC,L=Q/2+R-μ+ESC, EOHP,L=Q/2+R-μD/2-μ+ESP, EOHC,B=Q/2+R-μ, and EOHP,B=Q/2+R-μD/2-μ, depending on the case.

Step 4.

Repeat Steps 2 and 3 if TC improves from the previous value; Otherwise, stop.

4. Numerical Example

A numerical example that presents the developed method uses the following data: CP = 12.55 THB, CS = 4 THB/unit, CH = 0.012 THB/(unit-day), CO = 0.104 THB/(unit-day), W = 3,300 units, f(l) = {(1,0.365),(2,0.234),(3,0.257),(4,0.144)}, μL= 2.18 days, and μD = 360.27 units/day. Historical daily demand and lead time data are used to empirically construct f(d) and f(x). Note that these are real data of a fast moving item at a distribution center of a healthcare company in Thailand.

To evaluate the performance of the model, first the optimal solutions for all four cases are determined. Based on the data, the maximum possible demand during lead time is xmax = 8,116 units (demand data contain a few extreme outliers). For each policy, a complete enumeration of R and Q is performed, that is, evaluating (8,117)2 = 66,885,689 solutions. This takes approximately two hours of computing time, which is not a practical approach for an industrial user, and the distribution centers may have hundreds or thousands of items to manage.

Initial testing shows that the proposed method can generate optimal solutions for a continuous review policy and near optimal solutions for a periodic review policy. Table 2 contains iteration results of the proposed method. For periodic review cases, a detailed investigation of the iterations indicates that the search stops in one iteration, and a move from the first iteration solution yields worse total costs. Table 3 shows the results from using the developed method. For near optimal cases, the maximum off-optimal percentage of the total cost per period is only 0.0143% for a periodic review policy with backlog. The computational times for all four cases are also reported in Table 3.

Iteration steps to obtain the optimal R and Q.

Iteration C , B P , B C , L P , L
R Q R Q R Q R Q
Initial 868 868 868 868
( 1 ) 2,963 948 3,042 1,044 2,964 943 3,041 1,041
( 2 ) 2,912 995 2,915 989
( 3 ) 2,882 1,025 2,886 1,017
( 4 ) 2,864 1,042 2,870 1,032
( 5 ) 2,854 1,052 2,861 1,041
( 6 ) 2,848 1,058 2,856 1,046
( 7 ) 2,844 1,062 2,853 1,049
( 8 ) 2,841 1,065 2,851 1,051
( 9 ) 2,840 1,066 2,850 1,052
( 10 ) 2,849 1,053

Performance of the developed method, optimal and off-optimal.

 Continuous review with backlog (CPU time = 43.69 s) Continuous review with lost sales (CPU time = 48.27 s) Iter. R Q TC Iter. R Q TC 9 2,840 1,066 46.0347 10 2,849 1,053 46.1756 Opt. 2,840 1,066 46.0347 Opt. 2,849 1,053 46.1756 Off-opt 0% 0% 0% Off-opt 0% 0% 0% Periodic review with backlog (CPU time = 9.14 s) Periodic review with lost sales (CPU time = 9.13 s) Iter. R Q TC Iter. R Q TC 1 3,042 1,044 46.0413 1 3,041 1,041 46.1774 Opt. 3,020 1,066 46.0347 Opt. 3,029 1,053 46,1756 Off-opt 0.73% 2.06% 0.0143% Off-opt 0.40% 1.14% 0.0039%

From examining the optimal solutions, the differences between continuous review solutions and periodic review solutions are only in R, but not in Q, and it is not coincidental that the difference is equal to the expected undershoot amount of μD/2. Therefore, for a periodic review, one can simply solve the continuous review problem and add μD/2 to R, while Q remains unchanged, to obtain the optimal solution. Table 4 shows details of the optimal solutions.

Optimal solutions and their cost components.

Solution performance C , B P , B C , L P , L
R 2,840 3,020 2,849 3,029
Q 1,066 1,066 1,053 1,053

ES (units) 7.0931 7.0952 6.9565 6.9586
EO (units) 183.8970 183.8274 185.4028 185.3339
PS 0.0153 0.0153 0.0150 0.0151
PO 0.5077 0.6080 0.5102 0.6064

Fixed ordering cost (THB/cycle) 12.55 12.55 12.55 12.55
Shortage cost (THB/cycle) 28.3724 28.3808 27.8261 27.8343
Inventory holding cost (THB/cycle) 91.3130 91.3086 90.5150 90.5107
Over-ordering cost (THB/cycle) 4.8811 4.8774 4.9614 4.9577

Total cost (THB/cycle) 137.1165 137.1168 135.8525 135.8527

Cycle length (days) 2.9785 2.9786 2.9421 2.9421

Total cost per day (THB) 46.0347 46.0347 46.1756 46.1756

Annual cost (THB) 5,641.23 5,641.22 5,728.62 5,728.61

From the table, a comparison between continuous review and periodic review leads to the same findings for both backlog and lost sales cases. At the optimal solutions, (1) both policies have the same Q, (2) a periodic review has a higher R in order to compensate for the undershoot, but ES and shortage cost remain slightly higher, (3) holding cost and over-ordering cost are higher for a continuous review, and (4) a periodic review has a slightly higher total cost and longer cycle length, and together they lead to the same average total cost per day. In summary, there is no difference between continuous and periodic reviews, in terms of inventory management cost and performance. This is because of the inclusion of the undershoot amount to inventory OH in the calculation for the periodic policy; that is, a periodic review without considering undershoot would yield an inferior solution.

Another comparison between backlog and lost sales cases leads to the same findings for both a continuous review and a periodic review. Thus, it seems that there is no interaction between review policy and the nature of demands during stockout. At the optimal solutions, (1)R is higher for lost sales (as expected, because in our example, backlog cost and lost sales cost are the same, so it is only natural that the cycle service level would be set higher (less PS) to try to prevent shortages for lost sales), (2)Q is higher for backlog because it must satisfy demands during stockout, (3) higher R leads to higher EO and PO for lost sales, whereas higher Q leads to higher holding cost for backlog, and (4) although backlog has a higher total cost, it also has a longer cycle length, which leads to a lower average total cost per period. In summary, for the same cost components, backlog is more cost effective than lost sales.

5. Discussion

As previously mentioned, the problem under study is derived from a real inventory management problem at a distribution center of a healthcare company in Thailand. The operating conditions are as follows: The DC uses an (R,Q) policy under periodic (daily) review for a few hundred continuously stocked healthcare products, whose stocks are independently managed. For each item, storage space location and capacity are assigned. Replenishment lead time from the manufacturer varies between one and four days. This variation is caused by the difference between the day at which an order is placed and the delivery day. For example, an order placed on Monday and received by Friday has a four-day lead time, while an order placed on Tuesday and received by Friday has a three-day lead time, and so on. Items are classified by their movement (fast, medium, and slow). Most items’ demand data do not fit with known probability distributions. In a cycle where the order is received in a relatively short time and the demand during lead time is low, internal storage space capacity is exceeded. When this occurs, external storage space is rented to keep excessive inventory. In this system, shortages are mostly backlogged because the DC customers are retail stores of various types and sizes. The proposed mathematical model is applicable not only to the described system, but also to other systems with similar characteristics. Based on the authors’ experience, other industries in Thailand that face the same problems are distribution centers of retail businesses, authorized distributors of heavy equipment (agricultural tools and equipment, cooling systems, and compressors), spare part retailers and wholesalers, distribution networks for medicine, and medical supplies serving hospital systems.

6. Conclusions

This paper proposes a method to find optimal R and Q for both continuous and periodic reviews as well as for backlog and lost sales cases. The method is based on two newly developed mathematical expressions to find R and Q that minimize the average total cost per period and explicitly consider the over-ordering storage cost. Due to the complexity of the mathematical expressions, they must be solved iteratively. The method has a very promising performance based on the numerical test; that is, the solutions found are optimal for a continuous review and near optimal for a periodic review.

Three important fundamental insights are gained. First, one can solve the periodic review problem by simply solving the continuous review problem and adding the undershoot amount to the optimal R, while the optimal Q remains the same for both review policies. Continuous review and periodic review policies perform nearly the same, given that the undershoot amount is explicitly considered. For the same cost structure and cost components, cases of backlogged demand during stockout are more cost effective than those of lost demand. These findings are not subject to assumptions of the probability distribution of demand during varying lead time, because an empirical distribution fitted from real demand and lead time data is used.

Future research directions are as follows: (1) For cases of periodic review, our mathematical model can give only a near optimal solution. This calls for further investigation, which could lead to potential improvement of the method to handle periodic review cases. (2) The problem can be extended to multiple products sharing the same storage space.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank the industrial user, who provided the demand, cost, and operational data of their products to this study.

Mandal N. K. Roy T. K. Maiti M. Inventory model of deteriorated items with a constraint: a geometric programming approach European Journal of Operational Research 2006 173 1 199 210 MR2228121 10.1016/j.ejor.2004.12.002 2-s2.0-33746072112 Zhao X. Fan F. Liu X. Xie J. Storage-space capacitated inventory system with (r, Q) policies Operations Research 2007 55 5 854 865 10.1287/opre.1070.0394 MR2360952 Chou S.-Y. Julian P. C. Hung K.-C. A note on fuzzy inventory model with storage space and budget constraints Applied Mathematical Modelling 2009 33 11 4069 4077 2-s2.0-67649586300 10.1016/j.apm.2009.02.001 Zhao X. Qiu M. Xie J. He Q. Computing (r, Q) policy for an inventory system with limited sharable resource Computers & Operations Research 2012 39 10 2368 2379 10.1016/j.cor.2011.12.012 2-s2.0-84862813304 Zhong Y.-G. Zhou Y.-W. Improving the supply chain's performance through trade credit under inventory-dependent demand and limited storage capacity International Journal of Production Economics 2013 143 2 364 370 2-s2.0-84877577213 10.1016/j.ijpe.2012.07.013 Pan A. Hui C.-L. Ng F. An optimization of (Q, r) inventory policy based on health care apparel products with compound Poisson demands Mathematical Problems in Engineering 2014 2014 9 986498 10.1155/2014/986498 Ghosh S. K. Sarkar T. Chaudhuri K. A multi-item inventory model for deteriorating items in limited storage space with stock-dependent demand American Journal of Mathematical and Management Sciences 2015 34 2 147 161 2-s2.0-84939553972 10.1080/01966324.2014.980870 Beemsterboer B. Teunter R. Riezebos J. Two-product storage-capacitated inventory systems: A technical note International Journal of Production Economics 2016 176 92 97 2-s2.0-84962427386 10.1016/j.ijpe.2016.03.015 Zhang W. Rajaram K. Managing limited retail space for basic products: space sharing vs. space dedication European Journal of Operational Research 2017 263 3 768 781 MR3679230 10.1016/j.ejor.2017.05.045 Hariga M. A. A single-item continuous review inventory problem with space restriction International Journal of Production Economics 2010 128 1 153 158 2-s2.0-77957778506 10.1016/j.ijpe.2010.06.008 Huang Y.-F. An inventory model under two levels of trade credit and limited storage space derived without derivatives Applied Mathematical Modelling 2006 30 5 418 436 10.1016/j.apm.2005.05.009 Zbl1182.90007 2-s2.0-33144467467 Huang Y.-F. Lai C.-S. Shyu M.-L. Retailer's EOQ model with llimited storage space under partially permissible delay in payments Mathematical Problems in Engineering 2007 2007 2-s2.0-36248983064 10.1155/2007/90873 90873 Hariga M. A. Inventory models for multi-warehouse systems under fixed and flexible space leasing contracts Computers & Industrial Engineering 2011 61 3 744 751 2-s2.0-80053185357 10.1016/j.cie.2011.05.006 Zhou Y.-W. Zhong Y. Li J. An uncooperative order model for items with trade credit, inventory-dependent demand and limited displayed-shelf space European Journal of Operational Research 2012 223 1 76 85 MR2956535 10.1016/j.ejor.2012.06.012 2-s2.0-84864434959 Huang T.-S. The optimal buyer's ordering policy and payment policy for deteriorating items with limited storage capacity under supplier credit Journal of Information and Optimization Sciences 2013 34 6 417 454 MR3196361 10.1080/02522667.2013.867736 Ouyang L. Y. Ho C. H. Su C. H. Yang C. T. An integrated inventory model with capacity constraint and order-size dependent trade credit Computers & Industrial Engineering 2015 84 133 143 Sana S. S. An EOQ model for stochastic demand for limited capacity of own warehouse Annals of Operations Research 2015 233 383 399 10.1007/s10479-013-1510-5 MR3400081