A Two-Echelon Inventory Model with Lost Sales

Similar documents
A Two-Echelon Inventory Model for Single-Vender and Multi-Buyer System Through Common Replenishment Epochs

Evaluating Clustering Methods for Multi-Echelon (r,q) Policy Setting

SENSITIVITY ANALYSIS AND OPTIMIZATION OF DIFFERENT INVENTORY CONTROL POLICIES ALONG THE SUPPLY CHAIN

Essays on Joint Replenishment and Multi-Echelon Inventory Systems

Determination of a Desirable Inventory Policy in a three Echelon Multilayer Supply Chain with Normal Demand

Optimal Issuing Policies for Substitutable Fresh Agricultural Products under Equal Ordering Policy

Experiments with Protocols for Service Negotiation

Why do we have inventory? Inventory Decisions. Managing Economies of Scale in the Supply Chain: Cycle Inventory. 1. Understanding Inventory.

1 Basic concepts for quantitative policy analysis

Supplier selection and evaluation using multicriteria decision analysis

Ê«u U UO mm WF U WO UMB ««œù«re Ë W bmn «r * W??ö)«

A SIMULATION STUDY OF QUALITY INDEX IN MACHINE-COMPONF~T GROUPING

MULTIPLE FACILITY LOCATION ANALYSIS PROBLEM WITH WEIGHTED EUCLIDEAN DISTANCE. Dileep R. Sule and Anuj A. Davalbhakta Louisiana Tech University

Extended Abstract for WISE 2005: Workshop on Information Systems and Economics

A Scenario-Based Objective Function for an M/M/K Queuing Model with Priority (A Case Study in the Gear Box Production Factory)

Development of a Tool Management System for Energy Sector Company

Do Competing Suppliers Maximize Profits as Theory Suggests? An Empirical Evaluation

A Two-layer Time Window Assignment Vehicle Routing Problem

The Spatial Equilibrium Monopoly Models of the Steamcoal Market

A Multi-Product Reverse Logistics Model for Third Party Logistics

Appendix 6.1 The least-cost theorem and pollution control

STRATEGIES OF DELAY DELIVERY AND CONTROLLABLE LEAD TIME TO OPTIMIZE THE TOTAL COST IN SUPPLY CHAIN CONSIGNMENT INVENTORY

Planning of work schedules for toll booth collectors

Evaluating the statistical power of goodness-of-fit tests for health and medicine survey data

EVALUATING THE PERFORMANCE OF SUPPLY CHAIN SIMULATIONS WITH TRADEOFFS BETWEEN MULITPLE OBJECTIVES. Pattita Suwanruji S. T. Enns

Control of Yarn Inventory for a Cotton Spinning Plant: Part 1: Uncorrelated Demand. Russell E. King 1 Leigh Ann C. Brain 1 Kristin A.

A Group Decision Making Method for Determining the Importance of Customer Needs Based on Customer- Oriented Approach

Guidelines on Disclosure of CO 2 Emissions from Transportation & Distribution

Identifying Factors that Affect the Downtime of a Production Process

Aalborg Universitet. Positioning empty containers under dependent demand process. Dang, Vinh Quang; Yun, Won-Young; Kopfer, Herbert

Prediction algorithm for users Retweet Times

RECEIVING WATER HYDRAULICS ASSIGNMENT 2

6.4 PASSIVE TRACER DISPERSION OVER A REGULAR ARRAY OF CUBES USING CFD SIMULATIONS

FIN DESIGN FOR FIN-AND-TUBE HEAT EXCHANGER WITH MICROGROOVE SMALL DIAMETER TUBES FOR AIR CONDITIONER

Study on Productive Process Model Basic Oxygen Furnace Steelmaking Based on RBF Neural Network

THE IMPACTS OF DATA INACCURACY ON RETAILER S PERISHABLE INVENTORY

SIMULATION RESULTS ON BUFFER ALLOCATION IN A CONTINUOUS FLOW TRANSFER LINE WITH THREE UNRELIABLE MACHINES

Construction of Control Chart Based on Six Sigma Initiatives for Regression

Assignment II Design of Wastewater Discharge from the City of Gothenburg. Design of a Wastewater Discharge from the City of Gothenburg

Study on Multi-objective Optimization Model of Inventory Control and Supplier Selection Problem Under Uncertainty

Calculation and Prediction of Energy Consumption for Highway Transportation

On Effect of Delayed Differentiation on a Multiproduct Vendor Buyer Integrated Inventory System with Rework

Development and production of an Aggregated SPPI. Final Technical Implementation Report

Program Phase and Runtime Distribution-Aware Online DVFS for Combined Vdd/Vbb Scaling

Sporlan Valve Company

INCORPORATING WAITING TIME IN COMPETITIVE LOCATION MODELS: FORMULATIONS AND HEURISTICS 1

