Producer–retailer integrated EMQ system with machine breakdown, rework failures, and a discontinuous inventory issuing policy

This study develops two extended economic manufacturing quantity (EMQ)-based models with a discontinuous product issuing policy, random machine breakdown, and rework failures. Various real conditions in production processes, end-product delivery, and intra-supply chains such as a producer–retailer integrated scheme are examined. The first model incorporates a discontinuous multi-delivery policy into a prior work (Chiu et al. in Proc Inst Mech Eng B J Eng 223:183–194, 2009) in lieu of their continuous policy. Such an enhanced model can address situations in supply chain environments, where finished products are transported to outside retail stores (or customers). The second model further combines retailer’s stock holding costs into the first model. This extended EMQ model is applicable in situations in present-day manufacturing firms where finished products are distributed to company’s own retail stores (or regional sales offices) and stocked there for sale. Two aforementioned extended EMQ models are investigated, respectively. Mathematical modeling along with iterative algorithms are employed to derive the optimal production run times that minimize the expected total system costs, including the costs incurred in production units, transportation, and retail stores, for these integrated EMQ systems. Numerical examples are provided to demonstrate the practical application of the research results.

of perfect quality. However, in real manufacturing environments, due to process deterioration or other uncontrollable factors, both production of items of imperfect quality and machine breakdown are inevitable. Unsurprisingly, many studies have been carried out to enhance the classic EMQ model by addressing issues of imperfect product quality and random machine breakdowns (Barlow and Proschan 1965;Shih 1980;Bielecki and Kumar 1988;Grosfeld-Nir and Gerchak 2002;Inderfurth et al. 2006;Hishamuddin et al. 2014;Chiu and Chang 2014;Lin et al. 2014;Wu et al. 2014;Khedlekar et al. 2014;Pal et al. 2015;Ocampo 2015;Chiu et al. 2015a, b, c). Henig and Gerchak (1990) conducted a comprehensive analysis of a general periodic review production/inventory model with variable yield. Groenevelt et al. (1992) proposed two production control policies to cope with machine breakdown. The first policy considers that the production of an interrupted lot will not be resumed after a breakdown (i.e., the no resumption or NR policy), while the second assumes that the production of an interrupted lot resumed immediately after production is restored and if the current on-hand inventory is below a certain threshold level (i.e., the abort-resume or AR policy). Wee (1993) proposed an economic production policy for deteriorating items with partial back-ordering. Two numerical examples were used to illustrate this proposed theory, and the computational results indicated that the policy led to a lower cost. Gopalan and Kannan (1994) examined a two-stage transfer-line production system with inspection and rework. Transient state characteristics were analyzed for the system, subject to an initial buffer of infinite capacity, inspection at both inter-and end-stages, and rework. A stochastic model was developed to study the system, and as a result, some explicit analytical expressions of system characteristics were revealed. Moinzadeh and Aggarwal (1997) investigated a production-inventory system with random disruptions. They proposed an (s, S) policy for the system in which the time between breakdowns is exponential, restoration times are constant, and excess demand is backordered. A procedure for deriving optimal values from the policy was developed, and the system parameters that minimize the expected total cost per unit time were examined. Jabal Ameli et al. (2008) proposed a multi-objective integer linear programming approach for cell formation problems with alternative process routings and machine reliability issues. Their study aimed to simultaneously minimize total cost and maximize system reliability. Unlike the traditional reliability evaluation approaches, the approach used in their study was to model machine unreliability in terms of cost and time-based effects. Using the e-constraint method as an optimization tool for multi-objective programming, a numerical example was provided to demonstrate the capability of the proposed model to evaluate various effects of reliability issues. Chiu et al. (2009) studied the optimal production run time in an EMQ model with imperfect rework and Poisson machine breakdowns under the abort/resume (A/R) control policy. In their proposed system, a random defective rate is assumed and all defective items are reworked at the end of regular production, and there exists a certain percentage of rework failures. The system is subject to random breakdowns and the A/R inventory control policy is adopted when breakdowns occur. Mathematical modeling was used, and theorems related to conditional convexity and bounds of optimal production run times were proposed and proved in their study. A recursive searching algorithm was developed to locate the optimal run time that minimizes the expected production-inventory costs.
Unlike the implicit assumption of continuous inventory issuing policy in the classic EMQ model, in real-life supply chain systems, the most commonly adopted policy for distribution of finished products is the discontinuous (periodic) multiple delivery policy. Schwarz (1973) studied a deterministic one-warehouse N-retailer inventory problem with a continuous review policy. Schwarz's study aimed to determine the stocking policy that minimizes the average system cost. Goyal and Gupta (1989) reviewed and classified the vendor-buyer coordinated inventory models. They not only provided a scheme for models classification but also identified a few possible directions for future research. Sarker and Parija (1994) derived the optimal batch size for a production system operating under a fixed-quantity and periodic delivery policy. They studied a manufacturing system where raw materials were procured from vendors, and then processed and converted into finished items. An integrated inventory model was proposed and analyzed to simultaneously determine the optimal ordering policy for raw materials and the most economical production batch size so that total system costs could be minimized. Hill (1999) determined the optimal production and shipment policy for the single-vendor single-buyer integrated production-inventory problem. It was assumed that a vendor manufactures a product in batches at a finite rate and ships it to a customer who then consumes the product at a fixed rate. The objective was to determine a purchasing and production schedule that minimizes the overall system costs. As a result, a global optimal solution was derived. Viswanathan and Piplani (2001) proposed a model with a specific vendor-buyer coordinating discipline to analyze the benefit of supply chain inventories by using common replenishment epochs or time periods. A one-vendor, multi-buyer supply chain for a single product was analyzed under a specific strategy, whereby the vendor specifies common replenishment periods and requires all buyers to replenish only at those time periods. In return, the vendor offers a price discount to convince buyers to accept this strategy. As a result, they determined the optimal replenishment period and the price discount to be offered by the vendor for the proposed model. Giri and Maiti (2012) studied a supply chain model for a deteriorating product with time-varying demand and production rate. A single-vendor single-buyer two-echelon supply chain model was examined, whereby the buyer sells a seasonal product and its inventory is subject to deterioration, and the vendor's production rate is dependent on the buyer's demand rate, which is a linear function of time. Further, some nonconforming items may be randomly produced during a production run. Mathematical modeling was used to derive the average cost function of the proposed supply chain system. An algorithm for finding the optimal solution was developed. Additional studies (Kreng and Chen 2007;Gong and Chen 2012;Glock 2012;Wee and Widyadana 2013;Glock et al. 2014;Tseng et al. 2014;Safaei 2014;Rodger 2014) focused on various aspects of supply chain optimization.
To address real-life production-shipment situations, this study first extends Chiu et al. 's work (2009) by incorporating a discontinuous multi-delivery policy into their model in lieu of a continuous policy assumption to examine situations in real supply chain environments, where finished products are transported to retail stores (or customers) outside the production units. Next, the retailer's stock holding cost is incorporated into the proposed model in the second part of the study to explicitly deal with situations in present-day manufacturing firms where finished products are distributed to companies' own retail stores (or regional sales offices) and stocked there for sale. Accordingly, two separate extended EMQ models are developed, and with the aid of mathematical modeling and optimization techniques, we derive the optimal production run times that minimize the expected total system costs, including the costs incurred in production units, transportation, and retail stores, for these integrated EMQ systems. This study intends to fill the research gap in this specific area.

