A Sustainable Economic Recycle Quantity Model for Imperfect Production System with Shortages

: Recycling of products has a great impact on contemporary sustainable business strategies. In this study, a sustainable recycling process in a production-inventory model for an imperfect production system with a ﬁxed ratio of recyclable defective products is introduced. The piecewise constant demand rates of the non-defective items are considered under production run-time, production off-time with positive stock, and production off-time with shortages under varying conditions. Based on the production process, two cases are studied using this model. The ﬁrst case does not consider recycling processes, while the second case picks up all defective items before sending these items to recycling during the production off-time; the recycled items are added to the main inventory. The aim of this study is to minimize the total cost and identify the optimal order quantity. The manufacturing process with the recycling process provides a better result compared to without recycling in the ﬁrst case. Some theoretical derivations are developed to enunciate the objective function using the classical optimization technique. To validate the proposed study, sensitivity analysis is performed, and numerical examples are given. Finally, some managerial insights and the scope of future research are provided.


Introduction and Literature Review
Inventory management control monitors and determines the optimal inventory order to minimize the total inventory cost. In every business, retailers often implement numerous strategies to intensify the profitability and sustainability of the system (Mashud et al. 2021c;Popescu and Popescu 2019). The familiar economic order quantity (EOQ) model was first initiated by Harris (1915). Later, Alharkan et al. (2020) worked on a two-echelon supply chain model with lot-size-dependent lead-time, and Taft (1918) extended the EOQ model to consider incremental orders during the production process. Inventory of goods is essential to the production processes to meet customer demand. The standard economic production quantity (EPQ) models typically assume perfect product quality and perfect production processes. In reality, however, ideal production processes are rare, and defective items will occur. Many industries recycle these faulty items. Some examples of recyclable materials are jute, paper, glass, and cotton and other fabrics.
During the last two decades, many researchers have studied inventory models with imperfect production systems. Defective items are mainly generated due to faulty production process (Rosenblatt and Lee 1986;Hsu and Hsu 2014). Salameh and Jaber (2000) attempted to modify the traditional EPQ model to consider poor-quality items, while Chiu et al. (2014) investigated lot-sizing problems with arbitrary defective rates in finite manufacture systems under rework processes and shortages. Defective items also deteriorate in some cases and emit carbon into the environment (Daryanto and Wee 2019;Mashud et al. 2020a). Despite deterioration, some products are wasted during the production process in the whole supply chain. Oh and Hwang (2006) anticipated an inventory model considering recycled raw materials without shortage. They assumed unequal holding costs for the raw materials and serviceable objects. Benkherouf and Omar (2017) investigated optimal production batch sizing with rework procedures, while Mokhtari (2018) outlined an order lot size combination between producers' and suppliers' defective manufacturing systems considering production rework processes. Al-Salamah (2019) designed an EPQ model for defective manufacturing procedures using synchronous and asynchronous elastic rework rates. A model for a long-run single-stage manufacturing system with imperfect items was developed by Kang et al. (2018) who incorporated rework operation, the inspection process, and planned backordering. In that model, the unfulfilled customer demands during the production phase due to process imperfections were satisfied either at the end of the inspection process or after reworking the imperfect products. Hasan et al. (2020) provided a model for agricultural products where to prevent product deterioration, they used product separation and offered different discount policies, while Nobil et al. (2020) presented a model for imperfect production process with reorder point in an EOQ model. An imperfect multi-stage production system considering defective proportions in the production process and uncertain product demand, and incorporating lean philosophy, was introduced by Tayyab et al. (2019b) who tried to reduce system costs. Onwude et al. (2020) analyzed a food wastage supply chain for fresh agricultural products. However, none of these models reworked defective items at each production stage. This model assumes the reworking of defective items at each production stage to bridge the aforementioned gap, and consequently the defective items (scraped) are perfect after reworking. Very few studies have been carried out with consideration of the recycling process with shortages, so the mentioned model also tries to incorporate shortages in the proposed ERQ model.
Demand is a fundamental key factor in inventory models because it determines the perplexity of optimal results. Sometimes demand becomes uncertain (Mashud et al. 2020b) and sometimes it depends on the price of the products (Mashud et al. 2021a(Mashud et al. , 2021c. Fuzzy demand is also very popular among the practitioners of inventory research (Tayyab et al. 2019a). Bai and Varanasi (1996) formulated an optimal control problem with indefinite piecewise constant demand and developed an algorithm to determine the optimal production strategy using Pontryagin's minimum principle. Hsu and Hsu (2016) considered a model for defective items with stock-reliant demand with inflation for a Weibull distribution deterioration. Palanivel and Uthayakumar (2016) modified their model by adding payment delays. Viji and Karthikeyan (2016) developed a model for a three-level production process for Weibull distribution deterioration. Recently, Manna et al. (2017) offered an EPQ model for a production rate-reliant defective rate of items under advertisement with demand dependent on a defective production system.
In the collected works of imperfect production systems, most models were developed considering reworks, repair, remanufacture, and different production and supply policies. Some researchers considered two-level and three-level production rates, but production inventory models of three-level constant demand patterns are not available in the literature survey. The main contributions of the proposed model are: (i) The proposed model is developed focusing on the recycling of defective items collected from regular production after the proper screening with 100% recovery of raw materials used for production in the corresponding item category. (ii) An inventory model is developed for three-level piecewise constant demand, which varies under three different production time parts, from production run-time to production off-time with positive stock and production off-time with shortages.
The remainder of the paper is organized as follows: Notation and assumptions of the model are described in Section 2. Section 3 presents the mathematical formulation of the models. Theoretical results along with optimal solution uniqueness are shown in Section 4. Section 5 comprises numerical examples and convexity graphs. Sensitivity with cost-benefit is analyzed in Section 6. Managerial implementations are delivered in Section 7. Lastly, the conclusion and future research are discussed in Section 8.

Assumptions and Notation
In this section, some assumptions and notation necessary to the development of the paper are presented.

Assumptions
The following assumptions are worthy of mention in formulating the models discussed.
(i) The production rate is finite and fixed. (ii) The defective rate is known, and constant and defective items are produced randomly alongside the perfect product. (iii) The demand rate of the perfect product is a piecewise constant function (motivated from Bai and Varanasi (1996)): The demand during production-run and production off-times and stock-out period is as follows: during production − run time xD; during production − off time with positive stock yD; during production − off time with shortages (iv) The sum of the demand rate and defective rate is less than the rate of production.
(v) Lead-time is negligible, and the number of shortages is acceptable while the model is considered for a single items production system. (vi) Defective items are recyclable with 100% recovery of raw material usable for the production of same product in the next cycle time, and the value of the recycled materials is higher than that of purchased materials. (vii) The holding cost per defective item is same as that of fresh items, while the recycle cost per defective item is equal irrespective of its number.

Notation
The following notations (Table 1) are used in the development of the model:

Model Formulation
During the production run-time of cycle time T, a lot size of Q units is produced at a production rate P, and with it, W units of defective items are produced at a defective rate d. The production process runs during [0, K], and during the production-off time, the stock depletes due to customer demand. At time C, the stock becomes zero, and shortages start until G. Current demand is expressed at the demand rate "D" during production run-time, which is calculated again using a ratio coefficient to obtain a rate "xD" during production off-time while the stock is positive. Shortages also occur at the demand rate "yD" during the time t 3 while the stock is negative. Thus, several inventory flow situations exist for different values of the ratios "x" and "y", for which two examples are depicted in the following Figures 1 and 2

Model Formulation
During the production run-time of cycle time T, a lot size of Q units is produced at a production rate P, and with it, W units of defective items are produced at a defective rate d. The production process runs during [0, K], and during the production-off time, the stock depletes due to customer demand. At time C, the stock becomes zero, and shortages start until G. Current demand is expressed at the demand rate "D" during production runtime, which is calculated again using a ratio coefficient to obtain a rate " " during production off-time while the stock is positive. Shortages also occur at the demand rate " " during the time while the stock is negative. Thus, several inventory flow situations exist for different values of the ratios "x" and "y", for which two examples are depicted in the following Figures 1 and 2.   The stock is positive during the replenishment period, which increases by "(P − D − d)" units of non-defective items over time through the line OA of Figure 2 during the time period t1. Replenishment ends at point A with a maximum number of (non-defective) units. Therefore, Production stops at the end of the time period t1 at point K, and the inventory decreases at another demand rate of xD units per unit time, expressed in Figure 2 as line AC during t2. At the end of time period t2, the inventory reaches point C with zero on-hand stock. In this case, if the demand rate remains the same as "D", then the inventory decreases as indicated by the dashed line AB. Therefore,

=
(2) Therefore, At the end of time t2, thestock is 0 and then decreases at the demand rate yD units per unit time. During time t3, shortages reach a total of units.

= (4)
As soon as the shortage level reaches Qs at point E, production starts. The shortage amount Qs is served at the rate of (P − D − d) units with the current demand, which is fulfilled at the rate of D during timet4. After this period, the stock becomes zero again at the point F. Therefore, Hence, A total of Q units is produced in + time where W units are defective. Therefore, The following relation is obtained using Equation (7) in addition to Equations (1) and (5). The stock is positive during the replenishment period, which increases by "(P − D − d)" units of non-defective items over time through the line OA of Figure 2 during the time period t 1 . Replenishment ends at point A with a maximum number of Q d (non-defective) units. Therefore, Production stops at the end of the time period t 1 at point K, and the inventory decreases at another demand rate of xD units per unit time, expressed in Figure 2 as line AC during t 2 . At the end of time period t 2 , the inventory reaches point C with zero on-hand stock. In this case, if the demand rate remains the same as "D", then the inventory decreases as indicated by the dashed line AB. Therefore, Therefore, At the end of time t 2 , thestock is 0 and then decreases at the demand rate yD units per unit time. During time t 3 , shortages reach a total of Q s units.
As soon as the shortage level reaches Q s at point E, production starts. The shortage amount Q s is served at the rate of (P − D − d) units with the current demand, which is fulfilled at the rate of D during timet 4 . After this period, the stock becomes zero again at the point F. Therefore, Hence, A total of Q units is produced in t 1 + t 4 time where W units are defective. Therefore, The following relation is obtained using Equation (7) in addition to Equations (1) and (5).
If there are no defective items (i.e., d = 0) and no demand (i.e., D = 0), then the inventory is replenished through the line OH.
Given that the ratios of demand rates x and y modify the total demand D during production off-time with positive and negative stock respectively. The average demand rate in a cycle becomes 2+x+y 4m D and requires a total of (Q − W) units of fresh items for the cycle time T. Therefore, Thus and "m" is an arbitrary constant that has a single value for every set of system parameter values. The theoretical value "m" is explained in Section 4.3.
The foremost target of this work is to minimize the probable total cost of inventory where defective items are formed with perfect goods, but defective items are recyclable, and recycled raw materials are reusable for the production of the same category of new product. Firstly, the manufacturer must define all costs, production features, all abilities of the production procedure, and their position on recycling. The inventory holding cost of the inventory depends on average inventory. Average inventory also increases if defective items are stored properly during production run-time before they are sent to recycling. Recycling processes also incur a cost per defective item.
The following two cases are discussed regarding recycling of defective items. The total cost function of the inventory has the following components.
(1) Average setup cost or fixed cost (FC): (2) Average production cost (PC): (3) Average raw material cost (RMC): Each production lot size is a quantity for which raw materials are purchased for exact production in each time cycle. Therefore, (4) Average holding cost (HC): Defective items have no holding cost in this case. Therefore, Therefore, Hereafter, the total cost function is TCD = FC + PC + RMC + HC + SC The objective is to discover the optimal number of defective items and the shortage level in order to minimize total cost of the inventory per unit time. The corresponding optimal production lot size (EPQ) is also consequently determined.

Case 2 (ERQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)
The producer receives a fixed portion of recyclable defective items during regular production, and defective items are sent to be recycled during production off-time. Recyclable materials become raw materials again for the new production cycle, yet to meet the demand, supplementary raw materials must be purchased from the suppliers. The target of the study is to control raw material procurement and the recycle size and production procedure to fulfill demand while minimizing the total inventory cost. The recycling system is shown in Figure 3. J. Risk Financial Manag. 2021, 14, x FOR PEER REVIEW 7 of 21 Defective items have no holding cost in this case. Therefore, Total inventory = ( + ) Hereafter, the total cost function is TCD = FC + PC + RMC + HC + SC The objective is to discover the optimal number of defective items and the shortage level in order to minimize total cost of the inventory per unit time. The corresponding optimal production lot size (EPQ) is also consequently determined.

Case 2 (ERQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)
The producer receives a fixed portion of recyclable defective items during regular production, and defective items are sent to be recycled during production off-time. Recyclable materials become raw materials again for the new production cycle, yet to meet the demand, supplementary raw materials must be purchased from the suppliers. The target of the study is to control raw material procurement and the recycle size and production procedure to fulfill demand while minimizing the total inventory cost. The recycling system is shown in Figure 3. In this case, the production cycle starts with a shortage of amount in the inventory, which accumulates at the rate ( − − ) over the time . During this time, the current demand is also fulfilled at the demand rate . The inventory accumulates at the rate ( − − ) over time , and production stops at the end of . After that, the inventory level starts to decrease due to the demand rate xD over the time and due to the third level of demand rate yD, while shortages reach overtime . In the production-runtime ( + ), W units of defective items are produced and are recycled during production off-time ( + ) and added to the raw material of the next production cycle. Consequently, the inventory of the ERQ model procures raw materials for ( − ) units of items in each cycle time, and the process repeats. In this case, the production cycle starts with a shortage of amount Q s in the inventory, which accumulates at the rate (P − D − d) over the time t 4 . During this time, the current demand is also fulfilled at the demand rate D. The inventory accumulates at the rate (P − D − d) over time t 1 , and production stops at the end of t 1 . After that, the inventory level starts to decrease due to the demand rate xD over the time t 2 and due to the third level of demand rate yD, while shortages reach Q s overtime t 3 . In the production-runtime (t 1 + t 4 ), W units of defective items are produced and are recycled during production off-time (t 2 + t 3 ) and added to the raw material of the next production cycle. Consequently, the inventory of the ERQ model procures raw materials for (Q − W) units of items in each cycle time, and the process repeats.
The total cost function of the inventory system comprises the following cost components: (1) Average production cost (PC) is: (2) Average raw material cost (RMC): Each production lot size is Q, but raw materials for (Q − W) units of items in each time cycle are purchased. Therefore, (3) Average holding cost (HC): Holding cost of defective items is included during time (t 1 + t 4 ) before the items are sent for recycling. Therefore, Total inventory Hence, (4) Average shortage cost (SC): Hence, the total cost function is TCR = FC + PC + RMC + HC + SC + RC The objective is to find out the optimal number of defective items (ERQ) to be recycled and the optimal shortage level so as to minimize the total cost of the inventory per unit time. The corresponding optimal production lot size (EPQ) is also consequently determined.
In the next section, the uniqueness of optimal solution of the total cost functions TCD and TCR are explained and justified.

Theoretical Derivations
Convexity of total cost functions TCD and TCR relating to the decision variables was established using the Hessian matrix, allowing for the unique optimal solution to be derived. Special cases of the optimal values are also deliberated.

Case 1 (EPQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)
The convexity of TCD and uniqueness of optimal solution regarding the decision variables are established to prove the following theorem.
Theorem 1. The total cost function TCD gives minimum value regarding W and Q s simultaneously when (1 − x)D + d > 0, P − (1 − y)D − d > 0, and hence TCD emits a unique optimal solution W * and Q s * .
Proof. TCD is a function of W and Q s , and therefore the first order partial derivatives of TCD with regard to W and Q s are And The second order partial derivatives of TCD are Also And The Hessian matrix of TCD is Hence, the first principal minor It is assumed that,P − D − d > 0. Therefore, Once more the second principal minor Since the first and second principal minors of the Hessian matrix for TCD are positive, the Hessian matrix is positive definite. Hence TCD is a nonnegative, differentiable, and strictly convex function regarding the decision variables W and Q s simultaneously, and TCD is minimum at the unique optimal values W * and Q s * . Solving the necessary conditions  (14) and (15) the following results are obtained. And Therefore, the optimal solution is given by And Then the corresponding optimal production lot size is And the optimal cycle time is Also the maximum on-hand stock is Special Cases (i) If d = 0, then W * = 0.
Therefore, no defective items will be produced in the inventory if the defective rate is zero.
(ii) If x = y = 1 and d = 0, then we get This is the lot size in the standard EPQ model. Therefore, the optimal solution of the total cost function of this model conforms to the standard EPQ model.
(iii) If x = y, then, These are the results of the inventory model for two levels of piecewise constant demand.
(iv) If x = y = 1, then, These are the results of the inventory model of defective items for constant demand.

Case 2 (ERQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)
The convexity of TCR and hence the uniqueness of optimal values of the decision variables are determined to prove the following theorem.
Theorem 2. The total cost function TCR gives a minimum value regarding the decision variables W and Q s and hence a unique solution exists for W * and Q s * .
Proof. TCR is a function of W and Q s , and the relation between TCR and TCD is obtained as From the first order and second order partial derivatives, the following results are obtained from Equation (27).
The Hessian matrix of TCR is And the second principal minor Since the first and second principal minor of Hessian matrix for TCR are positive, the Hessian matrix is positive definite. Hence TCR is a nonnegative, differentiable, and strictly convex function regarding W and Q s concurrently and then TCR is minimum at the unique optimal solution W * and Q s * . The necessary conditions of minimization of TCR are the first order partial derivatives of TCR with respect to W and Q s are equivalent to zero, and hence the following results are obtained. And Therefore, the required economic recycle quantity (ERQ) is And shortage level is The corresponding production lot size of the inventory is And the optimal cycle time is Also, the maximum on-hand inventory is Therefore, no defective items will be produced in the inventory if the defective rate is zero.
(ii) If x = y = 1 and d = 0, then, This is the lot size in the standard EPQ model. Therefore, the optimal outcomes are conformable with the standard EPQ model.
(iii) If x = y then, These are the results of the ERQ model for two levels of piecewise constant demand.
(iv) If x = y = 1, then, These are the results of the ERQ model of defective items for constant demand.

Value of "m"
It is assumed that the arbitrary constant "m" determines the accurate average demand as well as total inventory cost but has no significance on the optimal solution. It is introduced to avoid the complexity of the theoretical results. The assumption is justified by the obtained optimal solutions of TCD and TCR and their special cases. The value of "m" is calculated in the following manner.
Addition of Equations (3) and (6) gives the cycle time Then, Equation (9) implies Putting the optimal values of Q s and Q d of both of the total cost function TCD and TCR, the same value of "m" is obtained as follows:

Numerical Illustrations and Case Study
The inventory model is analyzed considering a piecewise constant demand rate that may vary from one level to another. Certain products, for instance, such as bricks, mustard oil, and shopping mall products, have different demands during their production-run periods than during production-off shortage periods. The practical implications of the proposed model can be observed in a brick production company. For the production time of brick raw materials, some imperfect bricks are made before firing, for instance, during the setting and drying periods. These defective bricks are collected and stored properly in stock before they are sent for recycling during production off-time and added to the inventory for the next production cycle. Furthermore, the demand for bricks during the production off-time and shortage time may increase or decrease from the demand during their production run-time.
In order to study the applicability of the model, numerical illustrations are provided below considering the following static input values of parameters taken from MJB brick production company in Bangladesh where: P = 5000, D = 4500, d = 100, C o = 1000, C p = 50, C h = 10, C R = 50, C r = 5, C s = 3, x = 0.75, and y = 0.5. * + t 2 * + t 3 * + t 4 * = 0.3275 + 0.0388+ 0.1843+ 1.0370 = 1.5876 = T * Hence, the model is consistent. For the values of independent variables in the range 120 < W < 150 and 400 < Q s < 430, Figure 4 depicts the convexity of the total cost function TCD with respect to W and Q s . Cycle time verification for TCD: * + * + * + * = 0.3275 + 0.0388+ 0.1843+ 1.0370 = 1.5876 = * Hence, the model is consistent. For the values of independent variables in the range 120 < W < 150 and 400 < Qs < 430, Figure 4 depicts the convexity of the total cost function TCD with respect to W and Qs. For Qs = 414 and 120 < W < 150, Figure 5 depicts the convexity of the total cost function TCD vs. W. For Q s = 414 and 120 < W < 150, Figure 5 depicts the convexity of the total cost function TCD vs. W. For Qs = 414 and 120 < W < 150, Figure 5 depicts the convexity of the total cost function TCD vs. W. Fixing W = 136 and 250 < Qs < 350, the following Figure 6 shows that TCD is convex with respect to Qs, and the total inventory cost reaches minimum at "Qs = 414". Fixing W = 136 and 250 < Q s < 350, the following Figure 6 shows that TCD is convex with respect to Q s , and the total inventory cost reaches minimum at "Q s = 414". For the values of independent variables in the range 80 < W < 120 and 250 < Qs < 350, the following Figure 7 depicts the convexity of the total cost function TCR for W and Qs. Cycle time verification:t 1 * + t 2 * + t 3 * + t 4 * = 0.2357 + 0.0279 + 0.1326 + 0.7462 = 1.1425 = T * . Hence, the model is consistent.
For the values of independent variables in the range 80 < W < 120 and 250 < Q s < 350, the following Figure 7 depicts the convexity of the total cost function TCR for W and Q s .
For Q s = 300 and 80 < W < 120, the following figure (Figure 8) depicts the convexity of the total cost function with respect to "W".
For W = 98 and 250 < Q s < 350, Figure 9 shows that TCR is convex with respect to Q s and the total cost of the inventory is minimum at "Q s = 300".
For the values of independent variables in the range 80 < W < 120 and 250 < Qs < 350, the following Figure 7 depicts the convexity of the total cost function TCR for W and Qs. For Qs = 300 and 80 < W < 120, the following figure (Figure 8) depicts the convexity of the total cost function with respect to "W". For W = 98 and 250 < Qs < 350, Figure 9 shows that TCR is convex with respect to Qs and the total cost of the inventory is minimum at "Qs = 300".

Sensitivity Analysis
In this section, a sensitivity analysis of the anticipated model of various system parameters is examined. When the sensitivity of a parameter is analyzed, only changes to a For W = 98 and 250 < Qs < 350, Figure 9 shows that TCR is convex with respect to Qs and the total cost of the inventory is minimum at "Qs = 300".

Sensitivity Analysis
In this section, a sensitivity analysis of the anticipated model of various system parameters is examined. When the sensitivity of a parameter is analyzed, only changes to a parameter are considered keeping other parameters as fixed. The sensitivity analyses of the key parameters are given in Tables 2-9.

Sensitivity Analysis
In this section, a sensitivity analysis of the anticipated model of various system parameters is examined. When the sensitivity of a parameter is analyzed, only changes to a parameter are considered keeping other parameters as fixed. The sensitivity analyses of the key parameters are given in Tables 2-9. Table 2. Sensitivity of "x" to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, C o = 1000, C p = 50, C h = 10, C R = 50, C r = 5, C s = 3, y = 1.  Table 3. Sensitivity of "y" to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, C o = 1000, C p = 50, C h = 10, C R = 50, C r = 5, C s = 3, x = 1. Cost-benefit increases with the increase in "d". Cost-benefit decreases with increase in recycle cost. Table 6. Sensitivity of "holding cost" to cost-benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, C o = 1000, C p = 50, C R = 50, C r = 5, C s = 3, x = 1.5, y = 1.5. Cost-benefit decreases with increase in holding cost. Table 7. Sensitivity of "x" to cost-benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, C o = 1000, C p = 50, C h = 10, C R = 50, C r = 5, C s = 3, y = 1.5.

C
Explanations: From Table 3, it is clear that with the rise of "y", the total cost is amplified, while the total order quantity is also increased. In contrast, the total cycle length "T" declines. A noteworthy intensification in all the related costs, e.g., PC, HC, RC, can also be observed. Moreover, the maximum shortage Q s notably rises.

Cost-Benefitof ERQ Model
Cost-benefit is a tool used to gain insights into the decision-making process. In the cost-benefit analysis, we considered the comparative total cost of two inventory systems: the system without recycling and system with recycling of defective items. Thus, Average cost-benefit ACB is Percentage of cost-benefit (PCB) is Therefore, from Equation (38), we see that the cost-benefit increases when the cost of raw materials increases and the cost-benefit decreases when holding cost, recycling cost, or both costs increase.
While the effect of a parameter is analyzed, only changes to one parameter is considered while other parameters remain static.

Managerial Insights
The proposed model helps to control inventories and to determine the economic recycle quantity (ERQ), optimal production lot size, cycle time, total inventory cost, maximum on-hand stock, and maximum backlogs. From the numerical analysis, sensitivity analysis (Tables 2-9), and cost-benefit analysis, the following findings provide managerial insights to inventory managers.
Recycling larger or smaller quantities of defective items than the ERQ will increase the total inventory cost. Therefore, inventory managers must evaluate the ERQ of the item before recycling. When the company decides to recycle, production lot sizes should be set up in favor of the ERQ, and this strategy will provide greater facility to the manager to decline the total cost. If raw material costs increase significantly, recycling will benefit the company because it will diminish the load on raw-materials and use the wastage items. The total cost is more sensitive to demand increases during shortage time than it is to increases during production time with positive stock. If the defective rate is high, companies will receive better cost-benefit from recycling because it produces more defective items, which are usually treated as rejected, but the recycling process converts them to useable, which provides huge revenue for the managers.

Conclusions
In this study, we developed an imperfect production model with and without the recycling of defective items. The model is a generalization of models considering recycled raw materials. It is worth mentioning that recycled materials can be an alternative source of raw materials in the event of any uncertainties. Moreover, recycling defective production reduces waste and hence is also environmentally friendly. The study can be applied to bricks, textile, glass, paper, and jute production. It is demonstrated that the total cost function is strictly convex with a unique optimal solution. While this model considers a single item production system, multi-item inventory systems of the same recyclable raw material can be developed for further research. This paper outlines when to stop the production process and when to start it. Proper management of the recycling process is also illustrated in the model with some insights for managers. Furthermore, the model can be extended to consider time-dependent defective rates and selling prices (Mashud et al. 2020a). Considering the environmental emissions (Mashud et al. 2021b;Mishra et al. 2021) would be another interesting extension of this study. As the recycling process is shown in the model, a proper waste management (Tsai et al. 2021) could be another important discussion that is missing in this study. Some marketing strategies, for instance, discount on imperfect items (Mashud et al. 2020a), trade credit financing (Liao et al. 2018(Liao et al. , 2020Srivastava et al. 2018;Mashud et al. 2021aMashud et al. , 2021d can be implemented in the proposed model to make it more lucrative to the practitioners. Other imperfect process in both the EOQ and EPQ models (Lin and Srivastava 2015;Srivastava et al. 2021) can also be consider as an interesting extension. Funding: This work was completed at the authors' own cost. Funding from a third party was not used to carry out this research work.