ENHANCING OPERATIONAL EFFICIENCY OF A CONTAINER OPERATOR: A SIMULATION OPTIMIZATION APPROACH. Santanu Sinha Viswanath Kumar Ganesan

Trade Policies for Intermediate Goods under International Interdependence

Simulation and Optimization of a Pull-Strategy in the. Order-Picking Area of a Distribution Warehouse. Hertzstr. 16, Karlsruhe, Germany

Evaluating The Performance Of Refrigerant Flow Distributors

INTEGER PROGRAMMING 1.224J/ESD.204J TRANSPORTATION OPERATIONS, PLANNING AND CONTROL: CARRIER SYSTEMS

A Cost Model of Partial Postponement Strategy of the Single-Period Product under Stochastic Demand

The Credit Risk Assessment Model of Internet Supply Chain Finance: Multi-Criteria Decision-Making Model with the Principle of Variable Weight

The Credit Risk Assessment Model of Internet Supply Chain Finance: Multi-Criteria Decision-Making Model with the Principle of Variable Weight

Reprint from "MPT-Metallurgical P(ant and Technology International" issue No. 2/1990, pages Optimization of. Tempcore installations for

A Bi-Objective Green Closed Loop Supply Chain Design Problem with Uncertain Demand

Research on chaos PSO with associated logistics transportation scheduling under hard time windows

Optimization of the Replenishment Strategy for the Supplier Based on Genetic Algorithm

Transient analysis of collinear cracks under anti-plane dynamic loading

Application of Ant colony Algorithm in Cloud Resource Scheduling Based on Three Constraint Conditions

A Novel Gravitational Search Algorithm for Combined Economic and Emission Dispatch

DESIGNING TWO-ECHELON SUPPLY CHAIN USING SIMULATION AND PRICING STRATEGY

Non-stationary stochastic inventory lot-sizing with emission and service level constraints in a carbon cap-and-trade system

EVALUATION METHODOLOGY OF BUS RAPID TRANSIT (BRT) OPERATION

The Role of Price Floor in a Differentiated Product Retail Market

Numerical Analysis about Urban Climate Change by Urbanization in Shanghai

Modeling and Simulation for a Fossil Power Plant

TRAFFIC SIGNAL CONTROL FOR REDUCING VEHICLE CARBON DIOXIDE EMISSIONS ON AN URBAN ROAD NETWORK

How to Review the Performance/Adequacy of ECV Observations? - Science Perspective from an ECV Producer -

RIGOROUS MODELING OF A HIGH PRESSURE ETHYLENE-VINYL ACETATE (EVA) COPOLYMERIZATION AUTOCLAVE REACTOR. I-Lung Chien, Tze Wei Kan and Bo-Shuo Chen

Evaluation Method for Enterprises EPR Project Risks

QUANTITY DISCOUNTS, CAPACITY DECISIONS, AND CHANNEL STRUCTURE CHOICES IN SUPPLY CHAINS JONATHAN EUGENE JACKSON JR.

The Application of Uninorms in Importance-Performance Analysis

A TABU SEARCH FOR MULTIPLE MULTI-LEVEL REDUNDANCY ALLOCATION PROBLEM IN SERIES-PARALLEL SYSTEMS

An Analysis on Stability of Competitive Contractual Strategic Alliance Based on the Modified Lotka-Voterra Model

Multi-UAV Task Allocation using Team Theory

emissions in the Indonesian manufacturing sector Rislima F. Sitompul and Anthony D. Owen

Consumption capability analysis for Micro-blog users based on data mining

The Order Selection and Lot Sizing Problem in the Make-to-Order Environment

WORKING PAPER MASSACHUSETTS

Dynamic Task Assignment and Resource Management in Cloud Services Using Bargaining Solution

A Comparison of Unconstraining Methods to Improve Revenue Management Systems

Varunraj Valsaraj, Kara Kockelman, Jennifer Duthie, and Brenda Zhou University of Texas at Austin. Original Version: September 2007.

Honorable Kim Dunning Presiding Judge of the Superior Court 700 Civic Center Drive West Santa Ana, CA 92701

ON THE BULLWHIP AVOIDANCE PHASE: SUPPLY CHAIN COLLABORATION AND ORDER SMOOTHING

LIFE CYCLE ENVIRONMENTAL IMPACTS ASSESSMENT FOR RESIDENTIAL BUILDINGS IN CHINA

The research on modeling of coal supply chain based on objectoriented Petri net and optimization

An analysis on Kendall Lee queueing system with non-preemptive priority at BRI Ahmad Dahlan Yogyakarta

Economia Aziendale online 2000 Web

SIMULATION-BASED RISK MEASUREMENT IN SUPPLY CHAINS

Sources of information

Capacity allocation and flexibility in primary care

