On the Accuracy of Fluid Approximations to a Class of Inventory-level-dependent Eoq and Epq Models

Deterministic Economic Order Quantity EOQ models have been studied intensively in the literature, where the demand process is described by an ordinary differential equation, and the objective is to obtain an EOQ, which minimizes the total cost per unit time. The total cost per unit time consists of a " discrete " part, the setup cost, which is incurred at the time of ordering, and a " continuous " part, the holding cost, which is continuously accumulated over time. Quite formally, such deterministic EOQ models can be viewed as fluid approximations to the corresponding stochastic EOQ models, where the demand process is taken as a stochastic jump process. Suppose now an EOQ is obtained from a deterministic model. The question is how well does this quantity work in the corresponding stochastic model. In the present paper we justify a translation of EOQs obtained from deterministic models, under which the resulting order quantities are asymptotically optimal for the stochastic models, by showing that the difference between the performance measures and the optimal values converges to zero with respect to a scaling parameter. Moreover, we provide an estimate for the rate of convergence. The same issue regarding specific Economic Production Quantity EPQ models is studied, too.


Introduction
Consider an inventory item which is demanded.So the inventory level gradually decreases and is backed up by ordering new inventories from time to time.There are two costs to consider: a positive inventory level results in a holding cost, and every order induces a setup cost.The objective is to determine an order quantity that minimizes the total cost per unit time tcu .Such a minimizer is known as an Economic Order Quantity EOQ , and the model itself is known as an EOQ model.

Advances in Operations Research
Arguably the simplest EOQ model sometimes referred to as the classic EOQ model is based on the following assumptions: i the instantaneous holding cost rate is constant; ii the setup cost is constant; iii the demand comes in a deterministic and continuous process at a constant rate; iv no backlogging is allowed so that at inventory level zero all the arrived demand is rejected; v the inventory item is homogeneous and nonperishable so that only the demand reduces the inventory level; vi the inventory level is reviewed continuously so that it can be described by an ordinary differential equation; vii the replenishment takes place instantaneously after ordering.A rigorous description of the classic EOQ model is shortly given in Section 2.1.Amongst the efforts of generalizing the classic EOQ model, a great deal have been made on relaxing assumption iii .For instance, the EOQ model considered in 1 see also the references therein assumes the demand rate to be inventory level dependent, despite the dependence is of a specific form.This is a response to: "At times, the presence of inventory has a motivational effect on the people around it.It is a common belief that large piles of goods displayed in a supermarket will lead the customer to buy more 2 ."Generalizations regarding other assumptions include allowing backlogging and periodic reviews in 3 , accounting for perishable goods in 4-7 , and so on.A comprehensive review of the literature on EOQ models is available in 8 .Note that the aforementioned works mainly focus on EOQ models where the demand process is deterministic.Quite formally these deterministic models can be viewed as the fluid approximations to the corresponding stochastic EOQ models, where the demand comes in a stochastic jump process.For this reason, in what follows we call the deterministic models also fluid models, and EOQs derived from fluid models are referred to as fluid EOQs.Suppose now that a fluid EOQ is obtained and the corresponding stochastic EOQ model is appropriately formulated.Then the issue of interest is how to translate the fluid EOQ into an order quantity for the stochastic model, where the expected total cost per unit time TCU is nearly minimal.The fluid approximation would be justified if such a translation is obtained.The formal justification of the fluid approximations to various jump Markov optimization models with local transitions has been addressed by numerous authors, see 9-12 , all of which focus on queueing networks.In greater detail the optimization problems considered in 9, 12 are with a discounted criterion over an infinite horizon, the one in 11 is with an expected total cost criterion over a fixed finite time horizon, while the performance measure for the queueing network in 10 is the expected total cost up to the first moment the system gets empty.Here we emphasize that in fluid resp., stochastic EOQ models the performance measures are tcus resp., TCUs , which are long run averages.All of those works compare the optimal value for the fluid model with the performance measure for the scaled stochastic model under the translated policy and show that the difference converges to zero as the scaling parameter increases to infinity.However, none of them reveal the rate of convergence, which measures the accuracy of the fluid approximations and the efficiency of the underlying translations.Consequently the more recent development on this topic aims at obtaining the rate of convergence, see 13-15 .In particular, in 14 the author proposes a translation of the fluid EOQ and shows that it results in an asymptotically fluid optimal AFO and asymptotically optimal AO order quantity for the corresponding stochastic model.The accurate definitions of "AFO" and "AO" are postponed to Sections 2 and 3.However, fairly strong conditions on the system parameters are assumed there, restricting the applicability of the obtained results.
Therefore, the main contribution of this paper is to provide a refinement of the results obtained in 14 by relaxing the conditions assumed therein, and thus enlarging the applicability to cover a broader class of EOQ models.In greater detail, by taking the stochastic model as a continuous-time Markov chain, we justify the translation of fluid EOQ proposed in 14 and obtain its efficiency in the form of a rate of convergence.The Markov property of our stochastic model is a result of assuming exponentially distributed interarrival times in the demand process, which is standard in the current literature on inventory systems, see 7, 16-19 .In our models, the demand and holding cost rates are of a rather general inventory-level dependence, and thus with broad applicabilities.In particular, our results are applicable to the important case of discontinuous demand and holding cost rates compared to that globally Lipschitz continuous rates are assumed in 14 , see more discussions on this in Section 4.Moreover, in this paper results similar to those for EOQ models are derived for Economic Production Quantity EPQ models, too.
The rest of this paper is organized as follows.In Sections 2 and 3 we formulate EOQ and EPQ models and state the main results.In Section 4 some comments are given on the issues of possible applicabilities of the obtained results and thus illustrate the contribution of this paper.We finish this paper with conclusions.The proofs of the main statements are postponed to the appendix.