Model 1: Description, modeling, and solution process
The first proposed model in this study explores the optimal production run time for an economic manufacturing quantity model with a discontinuous inventory issuing policy, random machine breakdown, and rework failures. The first model extends the work of Chiu et al. (2009) by incorporating a discontinuous multi-delivery policy into their model in lieu of a continuous policy assumption. This enhanced model can be used to address the situations in real supply chain environments, where finished products are transported to the retail store or customer outside the production units.
Consider that a manufacturing system can produce a product at an annual production rate P 1 and the demand rate for this product is λ units per year. The production process may randomly produce x portion of defective items at a rate d 1 , where d 1 = P 1 x. In order to sustain regular operations (i.e., avoid the occurrence of shortage) (P 1 − d 1 − λ) > 0 must be satisfied. All products manufactured are screened and the unit inspection cost is included in the unit production cost C. All nonconforming items are reworked at a rate of P 2 , and the rework process starts immediately after the completion of the regular process. Since 100 % of reworks are not successful, a θ 1 portion (where 0 ≤ θ 1 < 1) of reworked items fails and becomes scraps.
It is assumed that finished products can only be delivered to the retail store if the entire production lot is quality assured after rework. A multi-delivery policy is used for transporting finished items. Under such a discontinuous product issuing policy, n fixed quantity installments of the finished lot are delivered at a fixed interval of time during delivery time t ′ 3 (see Fig. 1). Furthermore, during production time the machine is subject to a random breakdown which follows a Poisson distribution. The abort/resume (A/R) inventory control policy is used when a breakdown takes place. Under such a policy, the machine goes under repair immediately after a breakdown occurs, a constant machine repair time is assumed, and the interrupted lot will be resumed right after the restoration of machine (see Fig. 1).
In addition to the unit production cost C, other cost-related parameters used in the proposed modeling and analysis include the production setup cost K, a fixed cost M for machine repairing, unit holding cost h at the producer's end, unit rework cost C R , holding cost h 1 for each reworked item, disposal cost per scrap item C S , fixed delivery cost K 1 per shipment, unit shipping cost C T , and unit holding cost h 2 at the retailer's end. Additional notation used is listed below: the production time (i.e. uptime) to be determined in the proposed EMQ model, t: production time before a random breakdown occurs, β: number of breakdowns per year, a random variable which follows Poisson distribution, the maximum level of on-hand inventory in units when regular production ends (in the case of breakdown occurrence), level of on-hand perfect quality items at time t, I d (t): level of on-hand defective items at time t, I S (t): level of on-hand scrap items at time t, Q: production lot size per cycle, T ′ : cycle length (in the case of breakdown occurrence), h 3 : unit holding cost for safety stock, TC 1 (t 1 ): the total production-inventory-delivery costs per cycle (in the case of breakdown occurrence), E[TC 1 (t 1 )]: the expected production-inventory-delivery costs per cycle (in the case of breakdown occurrence), t 2 : time required to rework the defective items (in the case of no breakdown occurrence), t 3 : time required for transport all available perfect quality items (in the case of no breakdown occurrence), H 1 : the level of on-hand inventory in units when regular production ends (in the case of no breakdown occurrence), H: the maximum level of on-hand inventory in units when the rework process ends (in the case of no breakdown occurrence), TC 2 (t 1 ): total production-inventory-delivery costs per cycle when no breakdown occurrence, E[TC 2 (t 1 )]: the expected production-inventory-delivery costs per cycle (in the case of no breakdown occurrence), T: cycle length (in the case of no breakdown occurrence),