Two Measures of Production System Flexibility and their Application to Identify Optimal Capacity Investment

An Application of MILP-based Block Planning in the Chemical Industry

Simulation of Steady-State and Dynamic Behaviour of a Plate Heat Exchanger

USE OF ENTERPRISES DATA IN ECONOMIC PLANNING

COMPARISON OF DEMAND DRIVEN AND PRESSURE DEPENDENT HYDRAULIC APPROACHES FOR MODELLING WATER QUALITY IN DISTRIBUTION NETWORKS

PSO Approach for Dynamic Economic Load Dispatch Problem

On Advantages of Scheduling using Genetic Fuzzy Systems

Derived Willingness-To-Pay For Water: Effects Of Probabilistic Rationing And Price

Transcription:

DEPARTMENT OF OPERATIONS RESEARCH UNIVERSITY OF AARHUS Publcaton no. 99/4 A Two-Echelon Inventory Model wth Lost Sales Jonas Andersson Phlp Melchors ISSN 1398-8964 Department of Mathematcal Scences Telephone: +45 8942 3188 Buldng 530, Ny Munkegade DK-8000 Aarhus C, Denmark

A Two-Echelon Inventory Model wth Lost Sales Jonas Andersson Dvson of Producton Management, Department of Industral Engneerng, Unversty oflund, Box 118, SE-22100 Lund,Sweden Phlp Melchors æ Department of Operatons Research, Unversty ofaarhus, Ny Munkegade, Buldng 530, DK-8000 Aarhus C, Denmark. September 1999 Abstract Almost all mult-echelon nventory models assume that demand not satsæed mmedately can be backordered. In some stuatons ths assumpton s not realstc. For example, t may be more representatve to model stockouts as lost sales when the retalers are n a compettve market and customers can easly turn to another ærm when purchasng the good. Assumng lost sales at the retalers, we consder a one warehouse several retalers nventory system. Usng the well-known METRIC-approxmaton as a framework, we present a heurstc for ændng cost eæectve base-stock polces. In a numercal study we ænd that the cost of the polces suggested by the heurstc s on average 0.40è above the cost of the ès, 1;Sè-optmal polcy. Keywords : Inventory, Mult-echelon, Lost sales, METRIC 1 Introducton Consder a two-echelon nventory system wth one central warehouse and an arbtrary number of retalers. See Fgure 1. The retalers face customer demand and replensh ther stocks from the central warehouse. The warehouse, n turn, replenshes ts stock from an outsde suppler. Evaluaton and optmzaton of control polces for such nventory systems have attracted massve nterest n the lterature. See, for example, Axsíater ë3ë for an overvew. In the exstng lterature dealng wth mult-echelon nventory control the prevalent assumpton s that complete backloggng of orders s allowed n case of stockouts. For example, Axsíater ë4ë shows how to exactly evaluate the performance for dæerent èr; nqè-polces when the retalers face compound Posson demand and nventores æ Correspondng author. Tel.:+45-89 42 35 36; fax: +45-86 13 17 69; e-mal: phlp@mf.au.dk 1

Retalers Warehouse Fgure 1: Mult-echelon nventory system are contnuously revewed. Cachon ë5ë gves an exact method for the perodc revew case wth dentcal retalers. In some stuatons the assumpton of complete backloggng may not be so realstc. For example, t may be more representatve to model stockouts as lost sales when the retalers are n a compettve market and customers can easly turn to another ærm when purchasng the good. For some reason the research dealng wth mult-echelon nventory models has focused manly on the backorder case and the number of models dealng wth lost sales s rather lmted. Anupnd and Bassok ë1ë consder a perodc revew twoechelon nventory system where a part of the unsatsæed sales at the retalers are lost. Snce the transportaton tme between the manufacturer and the retalers s zero, the optmal order polcy at each retaler s a base-stock polcy. The manufacturer carres lnear producton cost and no holdng cost. The retalers can agree to centralze ther stocks and the problem consdered s whether or not ths wll lead to an ncrease n total expected sales at the manufacturer. Nahmas and Smth ë8ë also consder lost sales n a mult-echelon envronment n a paper more closely related to ths paper. However, ther model dæers from ours n several mportant aspects. Frst, they consder perodc revew batch order polces. The model s more general snce they deal wth partal lost sales. Ths means that, wth probablty u, demand not satsæed mmedately, s lost, and wth probablty 1, u, t s satsæed later by a specal order. Moreover, for the model to be tractable they assume nstantaneous delveres from the warehouse to the retalers. For sngle-echelon nventory models the lost sales assumpton s more common. The exact cost for a sngle level nventory system facng Posson demand and æxed leadtmes was ærst gven by Hadley and Whtn ë6ë. Smth ë11ë demonstrates how to evaluate and ænd optmal ès, 1;Sè-polces for an nventory system wth zero replenshment costs and generally dstrbuted stochastc leadtmes. Recently Hll ë7ë showed that for the lost sales case the ès, 1;Sè-polcy s not necessarly optmal. 2