Economic Order Quantity Models
In what follows, the trivial case of an order quantity taking zero is excluded from consideration, and the context should always make it clear when • stands for the function taking the largest integer part of its argument.

Description of Mathematical Models (a) Fluid Model
Suppose some order quantity y > 0 is fixed.Then let {x t , t ≥ 0} be the inventory level process subject to the dynamics dx/dt −μ x when 0 < x t ≤ y and x t order 0 y, where x t reaches state zero at t order .Here, μ x > 0 is the demand rate, and the impulsive jump of x t at t order reflects the instantaneous replenishment assumption.Let t cycle be the time duration between two consecutive jumps of the inventory level process, g x the instantaneous holding cost rate, and K > 0 the setup cost, incurred immediately whenever an order is made.We are interested in minimizing the long run average cost given by tcu y lim Let us call y * the EOQ for the fluid model, so that tcu y * inf y>0 tcu y .In particular, if μ x μ > 0 and g x hx i.e., the demand rate and the holding cost rate per unit of inventory are constants μ > 0, h > 0 , then we have the so called classic EOQ model.

(b) (Scaled) Stochastic Model
In the corresponding stochastic model, the inventory item is measured in small units, so that a scaling parameter n 1, 2, . . . is present to indicate the units.The intuitive meaning of this scaling parameter is explained in Remark 2.1 below.In greater detail, fixing some integer Advances in Operations Research order quantity n Y > 0, the inventory level process { n X t , t ≥ 0} is modelled as a continuoustime Markov chain with the state space {0, 1, 2, . . ., n Y } and the transition rates where nμ i/n is the instantaneous demand rate and nκ is the parameter of the exponentially distributed lead time between the ordering and the corresponding replenishment.In other words, the time between two consecutive demand arrivals is exponentially distributed with mean 1/nμ i/n when the current inventory level is i, and new inventories are ordered when the inventory level hits zero.
Remark 2.1.The scaled stochastic model can be linked to the fluid model by taking n Y ny , where y > 0 is the order quantity for the fluid model Here ny > 0 for big enugh n.So below we assume that ny > 0, that is, we consider large enough n. .The above-described scaling is often referred to as a fluid scaling, and its intuitive meaning can be understood as follows.Clearly, when n 1, the stochastic model is a corresponding version of the fluid model.As we increase n, the demand comes in smaller units, and inventories are measured more accurately.Take μ i/n M as an example.Suppose that n 1 corresponds to the unit of kg so that on average 1 • μ i/1 M kg units of demand come per time unit.In the case of n 1000, the unit will be g : on average nμ i/n 1000M g units of demand come per time unit.Meanwhile, the cost rate is not amplified by multiplying n because it costs the same to hold either 1000M g or M kg of inventories.
We are interested in minimizing the performance measure given by In what follows, the EOQs for both fluid and stochastic models are assumed to be unique, and the similar assumption applies to the EPQ models.