Fig. 1
On-hand inventory level of perfect quality products in the extended EMQ model 1 with rework failures, machine breakdown, and multi-delivery policy TCU(t 1 ): the total production-inventory-delivery costs per unit time whether a breakdown takes place or not, E[TCU(t 1 )]: the long-run expected production-inventory-delivery costs per unit time whether or not a breakdown takes place, T: the cycle length whether or not a machine breakdown takes place.
Suppose t denotes time before a random breakdown taking place during production uptime t 1 , then the following two different cases must be analyzed, respectively.

Case 1: A breakdown occurs during uptime t 1
In this case, t < t 1 . That is the time before a random breakdown taking place is shorter than the production uptime. In other words, a machine breakdown takes place during the production process. Under the AR inventory policy, the machine goes under repair immediately, and, once it is fixed and restored, the production of interrupted lot is resumed right away. The on-hand inventory level of perfect quality products at the time a breakdown occurs is H ′ 2 (refer to Fig. 1), and the level of on-hand inventory remains at H ′ 2 until the machine is repaired. While H ′ 1 denotes the level of on-hand inventory in units when regular production of the remaining of interrupted lot is completed. Then, the reworking of defective products begins at a rate of P 2 per year (see Fig. 2) and the maximum number of defective items per cycle is given in Eq. (1).
It is assumed that 100 % of reworks are not successful, a θ 1 portion of reworked items fails and becomes scrap during t ′ 2 (see Fig. 3) and Eq. (2) shows maximum number of scrap items.
Upon the completion of rework, the on-hand inventory level perfect quality products is H ′ (Fig. 1). Next, in product delivery time t ′ 3 , fixed quantity of n installments of finished batch is delivered to the customer at a fixed interval of time. It can be seen from Fig. 1 Fig. 2 On-hand inventory level of defective products in the extended EMQ model 1 In the case of breakdown occurrence, the total production-inventory-delivery costs per cycle, TC 1 (t 1 ) consists of (1) the variable production cost; (2) production setup cost; (3) machine repair cost; (4) rework cost; (5) disposal cost for scraps; (6) fixed and the variable product transporting costs; (7) inventory holding cost during rework; (8) holding cost for safety stocks (i.e., stocks to prevent shortage occurrence due to a machine breakdown); and (9) inventory holding costs during the entire production cycle. Thus, TC 1 (t 1 ) can be expressed as In this study, it is assumed that x is a random variable with a known probability density function, thus the expected values of x are used in our analysis to take the randomness of x into account. By substituting all related system parameters into Eq. (4) (Chiu et al. 2009) and with further derivations, the expected system costs per cycle for the case of breakdown occurrence, E[TC 1 (t 1 )] can be obtained as