In ths paper we analyze a model for a one warehouse, multple retaler nventory system. Demand occurs only at the retalers and follows ndependent Posson processes. All leadtmes are assumed to be constant. All nstallatons use ès, 1;Sè-polces wth contnuous revew. It s assumed that backloggng of customer demand s not allowed. The analyss departs n one of the most wdely known mult-echelon nventory models, the METRIC-model developed by Sherbrooke ë10ë. In ts orgnal settng, t s assumed that stockouts at the retalers are completely backlogged. We demonstrate how the METRICmodel can be modæed to handle the lost sales case. Our approach gves an approxmate model whch s qute smple and eæcent from a computatonal pont of vew. Smulaton experments ndcate that the performance s very good. The outlne of ths paper s as follows: In Secton 2 we gve a detaled problem formulaton and pose all assumptons. Secton 3 gves the soluton procedure. The numercal results are gven n Secton 4, and n Secton 5 we gve some conclusons and pont out some possble drectons for future research. 2 Problem Formulaton The nventory system under consderaton conssts of one central warehouse and an arbtrary number of retalers. The retalers face Posson customer demand. No backloggng s allowed at the retalers. Consequently, the customers that arrve to a retaler that s out of stock wll become lost sales for the retaler. When stockouts occur at the warehouse, all demands from the retalers are fully backlogged and the backorders are ælled accordng to a FIFO-polcy. The transportaton tme between the warehouse and a gven retaler s assumed to be constant as well as the transportaton tme from the external suppler to the warehouse. The cost of a replenshment s assumed to be zero or neglgble compared to the holdng and stockout costs. The external suppler s assumed to have nænte capacty, whch means that the replenshment leadtme for the central warehouse s constant. All nstallatons use ès, 1;Sè-polces wth contnuous revew. Unts held n stock both at the warehouse and at the retalers ncur holdng costs per unt and tme unt. Moreover, a æxed penalty cost per lost customer s ncurred at the retalers. In ths paper we present a model for the consdered nventory system, whch can be used to evaluate the long-run average cost for dæerent polces wthn the class of ès, 1;Sè-polces. The objectve s to ænd the polcy that mnmzes the long-run average cost for the nventory system. Let us ntroduce the followng notaton: N = the number of retalers, ç = demand ntensty at retaler, =1; 2;:::;N, L = transportaton tme for the delveres from the warehouse to retaler, =1; 2;:::;N, L 0 =transportaton tme for the delveres from the external suppler to the warehouse, 3

S 0 = order-up-to level at the warehouse, S = order-up-to level at retaler, =1; 2;:::;N, h 0 = holdng cost rate at the warehouse, h = holdng cost rate at retaler, =1; 2;:::;N, ç = penalty cost for a lost sale at retaler, =1; 2;:::;N. We want to determne the total cost for the nventory system n steady state. Deæne TC = total cost for the nventory system per tme unt n steady state, C 0 = cost per tme unt for the warehouse n steady state, C = cost per tme unt for retaler n steady state, =1; 2;:::;N. Obvously, TC = C 0 + Our objectve s to determne a control polcy, S 0 ;S 1 ; ::; S N that mnmzes the total cost, TC. NX =1 C : è1è 3 Soluton Procedure In ths secton we ærst demonstrate how the total cost for dæerent control polces can be evaluated. For the backorder case the exact cost of the system can be derved by observng that any unt ordered by retaler s used to fulæll the S th demand. The cost can then be derved by condtonng on the arrval tme of the S th demand èwhch s Erlang dstrbutedè and the arrval of the ordered unt èsee Axsíater ë3ëè. In a lost sales envronment the correspondng observaton s that any unt ordered by retaler s used to fulæll the S + X th demand, where X s a random varable denotng the number of lost sales ncurred at the retaler durng the replenshment lead tme. X s obvously very hard to characterze and we have therefore chosen to focus on a heurstc rather than on the exact soluton. The analyss has many smlartes wth the analyss n Sherbrooke ë10ë. However, our assumpton of lost sales at the retalers destroys some of the nce propertes vald for the backorder model. The analyss of the warehouse, e.g., becomes more complex for the lost sales case. In the backorder case, all customers arrvng at the retalers generate demands at the warehouse mmedately at the arrval epoch, snce all retalers use contnuous revew ès, 1;Sè-polces. Consequently, the warehouse faces a Posson process wth ntensty ç 0 = ç 1 + ç 2 + æææ+ ç N. For the lost sales case ths s not true. When backorderng s not allowed, customer demands can be lost due to stockouts at the retalers. Therefore the demand at the warehouse s not Posson process anymore. 4