Main Results
Condition 1. a There exist constants d 1 ≥ 0, k 1 > 0 and δ > 0 such that δ ≤ μ x ≤ k 1 and |g x | ≤ d 1 ; here functions μ • and g • are measurable and both defined on 0, ∞ .b There exist finite intervals 0, x 1 , x 1 , x 2 , . . .with lim j → ∞ x j ∞ such that on each of them g x and μ x are Lipschitz continuous with a common Lipschitz constant d g and d μ , respectively.
Note that Condition 1 implies that for any fixed y > 0, there exist L possibly ydependent finite intervals 0, x 1 , x 1 , x 2 , . . ., x L , 3 y 1 such that on each interval 1/μ x and g x /μ x are Lipschitz continuous with Lipschitz constants d μ /δ 2 and d 1 d μ k 1 d g /δ 2 , respectively.For simplicity, we define

2.7
In particular, Proposition 2.2 implies that lim n → ∞ | n TCU ny * − tcu y * | 0, that is, ny * is AFO.The same calculations as in 14, page 406 result in the next lemma.
Lemma 2.3.Suppose that Condition 1 is satisfied and g x ≥ 0, then 3), then one has Corollary 2.5.Under Condition 1, suppose in addition that L is y-independent and g x ≥ 0 on 0, ∞ .Then the following statements hold.

2.9
Here the y-independent E n goes to zero as fast as 1/n in the sense of

2.10
Here we recall that n Y * and y * are the EOQ for the (scaled) stochastic model and fluid model, respectively.
Corollary 2.5 refines Corollary 2.4 for certain cases by providing an estimate for the rate of convergence.
Fix some order quantity for the fluid model y and scaling parameter n, and let n Π i , i 0, 1, . . ., ny be the stationary distribution of the inventory process { n X t , t ≥ 0} and π x the invariant density in the fluid model of the underlying dynamics of x t .Then the following proposition shows that the fluid model can also be used to provide approximations to stationary distributions of the inventory level process in the scaled stochastic model.Proposition 2.6.Suppose that Condition 1 is satisfied and some order quantity for the fluid model y > 0 is fixed.Then one has and for i 1, 2, . . ., ny

A Comparison with Section 4 of [14]
The fluid approximations of the EOQ models are also briefly considered in 14, Section 4 .Therefore, we mention the main difference of this paper from that one in this subsection.Firstly, the present paper is based on a weaker condition.Indeed, instead of Condition 1, the following stronger condition is assumed in 14 .The global Lipschitz property is essential to the corresponding proof in 14 .
Secondly, the approach in 14 is based on the explicit expression of n TCU n Y obtained by solving the associated Poisson equations.Instead of doing that, the present paper employs the recent results of 15 .One advantage of this approach lies in the weaker condition required only piecewise Lipschitz continuity is needed .Another advantage is that it allows one to study EPQ models, because the Poisson equations for EPQ models are much more difficult to solve compared to those for EOQ models.That is why EPQ models are not considered in 14 .

Economic Production Quantity Models
In EOQ models the inventory is backed up at once by ordering new inventory items from external suppliers.In this section we consider the situation where the inventory is gradually backed up by producing new items.In greater detail, the inventory level decreases gradually to meet the demand, and when it hits zero, the production is switched on and new inventory items are being produced to back up the inventory.The production is switched off as soon as the inventory is backed up to a predetermined level.Here we have to account for the cost incurred from switching on the production as well as from holding the inventory items.The aim is to obtain an Economic Inventory Backup Level EIBL that minimizes TCU, and the resulting model is called an EPQ model.Similar to the previous section, below we justify the fluid approximations to stochastic EPQ models, whose rigorous description is shortly given.