Case 2: No breakdown occurs during uptime t 1
In this case, t > t 1 . That is the time before a random breakdown taking place is longer than the production uptime. In other words, no breakdown occurs during the production process (Chiu et al. 2009).
In case 2 the production cycle length is T = t 1 + t 2 + t 3 , and total production-inventory-delivery costs per cycle (in no breakdown occurrence case) TC 2 (t 1 ) is Fig. 3 On-hand inventory level of scrap items in the extended EMQ model 1 To take randomness of defective items into account and substitute all related parameters in Eq. (6), and with further derivations E[TC 2 (t 1 )] can be obtained as follows:

Integration of the proposed EMQ models with and without breakdown
This study assumes a breakdown can occur randomly and it follows a Poisson distribution with mean equals to β per year. Let f(t) represent the probability density function of random time t before a breakdown occurs and F(t) denote the cumulative density function of t. Then, the expected production-inventory-delivery costs per unit time, E[TCU(t 1 )] is

The expected cycle length E[T] is
Since the number of breakdown per unit time is a random variable that follows a Poisson distribution with mean equal to β. The time between breakdowns obeys the Exponential distribution with density function f(t) = βe −βt , and the cumulative density function F(t) = 1 − e −βt .
By substituting E[TC 1 (t 1 )], E[TC 2 (t 1 )], and E[T] into Eq. (8) and solving the integration of the mean time to breakdown in E[TCU(t 1 )], we obtain where (6) TC 2 (t 1 ) = C(P 1 t 1 ) + K + C R (xt 1 P 1 ) + C S (θ 1 xt 1 P 1 ) and Determining the optimal production run time With the aim of determining the optimal production run time t * 1 , we first must prove that E[TCU(t 1 )] is convex. Let w(t 1 ) represent the following: Equation (14) shows the second derivative of E[TCU(t 1 )] with respect to t 1 .
The first term in the right-hand side (RHS) of Eq. (14) is positive because annual demand λ > 0. Therefore, we obtain With further derivations, the left-hand size (LHS) in Eq. (15) Therefore, Eq. (15) becomes or If E[TCU(t 1 )] is convex, then to search for the optimal production run time t * 1 , we can set the first derivative of E[TCU(t 1 )] = 0 as follows: Again, the first term in the RHS of Eq. (21) is positive, so the second term equals to zero. In order to find the bounds for t * 1 , let To prove t * 1 falls within bounds, we first multiply the second term of Eq. (21) by 2P 1 t 2 1 β and obtain In order to search for the optimal t * 1 , Eq. (24) can be rearranged as or Because e −βt 1 is the complement of the cumulative density function F (t 1 ) = 1 − e −βt 1 . As 0 ≤ F (t 1 ) ≤ 1, so 0 ≤ e −βt 1 ≤ 1. Let e −βt 1 = 0 and e −βt 1 = 1 represent bounds for e −βt 1 , respectively. By substituting them into Eq. (25), we obtain and t * 1L < t * 1 < t * 1U . It can be seen that although the optimal t * 1 cannot be expressed in a closed form, it falls within the aforementioned bounds, and it can be located with the use of a proposed recursive searching algorithm as follows.
Since e −βt 1 is the complement of cumulative density function, thus, 0 ≤ e −βt 1 ≤ 1. Let The proposed recursive searching algorithm to find t * 1 is listed below: 1. Let u(t 1 ) = 0 and u(t 1 ) = 1 initially and compute the upper and lower bounds for t * 1 , respectively (i.e., to obtain the initial values of [t * 1L , t * 1U ]).