Another mportant dæerence compared wth the backorder case s that the order-upto level S at retaler, aæects the costs at all retalers and at the warehouse. In the backorder case S only aæects the local cost at retaler, snce the warehouse demand process s unaæected by the order-up-to levels at the retalers. For the lost sales case the order-up-to level aæects the number of lost sales and consequently, the demand process at the warehouse s not ndependent of the polces at the retalers. Therefore the order-up-to level at a certan retaler aæects the costs at all nstallatons n the nventory system. We wll ærst show how to evaluate the costs at the retalers gven a certan replenshment leadtme provded by the warehouse. We then show how to calculate the cost at the warehouse gven the demand ntensty from the retalers. Fnally we ntroduce an teratve procedure from whch we obtan the total cost for the nventory system. 3.1 Approxmate retaler cost As Sherbrooke ë10ë we use a queueng system analogy when evaluatng the costs for the retalers. For a retaler where backloggng s allowed, the number of outstandng orders towards the central warehouse s the same as the occupancy level n an M=G=1 queue. Recall that the customer demand s Posson and the replenshment leadtmes are stochastc, snce orders can be delayed due to stockouts at the central warehouse. For ths type of queue a famous theorem by Palm ë9ë states that the steady state occupancy level s Posson dstrbuted wth mean çl, where ç s the arrval rate and L s the mean servce tme. Palm's theorem holds for..d. servce tmes. The stochastc leadtmes n our case are evdently not ndependent, but f we dsregard ths correlaton we can approxmate the number of outstandng orders wth a Posson dstrbuton. Ths s the dea behnd the METRIC-approxmaton. When demand s lost, the queueng system of nterest s an M=G=S=S queue, wth S servers, each wth generally dstrbuted servce tmes and no queueng allowed. If the servce tmes are ndependent random varables wth mean ç L, Erlang's loss formula states the steady-state dstrbuton for the occupancy level as q S èjè= èçlè ç j =j! P S n=0 èç Lèn ç =n! for 0 ç j ç S where q S èjè = the probablty that j servers èout of Sè are occuped n steady state. Followng METRIC we approxmate the number of outstandng orders wth ths dstrbuton. Suppose that the mean leadtme for retaler s L ç and let q S èjè be the steady state probablty of j outstandng orders gven a desred base-stock level S. The expected number of lost sales per tme unt s clearly ç q S ès è and the expected number of unts 5

n stock s XS j=0 ès, jèq S èjè=s, ë1, q S ès èëç ç L : è2è The total relevant cost for retaler s therefore C ès ; L ç è=ç ç q S ès è+h çs, ë1, q S ès èëç L ç ç and the rate of demand from retaler whch s not lost s è1, q S ès èèç. The dervaton of the exact cost of a ès, 1;Sè lost sales sngle stage nventory system wth generally dstrbuted leadtmes was ærst presented by Smth ë11ë. He also proves that C ès ; ç L èsconvex n S for æxed ç L, whch means that the optmal value can be found by a local search routne. 3.2 Approxmate warehouse cost In the backorder case the demand process at the warehouse s a Posson process. In the lost sales case ths s not the case. If, for example, the base-stock level at a retaler s one, the smallest nterval between two successve demands from that retaler wll be the retaler leadtme. We wll gnore ths and approxmate the demand process at the warehouse wth aposson process wth mean æ. æ depends on how much demand s lost at the retalers and s determned as æ= NX =1 ç è1, q S ès èè è3è Snce we have a æxed determnstc leadtme L 0,we can ænd the average holdng cost ncurred at the warehouse as a functon of æ and S 0. X S 0 C 0 ès 0 ; æè = h 0 ès 0, jè èæl 0è j expè,æl 0 è j! j=0 We can also derve the mean delay due to stockouts at the warehouse by ærst calculatng B 0, the average number of backorders at the warehouse. B 0 = 1X j=s 0 +1 èj, S 0 è èæl 0è j expè,æl 0 è; j! We then apply Lttle's formula to obtan the average delvery delay, B 0 =æ. The mean leadtme for retaler s then è4è çl = L + B 0 =æ è5è Fnally, we obtan the total cost from è1è. 6