Description of Mathematical Models (a) Fluid Model
Suppose that we fix some real inventory backup level y > 0, meaning that the production is always on until the inventory reaches the level y.Let {x t , t ≥ 0} represent the inventory level process in the fluid model with state space 0, y , and instantaneous demand and production rates μ x > 0 and λ x > 0, respectively.Then the inventory level process is subject to the dynamics dx dt where the production-off phase and production-on phase, superseding each other, are triggered by x t y and x t 0, respectively.In words, without any delay, once the inventory level reaches zero, production is switched on till it reaches the inventory backup level y.Let g x be the holding cost rate, K > 0 the setup cost incurred with switching on the production, and t cycle the time duration between two consecutive production switching-offs.So tcu y lim Let us denote by y * the EIBL for the fluid model, so that tcu y * inf y>0 tcu y .

(b) (Scaled) Stochastic Model
Suppose that we fix some inventory backup level positive integer n Y , meaning that the production is always on until the inventory reaches level n Y .Let { n X t , t ≥ 0} represent the inventory level process.We model it as a continuous-time Markov chain with the state space { n Y, off , n Y − 1, off , . . ., 0, off , 0, on , . . ., n Y − 1, on }, where i, off indicates that the inventory level is i and the production is off and the denotation of i, on can be understood in the same way.Its transition rates are given by for all i 1, . . ., n Y : and finally for all i 1, . . ., n Y −2: where nλ i/n and nμ i/n stand for the instantaneous production and demand rates, nκ is the parameter of the exponentially distributed lead time between the switching and the Advances in Operations Research 9 actual production-on, and we have ignored all the cases when the transition rates take zero.So we have
Note that Condition 3 implies that for any fixed y > 0, there exists an integer L possibly y-dependent and L 1 intervals 0, x 1 , x 1 , x 2 , . . ., x L , 3 y 1 such that on each interval 1/μ x , g x /μ x , and where one puts δ max{δ, δ μκ } and δ min{δ, δ μκ } and the n-independent terms B 1 , B 2 are given by A.2 in the appendix.
As in the case of EOQ models, we observe from Proposition 3.1 that ny * is AFO.

Advances in Operations
c If L is y-independent and g x ≥ 0, then the following two substatements hold.

3.8
Here F n is y-independent, and

3.9
One has c2

Example and Comments
In this section we firstly verify our results by considering a specific EPQ model, where tcu • and n TCU • can be analytically computed.Then we comment on the applications of our results.
Example 4.1.As for EOQ models, one may refer to 14 for an example.Hence we study the following classic setting for the scaled stochastic EPQ model: assume constant demand and production rates μ x D > 0 and λ x R > 0, linear holding cost g x hx with a constant h > 0, constant setup cost K > 0, and finally no lead time between "switching" and "actual production-on," corresponding to if we take κ → ∞.Therefore, if we consider the underlying continuous time Markov chain { n X t , t ≥ 0}, state 0, off will be excluded.The transition rates are modified accordingly.Clearly, Condition 3 is satisfied with this classic setting.The following lemma gives the explicit formula for n TCU Z , where for simplicity we have put Z instead of n Y for the inventory backup level.

Proposition 4.2. For the stochastic EPQ model described above,
The proof of this proposition is based on solving quite tediously the associated Poisson equation for n TCU Z and is omitted here.
The  1.Note that when n 100, the curves of tcu y and n TCU ny nearly coincide.Secondly, if n increases, then by inspecting the numerator especially the first two terms and the denominator of the expression for n TCU Z , we see that if Z does not increase as fast as n, n TCU Z will blow up to ∞ it can be easily checked that the expression