Substitute the current values of [t *
1L , t * 1U ] into e −βt 1 and calculate the new bounds (expressed as u L and u U ) for e −βt 1 . Hence, u L < u(t 1 ) < u U . 3. Let u(t 1 ) = u L and u(t 1 ) = u U , and re-compute the new upper and lower bounds for t * 1 , respectively (i.e. to update the current values of [t * 1L , t * 1U ]). 4. Repeat steps 2 and 3, until there is no significant difference between t * 1L and t * 1U (or there is no significant difference in terms of their effects on E[TCU(t * 1 )]). 5. Stop. The optimal production run time t * 1 is obtained.

Problem description and modeling
In the manufacturing sector, some producers of consumer goods may have their own retail stores or regional sales offices to promote and sell their finished products to (Fig. 4). In order to address such a real life intra-supply chains situation, the second model of this study enhances the first model to further explore the optimal run time for a producer-retailer integrated EMQ model, wherein, the retailer's stock holding cost is incorporated into the first model. According to the proposed multi-delivery policy, Fig. 5 depicts the stock holding status at the retailer's side for the case of breakdown taking place. Additional notations used at the modeling and analysis of this extended producer-retailer integrated EMQ model 2 are listed as follows.

D:
number of finished items (a fixed quantity) distributed to retail store per delivery, I: number of left over items in each interval of time t n , after satisfying the demand in t n , Fig. 4 Extension to a producer-retailer integrated EMQ system (model 2)

Fig. 5
On-hand Inventory level of finished products at the retailer's side in the extended EMQ model 2 with breakdown occurrence I c (t): on-hand inventory level of finished products at the retailer's side at time t, h 2 : unit holding cost for products stored at the retailer's side, TC 3 (t 1 ): total production-inventory-delivery costs per cycle in this enhanced model (in the case of breakdown occurrence), TC 4 (t 1 ): total production-inventory-delivery costs per cycle in this enhanced model (in the case of no breakdown occurrence), E[TC 3 (t 1 )]: the expected production-inventory-delivery costs per cycle in this enhanced model (in the case of breakdown occurrence), E[TC 4 (t 1 )]: the expected production-inventory-delivery costs per cycle in this enhanced model (in the case of no breakdown occurrence), E[TCU 2 (t 1 )]: the long-run expected production-inventory-delivery costs per unit time in this enhanced model, whether a breakdown occurring or not.
It can be seen from Fig. 5, at the retailer's side the demand between shipments is λt n , the number of left over items in each interval of time t n (after satisfying the demand in t n ) is Total inventory holding costs on the retailer's side with and without breakdown are expressed in Eqs. (29) and (30), respectively (also refer to Fig. 5): By incorporating the retailer's holding costs into the original models with and without breakdown, respectively (i.e., into Eqs. (4) and (6) in "Model 1: Description, modeling, and solution process" section), we obtain To take the randomness of defective items into account and substitute all related parameters into Eqs. (31) and (32), and with further derivations, we obtain E[TC 3 (t 1 )] and E[TC 4 (t 1 )] as