3.3 Overall soluton procedure We can now establsh the soluton procedure. The procedure enumerates over S 0. It can be shown that for a cost mnmzng soluton, S 0 can not be negatve. See, for example, Axsíater ë2ë. Consequently our procedure starts wth S 0 =0. Moreover, S 0 s bounded from above by an aborton crtera. We need the followng new notaton: C mn = mn S C ès ;L è = mnmum cost per tme unt for retaler n steady state when the leadtme, ç L, s equal to the transportaton tme, L, =1; 2;:::;N. S èkè = order-up-to level at retaler n teraton k. TC æ ès 0 è = mnmum value of TC gven a æxed value of S 0. Let us ærst consder two smple lemmas. The proofs can be found n the Appendx. The ærst lemma gves a lower bound for the retaler costs, and the second establshes two mportant propertes for the warehouse cost. Lemma 1. C mn s a lower bound for the retaler cost, C ès ; ç L è for all S and any çl él. Lemma 2. C 0 ès 0 ;ç 0 è ç C 0 ès 0 ; æè, for all S 0 and all æ ç ç 0. Moreover, C 0 ès 0 ;ç 0 è s convex n S 0. To construct an aborton crtera for the procedure, consder the cost functon TC lb ès 0 è=c 0 ès 0 ;ç 0 è+ NX =0 C mn : By Lemma 1 and Lemma 2, TC lb ès 0 èsalower bound for the cost functon, TC æ ès 0 è. Moreover, snce the cost functon TC lb ès 0 è s convex n S 0 the search over S 0 can be aborted when S 0 satsæes mn xçs 0 TC æ èxè ç TC lb ès 0 è: The aborton crtera s llustrated n Fgure 2. 7

TC æ ès 0 è TC lb ès 0 è mn xçs0 TC æ èxè S max 0 S 0 Fgure 2: Illustraton of the aborton crtera. The search for the optmal S 0 s aborted at S max 0. The soluton procedure can now be establshed as: STEP 0: Set S 0 =0andTC mn = 1. STEP 1: Set k = 0 and æ = ç 0. STEP 2: For each =1; 2;:::;N calculate ç L by è4è and è5è Let z æ = mn S C ès; ç L è and set S èkè =mnfsjc ès; ç L è=z æ g STEP 3: If ké0ands èkè =S èk, 1è for all =1; 2;:::;N then goto STEP 4, STEP 4: else calculate æ by è3è, set k P := k + 1 and goto STEP 2. Set TC æ n ès 0 è=c 0 ès 0 ; æè + =1 C ès èkè; L ç è. If TC æ ès 0 è étc mn then set TC mn = TC æ ès 0 è and let S opt 0 = S 0 and S opt = S èkè for =1; 2;:::;N. If TC mn étc lb ès 0 è then STOP, else set S 0 = S 0 + 1 and goto STEP 1. 4 Numercal Results In order to examne the eæectveness of the presented methodology we have performed a small numercal study. In total we consder 36 dæerent test problems wth æve dentcal retalers. For each test problem we ænd the best order-up-to levels accordng to our method. We also obtan the approxmate total holdng and stock out costs for the nventory system. The accuracy of these results are then evaluated by smulaton. Each smulaton conssts of 10 runs, each wth a run length of 100 000 tme unts. The result s a conædence nterval for the exact cost. We express the conædence lmts on a 95è sgnæcance level. A comparson between the total cost gven by our method and the total cost for the smulaton gves an ndcaton of how accurate our method s when estmatng the total cost for the nventory system. We also use smulaton to determne the optmal polcy for the system. The cost for ths polcy can then be compared wth the cost for the polcy determned by our technque, 8

to obtan an estmate for the performance of the method when optmzng the orderng polces. The polcy that we report as the optmal polcy s the polcy wth the lowest average cost. However, ths polcy does not necessarly domnate all the other polces when takng conædence ntervals nto consderaton. Moreover, we only search wthn polces where the order-up-to levels are dentcal for the retalers. The problem data and results can be found n Table 1. We only report the optmal polcy when t s dæerent from the one obtaned from our algorthm. From Table 1 we can see that our method performs rather well for all the consdered problems. It seems that we mostly tend to underestmate the total cost, especally n the problems wth hgh stockout costs at the retalers. Ths s due to the METRIC-approxmaton, where the stochastc leadtmes are replaced by ther averages when evaluatng the costs for the retalers. On average the method underestmates the costs wth 1.1 è. In 13 problems we can observe èon a sgnæcance level of 95èè that the method fals to ænd the optmal polcy. In 9 more problems the polcy suggested by our method does not have the lowest average cost accordng to the smulaton runs. However, n these cases the devatons are not sgnæcant on a 95è conædence level. In comparson to the optmal polces obtaned by smulaton, the ncrease n costs by usng the polces obtaned by our method s only 0.40 è, on average. In 16 of the 22 problems where we fal to ænd the true optmal polcy, the method merely underestmates the order-up-to level at the warehouse by a sngle unt. In one problem the warehouse order-up-to level s underestmated by two unts. In the other 5 problems where the optmal polcy s not found, the method tends to allocate more stock to the retalers and less stock to the warehouse than what s optmal from a cost perspectve. Fnally t seems that our methodology performs better f the warehouse leadtme s small compared to the transportaton tme from the warehouse to the retalers. In the 12 problem nstances wth L =1:5 the average cost ncrease, SCèCC s only 0.07è, whereas n the problems wth L =0:5, the correspondng ægure s consderable hgher, 0.67è. Ths behavor s due to the METRIC-approxmaton, where the stochastc replenshment leadtme facng a retaler s replaced by ts mean value. If the constant transportaton tme to the retalers s large compared to the warehouse leadtme, the stochastc delvery delays tend to have less relatve varaton and consequently the mpact of the METRICapproxmaton wll be smaller. 5 Conclusons and drectons for future research Ths paper presents a heurstc method for evaluaton and optmzaton of ès,1;sè-polces for a one warehouse, multple retalers nventory system. The evaluaton technque uses the well-known METRIC-approxmaton as a framework. From a computatonal pont of 9