Comments
Let us comment on the applicability issues of our main results Propositions 2.2 and 3.1 .We mainly focus on EOQ models, as absolutely similar comments can be made on EPQ models in the same manner.
Although we assume the ordering point to be always zero, our results are still applicable when it is set to be some fixed positive level, because Lemma A.1 holds if we put another absorbing state instead of zero.In particular, if one allows the state taking negative values, by putting some negative state absorbing, our results also impound the case of backlogging.This flexibility regarding the ordering point together with the fact that g • is unrestricted in signs enriches the applications of our results, in that although we require μ to be separated from zero, when profit rather than solely operational cost is counted, the ordering point is most likely positive, meaning in cases of μ x αx β , α > 0, 0 < β < 1 as in 20 and μ x αx −β , α > 0, β ≥ 1 as in 1 , μ will be essentially separated from zero, validating our results.
The state-dependence given in Conditions 1 and 3 is fairly general.In particular, that functions λ x , μ x and g x being bounded is not restrictive, because once some EOQ for the fluid model y > 0 is fixed, to validate Propositions 2.2 and 3.1, they are only required to be bounded on bounded intervals.Note that in addition to the demand rate, some authors such as those of 5, 6 also include a state-dependent deteriorating rate, to indicate that the underlying goods are perishable.Our results are also applicable to such cases: one only needs explain μ • as the total reduction rate of the inventory level.
Finally, Propositions 2.2 and 3.1 are significant extensions of the relevant results in 14 , where the author only focuses on EOQ models and requires global Lipschitz continuity of μ • and g • .However, from the modelling point of view, the case of discontinuous functions is interesting and important as demonstrated by 1, 3, 5, 21, 22 , where 1 considers a piecewise constant function μ • and the others consider discontinuous μ • taking either a constant value or according to μ x αx β , α > 0, 0 < β < 1.The results in 14 are derived based on the closed form of the solution to a Poisson equation satisfied by n TCU • , which is tremendously difficult to get explicitly in the case of stochastic EPQ models.

Conclusions
To sum up, in this work we formally justified a general class of inventory level-dependent deterministic EOQ and EPQ models, regarded as the fluid approximations to their stochastic versions, by showing a translation of the fluid EOQ EIBL to provide an order quantity inventory backup level asymptotically achieving some optimality for the stochastic model.The efficiency of the translation was obtained, as distinguished from the majority of the works on fluid approximations.The class of inventory models are quite broad so that to various extent, the obtained results are directly applicable to the existing works such as 1, 5, 6, 8, 20-22 .The present work is a significant extension of the relevant results in 14 .

Appendix
To aid our proof, firstly, let us consider the following one-dimensional birth-and-death process { n Z t , t ≥ 0} with state space {0, 1, . ..} and birth and death rates nα i/n and nβ i/n , respectively, where nonnegative measurable functions α and β are defined on 0, ∞ and i indicates the current state of the process.In addition, α 0 β 0 0, where the equality holds only at 0, meaning that state zero is absorbing.Let E i denote the expectation of any underlying functional of the process with the initial state n Z 0 i.Let a real measurable function γ • defined on 0, ∞ be fixed with γ 0 0. Now we are in the position to state the following condition.
b There exist finite intervals z 0 0, z 1 , z 1 , z 2 , . . .with lim j → ∞ z j ∞ such that on each of them, γ z / β z − α z is a Lipschitz continuous function with a common Lipschitz constant d 2 .
Note that Condition 1 b implies that for any fixed y > 0 there exists an integer L possibly z-dependent and L 1 finite intervals 0, z 1 , z 1 , z 2 , . . ., z L , 3 y 1 such that on each interval, function γ • / β − α is Lipschitz continuous with a common Lipschitz constant d 2 .
The following lemma is a slightly stronger version of 15, Theorem 2 and will play an important role in our proof.
where regarding the second integral the underlying dynamics is given by dz/ds α z − β z ,

Advances in Operations Research
z 0 i/n and B 1 and B 2 are given by A.2 Proof.It can be easily checked in the proof of 14, Theorem 2 that our Condition 1, weaker than the original conditions imposed therein, is sufficient for the statement.See also 23 .