Integration of producer-retailer integrated models with and without breakdown
As stated in section "Integration of the proposed EMQ models with and without breakdown", in the proposed study a breakdown can occur randomly and the mean time to machine breakdowns obeys an exponential distribution with density function f (t) = βe −βt , and its cumulative density function F(t) = 1 − e −βt . Hence, the expected production-inventory-delivery costs per unit time, E[TCU 2 (t 1 )] is By substituting E[TC 3 (t 1 )], E[TC 4 (t 1 )], and E[T] into Eq. (35) and solving the integration of the mean time to breakdown in E[TCU 2 (t 1 )], we obtain where δ 1 and δ 2 are provided in Eqs. (11) and (12), and let δ 5 is denote the following Determining the optimal production run time Before determining the optimal production run time t * 1 , we must first prove that E[TCU 2 (t 1 )] is convex. Let π(t 1 ) represent the following: Equation (39) shows the second derivative of E[TCU 2 (t 1 )] with respect to t 1 .
The first term in the RHS of Eq. (39) is positive because annual demand λ > 0. Therefore, we have With further derivations, the left-hand side (LHS) of Eq. (40) becomes Let Then Eq. (40) becomes or If E[TCU 2 (t 1 )] is proven to be convex, then the optimal production run time t * 1 can be solved by setting the first derivative of E[TCU 2 (t 1 )] = 0.
Since the first term in the RHS of Eq. (44) is positive, and so the second term is equal to zero. In order to find the bounds for t * 1 , let Please refer to the similar proof as that for Theorem 2 (presented in "Model 1: Description, modeling, and solution process" section).
Upon showing that t * 1 falls within the upper and lower bounds, by multiplying the second term of Eq. (44) by 2P 1 t 2 1 β one has With further rearrangement, one obtains Because e −βt 1 is the complement of the cumulative density function F (t 1 ) = 1 − e −βt 1 . As 0 ≤ F (t 1 ) ≤ 1, so 0 ≤ e −βt 1 ≤ 1. Let e −βt 1 = 0 and e −βt 1 = 1 represent the upper and lower bounds of e −βt 1 , respectively. Then, by applying the proposed recursive searching algorithm provided at the end of section "Model 1: Description, modeling, and solution process", we can find the optimal production run time t * 1 .

Numerical demonstration of proposed EMQ model 1
This section uses the same numerical example as in Chiu et al. (2009) to provide a comparison for readers. The following values are used for the corresponding system variables: β Poisson breakdown rate with a mean of 0.5 per year, M fixed machine repair cost of $500 per breakdown, g machine repair time (t r ) of 0.018 year, λ demand rate of 4000 items per year, P 1 production rate of 10,000 items per year, x uniform distribution nonconforming rate over the range [0, 0.2], K setup cost of $450 per cycle, C unit manufacturing cost of $2, h holding cost of $0.6 per item per unit time, P 2 rework rate of 5000 items per year, C R unit rework cost of $0.50 per product per year, h 1 holding cost of $0.80 for each reworked item, θ 1 rework failure rate of 0.1, C S unit disposal cost of $0.3 for scrap items, K 1 fixed transportation cost of $90 per delivery, n number of deliveries per cycle: 4, C T variable delivery cost of $0.001 per product shipped To test for the convexity of E[TCU(t 1 )] (Theorem 1), we can apply the upper and lower bounds of t * 1 (i.e., Eqs. (22) and (23)) to Eq. (13), and obtain that t * 1U = 0.4595 < w(t * 1U ) = 2.87230, and t * 1L = 0.3056 < w(t * 1L ) = 2.62608. Therefore, the expected costs function [TCU(t 1 )] is convex.
Next, to find the optimal production run time t * 1 , we substitute t * 1U and t * 1L into Eq. (10), and find that E[TCU(t * 1U )] = $10,417.21 and E[TCU(t * 1L )] = $10,323.69, respectively. Since the optimal run time t * 1 falls within the bounds of t * 1L and t * 1U , by applying the proposed searching algorithm, we find that t * 1 = 0.3320 years and that the optimal system cost E[TCU(t * 1 )] = $10,317.27. Figure 6 depicts the effect of variations in production run time t 1 on the expected system cost E[TCU(t 1 )]. Table 1 exhibits the step-by-step iterations of the proposed searching algorithm.
The proposed EMQ-based model 1 is intended to address the effect of discontinuous multi-delivery policy on the EMQ model with machine breakdown and rework failures. For this reason, sensitivity analyses of variations in the fixed transportation cost and their effects on the expected system cost E[TCU(t * 1 )] and on the production run time t * 1 are carried out.
The results indicate that, as the fixed transportation costs K 1 (or the ratio of K 1 /K) increase, the production run time t * 1 increases significantly (as illustrated in Fig. 7) along with the expected system cost E[TCU(t * 1 )] (see Fig. 8).