vew the presented technque s very eæcent and smple. Numercal results also ndcate that the performance s qute good. Up to our knowledge, no paper s yet publshed, whch deals wth lost sales n a contnuous revew mult-echelon nventory settng. Moreover, the orgnal backorder METRICmodel ë10ë s one of the most wdely used mult-echelon nventory models. Our lost sales generalzaton makes the polcy evaluaton a bt more complex, snce we have to use an teratve procedure to obtan the cost. Stll, the model s rather smple and easy to mplement. Moreover, n many practcal stuatons lost sales s a reasonable way to model stockouts. Therefore our technque s also relevant for practtoners. In a research perspectve our model can form a framework n whch dæerent generalzatons can be consdered as optons for future research. For example, batch orderng polces and more general demand processes may be analyzed, stll usng the deas presented n ths paper. Generalzatons to perodc revew polces s also mportant. The dervaton of an exact evaluaton of costs seems to be a very dæcult problem to solve. Ths s a real challenge for future research. Appendx Proof for Lemma 1 We need to show that mn S C ès ;L è ç mn C ès ; L ç èforl ç L ç : è6è S Let l be an arbtrarly chosen leadtme, where L ç l ç ç L. Consder the cost C ès ;l è, where S s set to ts optmal value for each l. Obvously, C ès ;l è ç C ès, 1;l è for each l such that L ç l ç ç L. Start wth l = ç L and let l be contnuously lowered untl we reach l =L, whle S s set to ts optmal value for each l. Snce C ès ;l ès acontnuous functon of l for æxed S, t also s a contnuous functon of l when S s optmally chosen. Moreover, the fact that S mnmzes the cost C ès ;l è, mples that C ès ;l è ç C ès, 1;l è. Consequently, è6è follows f C ès ;l è ç C ès, 1;l è è @C ès ;l è @l ç 0 for L ç l ç ç L : è7è For notatonal reasons we omt the ndex from all varables. It can be shown that @CèS; lè @l =,hçè1, q S èsèè + çèhçl + ççèèq S ès, 1è, q S èsè+q S èsè 2 è: è8è Moreover, CèS; lè ç CèS, 1;lè mples that h hçl + çç ç qs,1 ès, 1è, q S èsè è9è 10

Let From è8è we have that A = è9è and è11è now gve h hçl + çç A = 1 çèhçl + ççè @C ès; lè : è10è @l, q S èsè, 1 æ + q S ès, 1è, q S èsè+q S èsè 2 è: è11è A ç èq S,1 ès, 1è, q S èsèèèq S èsè, 1è + q S ès, 1è, q S èsè+, q S èsè æ 2 = q S,1 ès, 1èq S èsè, q S,1 ès, 1è + q S ès, 1è = q S,1 ès, 1èq S èsè, èq S,1 ès, 1èq S èsèè = 0: Note that q S èsè ç 1. complete. Consequently, A ç 0 and therefore è7è holds and the proof s Proof for Lemma 2 Snce æ ç ç 0,we only need to show that @C 0èS 0 ;æè ç 0. The convexty ofc @æ 0 ès 0 ;ç 0 èns 0 follows, for example, from Axsíater ë3ë. @C 0 ès 0 ; æè @æ 0 =,h 0 L 0 expè,æl 0 è@ S0 + X S 0,1 =,h 0 L 0 expè,æl 0 è ç 0: j=0 S X 0 j=1 èæl 0 è j j! ès 0, jè æ èèæl 0 è j j!, èæl 0è èj,1è èj, 1è!!1 A References ë1ë R. Anupnd and Y. Bassok. Centralzaton of stocks: Retalers vs. manufacturer. Management Scence, 45:178í191, 1999. ë2ë S. Axsíater. Smple soluton procedures for a class of twoíechelon nventory problems. Operatons Research, 38:64í69, 1990. ë3ë S. Axsíater. Contnuous revew polces for multílevel nventory systems wth stochastc demand. In S.C. Graves, A.H.G. Rnnooy Kan, and P. Zpkn, edtors, Handbooks n OR & MS, Vol. 4, pages 175í197. Elsever Scence Publshers B.V., NorthíHolland, 1993. 11