Proofs of Proposition 2.2, Corollaries 2.4 and 2.5, and Proposition 2.6
For both the fluid model and scaled stochastic model let us call the time duration between two consecutive replenishments a cycle and denote them by t cycle and n T cycle , respectively.
Here for simplicity, we do not explicitly indicate the y-dependence resp., n Y -dependence of t cycle resp.,n T cycle .Clearly { n X t , t ≥ 0} is a regenerative process 24, page 425 in that it probabilistically repeats itself from one cycle to the next.It then follows from 25, Theorem 1.1, Proposition 131 see also 24, Proposition 7.3 that as far as the long-run average n TCU ny is concerned, it suffices to consider the inventory level process and the cost incurred with it over only one cycle.For simplicity, we always consider the cycle starting at time t 0 with the initial position n X 0 ny .Let us denote by n TC and tc the total cost incurred over the cycle in the stochastic and fluid model, respectively.Then the following lemma indicates that the difference between E ny n T cycle and t cycle and the one between E ny n TC and tc cannot be too big.
Lemma A.2.Under Condition 1 the following two inequalities hold with nonnegative B 1 and B 2 given by A.2 : Proof. a Let us denote by n T absorbing the time duration from the starting point t 0 up to the point when n X t firstly reaches state zero.Obviously we have Then E ny n T cycle E ny n T absorbing 1/nκ, where the second term on the right hand side is the expected lead time.Now observe firstly that Condition 1 is a specific version of Condition 1: one can take η > 1 to be arbitrary, and put functions Here it does no matter to put μ 0 0. α x ≡ 0, β x μ x and I{x > 0} g x γ x ; and secondly that the inventory level process from t 0 up to n T absorbing is a pure death process.Therefore, one can refer to Lemma A.
b Let us denote by n TC absorbing the cost incurred during the interval 0, n T absorbing , so that E ny n TC E ny n TC absorbing g 0 /nκ, where the second term on the right hand side corresponds to the cost incurred over the lead time.In the same way as in part a , comparing E ny n TC absorbing with tc first, and then adding g 0 /nκ results in the statement.Remember that the setup cost cancels out.

Proof of
where the last inequality follows from the facts that y/k 1 ≤ t cycle ≤ y/δ, tc ≤ d 1 y/δ, E ny n T cycle ≥ ny /nk 1 1/nκ and Lemma A.2. Now let us easily observe that η 1 / η−1 and η/ η−1 both decrease with η ∈ 1, ∞ .It follows that B 1 /n, B 2 η −2n y 1 , and thus the above-derived expression all decrease with η, where we recall that η can be an arbitrary number on the interval 1, ∞ , see Condition 1.This implies that A.5 Proof of Corollary 2.4.For any fixed n, let us denote n Y * n y n .We do the proof in two parts.
Part 1.We consider the case of a convergent sequence y n .Suppose now that as n → ∞, y n does not go to y * but lim n → ∞ y n y > 0; here we allow y to be from the extended real line.In particular, for big enough n, y n is separated from zero.According to Proposition 2. Part 2. Now consider the case of a divergent sequence y n .One only needs consider the following two situations: either it has a bounded subsequence, which by Bolzano-Weierstrass theorem further has a convergent subsequence; or it does not have a bounded subsequence, which means that it has a subsequence blowing up to ∞.However, by taking the corresponding subsequences, we find that both situations have been essentially covered in Part 1. Part 2 is thus proved.
Proof of Corollary 2.5.a Under the conditions of the statement we have RHS of 2.7 Recall, here ny − 1 > 0 A.6 Here we use the fact that y/ ny − 1 decreases with y and y ≥ δ 2 K/k A.9 Here we put π 0 0 for convenience.Then A.10

Advances in Operations Research
Here we recall a of Lemma A.2. Recall that in the above derived expression, η can be any number from 1, ∞ .After passing to the limit η → ∞, we eventually end up with as required.

