Technical Note on Q , r , L Inventory Model with Defective Items

and Applied Analysis 3 Wu and Ouyang 1 stated that the optimal solution can be obtained by the iterative method. Nevertheless, we point out that two sequences, generated by Wu and Ouyang 1 , are not guaranteed to converge. On the other hand, even when the sequence converges, the reason why the limit is the optimal solution is not discussed or verified. In Section 3, we will develop an approach to show that there is an optimal order quantity under a reasonable assumption, based on which the feasible domain with an upper bound and a lower bound can be derived. 3. Our Revision We simplify the expression in 2.2 and 2.3 as follows: Q2 α1 α2 (√ 1 k2 − k ) , 3.1 √ 1 k2 − k √ 1 k2 2α3Q ( 1 − βα3Q α4 , 3.2 where α1 2D/hδ A ci Li−1 − L ∑i−1 j 1 cj bj − aj , α2 D/hδ π π0 1 − β σ √ L, α3 h 1 − E p , and α4 D π π0 1 − β . Owing to 3.2 , we compute 1 − √ 1 k2 − k √ 1 k2 1 − 2α3Q ( 1 − βα3Q α4 , 3.3


Introduction
Wu and Ouyang 1 considered the imperfect production of the supplier and/or damage in delivery so that an arrival order lot may contain defective items.They extended the inventory model presented by Paknejad et al. 2 with constant lead time and fixed defective rate in an order lot to stochastic inventory models with crashable lead time and a random number of defective items.Wu and Ouyang 1 assumed that the purchasers inspect all the items they have ordered.Inspection is proposed to be nondestructive and error-free.All defective items are detected and will be returned to the vendor at the time of delivery of the next lot.The inventory model is continuously reviewed and an order of size, Q, is made whenever the inventory level falls to the reorder point r.In essence, Wu and Ouyang 1 applied the minimax approach for the stochastic inventory models with distribution-free demand to derive a mixture inventory model with back orders and lost sales in which an arrival order lot may contain defective items and the number of defective items is a random variable.The decision variables include the order quantity, the reorder point, and the lead time.Wu and Ouyang 1 developed two inventory models: for the first model the lead time demand follows a normal distribution, and for the second model, the distribution of the lead time demand is unknown except for the finite first and second moments.Because the information about the form of the probability distribution of lead time demand is often limited in practice, we consider in this paper only the second model of Wu and Ouyang 1 in which the minimax approach proposed by Moon and Gallego 3 is applied.Wu and Ouyang 1 claimed that they proposed an algorithm procedure to obtain the optimal ordering strategy.However, we have found and shown that the iterative method they employed in the attempt to minimize the average cost may not be able to locate or guarantee the optimal solution.For rectification and improvement, the aim of this study is to develop an analytical approach to find an upper bound and a lower bound for the order quantity.Our approach will ensure the uniqueness for the optimal solution.The same problem as Wu and Ouyang 1 is examined to demonstrate our proposed approach.

Review of Previous Results
To be compatible with the results of Wu and Ouyang 1 , we use the same notation and assumptions as them.For the distribution-free model, we directly quote the objective function of Wu and Ouyang 1 where Abstract and Applied Analysis 3 Wu and Ouyang 1 stated that the optimal solution can be obtained by the iterative method.Nevertheless, we point out that two sequences, generated by Wu and Ouyang 1 , are not guaranteed to converge.On the other hand, even when the sequence converges, the reason why the limit is the optimal solution is not discussed or verified.In Section 3, we will develop an approach to show that there is an optimal order quantity under a reasonable assumption, based on which the feasible domain with an upper bound and a lower bound can be derived.

Our Revision
We simplify the expression in 2.2 and 2.3 as follows: where If we observe the left hand side of 3.4 , then k is a safe factor with k ≥ 0. We try to locate interior solution for the first partial derivatives system, so we only consider k > 0. It turns out that an upper bound for the lot size, Q, is derived If we flip over and take square and then minus one on both side of 3.4 , it yields that

3.6
According to 3.5 , the relation, α 4 − βα 3 Q > 0, holds such that 3.6 is well defined and based on which we find a relation to express k as a function in Q so that Based on 3.7 , we derive that Using 3.7 and 3.8 , equation 3.1 turns into From 3.9 , we obtain the lower bound 3.10 Equation 3.5 yields α 3 Q/ α 4 − βα 3 Q < 1, and therefore we have 3.11 owing to 3.9

3.11
Next, we will compare the two upper bounds α 4 / 1 β α 3 from 3.5 and √ α 1 α 2 from 3.11 .We use the same data as that in Wu and Ouyang 1 .We compute the two upper bounds for different values of β and L i , and then list the ratio of α 4 / 1 β α 3 over Hence, by 3.5 , 3.10 and 3.11 , and our comparison in Table 1, we conclude that the feasible domain for the lot size, Q, is

3.12
We now solve 3.9 under the restriction of 3.12 by assuming the auxiliary function, 3.14 From 3.5 , we see that

3.15
Together with 3.12 , it yields that

3.16
By combining the results of 3.14 through 3.16 , we obtain Therefore, f Q is a convex function.The auxiliary function can be rewritten as Next, we will show that f 1, the following expression can be derived, We may separate the graph of f Q into two parts.On the left wing, f Q decreases to its minimum.On the right wing, f Q increases from its minimum.Accordingly, we divide the problem into two cases.For the first case, √ α 1 is on the left wing, so that from f √ α 1 < 0, f Q continuously decreases to its minimum, say f Q min , with f Q min < 0 and then f Q changes, to increase.Therefore, there is a unique point, say Q * , that satisfies For the second case, √ α 1 is on the right wing, with that there is a unique point, we still say Q * , that satisfies We are handling a minimum problem that is bounded below by zero so that the existence of the optimal solution is trivial.
Under our assumption, α 4 > √ α 1 α 2 1 β α 3 , we have shown that there is a unique solution, Q * from 3.21 or 3.22 , respectively, for the first derivative system such that Q * and k Q * , derived by 3.7 , are the optimal quantity and the minimum solution for the stochastic inventory model, respectively.

Numerical Examples
In order to demonstrate that our approach can derive the minimum solution, we refer to the previous example and list our computed results in Table 2. Q i is the unique solution for f Q 0 of 3.13 .k i and EAC u Q i , k i , L i are derived by 3.7 and 2.1 , respectively.From Table 2, we know the minimum for each backordered rate, and then we list the comparison between our findings and that of Wu and Ouyang 1 in Table 3.
By comparing the third and the fifth column of Table 3, it shows that our findings are slightly better than those of Wu and Ouyang 1 .Moreover, we try to consider the iterative method, which is the solution algorithm of Wu and Ouyang 1 .In 2.2 , Q is an explicit function of k.By plugging, into a value of k, we derive the corresponding value of Q.
However, in 2.3 , k is expressed as an implicit function of Q.When we plug into a value of Q, then 2.3 only derives the value of √ 1 k 2 / √ 1 k 2 − k .Hence, we will not use 2.3 .Instead, we apply the equivalent relation in 3.7 .Consequently, we consider 2.2 and 3.7 in studying the iterative algorithm of Wu and Ouyang 1 .
With k 0 0 proposed by them , and for example, β 0 and i 3, it yields that with θ 1 2D/hδ A 3 j 1 c j b j −a j , θ 2 Dσ/hδ π π 0 L 3 , and θ 3 h 1−E p /D π π 0 .The computation results for two sequences Q n and k n are listed in Table 4.
From Table 4, we see that by iterative algorithm proposed by Wu and Ouyang 1 , for β 0 and lead time, L 3 , Wu and Ouyang 1 should derive that k * 3.15 and Q * 179.74.However, they derived that Q * 183.It may indicate that the iterative method with two generated sequences may be too complex to execute, such that our approach provides an improvement to locate the optimal solution.

Table 3 :
The comparison between Wu and Ouyang 1 and ours.

Table 4 :
The iterative method of Wu and Ouyang 1 .