Numerical demonstration of proposed EMQ model 2
This subsection provides a numerical example to explain the proposed producer-retailer integrated EMQ model 2. First, we assume an extra parameter-unit holding cost on the retailer's side h 2 = $1.50. By applying Eqs. (45) and (46), we obtain that t * 1U = 0.2848 and t * 1L = 0.1951, respectively. Then, we use them to test for the convexity of the system cost E[TCU(t * 1 )] (Theorem 3) and find that t * 1U < π(t * 1U ) = 2.5894 and t * 1L < π(t * 1L ) = 2.4511. Hence, the system cost function [TCU(t 1 )] is convex.
Once the convexity of [TCU(t 1 )] is proven, we apply the proposed recursive algorithm (Theorem 4) to search for the optimal run time t * 1 over the interval [t * 1L , t * 1U ]. By substituting t * 1U and t * 1L in Eq. (36), we obtain the initial values of E[TCU(t * 1U )] = $11,642.46 and E[TCU(t * 1L )] = $11,481.07. Then, by further applying the algorithm (as provided at the end of "Model 1: Description, modeling, and solution process" section), we find that t * 1 = 0.2051 years and E[TCU(t * 1 )] = $11,477.27. The step-by-step iterative results are shown in Table 2.
A further analysis on the behavior of the system costs E[TCU(t 1 )] related to the production run time t 1 is depicted in Fig. 9. It shows that, by simply applying the result derived from the second EMQ model, the management of such a producer-retailer Fig. 7 The effects of variations in the K 1 /K ratios on the optimal t * 1 in the extended EMQ model 1 Fig. 8 The effects of variations in the K 1 /K ratios on the optimal E[TCU(t * 1 )] in the extended EMQ model 1 integrated system can help realize cost savings of $359.00 (or 3.13 % over the total system costs) as compared to the result from the EMQ model 1.
The purpose of the proposed EMQ-based model 2 is to explore the optimal run time t * 1 for a producer-retailer integrated EMQ model and to study the effect of the retailer's holding costs on the replenishment decision. For this reason, sensitivity analyses of variations in h 2 (or the ratio of h 2 /h) and their effects on t * 1 and the expected system cost E[TCU(t * 1 )] are performed (see Figs. 10, 11). The figures show that, as the retailer's holding costs h 2 (or the ratio of h 2 /h) increase, the optimal production run time t * 1 decreases and the expected cost E[TCU(t * 1 )]  . 9 The behavior of E[TCU 2 (t 1 )] and t 1 in the extended EMQ model 2 Fig. 10 The effects of variations in the h 2 /h ratios on the optimal t * 1 in the extended EMQ model 2 increases. For management of the producer-retailer integrated system, the results of these sensitivity analyses can provide helpful information during the decision-making process, by offering insights into the effects of various inventory-holding costs in different retail stores.

Conclusions
In this study, two extended EMQ-based models with a discontinuous product issuing policy, random machine breakdown, and rework failures are developed. Various real conditions in the production process, end-product delivery, and intra-supply chains such as a producer-retailer integrated scheme are examined. With the aid of mathematical modeling and optimization techniques, we derive the optimal replenishment run time decision and reveal various important factors for the system parameters of the studied models (refer to Figs. 6,7,8,9,10,11 as examples). The obtained results can help production planners determine the optimal production run time in an intra-supply chain situations where finished products are distributed to companies' own retail stores or regional sales offices and stocked there for sale. An interesting area for future study would be the effect of variable production rates on these models. Fig. 11 The effects of variations in the h 2 /h ratios on the optimal E[TCU 2 (t * 1 )] in the extended EMQ model 2