Optimal stock allocation in single echelon inventory systems subject to a service constraint
|
|
- Dale Fisher
- 6 years ago
- Views:
Transcription
1 Abstract number: Optimal stock allocation in single echelon inventory systems subject to a service constraint Annalisa Cesaro Dario Pacciarelli Dipartimento di Informatica e Automazione, Università Roma Tre, Via della Vasca Navale, 79, Rome I-00146, Italy cesaro@dia.uniroma3.it pacciarelli@dia.uniroma3.it phone: POMS 20th Annual Conference Orlando, Florida U.S.A. May 1 to May 4, 2009 This paper deals with spare optimization in multi-location inventory systems with single item and repairable spare parts. Lateral and emergency shipments occur in response to stockouts. A continuous review basestock policy is assumed for the inventory control of the spare parts. The resulting model is a queueing network with blocking, which can be studied using a Markov chain modeling approach. The objective is to minimize the total costs for inventory holding, lateral transshipments and emergency shipments subject to a target level of operational availability of the whole system. We model the stock allocation problem as a non linear integer program and exploit its special structure to design a primal-based and dual-based algorithm. Our solution procedure is based on a convergent Lagrangian and objective level cut method. Computational experiments, carried on practical data from an airport equipment maintenance context show that this method finds an optimal allocation for all practical instances for which the optimum is known. 1. Introduction Single echelon inventory systems are experiencing an increasing interest in practice, in particular for the management of expensive spare parts. In such a context, the supply chain involves at least three actors: equipment users, logistics companies and equipment suppliers. The users need spare parts to carry on their business without interruptions. Intermediate logistic companies are in charge of replenishing spare parts in the short term, by granting 1
2 the contractual service level to the users at minimum cost. The suppliers are responsible for supplying new components and/or repaired items to the logistic companies. As observed by several authors, see e.g. by [10], the logistics of spare parts differs from those of other materials in several ways. Equipments may have remarkable costs, long repairing times and sporadic failures. The latter are difficult to forecast and may cause relevant financial effects, due to the economical implications of a lack of equipment at the operational sites. In such cases, continuous review policies are typically adopted to reduce both reaction time to stockouts and inventory levels [2, 9]. Several heuristic procedures can be found in the literature for allocating spares to warehouses in a single echelon context with complete pooling [3, 12, 21]. This paper proposes and evaluates a new solution procedure for spare parts allocation in this context. Our work is motivated by a practical problem faced by a large Italian logistics company. The company handles 17 warehouses supporting the daily activity of 38 airports spread over the Italian territory. Stock levels are currently determined with the VARIMETRIC algorithm of Sherbrooke (Sherbrooke 2004), based on a stiff hierarchic structure. However, in operation lateral transshipments take place between stocking points whenever there is an emergency request for a spare part, using couriers and overnight carriers to rapidly transfer the required items. The company is therefore interested in determining the potential savings deriving from explicit inclusion of lateral transshipments in the model. To this aim, a new method for stock level definition and spare parts allocation is necessary, and this need motivates the present paper. Specifically, this paper develops a new method for finding the optimal stock allocation in a single echelon inventory system with continuous review, lateral and emergency transshipments, non-negligible transshipment times, subject to a minimum level of operational availability (e.g. the fraction of time during which all operational sites are working). Our optimization method is independent from the specific method chosen for assessing the operational availability of a given allocation, and therefore it applies as well to different contexts with different performance indicators. This paper is organized as follows. Relevant notation is described in Section 2. Section 3 describes the single echelon one-for-one ordering model with complete pooling. In Section 4 the spare parts allocation problem is formulated as a non-linear integer program, and a convergent Lagrangian and objective level cut method is described for finding an optimal allocation. Computational experiments are presented in Section 5, based on practical data from the airport maintenance context. Some conclusions follow in Section 6. 2
3 2. Notation In order to formally define the problem, let us introduce the following notation. Let A = {1, 2,..., a} be a set of operational sites (e.g., airports) where working equipments are located. We assume that operational sites are grouped on a regional basis, with a warehouse of spare parts for each regions. Let W = {1, 2,..., w} be the set of regional warehouses. Let s i be the number of spare parts to allocate to each warehouse i W, and with s = (s 1,..., s w ) an allocation of spares to warehouses, i.e., the vector of decision variables. Let µ be the service rate of a server at warehouse h, T hi be the transfer time for a spare from warehouse h to warehouse i and T s (j, h) be the substitution time, i.e., the time needed to transfer a spare part to the site j A from the warehouse h W and to physically replace the failed item. Let λ jh be the rate of requests for a spare part from site j to warehouse h, λ h = j A λ jh be the aggregated arrival rate of failure processes at warehouse h and λ h the effective arrival rate at warehouse h including lateral transshipments. Let P B (s) be the network blocking probability for a given allocation s, i.e., the probability that a failure occurs at some site and no warehouse can satisfy the spare demand. Given an allocation s, let π hi (s) be the probability of the event: there are no spares in warehouse h W and i W is the closest warehouse with available spares (i.e., every warehouse l such that T hl < T hi, including the case l = h, is in stockout condition). Let n = (n 1,..., n w, n w+1 ) be a vector representing the state of the network, in which n i is the number of outstanding requests at warehouse i W, and n w+1 is the number of outstanding emergency requests to the external supplier. In our system, n w+1 > 0 only if n i = s i for all i W. Let p(n) be the probability of being in state n for the whole warehouses network. Denote p i (s i ) as the marginal probability of having s i outstanding orders at warehouse i, i.e., the probability of having no stock available at warehouse i. We also denote with MT T R the mean time to repair, i.e. the average replenishment time of the external supplier, with M T BF the mean time between failures, with OS the order and ship time, with MCMT the mean corrective maintenance time and with OA the operational availability of all the a sites. Let c h be the inventory holding cost for warehouse h, c tr ij be the cost for a lateral transshipment from warehouse j to warehouse i, in stockout condition, and c em be the emergency transshipment cost. 3
4 3. Problem description In our model, a logistic company wants to find the stock level s i of each warehouse i W such that a minimum level of service is provided at the operational sites and the overall cost is minimum. Costs are related to inventory holding, transshipments and emergency transshipments. Given an allocation s of spares to warehouses, the model we use to compute the level of service is a single item, single echelon, w-locations, continuous review, one-for-one replenishment policy inventory system, which allows for lateral and emergency transshipments with complete pooling and non-negligible transshipment times in operational availability computations. The resulting model is a queueing network with blocking, in which each warehouse i W is viewed as a queue with s i servers and no buffer and the external supplier is viewed as a queue with infinitely many servers. 3.1 Performance indicator The contractual service level to grant is the operational availability (OA) of all operational sites for each item, computed as in [17]. OA = MT BF MT BF + MCMT. (1) In equation (1), MCMT is the average time occurring from the failure of an item to its physical substitution. This is the substitution time if the spare is available at the regional warehouse. If no spares are locally available, the request is forwarded to the closest warehouse with available spares and MCMT increases by the deterministic transfer time between the two warehouses. When no warehouse has spares available, MCMT equals the substitution time plus the replenishment time from the external supplier. The MCMT can be therefore computed as follows: MCMT = h W j A λ jh T s (j, h)+ h W (λ h i W π hi(s) T ih )+ ( h W λ h) P B (s) (MT T R + OS) We observe that the first term h W A j=1 λ jh T s (j, h) of Equation (2) only depends on the failure process and on the distance between the sites and their respective regional warehouses. In other words, it does not depend on the specific spare parts management (2) 4
5 policy being used. Therefore, for sake of simplicity we assume it negligible in our model and omit its computation in the rest of this paper. As for the quantity π hi (s) defined in Section 2, we assume that a strict deterministic nearest chosen neighbor rule is adopted for sourcing a lateral transshipment, as in Kukreja [12]. Differently form [12], we use equation (3) to compute this value, which directly follows from the definition of π hi (s). π hi (s) = (1 p i (s i )) (p l (s l )). (3) l:t hl <T hi 3.2 Assumptions We use the Poisson distribution for the demand process, which is a typical assumption for modeling low demand processes [18]. The MTBF of an equipment depends on its age and on other exogenous agents, such as the damp, the temperature and other operational conditions. Therefore, we use specific values for each site. The replenishment time of the external supplier is a random variable, exponentially distributed, with known mean value MTTR. The capacity of the supplier repair shop is assumed to be infinite. It follows that also the number of replenishments from the external supplier follows the Poissonian distribution. These common assumptions make possible to use the Markovian analysis for modeling the multi-dimensional inventory system. However, as shown by [12], the results obtained with these assumptions holds under less restrictive hypothesis on the replenishment distribution. 3.3 Multi-dimensional markovian model In order to compute the operational availability of the system under study for a given allocation s, we must compute the probability p(n) of each state n = (n 1,..., n w, n w+1 ) of the associated Markov chain, defined in Section 2. To this aim we use a model similar to that of Wong et al. [19]. The main difference is that we include the external supplier explicitly in the model, while in [19] a failure of a part occurring when all warehouses are in stockout condition is lost. In case of blocked network, the first repaired item returned by the external supplier is used for replacing a failed item in some operative site, if any. There are direct transitions among states just in case of a single arrival event (i.e. a request for a spare at some warehouse) or a single departure event (i.e. the replenishment of a repaired item by the external supplier). 5
6 Figure 1: The infinite state space Markov chain for a system of two warehouses (a) and finite space state model (b) and the aggregated birth death model (c) Let e i be a vector with w + 1 elements, all equal to 0 but the element in position i that is equal to 1, and let ψ(h, i) be equal to 1 if (n i < s i and n l = s l for each l W such that T lh < T ih, included l = h) and be equal to 0 otherwise. With this notation, n+e i is the state of the Markov chain representing an arrival at the i-th warehouse (if i W and n i < s i ), due either to a failure in the i-th service region or to a re-forwarded request from some other warehouse h in stockout conditions, i.e., such that ψ(h, i) = 1. Similarly, n e i is the state with a departure from the i-th warehouse (if i W and n i > 0). For the external supplier, n + e w+1 represents a new emergency request (if n i = s i for each i W ) and n e w+1 represents the fulfillment of an emergency request (if n w+1 > 0). The transition rate q(n, m) from state n towards state m = n ± e i and n ± e w+1 is as follows. q(n, n + e i ) = λ i + h W {i} ψ(h, i) λ h, for i W and n i = 0, 1,..., s i 1; q(n, n + e w+1 ) = i W λ i, if n i = s i i W ; q(n, n e i ) i µ, for i W and n i > 0 and n w+1 = 0; q(n, n e w+1 ) = w+1 i=1 n i µ, for n i = s i i W and n w
7 Figure 1(a) shows the Markov chain for two warehouses, the first having two spares and the second having three available spares. Steady state probabilities can be computed for each state in the Markov chain by solving a linear system. Note that, to this aim the proposed Markov chain is equivalent to the one shown in figure 1(b), having a finite number of states. In the latter Markov chain, all states in which all warehouses are in stockout condition are grouped in one single state with different departure transition rates. Let n B be the state in the original model where i W, (n B ) i = s i, (n B ) w+1 = 0 and p(n) be equal to the probability of being in state n. Denote P B (s) as in Section 2. In the equivalent model the modified departure rates of the single state corresponding to n B and all other blocked states in the original model become: q(n B, n B e i ) = (n B ) i µ F for i W and (n B ) i = s i and (n B ) w+1 = 0 and by using F = p(n B) P B (s). p(n B) and P B (s) may be easily computed, as stated in Theorem 1 in [6],by using the equivalence with the proposed Markov chain model with a simple birth death model, as shown in figure 1(c). Specifically let S = i W s i W i be the total stock level in the network, and let ρ = λ i. Given a set W of µ warehouses, with total stock level S, in which the service process is exponentially distributed with average rate µ for each server and the demand flow to warehouse i W is Poissonian with average rate λ i, the blocking probability of all warehouses is P B (s) = 1 S 1 k=0 ρ k k! e ρ. Finally p(n B ) is equal to ρn B N B! e ρ, where N B = w+1 i=1 (n B) i. By using the state probabilities computed solving the above finite state Markov chain model, we are able to compute p i (s i ) for each i W. p i (s i ) is the marginal blocking probability for the i-th warehouse, as defined in Section 2. It may computed as follows: p i (s i ) = p(n) (4) n N:n i =s i where N is the Markov chain state space. 4. Solution procedure In this section we face the main problem of defining an optimal allocation of spares to warehouses. We model the spare part optimization problem as a non linear program and propose a primal-based, dual-based algorithm for solving it at optimality. To this aim we study a Lagrangian relaxation of our model and use it within an iterative procedure known as convergent Lagrangian and objective level cut method [15]. 7
8 4.1 Problem formulation The Spares Allocation Problem is the problem of finding an allocation s which minimizes the overall cost for inventory holding, lateral and emergency transshipments, subject to a constraint on the minimum operational availability of the system, defined by equation (2) with the assumption A h W j=1 λ jh T s (j, h) = 0. Letting L be the minimum operational availability level that must be satisfied by a feasible allocation, the Spares Allocation Problem P 0 can be formulated as an integer program with non-linear objective function and a single non-linear constraint. With the notation defined in Section 2 we have: Problem P 0 : min s.t. : w i=1 ch s i + λ i j W π ij(s) c tr ji + λ i P B (s) c em w i=1 [λ i j W π ij(s) T ji + λ i P B (s) (MT T R + OS)] > (1 L) MT BF L (5) In the following subsections we focus on the computation of an optimal allocation s. The solution method consists of an iterative procedure. At the k-th step, the Lagrangian relaxation of a revised problem P k is computed. Problem P k is obtained from the original problem P 0 by adding constraints on the objective value, such that the optimum still satisfies all the constraints and the difference between the Lagrangian lower bound of problem P k and the current best upper bound is reduced. This method is called convergent Lagrangian and objective level cut method [15]. Let us briefly recall the main concepts of the method. Given a set X R n and r+1 functions t : R n R and g j : R n R, j = 1,..., r, the primal problem P is as follows: In our case: g(s) = t(s) = min{t(x) : x X, g j (x) 0 for j = 1,..., r}. (6) w c h s i + λ i i=1 (1 L) MT BF L j W π ij (s) c tr ji + λ i P B (s) c em, (7) [ w λ i ] π ij (s) T ji + λ i P B (s) (MT T R + OS), (8) i=1 j W X = {s i 0, s i integer i W }. (9) 8
9 The Lagrangian dual D is defined as: max{q(γ) : γ 0}, (10) where q(γ) = inf{l(x, γ) : x X} and L(x, γ) = t(x) + γ g(x). Function q(γ) is called the Lagrangian relaxation of the primal problem P, L(x, γ) is the Lagrangian function and γ is the transpose of vector γ. The following properties hold [5]. The Lagrangian relaxation q(γ) is a concave function. Given any γ 0 and any feasible solution x for the primal problem, q(γ) t(x) holds. This property is called weak duality. Let x be a solution to the Lagrangian relaxation, i.e., such that q( γ) = L( x, γ). If x is feasible for the primal problem and γ g( x) = 0, then x is an optimal solution for the primal problem. In this case strong duality holds, i.e., q( γ) = t( x). The conditions γ g( x) = 0 are called complementary slackness conditions. Let x be an optimal solution for the primal problem and γ an optimal solution for the Lagrangian dual. The difference t(x ) q(γ ) (non-negative by weak duality) is called the duality gap between the primal problem and the Lagrangian dual. Let UB be an upper bound to t(x ). We denote UB q(γ ) a duality bound between the primal problem and the Lagrangian dual. It is clear than a duality bound is always larger than or equal to the duality gap. If t = q(γ ), then strong duality holds. Unfortunately, strong duality rarely occurs in integer programming. In this paper we use a method to strengthen the primal formulation in order to achieve strong duality. In subsection we deal with the exact computation of the Lagrangian relaxation q(γ) of problem P 0. Our approach is similar to that of Wong [19]. Then, in subsection 4.1.2, we use the sub-gradient method to find the maximum q(γ ) of the Lagrangian relaxation. A simple upper bound to the optimal primal solution is obtained heuristically in subsection Finally, in subsection 4.2, we deal with the computation of an optimal solution s to Problem P The Lagrangian relaxation For a given γ, the Lagrangian relaxation q P0 (γ) of problem P 0 is obtained by relaxing the operational availability constraint: 9
10 { w q P0 (γ) = min s {L(s, γ)} = min s i=1 ch s i + λ i j W π ij(s) c tr ji + λ i P B (s) c em + [ +γ (1 L) MT BF w L i=1 (λ i )]} j W π ij(s) T ji + λ i P B (s) (MT T R + OS) (11) Let ZP 0 be the optimum of problem P 0. Given an allocation s, let Z P0 (s) be the value of the objective function of problem P 0 for the given s, let OA(s) be the operational availability of the network under s, and MCMT (s) be the mean corrective maintenance time under s. From the properties of the Lagrangian relaxation shown in the previous section, it follows that [19]: Z P 0 q P0 (γ) for any γ 0 Z P 0 max γ 0 q P0 (γ) If for some γ 0 the optimal solution to problem min s {L(s, γ)} is s and s is a feasible allocation for problem P 0, then (Z P0 (s ) ZP 0 ) γ (MCMT (s (1 L) MT BF ) ( )) L If for some γ 0 the optimal solution s to problem min s {L(s, γ)} is feasible for P 0 and γ ( (1 L) MT BF L then s is an optimal solution to problem P 0. ) MCMT (s ) = 0, Let us now focus on the computation of q P0 (γ) = min s {L(s, γ)} for a given γ. Problem (11) can be rewritten as: { w q P0 (γ) = min s i=1 [c h s i + λ i j W π ij(s) (c tr ji γ T ji )+ ] } (12) +λ i P B (s) (c em γ (MT T R + OS)) + γ (1 L) MT BF L Note that the last term does not depend on s and can be omitted for the purpose of finding s. Moreover, the quantity [ ] w i=1 c h s i + λ i P B (s) (c em γ (MT T R + OS)) only depends on the overall stock level S = w i=1 s i (see Section 3.3). We can therefore define the functions f 1 (S) = w i=1 [ ] c h s i + λ i P B (s) (c em γ (MT T R + OS)) f 2 (s) = w λ i i=1 j W π ij (s) (c tr ji γ T ji ) 10
11 and rewrite the problem as: min s {f 1 (S) + f 2 (s)} (13) This problem can be solved at optimality by using a procedure by Wong [19]. The idea of the procedure is to introduce the quantities f 3 (S) = min s {f 2 (s) : j W s i = S} and to solve the problem min S {f 1 (S) + f 3 (S)} for increasing values of S, starting from S = The sub-gradient optimization procedure When problem (13) is solved for a given γ, we can easily compute q P0 (γ), a lower bound on Z P 0. In order to solve the Lagrangian dual (10) we use the sub-gradient optimization method. Note that, letting sˆγ = arg min s {f 1 (S) + f 2 (s)}, the function g(sˆγ ) defined by equation (8) provides a sub-gradient of function q(γ) for γ = ˆγ. The sub-gradient method generates a sequence of dual feasible points obtained from an initial value γ 0, for k = 0, as follows: γ k+1 = [γ k + o k g k ] + (14) Here, g k = g(s γ k) denotes the sub-gradient at point γ k, [ ] + denotes projection on the set {γ : γ 0, q P0 (γ) > }, and o k is a positive scalar stepsize computed as in [19]. o k = σ k q(γk ) Ẑ g(s γ k) 2 2 In this formula, Ẑ is the the best upper bound to problem P 0 so far. It is initially computed with the algorithm described in subsection and then updated during the procedure execution. The value σ k is a scalar chosen between 0 and 2. A sketch of the sub-gradient method is reported in algorithm Upper-bound computation A simple upper bound to Z P 0 is computed by allocating spare parts only to those warehouses i W with positive demand λ i > 0 and by avoiding concentration in just few warehouses. Simulation experiments [7] show that this allocation policy is effective in many practical contexts. Therefore its aim is that of allocating parts giving preference to warehouses with larger demand. 11
12 Algorithm 1 Sub-gradient Algorithm 1: k 0, γ 0 0, Ẑ solution of heuristics described in subsection : while sub-gradient direction NOT close to 0 and lower bound value not respecting any Lagrangian value constraint (refer to section 4.2) do 3: for all i W do 4: s i = 0, s i = 0 5: end for 6: min t(0), refer to equation 7; 7: go = true 8: while go == true do 9: S = S + 1, compute f(s), referring to 13; 10: if f(s) min then 11: go = false 12: else 13: for all possible stock allocations ŝ with ŝ ŝ w = S and respecting any objective level constraint (refer to section 4.2) do 14: calculate f(s) and g(ŝ) 15: if f(s) + g(ŝ) < min then 16: min = f(s) + g(ŝ), s i = ŝ i i = 1,..., w 17: end if 18: if ŝ feasible for problem P 0 5 and t(ŝ) < Ẑ, by using equation 7 then 19: Ẑ t(ŝ) 20: end if 21: end for 22: end if 23: end while 24: k = k + 1; 25: compute γ k and the new sub-gradient direction, refer to 14, by using s ; 26: end while 12
13 Our heuristic procedure greedily allocates one spare at a time and checks the operational availability of system. The warehouses i W with λ i > 0 are ordered for decreasing value of λ i, and a spare is allocated to each warehouse in this order until the contractual level of the Operational Availability OA is met. 4.2 Convergent Lagrangian and objective level cut method In order to achieve an optimal allocation s for problem P 0 we follow the convergent Lagrangian and objective level cut method described in [15] with slight modifications to adapt the method to our single echelon spare part optimization problem. The objective level cut is used to eliminate the duality gap and thus to guarantee the convergence of the Lagrangian method on a revised domain Objective level cut method for the spare optimization problem In this section, we describe the solution scheme of the convergent Lagrangian and objective level cut method for solving problem P 0. This method has been proposed and applied in [15] for separable non linear integer programming and polynomial 0-1 programming. In [15] an integer objective value is required for optimality, whereas the objective function of problem P 0 may not be integer. However, we show that the method can be adapted to deal with problem P 0. The key ideas consist of replacing exact lower and upper bound values with their respective integer rounded up values and adding Lagrangian level cuts. We can prove that the properties of the Lagrangian and objective level cut method still applies, and that the proposed algorithm, sketched in algorithm 2, finds an optimal solution to Problem P 0, or proves infeasibility, in at most UB 1 LB iterations (see theorem 3). UB 1 and LB 1 are respectively the upper and lower bound values computed by solving the conventional Lagrangian dual problem (10). Therefore, the efficiency of the algorithm depends on the solution algorithm chosen for solving the Lagrangian dual in the revised form, where objective and Lagrangian level cuts have been added. The algorithm starts solving the conventional Lagrangian dual problem. In case of a positive duality gap, a revised dual problem is iteratively formulated and solved until the duality gap reduces down to zero. Specifically, a lower objective level cut is added to problem P 0, thus reducing the duality gap (see theorems 1 and 2), and a lower Lagrangian cut is imposed to the Lagrangian dual problem (10), so that the Lagrangian dual value is increased with respect to the current lower bound and the duality gap is further reduced (see theorem 2). 13
14 Note that the optimal solution is always valid for both cuts. Let l and u be a lower and an upper bounds to problem (5), respectively. Then the revised problem is as follows: Revised problem P 0 (l, u): min t(s) = w i=1 [ch s i + λ i j W π ij(s) c tr ji + λ i P B (s) c em ] s.t. : w i=1 [λ i j W π ij(s) T ji + λ i P B (s) (MT T R + OS)] > (1 L) MT BF L s X(l, u) = {s : s i 0, s i integer, i W, l t(s) u } Clearly, P 0 (l, u) is equivalent to P 0 when l t u. The Lagrangian relaxation of P 0 (l, u) is: Revised problem q (l,u) P 0 (γ): min s {L(s, γ)} = min s { w i=1 ch s i + λ i j W π ij(s) c tr ji + λ i P B (s) c em + +γ [ ( w i=1 [λ i j W π ij(s) T ji + λ i P B (s) (MT T R + OS)])+ (1 L) MT BF + } L (15) s X(l, u) {s : s i 0, s i integer, i W, L(s, γ) > l } The Lagrangian dual problem of P 0 (l, u) is then given as: Revised Lagrangian dual problem D(l, u): max q (l,u) P 0 (γ) s.t. : γ 0 (16) (17) In case of duality gap the algorithm starts solving the revised Lagrangian dual problems by using as initial lower and upper bound values the ones computed by solving the conventional Lagrangian dual problem (10). Specifically, u = q(γ) and l = t(s), computed according to equations (10) and (7), the latter firstly generated heuristically (as described in subsection 4.1.3) and updated eventually during the subgradient search. Then, the method updates these lower and upper bound values iteratively, by adding new objective and Lagrangian cuts, computed by solving problem (17). Given a Lagrange multiplier ˆγ, and since the dual function q (l,u) P 0 (γ) described in (17) is a concave function, its maximum can be computed through sub-gradient optimization, as described in subsection A sketch of the overall 14
15 Algorithm 2 Convergent Lagrangian and objective level cuts algorithm 1: LB 0 and UB 0 Ẑ (Ẑ feasible solution for P 0, found by applying heuristics described in subsection 4.1.3) 2: k = 0, go = true 3: while go do 4: k = k + 1, solve D( LB k, UB k ), for computing the best γ k ; 5: LB q ( LB k, UB k ) (γ k ) and UB k Ẑ (Ẑ best feasible solution for P 0 ( LB k, UB k ), found during the dual search) 6: if LB k = UB k then 7: go = false 8: else 9: go = true 10: end if 11: if UB k < LB k then 12: go = f alse, problem infeasible 13: end if 14: end while method is provided by Algorithm 2. It can be proved [15] that at least one infeasible solution will be cut when adding a new objective cut. Theorem 1 The following properties hold: 1. Let γ (l, u) denote the optimal solution to the Lagrangian dual (10). The optimal dual value q P0 (γ, l, u) is a non decreasing function of l 2. If l t u then q P0 (γ ) q P0 (γ, l, u) t. Moreover let σ = max{t(s) t(s) < t, x X(l, u)}. If σ < l t, then γ = 0 and q P0 (γ, l, u) = t 3. For l < t, we have q P0 (γ, l, u) l Proof. The proof is similar to that provided in [14] for an analogous theorem. The only difference is the introduction of rounded up values for the upper and lower bounds, which does not affect the proof. Moreover we have: Theorem 2 Let Q P0 (γ ) be the set of feasible solutions of q P0 (γ ). If q P0 (γ, l, u) < t, then min{f(s) s Q P0 (γ ), g(s) > 0} q P0 (γ, l, u) holds. 15
16 Proof. The proof is similar to that provided in [14] for an analogous theorem. The only difference is the introduction of rounded up values for the upper and lower bounds, which does not affect the proof. Theorem 1 implies that the revised dual problem provides a value which is greater or equal to the current lower bound. Theorem 2 guarantees that at least one infeasible solution is cut by an objective cut higher than the current lower and upper bounds l and u. The introduction of rounded up values for lower and upper bounds allows us to prove the following theorem: Theorem 3 Algorithm 2 either finds an optimal solution of problem 5 or reports an infeasibility of problem 5 in at most UB 1 LB iterations. Proof. From Algorithm 2, it follows that LB k t holds. If the algorithm stops at step 12 problem 5 is infeasible. Optimality holds when Algorithm 2 stops at step 7. If the algorithm does not stop at iteration k, then LB k+1 LB k + 1 holds at each step. Notice that for any k, LB k t UB k + 1 holds. Therefore, in at most UB 1 LB 1 iterations UB k = LB k must hold. If the algorithm does not stop before UB 1 LB iterations, then it must stop at UB 1 LB th iteration either with an infeasibility or reaching the optimal solution. 5. Case study from the corrective airport maintenance context In this section we evaluate the effectiveness of our method for a practical problem arising in the airport maintenance context. We analyze 30 practical instancesprovided by an Italian logistics company supporting the activity of the 38 civil airports spread over the Italian territory. The company handles 17 warehouses and manages the overall processes of purchasing, holding, ensuring that the overall reliability of safety equipments is always within contractual limits. The aim of the company is therefore to grant the prescribed quality of service at minimum cost. While the company currently follows a two echelon policy for spare part management, in this paper we evaluate the potential of the single echelon policy, which is generally acknowledged to achieve better performance in this context. 16
17 The instances differ in the mean demand at each warehouse and in the number of warehouses with positive demand. The replenishment time of an item has been set equal to three months for all items and warehouses. Table 1 reports the main data for our instances. Item Cost Replenishment time MTBF number airports number items months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h months h Table 1: Characteristic data for 30 items A proven optimum has been computed for each instance by using a complete enumeration technique. Starting from S = 1, and then increasing S by one at each step, the costs associated to all possible stock allocations s with s s w = S are evaluated until a feasible stock allocation is found for some S. We can conclude at this point that no better solutions can be obtained for larger values of S, since purchasing and holding costs are increasing 17
18 with S and dominate the transshipment costs. For each instance, we compute the same stock allocation either by using our algorithm or by using the enumerative technique. In all cases, the algorithm finds a solution with the same proven minimum cost of the enumeration technique. Each optimal solution has been found within one minute of computation time, whereas the computation time required by the enumeration technique is sometimes larger than 30 minutes. 6. Conclusions In this paper we propose and evaluate a solution methodology for optimizing inventory stock allocation of repairable spare parts in a single echelon, w-locations system, where lateral and emergency shipments occur in response to stockouts. We model our problem as a nonlinear integer program and apply a primal based, dual based algorithm for allocating the spare parts optimally. We evaluate this technique by using practical data from airport corrective maintenance context, by determining the stocking policy for items supporting the daily activities of 38 Italian airports. The technique finds optimal stock levels for all test instance for which the proven optimum can be found via total enumeration, thus demonstrating its validity. References [1] Airports Council International (ACI). Worldwide and Regional Forecasts, Airport Traffic , (2005). [2] Alfredsson, P. and J. Verrijdt. Modeling emergency supply flexibility in a two echelon inventory system. Management Science 45, (1999) [3] Axsater,S.. Modelling emergency lateral transshipments in inventory systems. Management Science 6, (1990) [4] Axsater,S.. Inventory control. Springer Verlang (2006). [5] Bertsekas, D.P.. Nonlinear programming. Athena Scientific second edition (2003). [6] Cesaro, A., Pacciarelli, D.. Performance assessment for single echelon airport spare part management. Technical Report University Roma Tre (2009). 18
19 [7] Cesaro, A., Pacciarelli, D.. Evaluation of peaked lateral transshipment in inventory system subject to a service constraint. MSOM 2007 conference (June 18 19,2007). [8] Dada, M.. A two echelon inventory system with priority shipments. Management Science 38, (1992) [9] Grahovac, J., A. Chakravarty. Sharing and Lateral Transshipment of Inventory in a Supply Chain with Expensive Low-Demand Items. Management Science 47, (2001) [10] Huiskonen, J.. Maintenance spare parts logistics: special characteristics and strategic choices. Int. J. Production Economics 71, (2001) [11] Kennedy, W.J., J. Wayne Patterson, and L. D. Fredendall. An overview of recent literature on spare parts inventories. Int. J. Production Economics 76, (2002) [12] Kukreja, A., C.P. Schmidt, D.M. Miller. Stocking decisions for low-usage items in a multilocation inventory system Management Science 47, (2001) [13] Kukreja, A., C.P. Schmidt A model for lumpy demand parts in a multi-location inventory system with transshipments. Computers and Operations Research 32, (2005) [14] Li, D., X.L. Sun. Towards strong duality in integer programming. J. Glob. Optim 35, (2006) [15] Li, D., J. Wang, X.L. Sun. Computing exact solution to nonlinear integer programming: Convergent Lagrangian and objective level cut method. J. Glob. Optim 39, (2007) [16] Li, D., X.L. Sun. Nonlinear integer programming. Springer s INTERNATIONAL SE- RIES (2006) [17] Sherbrooke, C.C.. Optimal inventory modeling of systems: multi-echelon techniques. Wiley (2004). [18] Taragas, G., D. Vlachos. Effectiveness of stock transshipment under various demand distributions and nonnegligible transshipment times. Production And operations Management 11, (2002)
20 [19] Wong, H.,D. Catrysse,D. Van Houdeusden, Inventory pooling of repairable spare parts with non-zero lateral transshipment time and delayed lateral transshipments. European Journal of Operational Research 165 (2005) [20] Wong, H.,D. Catrysse,D. Van Houdeusden. Stocking decisions for repairable spare parts pooling in a multi-hub system. Int. J. Production Economics 93 (2005) [21] Wong, H.,G.J. Van Houtum, D. Catrysse,D. Van Houdeusden. Simple, efficient heuristics for multi-item multi-location spare parts systems with lateral transshipments and waiting time constraints. Journal of the Operational Research Society 56 (2005),
A New Approximate Evaluation Method for Two-Echelon Inventory Systems with Emergency Shipments
A New Approximate Evaluation Method for Two-Echelon Inventory Systems with Emergency Shipments Erhun Özkan, Geert-Jan van Houtum, Yasemin Serin Beta Working Paper series 363 BETA publicatie WP 363 (working
More informationOnline Appendices: Inventory Control in a Spare Parts Distribution System with Emergency Stocks and Pipeline Information
Online Appendices: Inventory Control in a Spare Parts Distribution System with Emergency Stocks and Pipeline Information Christian Howard christian.howard@vti.se VTI - The Swedish National Road and Transport
More informationSolving Dual Problems
Lecture 20 Solving Dual Problems We consider a constrained problem where, in addition to the constraint set X, there are also inequality and linear equality constraints. Specifically the minimization problem
More informationTutorial: System-Oriented Inventory Models for Spare Parts Geert-Jan van Houtum
YEQT VIII, 3-5 November, 2014 Eindhoven Tutorial: System-Oriented Inventory Models for Spare Parts Geert-Jan van Houtum Literature Seminal paper: Sherbrooke, C.C., METRIC: A multi-echelon technique for
More informationOn service level measures in stochastic inventory control
On service level measures in stochastic inventory control Dr. Roberto Rossi The University of Edinburgh Business School, The University of Edinburgh, UK roberto.rossi@ed.ac.uk Friday, June the 21th, 2013
More informationA TWO-ITEM TWO-WAREHOUSE PERIODIC REVIEW INVENTORY MODEL WITH TRANSSHIPMENT
Annals of Operations Research manuscript No. (will be inserted by the editor) A TWO-ITEM TWO-WAREHOUSE PERIODIC REVIEW INVENTORY MODEL WITH TRANSSHIPMENT Koushik S. Ramakrishna Moosa Sharafali Yun Fong
More informationIn the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight.
In the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight. In the multi-dimensional knapsack problem, additional
More informationExpected Time Delay in Multi-Item Inventory Systems with Correlated Demands
Expected Time Delay in Multi-Item Inventory Systems with Correlated Demands Rachel Q. Zhang Department of Industrial and Operations Engineering, University of Michigan, Ann Arbor, Michigan 48109 Received
More informationImproved Algorithms for Machine Allocation in Manufacturing Systems
Improved Algorithms for Machine Allocation in Manufacturing Systems Hans Frenk Martine Labbé Mario van Vliet Shuzhong Zhang October, 1992 Econometric Institute, Erasmus University Rotterdam, the Netherlands.
More informationA A Multi-Echelon Inventory Model for a Reparable Item with one-for. for- one Replenishment
A A Multi-Echelon Inventory Model for a Reparable Item with one-for for- one Replenishment Steve Graves, 1985 Management Science, 31(10) Presented by Hongmin Li This summary presentation is based on: Graves,
More informationTHIELE CENTRE. The M/M/1 queue with inventory, lost sale and general lead times. Mohammad Saffari, Søren Asmussen and Rasoul Haji
THIELE CENTRE for applied mathematics in natural science The M/M/1 queue with inventory, lost sale and general lead times Mohammad Saffari, Søren Asmussen and Rasoul Haji Research Report No. 11 September
More informationStochastic Optimization
Chapter 27 Page 1 Stochastic Optimization Operations research has been particularly successful in two areas of decision analysis: (i) optimization of problems involving many variables when the outcome
More information3E4: Modelling Choice. Introduction to nonlinear programming. Announcements
3E4: Modelling Choice Lecture 7 Introduction to nonlinear programming 1 Announcements Solutions to Lecture 4-6 Homework will be available from http://www.eng.cam.ac.uk/~dr241/3e4 Looking ahead to Lecture
More informationInventory optimization of distribution networks with discrete-event processes by vendor-managed policies
Inventory optimization of distribution networks with discrete-event processes by vendor-managed policies Simona Sacone and Silvia Siri Department of Communications, Computer and Systems Science University
More informationDesigning the Distribution Network for an Integrated Supply Chain
Designing the Distribution Network for an Integrated Supply Chain Jia Shu and Jie Sun Abstract We consider an integrated distribution network design problem in which all the retailers face uncertain demand.
More informationQueueing Theory I Summary! Little s Law! Queueing System Notation! Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K "
Queueing Theory I Summary Little s Law Queueing System Notation Stationary Analysis of Elementary Queueing Systems " M/M/1 " M/M/m " M/M/1/K " Little s Law a(t): the process that counts the number of arrivals
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationPBW 654 Applied Statistics - I Urban Operations Research
PBW 654 Applied Statistics - I Urban Operations Research Lecture 2.I Queuing Systems An Introduction Operations Research Models Deterministic Models Linear Programming Integer Programming Network Optimization
More informationRouting. Topics: 6.976/ESD.937 1
Routing Topics: Definition Architecture for routing data plane algorithm Current routing algorithm control plane algorithm Optimal routing algorithm known algorithms and implementation issues new solution
More informationCoordinated Replenishments at a Single Stocking Point
Chapter 11 Coordinated Replenishments at a Single Stocking Point 11.1 Advantages and Disadvantages of Coordination Advantages of Coordination 1. Savings on unit purchase costs.. Savings on unit transportation
More informationLinear Programming: Chapter 5 Duality
Linear Programming: Chapter 5 Duality Robert J. Vanderbei September 30, 2010 Slides last edited on October 5, 2010 Operations Research and Financial Engineering Princeton University Princeton, NJ 08544
More information5/15/18. Operations Research: An Introduction Hamdy A. Taha. Copyright 2011, 2007 by Pearson Education, Inc. All rights reserved.
The objective of queuing analysis is to offer a reasonably satisfactory service to waiting customers. Unlike the other tools of OR, queuing theory is not an optimization technique. Rather, it determines
More informationA TANDEM QUEUEING SYSTEM WITH APPLICATIONS TO PRICING STRATEGY. Wai-Ki Ching. Tang Li. Sin-Man Choi. Issic K.C. Leung
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 00X Website: http://aimsciences.org pp. X XX A TANDEM QUEUEING SYSTEM WITH APPLICATIONS TO PRICING STRATEGY WAI-KI CHING SIN-MAN CHOI TANG
More informationThe Single and Multi-Item Transshipment Problem with Fixed Transshipment Costs
The Single and Multi-Item Transshipment Problem with Fixed Transshipment Costs Reut Noham, Michal Tzur Department of Industrial Engineering, Tel-Aviv University, Tel Aviv, Israel Received 1 October 2013;
More informationMULTI-PERIOD MULTI-DIMENSIONAL KNAPSACK PROBLEM AND ITS APPLICATION TO AVAILABLE-TO-PROMISE
MULTI-PERIOD MULTI-DIMENSIONAL KNAPSACK PROBLEM AND ITS APPLICATION TO AVAILABLE-TO-PROMISE Hoong Chuin LAU and Min Kwang LIM School of Computing National University of Singapore, 3 Science Drive 2, Singapore
More informationQueueing Theory. VK Room: M Last updated: October 17, 2013.
Queueing Theory VK Room: M1.30 knightva@cf.ac.uk www.vincent-knight.com Last updated: October 17, 2013. 1 / 63 Overview Description of Queueing Processes The Single Server Markovian Queue Multi Server
More informationPerformance Evaluation of Queuing Systems
Performance Evaluation of Queuing Systems Introduction to Queuing Systems System Performance Measures & Little s Law Equilibrium Solution of Birth-Death Processes Analysis of Single-Station Queuing Systems
More informationOn the Approximate Linear Programming Approach for Network Revenue Management Problems
On the Approximate Linear Programming Approach for Network Revenue Management Problems Chaoxu Tong School of Operations Research and Information Engineering, Cornell University, Ithaca, New York 14853,
More informationMotivation. Lecture 2 Topics from Optimization and Duality. network utility maximization (NUM) problem:
CDS270 Maryam Fazel Lecture 2 Topics from Optimization and Duality Motivation network utility maximization (NUM) problem: consider a network with S sources (users), each sending one flow at rate x s, through
More informationA PARAMETRIC DECOMPOSITION BASED APPROACH FOR MULTI-CLASS CLOSED QUEUING NETWORKS WITH SYNCHRONIZATION STATIONS
A PARAMETRIC DECOMPOSITION BASED APPROACH FOR MULTI-CLASS CLOSED QUEUING NETWORKS WITH SYNCHRONIZATION STATIONS Kumar Satyam and Ananth Krishnamurthy Department of Decision Sciences and Engineering Systems,
More information56:171 Operations Research Fall 1998
56:171 Operations Research Fall 1998 Quiz Solutions D.L.Bricker Dept of Mechanical & Industrial Engineering University of Iowa 56:171 Operations Research Quiz
More informationAppendix A.0: Approximating other performance measures
Appendix A.0: Approximating other performance measures Alternative definition of service level and approximation. The waiting time is defined as the minimum of virtual waiting time and patience. Define
More informationExercises Stochastic Performance Modelling. Hamilton Institute, Summer 2010
Exercises Stochastic Performance Modelling Hamilton Institute, Summer Instruction Exercise Let X be a non-negative random variable with E[X ]
More informationA Duality-Based Relaxation and Decomposition Approach for Inventory Distribution Systems
A Duality-Based Relaxation and Decomposition Approach for Inventory Distribution Systems Sumit Kunnumkal, 1 Huseyin Topaloglu 2 1 Indian School of Business, Gachibowli, Hyderabad 500032, India 2 School
More informationThe Dual Simplex Algorithm
p. 1 The Dual Simplex Algorithm Primal optimal (dual feasible) and primal feasible (dual optimal) bases The dual simplex tableau, dual optimality and the dual pivot rules Classical applications of linear
More informationElectronic Companion to Initial Shipment Decisions for New Products at Zara
e-companion to Gallien et al: Initial Shipment Decisions for New Products at Zara ec1 Electronic Companion to Initial Shipment Decisions for New Products at Zara This document includes appendices to Initial
More informationWe consider a joint location-inventory problem involving a single supplier and multiple
Zuo-Jun Max Shen Collette Coullard Mark S. Daskin Department of Industrial & Systems Engineering, University of Florida, Gainesville, Florida 32611 Department of Industrial Engineering and Management Sciences,
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/39637 holds various files of this Leiden University dissertation Author: Smit, Laurens Title: Steady-state analysis of large scale systems : the successive
More informationCPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Motivating Quote for Queueing Models Good things come to those who wait - poet/writer
More informationSession-Based Queueing Systems
Session-Based Queueing Systems Modelling, Simulation, and Approximation Jeroen Horters Supervisor VU: Sandjai Bhulai Executive Summary Companies often offer services that require multiple steps on the
More informationA STAFFING ALGORITHM FOR CALL CENTERS WITH SKILL-BASED ROUTING: SUPPLEMENTARY MATERIAL
A STAFFING ALGORITHM FOR CALL CENTERS WITH SKILL-BASED ROUTING: SUPPLEMENTARY MATERIAL by Rodney B. Wallace IBM and The George Washington University rodney.wallace@us.ibm.com Ward Whitt Columbia University
More informationF 1 F 2 Daily Requirement Cost N N N
Chapter 5 DUALITY 5. The Dual Problems Every linear programming problem has associated with it another linear programming problem and that the two problems have such a close relationship that whenever
More informationAppendix A Taylor Approximations and Definite Matrices
Appendix A Taylor Approximations and Definite Matrices Taylor approximations provide an easy way to approximate a function as a polynomial, using the derivatives of the function. We know, from elementary
More informationOn Two Class-Constrained Versions of the Multiple Knapsack Problem
On Two Class-Constrained Versions of the Multiple Knapsack Problem Hadas Shachnai Tami Tamir Department of Computer Science The Technion, Haifa 32000, Israel Abstract We study two variants of the classic
More informationLinear Programming. H. R. Alvarez A., Ph. D. 1
Linear Programming H. R. Alvarez A., Ph. D. 1 Introduction It is a mathematical technique that allows the selection of the best course of action defining a program of feasible actions. The objective of
More informationM/M/1 Queueing systems with inventory
Queueing Syst 2006 54:55 78 DOI 101007/s11134-006-8710-5 M/M/1 Queueing systems with inventory Maike Schwarz Cornelia Sauer Hans Daduna Rafal Kulik Ryszard Szekli Received: 11 August 2004 / Revised: 6
More informationA Study on Performance Analysis of Queuing System with Multiple Heterogeneous Servers
UNIVERSITY OF OKLAHOMA GENERAL EXAM REPORT A Study on Performance Analysis of Queuing System with Multiple Heterogeneous Servers Prepared by HUSNU SANER NARMAN husnu@ou.edu based on the papers 1) F. S.
More informationFacility Location and Distribution System Planning. Thomas L. Magnanti
Facility Location and Distribution System Planning Thomas L. Magnanti Today s Agenda Why study facility location? Issues to be modeled Basic models Fixed charge problems Core uncapacitated and capacitated
More informationAn Uncertain Bilevel Newsboy Model with a Budget Constraint
Journal of Uncertain Systems Vol.12, No.2, pp.83-9, 218 Online at: www.jus.org.uk An Uncertain Bilevel Newsboy Model with a Budget Constraint Chunliu Zhu, Faquan Qi, Jinwu Gao School of Information, Renmin
More informationLecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem
Lecture 3: Lagrangian duality and algorithms for the Lagrangian dual problem Michael Patriksson 0-0 The Relaxation Theorem 1 Problem: find f := infimum f(x), x subject to x S, (1a) (1b) where f : R n R
More informationQueueing systems. Renato Lo Cigno. Simulation and Performance Evaluation Queueing systems - Renato Lo Cigno 1
Queueing systems Renato Lo Cigno Simulation and Performance Evaluation 2014-15 Queueing systems - Renato Lo Cigno 1 Queues A Birth-Death process is well modeled by a queue Indeed queues can be used to
More informationQueues and Queueing Networks
Queues and Queueing Networks Sanjay K. Bose Dept. of EEE, IITG Copyright 2015, Sanjay K. Bose 1 Introduction to Queueing Models and Queueing Analysis Copyright 2015, Sanjay K. Bose 2 Model of a Queue Arrivals
More informationThe Multi-Path Utility Maximization Problem
The Multi-Path Utility Maximization Problem Xiaojun Lin and Ness B. Shroff School of Electrical and Computer Engineering Purdue University, West Lafayette, IN 47906 {linx,shroff}@ecn.purdue.edu Abstract
More information- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs
LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs
More informationHow to Take the Dual of a Linear Program
How to Take the Dual of a Linear Program Sébastien Lahaie January 12, 2015 This is a revised version of notes I wrote several years ago on taking the dual of a linear program (LP), with some bug and typo
More informationTwo Heterogeneous Servers Queueing-Inventory System with Sharing Finite Buffer and a Flexible Server
Two Heterogeneous Servers Queueing-Inventory System with Sharing Finite Buffer and a Flexible Server S. Jehoashan Kingsly 1, S. Padmasekaran and K. Jeganathan 3 1 Department of Mathematics, Adhiyamaan
More informationCENTER FOR MULTIMODAL SOLUTIONS FOR CONGESTION MITIGATION (CMS)
Final Report to the CENTER FOR MULTIMODAL SOLUTIONS FOR CONGESTION MITIGATION (CMS) CMS Project Number: 2011-023 Project Title: The Impacts of Freight Mode Splitting on Congestion, Risk, and Delivery Reliability
More informationT. Liggett Mathematics 171 Final Exam June 8, 2011
T. Liggett Mathematics 171 Final Exam June 8, 2011 1. The continuous time renewal chain X t has state space S = {0, 1, 2,...} and transition rates (i.e., Q matrix) given by q(n, n 1) = δ n and q(0, n)
More informationAvailability. M(t) = 1 - e -mt
Availability Availability - A(t) the probability that the system is operating correctly and is available to perform its functions at the instant of time t More general concept than reliability: failure
More informationOptimisation. 3/10/2010 Tibor Illés Optimisation
Optimisation Lectures 3 & 4: Linear Programming Problem Formulation Different forms of problems, elements of the simplex algorithm and sensitivity analysis Lecturer: Tibor Illés tibor.illes@strath.ac.uk
More informationLecture 7: Simulation of Markov Processes. Pasi Lassila Department of Communications and Networking
Lecture 7: Simulation of Markov Processes Pasi Lassila Department of Communications and Networking Contents Markov processes theory recap Elementary queuing models for data networks Simulation of Markov
More informationImproving Supply Chain Performance: Real-Time Demand Information and Flexible Deliveries
Improving Supply Chain Performance: Real-Time Demand Information and Flexible Deliveries Kevin H. Shang Sean X. Zhou Geert-Jan van Houtum The Fuqua School of Business, Duke University, Durham, North Carolina
More informationDistributed Optimization. Song Chong EE, KAIST
Distributed Optimization Song Chong EE, KAIST songchong@kaist.edu Dynamic Programming for Path Planning A path-planning problem consists of a weighted directed graph with a set of n nodes N, directed links
More informationA Queueing System with Queue Length Dependent Service Times, with Applications to Cell Discarding in ATM Networks
A Queueing System with Queue Length Dependent Service Times, with Applications to Cell Discarding in ATM Networks by Doo Il Choi, Charles Knessl and Charles Tier University of Illinois at Chicago 85 South
More informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More informationOPTIMALITY OF RANDOMIZED TRUNK RESERVATION FOR A PROBLEM WITH MULTIPLE CONSTRAINTS
OPTIMALITY OF RANDOMIZED TRUNK RESERVATION FOR A PROBLEM WITH MULTIPLE CONSTRAINTS Xiaofei Fan-Orzechowski Department of Applied Mathematics and Statistics State University of New York at Stony Brook Stony
More informationA Market Mechanism to Assign Air Traffic Flow Management Slots
A Market Mechanism to Assign Air Traffic Flow Management Slots Andrea Ranieri, Lorenzo Castelli Università degli Studi di Trieste Dipartimento di Elettrotecnica, Elettronica ed Informatica 8 th USA/Europe
More informationThe Effect of Supply Disruptions on Supply Chain. Design Decisions
The Effect of Supply Disruptions on Supply Chain Design Decisions Lian Qi Department of Supply Chain Management & Marketing Sciences Rutgers Business School, Rutgers University, Newark, NJ Zuo-Jun Max
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationOne billion+ terminals in voice network alone
Traffic Engineering Traffic Engineering One billion+ terminals in voice networ alone Plus data, video, fax, finance, etc. Imagine all users want service simultaneously its not even nearly possible (despite
More informationOPTIMAL CONTROL OF A FLEXIBLE SERVER
Adv. Appl. Prob. 36, 139 170 (2004) Printed in Northern Ireland Applied Probability Trust 2004 OPTIMAL CONTROL OF A FLEXIBLE SERVER HYUN-SOO AHN, University of California, Berkeley IZAK DUENYAS, University
More informationLagrangian Duality. Richard Lusby. Department of Management Engineering Technical University of Denmark
Lagrangian Duality Richard Lusby Department of Management Engineering Technical University of Denmark Today s Topics (jg Lagrange Multipliers Lagrangian Relaxation Lagrangian Duality R Lusby (42111) Lagrangian
More informationPart A. Ch (a) Thus, order quantity is 39-12=27. (b) Now b=5. Thus, order quantity is 29-12=17
Homework2Solution Part A. Ch 2 12 (a) b = 65 40 = 25 from normal distribution table Thus, order quantity is 39-12=27 (b) Now b=5 from normal distribution table Thus, order quantity is 29-12=17 It is interesting
More information56:270 Final Exam - May
@ @ 56:270 Linear Programming @ @ Final Exam - May 4, 1989 @ @ @ @ @ @ @ @ @ @ @ @ @ @ Select any 7 of the 9 problems below: (1.) ANALYSIS OF MPSX OUTPUT: Please refer to the attached materials on the
More information56:171 Operations Research Final Exam December 12, 1994
56:171 Operations Research Final Exam December 12, 1994 Write your name on the first page, and initial the other pages. The response "NOTA " = "None of the above" Answer both parts A & B, and five sections
More informationMIT Manufacturing Systems Analysis Lectures 6 9: Flow Lines
2.852 Manufacturing Systems Analysis 1/165 Copyright 2010 c Stanley B. Gershwin. MIT 2.852 Manufacturing Systems Analysis Lectures 6 9: Flow Lines Models That Can Be Analyzed Exactly Stanley B. Gershwin
More informationA Shadow Simplex Method for Infinite Linear Programs
A Shadow Simplex Method for Infinite Linear Programs Archis Ghate The University of Washington Seattle, WA 98195 Dushyant Sharma The University of Michigan Ann Arbor, MI 48109 May 25, 2009 Robert L. Smith
More informationSystems Simulation Chapter 6: Queuing Models
Systems Simulation Chapter 6: Queuing Models Fatih Cavdur fatihcavdur@uludag.edu.tr April 2, 2014 Introduction Introduction Simulation is often used in the analysis of queuing models. A simple but typical
More informationOn the Partitioning of Servers in Queueing Systems during Rush Hour
On the Partitioning of Servers in Queueing Systems during Rush Hour Bin Hu Saif Benjaafar Department of Operations and Management Science, Ross School of Business, University of Michigan at Ann Arbor,
More informationmaxz = 3x 1 +4x 2 2x 1 +x 2 6 2x 1 +3x 2 9 x 1,x 2
ex-5.-5. Foundations of Operations Research Prof. E. Amaldi 5. Branch-and-Bound Given the integer linear program maxz = x +x x +x 6 x +x 9 x,x integer solve it via the Branch-and-Bound method (solving
More informationA Decomposition Approach for a Class of Capacitated Serial Systems 1
A Decomposition Approach for a Class of Capacitated Serial Systems 1 Ganesh Janakiraman 2 IOMS-OM Group Stern School of Business New York University 44 W. 4th Street, Room 8-71 New York, NY 10012-1126.
More informationQUEUING MODELS AND MARKOV PROCESSES
QUEUING MODELS AND MARKOV ROCESSES Queues form when customer demand for a service cannot be met immediately. They occur because of fluctuations in demand levels so that models of queuing are intrinsically
More information57:022 Principles of Design II Final Exam Solutions - Spring 1997
57:022 Principles of Design II Final Exam Solutions - Spring 1997 Part: I II III IV V VI Total Possible Pts: 52 10 12 16 13 12 115 PART ONE Indicate "+" if True and "o" if False: + a. If a component's
More information0utline. 1. Tools from Operations Research. 2. Applications
0utline 1. Tools from Operations Research Little s Law (average values) Unreliable Machine(s) (operation dependent) Buffers (zero buffers & infinite buffers) M/M/1 Queue (effects of variation) 2. Applications
More informationA Customer-Item Decomposition Approach to Stochastic Inventory Systems with Correlation
A Customer-Item Decomposition Approach to Stochastic Inventory Systems with Correlation Yimin Yu Saif Benjaafar Program in Industrial and Systems Engineering University of Minnesota, Minneapolis, MN 55455
More informationCHAPTER-3 MULTI-OBJECTIVE SUPPLY CHAIN NETWORK PROBLEM
CHAPTER-3 MULTI-OBJECTIVE SUPPLY CHAIN NETWORK PROBLEM 3.1 Introduction A supply chain consists of parties involved, directly or indirectly, in fulfilling customer s request. The supply chain includes
More informationMethodology for Computer Science Research Lecture 4: Mathematical Modeling
Methodology for Computer Science Research Andrey Lukyanenko Department of Computer Science and Engineering Aalto University, School of Science and Technology andrey.lukyanenko@tkk.fi Definitions and Goals
More information56:171 Operations Research Final Examination December 15, 1998
56:171 Operations Research Final Examination December 15, 1998 Write your name on the first page, and initial the other pages. Answer both Parts A and B, and 4 (out of 5) problems from Part C. Possible
More informationSlides 9: Queuing Models
Slides 9: Queuing Models Purpose Simulation is often used in the analysis of queuing models. A simple but typical queuing model is: Queuing models provide the analyst with a powerful tool for designing
More informationSection Notes 9. IP: Cutting Planes. Applied Math 121. Week of April 12, 2010
Section Notes 9 IP: Cutting Planes Applied Math 121 Week of April 12, 2010 Goals for the week understand what a strong formulations is. be familiar with the cutting planes algorithm and the types of cuts
More informationExercises Solutions. Automation IEA, LTH. Chapter 2 Manufacturing and process systems. Chapter 5 Discrete manufacturing problems
Exercises Solutions Note, that we have not formulated the answers for all the review questions. You will find the answers for many questions by reading and reflecting about the text in the book. Chapter
More informationNo-Holdback Allocation Rules for Continuous-Time Assemble-to-Order Systems
OPERATIONS RESEARCH Vol. 58, No. 3, May June 2010, pp. 691 705 issn 0030-364X eissn 1526-5463 10 5803 0691 informs doi 10.1287/opre.1090.0785 2010 INFORMS No-Holdback Allocation Rules for Continuous-Time
More informationSection 3.3: Discrete-Event Simulation Examples
Section 33: Discrete-Event Simulation Examples Discrete-Event Simulation: A First Course c 2006 Pearson Ed, Inc 0-13-142917-5 Discrete-Event Simulation: A First Course Section 33: Discrete-Event Simulation
More informationIntroduction to queuing theory
Introduction to queuing theory Claude Rigault ENST claude.rigault@enst.fr Introduction to Queuing theory 1 Outline The problem The number of clients in a system The client process Delay processes Loss
More informationQueuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem. Wade Trappe
Queuing Networks: Burke s Theorem, Kleinrock s Approximation, and Jackson s Theorem Wade Trappe Lecture Overview Network of Queues Introduction Queues in Tandem roduct Form Solutions Burke s Theorem What
More informationBIRTH DEATH PROCESSES AND QUEUEING SYSTEMS
BIRTH DEATH PROCESSES AND QUEUEING SYSTEMS Andrea Bobbio Anno Accademico 999-2000 Queueing Systems 2 Notation for Queueing Systems /λ mean time between arrivals S = /µ ρ = λ/µ N mean service time traffic
More informationInventory Ordering Control for a Retrial Service Facility System Semi- MDP
International Journal of Engineering Science Invention (IJESI) ISS (Online): 239 6734, ISS (Print): 239 6726 Volume 7 Issue 6 Ver I June 208 PP 4-20 Inventory Ordering Control for a Retrial Service Facility
More information5 Handling Constraints
5 Handling Constraints Engineering design optimization problems are very rarely unconstrained. Moreover, the constraints that appear in these problems are typically nonlinear. This motivates our interest
More informationThe Kuhn-Tucker Problem
Natalia Lazzati Mathematics for Economics (Part I) Note 8: Nonlinear Programming - The Kuhn-Tucker Problem Note 8 is based on de la Fuente (2000, Ch. 7) and Simon and Blume (1994, Ch. 18 and 19). The Kuhn-Tucker
More informationII BSc(Information Technology)-[ ] Semester-III Allied:Computer Based Optimization Techniques-312C Multiple Choice Questions.
Dr.G.R.Damodaran College of Science (Autonomous, affiliated to the Bharathiar University, recognized by the UGC)Re-accredited at the 'A' Grade Level by the NAAC and ISO 9001:2008 Certified CRISL rated
More information