ë4ë S Axsíater. Exact analyss of contnuous revew èr; Qèípolces n twoíechelon nventory systems wth compound Posson demand. Lund Unversty, Sweden, 1997. ë5ë G.P. Cachon. Exact evaluaton of batch-orderng nventory polces n twoíechelon supply chans wth perodc revew. Fuqua School of Busness, Duke Unversty, 1995. ë6ë G. Hadley and T.M. Whtn. Analyss of Inventory Systems. PrentceíHall, Englewood Clæs, NJ, 1963. ë7ë R.M. Hll. On the suboptmalty of ès,1;sè lost sales nventory polces. Internatonal Journal of Producton Economcs, 59è1-3è:377í385, 1999. ë8ë S. Nahmas and S.A. Smth. Optmzng nventory levels n a two-echelon retaler system wth partal lost sales. Management Scence, 40:582í596, 1994. ë9ë C. Palm. Analyss of the Erlang traæc formula for busy sgnal assgnment. Ercson Techncs, 5:39í58, 1938. ë10ë C.C. Sherbrooke. METRIC: A mult-echelon technque for recoverable tem control. Operatons Research, 16:122í141, 1968. ë11ë S.A. Smth. Optmal nventores for an ès, 1;Sè system wth no backorders. Management Scence, 23:522í528, 1977. 12

ç ç L S Calc cost Sm cost Spread Opt pol Cost Spread Cost dev 1 1 5 0.5 4,2 10.82 10.74 0.01 2 1 5 1.0 2,3 11.99 12.03 0.02 3 1 5 1.5 3,3 12.51 12.46 0.01 4 1 25 0.5 5,3 15.86 16.15 0.03 5 1 25 1.0 4,4 18.18 18.43 0.04 5,4 18.41 0.03 0.1è 6 1 25 1.5 3,5 19.94 20.16 0.06 4,5 20.12 0.04 0.2è 7 1 125 0.5 5,4 20.27 21.27 0.05 6,4 20.90 0.06 1.7è 8 1 125 1.0 5,5 23.64 24.16 0.06 6,5 24.09 0.06 0.3è 9 1 125 1.5 5,6 26.15 26.45 0.08 10 2 5 0.5 8,3 15.43 15.41 0.02 9,3 15.33 0.02 0.5è 11 2 5 1.0 6,5 17.11 17.27 0.02 8,4 17.16 0.02 0.6è 12 2 5 1.5 5,6 18.23 18.34 0.02 7,5 18.29 0.03 0.3è 13 2 25 0.5 8,5 21.56 22.52 0.04 9,5 22.30 0.04 1.0è 14 2 25 1.0 9,6 24.96 25.35 0.06 10,6 25.33 0.05 0.1è 15 2 25 1.5 7,8 27.49 27.92 0.08 8,9 27.90 0.04 0.1è 16 2 125 0.5 9,6 26.82 28.56 0.12 10,6 28.15 0.07 1.5è 17 2 125 1.0 9,8 31.84 32.75 0.08 18 2 125 1.5 10,9 35.49 36.07 0.11 19 1 5 0.5 4,1 16.96 16.41 0.02 20 1 5 1.0 2,2 17.55 17.51 0.02 21 1 5 1.5 2,2 18.14 18.04 0.03 22 1 25 0.5 4,3 27.50 27.99 0.03 5,3 27.96 0.02 0.1è 23 1 25 1.0 5,3 30.81 30.72 0.06 6,3 30.65 0.05 0.2è 24 1 25 1.5 4,4 33.07 33.11 0.05 25 1 125 0.5 7,3 36.86 37.21 0.12 8,3 36.92 0.08 0.8è 26 1 125 1.0 4,5 42.75 43.64 0.14 7,4 43.19 0.17 1.1è 27 1 125 1.5 6,5 46.83 46.92 0.10 28 2 5 0.5 6,3 24.21 24.42 0.02 8,2 24.14 0.03 1.2è 29 2 5 1.0 5,4 26.48 26.60 0.02 7,3 26.18 0.03 1.6è 30 2 5 1.5 6,4 27.61 27.43 0.02 31 2 25 0.5 10,4 36.44 37.10 0.06 11,4 36.95 0.05 0.4è 32 2 25 1.0 10,5 42.85 42.91 0.12 11,5 42.85 0.07 0.1è 33 2 25 1.5 8,7 46.74 47.11 0.08 34 2 125 0.5 11,5 47.24 48.92 0.08 13,5 48.21 0.12 1.5è 35 2 125 1.0 10,7 56.53 57.73 0.16 11,7 57.37 0.11 0.6è 36 2 125 1.5 11,8 63.77 64.19 0.19 12.8 64.05 0.10 0.2è Table 1: Numercal results. The optmal polcy s only reported when t s dæerent than the polcy suggested by our algorthm. 13