Proofs of Proposition 3.1 and Corollary 3.2
Let us call a cycle the time duration between two consecutive moments when the inventory is fully backed up.Arguing similarly as for EOQ models, it suffices to consider the inventory level process { n X t , t ≥ 0} and the cost incurred over one complete cycle, for which we put the starting time of t 0. Let us denote by t cycle , n T cycle and tc, n TC the duration of a cycle and the cost incurred over a cycle in the fluid and scaled stochastic model, respectively.Notice additionally that a cycle is always constituted to by two phases corresponding to the on and off of the production.This raises another set of denotations: let t on , n T on t off , n T off , and tc on , n TC on tc off , n TC off be the total cost incurred during the production-on off phases in the fluid and scaled stochastic model, respectively.We agree on that in both fluid and scaled stochastic model the setup cost is accounted for in tc off and n TC off .
Proof. a Let us concentrate on the inventory level process over the production-on phase.
In the fluid model, it appears convenient to reflect the trajectory {x t , t ∈ t on , t cycle } corresponding to the solid curve in Figure 2 about the horizontal t-axis first, and then shift the resulting trajectory corresponding to the curve of crosses in Figure 2 upwards by y units, and finally further shift the resulting trajectory to the left by shifting the time by t off units to the right to get { on x t , t ∈ 0, t on } corresponding to the curve of solid boxes in Figure 2 .Note now, for { on x t , t ∈ 0, t on } with on x 0 y the roles of production and demand have switched over: each produced unit reduces on x by one unit, and each demanded unit increases on x by one unit.More precisely, let us define the following functions: on μ y λ 0 , on μ 0 0, on μ x λ y − x , x ∈ 0, y , on λ y μ 0 0, on λ 0 0, on λ μ y − x , x ∈ 0, y , on g y g 0 , on g 0 0, on g x g y − x , x ∈ 0, y .

A.13
Then the dynamics of on x 0 y with d on x t /dt| on x t x on λ x − on μ x for x ∈ 0, y is of our interest, because we can write tc on ∞ 0 g on x t dt.Absolutely similar arguments are applicable to the scaled stochastic model.Consequently, we can consider the inventory level process during a production-on phase as a birth-and-death process {

A.15
Proof of Corollary 3.2.a Suppose that the statement does not hold.That is, for some subsequence {n j , j 1, 2, . ..} with n j → ∞ as j → ∞, n j Y * o n j in that lim j → ∞ n j Y * /n j 0. Under n j Y * we have where by 26, Theorem 1, page 175 the term inside the first curry bracket corresponds to E n j T off and the second last sum corresponds to E n j T on .Here we agree on that when i 0, the term in the second curry bracket reduces to 1/n j λ 0 .This gives

Condition 2 .
a Condition 1 a holds.b Functions g • and μ • are globally Lipschitz continuous.

Figure 2 :
Figure 2: The illustrative graph of on x t .

k
0 μ i/n j μ i − 1 /n j • • • μ i − k /n j /n j λ i/n j λ i − 1 /n j • • • λ i − k − 1 /n j ,where the last inequality follows from the fact g x ≥ 0. Clearly, the right hand side expression of the above inequality goes to infinity as n j → ∞, because λ x and μ x are dt , 3.6 with the holding cost g i/n and setup cost K > 0. Let us denote by n Y * the EIBL for the scaled stochastic model, so that n TCU n Y * infn Y 1,2,... n TCU n Y .The concept of AFO and AO inventory backup level ny can be understood in the same manner as introduced at the end of Section 2.1.
corresponding deterministic EPQ model can be solved easily, and we have which agree with Proposition 3.1.For the means of illustration, we put K 5, h 1, D 5, R 10. Then the graphs of tcu • and n TCU • are plotted in Figure Proposition 2.2.Under Condition 1 we have is continuous in y.This gives n TCU n y n → tcu y > tcu y * .However, it follows from Proposition 2.2 that n TCU ny * → tcu y * .This indicates that at least for big enough n, n TCU ny * < n TCU n y n n TCU n Y * , which is a desired contradiction.Hence lim n → ∞ y n y * , and consequently, lim n → ∞ | n TCU ny * − n TCU n Y * | 0, as required.
1 d 1 − δ/nκ.a is now clear.b According to Lemma 2.3, y * ≥ δ 2 K/k 1 d 1 , and n Y * satisfies n Y * /n ≥ δ 2 K/k 1 d 1 − δ/nκ.Therefore, according to part a , for n ≥ N we have Proof of Proposition 2.6.Now n Π i and π x can be easily computed as done in 14 .So we have n TCU ny * ≤ tcu y * E n ≤ tcu n Y * n E n ≤ n TCU n Y * n TCU ny * ≥ n TCU n Y * ≥ tcu n Y * cycle , 0 < x ≤ ny .