key: cord-277918-a2iisve0 authors: Malik, Asif Iqbal; Sarkar, Biswajit title: Disruption management in a constrained multi-product imperfect production system date: 2020-06-26 journal: J Manuf Syst DOI: 10.1016/j.jmsy.2020.05.015 sha: doc_id: 277918 cord_uid: a2iisve0 Over several decades, production and inventory systems have been widely studied in different aspects, but only a few studies have considered the production disruption problem. In production systems, the production may be disrupted by priorly unknown disturbance and the entire manufacturing plan can be distorted. This research introduces a production-disruption model for a multi-product single-stage production-inventory system. First, a mathematical model for the multi-item production-inventory system is developed to maximize the total profit for a single-disruption recovery-time window. The main objective of the proposed model is to obtain the optimal manufacturing batch size for multi-item in the recovery time window so that the total profit is maximized. To maintain the matter of multi-product, budget and space constraints are used. A genetic algorithm and pattern search techniques are employed to solve this model and all randomly generated test results are compared. Some numerical examples and sensitivity analysis are given to explain the effectiveness and advantages of the proposed model. This proposed model offers a recovery plan for managers and decision-makers to make accurate and effective decisions in real time during the production disruption problems. A supply chain comprises of facilities and entities that are involved in transforming raw materials into completed products and later delivered to end customers through the supplier. Mostly, conventional supply chains are designed for problem-free and smoothly-operating environments. However, in reality, some unexpected events may occur during the production process such as labor strikes, natural disasters, machine breakdowns, raw materials shortage, and transportation problems. These events are almost unavoidable and cause disruptions at different stages of a supply chain system. In real-life supply chain environments, production disruption is one of very familiar disruptions/ interruptions. In literature, production disruption is defined as any form of disturbance during a production process, including power cut, tool failure, machine breakdown, material shortage, or any type of manmade or accidental interruption [39, 40] . In 1995, a disastrous earthquake hit Kobe and destroyed all available transportation links to the site. This earthquake affected the Toyota production setup negatively and a large amount of revenue was lost due to this unavoidable event [50] . Recovering from production interruptions due to these kinds of events without having a proper response results in an additional cost for the organization, which is considerably high [31] . The most recent example of this type of production disruptions is the COVID-19 epidemic and it caused manufacturers like Hyundai and Fiat Chrysler Automobiles NV to halt their productions [16] . Batch production systems are a well-known and very popular technique in advanced manufacturing setups [5, 46] . In these manufacturing setups, both single and multiple products depending on the type of products are produced and delivered in batches. Batch production helps to reduce cost and increase profitability, the production lot size is determined to minimize the cost of the manufacturing system. There are numerous industries that produce single or multiple products using the batch production technique [27, 48, 53] . However, in real-life cases, there are lots of risk factors involved and these should be taken into consideration when a system is analyzed [37] . The system may face production interruptions due to the above mentioned inevitable events or any other type of production system failure. Because of the imperfect production, the process reliability (which is usually less than 100%) is also an important factor in real-life production systems. It is very difficult to obtain a system which produces all perfect products. Production disruption may affect the organization financially and the reputation of the organization in the market may suffer due to the shortage of goods and customer demand that is not fulfilled during this period. To avoid financial and reputation losses, the organization should have a recovery plan for these types of inevitable events [38] . The organization can use such a plan in a time of need and quickly alleviate the disruption effect. Another aspect of production systems is the impact they have on the environment and human lives. This is why production disruption has recently gathered the interest of researchers. Moreover, reliability of the system is highly related to production disruption, and if the system is more reliable, the chances of moving from in-control to out-of-control are rare, resulting in a less imperfect production [19] . Usually, rejected imperfect products are either reworked to obtain perfect quality or disposed. Disposing of rejected products is more dangerous for the environment and should be avoided. One of the most widely studied research topics is the analysis of production-inventory systems. In early studies, researchers developed models under ideal conditions with simple assumptions. Later, a singleitem imperfect production inventory control model with reworking was considered by Sana [43] . Pal et al. [30] incorporated price dependent production rate in an inventory model under imperfect conditions. Sarkar et al. [45] studied a single-stage production-inventory model for rework at the same production setup with random defective rates and backorders. Cárdenas-Barrón and Sana [8] investigated a two-echelon supply chain model for multiple items under promotion dependent demand and delayed payments. Moreover, Cárdenas-Barrón et al. [6] proposed a multi-period supplier selection problem for multiple items in an inventory model. An improved multi-stage imperfect production model was analyzed by Kim and Sarkar [20] under controllable lead time. Further, Cheikhrouhou et al. [9] developed an inventory inspection model for order size and sample size with the return of imperfect products. Recently, Malik and Sarkar [26] studied an inventory control model for a coordination supply chain management to reduce the leadtime by considering two different modes of transportation. Malik and Kim [24] introduced a flexible production model with production rate and economic lot size optimization considering the relation between production rates and carbon emissions. Multi-product production inventory models have been considered by several researchers for previous two decades. A multi-product optimal production quantity model was developed by Pasandideh et al. [32] with imperfect production, permissible rework, and a space constraint. Taleizadeh et al. [51] analyzed a multi-item production model with limited production capacity and service level constraint. Further, Taleizadeh et al. [52] updated the imperfect production model with a random defective rate under limited production capacity. They minimized total cost on optimal values of back ordered quantity, cycle length, and production quantity. Pirayesh and Poormoaied [41] solved a multi-item imperfect production model under limited production capacity by using the genetic algorithm and swarm optimization. Pasandideh et al. [33] developed a multi-item lot size inventory model considering budget and space constraints. Recently, Malik and Sarkar [25] analyzed a multi-product model with uncertain demand, service level and storage space constraints to optimize the order quantity and process quality. Recently, researchers focused on developing different types of disruption recovery models for supply chain and production-inventory systems. They investigated supply chain systems for several types of disruptions such as supply disruption [21, 29] , demand disruption [14, 1, 17] , transportation disruption [57, 13, 34] , machine breakdowns [22, 2, 59] and production disruption [39] . One can find a recent review on rescheduling of production systems in Uhlmann and Frazzon [55] . Studies for optimal recovery strategies in production disruptions are scarce and limited. Xia et al. [58] introduced a disruption management general approach for the two-stage production-inventory system. They incorporated a penalty cost of deviations from the original plan to the new plan and a recovery time window for production disruption. The work of Xia et al. [58] was extended for real-time recovery strategy in a single-stage perfect production-inventory system, by Hishamuddin et al. [12] . The model considered constrained non-linear optimization problem for a known disruption period. Two different solution methodologies, heuristic and evolutionary algorithms were used and compared the results. Paul et al. [37] suggested a real-time production disruption management plan for a two-stage batch production system under the reliability considerations. Initially, they developed a single (independent) production disruption recovery model and then extended it to multiple (dependent) disruptions scenarios. Further, Paul et al. [39] extended model of Hishamuddin et al. [12] to imperfect production system with random production disruptions and a uniform random distribution was followed by disruption occurrences. They introduced a mathematical model that deals with both dependent and independent (single and multiple) disruptions on a realtime basis. Moreover, a three-stage production-inventory model with single and multiple disruptions was proposed by Paul et al. [38] . They considered perfect production and presented a heuristic solution and a standard search algorithm to solve the model. Paul et al. [40] proposed a three-stage supply chain model for mitigation of production disruptions. They developed three different models for different scenarios, an ideal plan with infinite planning horizon, a predictive disruption mitigation plan, and a reactive disruption mitigation plan. The model was solved by a heuristic approach and optimization with LINGO software and compared the results of both. Another important factor that should be considered when developing a production-inventory model is the process reliability, which has a significant impact on system cost and profit. Lulu and Black [23] studied the multi-component manufacturing and assembly environment and analyzed the impact of process unreliability on the performance of the manufacturing production system. They proved that process unreliability results in lower system utilization. Cheng [11] defined the production process reliability as the percentage of non-defective products produced in a system. He considered process reliability for a single-period inventory system and formulated it as a geometric problem. A production model with production process reliability was developed by Bag et al. [3] . They considered the set-up cost, production period, and process reliability as the decision variables for that model. Sarkar [44] analyzed an economic manufacturing quantity model for an imperfect production process under the effects of inflation while considering reliability as the decision variable. Paul et al. [36] analyzed a production-inventory system using process reliability and uncertain demand. Furthermore, Paul et al. [35] introduced a two-stage batch production-inventory system for real-time disruption management under consideration of the process reliability. They formulated the model for both kinds of disruptions, single and multiple (dependent and independent), and revised the solution with changed parameters. See Table 1 for the contributions from different authors, only the most relevant papers are included. The research contributions regarding the production type, disruption type, reliability, and constraints are given. In previous studies, no one has considered the production-disruption problem for multi-item production-inventory systems. In this study, the production disruption model for a single-stage single-product production-inventory system introduced by Hishamuddin et al. [12] is considered. They considered a perfect production-inventory system facing production disruptions during the production cycle. In their model, they assumed all the produced products are perfect. However, in real life situations production systems also produce some defective products due to system limitations. Production disruptions due to inevitable events such as natural disasters or machine breakdowns cannot be predicted in advance. However, when the production system experiences disruption, the original plan is revised and updated quickly to come to alleviate the effects of the disruption. Till now in literature, no study is available regarding multi-product imperfect production systems under disruption. This research considered a single-stage multi-product imperfect production-inventory system that faces random production disruptions at different stages of the production cycle. In this paper, the authors first developed a multi-constrained non-linear mathematical model to deal with the production disruption in an imperfect production system where multiple products are produced. The developed mathematical model is solved by two standard search algorithms, genetic algorithm (GA) and pattern search (PS). Authors have generated number of test problems for the production disruption by using a uniform distribution. Numerical examples are presented to demonstrate the effectiveness of the proposed research and the obtained results from both search techniques, GA and PS, are compared. The production system is faced with eight types of costs: setup cost, production cost, holding cost, inspection cost, imperfect products rejection cost, interest and depreciation cost, backorder cost, and lost sales cost. The research contributions regarding the production type, disruption type, reliability, and constraints are given in Table 1 . Based on the research gap from the literature review and authors' contribution table, the overall contribution of this paper is twofold 1 First, a multi-item single-stage production system is studied under the consideration of unavoidable production disruption. To the best of authors' knowledge, in existing literature, only a very few studies exist with production disruption scenario while none of them has considered multi-item production system. In practical environments, multi-item production systems are getting more famous among manufacturers. Therefore, it is need of time for managers and manufacturers to have multi-item production recovery models with disruptions during production cycles. 2 Second, most of the existing models do not consider budget and space limitations for economic and production decisions while these two factors are the most important for multi-item production-inventory systems. This model considers limitations over both, budget and space, to make the model more practical and real-life problembased. The remainder of the paper is formulated as follows. In Section 2, the problem definition, notation, and assumptions are given. A mathematical model is given in Section 3. Section 4, presents the solution approach, and numerical examples are given in Section 5. Section 6 consists of a sensitivity analysis. Finally, conclusions are given in Section 7. In this section, first the problem definition is given, which indicates the major motivation for this research. After that, notation and assumptions are given for the mathematical model. In manufacturing systems, it is no surprise that disruption problem can emerge at any time within the manufacturing uptime. In this particular section, the disruption problem in a multi-item production system is described. Almost every manufacturing system may face disruption during the production run-time and recover from disruption; it is a difficult and formidable task for the manufacturer. Usually, it is very difficult to devise a recovery plan after the occurrence of the disruption. Therefore, manufacturers should have a recovery plan prior to disruption. When disruption happens, it takes time to get to a point, where the manufacturer can start production again. Thus, the system loses time due to this recovery process. Occasionally the system may even face another disruption. However, if the system faces another disruption after recovering from the previous one, researchers will consider it as a single-disruption case. This is because the second disruption will not be affected by the previous disruption and recovery plan, and the system recovers from the second disruption on the basis of the original recovery plan. Paul et al. [39] developed a single item production disruption model, but this model does not consider any constraint for the budget or space. Most of the studies have been for the single-item production system because of the complexity of the multiitem production system. This paper presents the first attempt at a multiitem production system recovery plan for a disrupted system. This model considers a multi-item production system for the manufacturer, with limits on the available budget and storage space for the inventory. Fig. 1 (a) shows an ideal plan for the single-stage batch productioninventory system where multiple products are produced. These products are produced in batches. After the completion of each batch, there is a production down time which is the summation of idle time and setup. In an ideal production-inventory system, the production quantity for each cycle i is X in,0 (i = 1, 2, .. ., M), and summation of these quantities is equal to Q n . The recovery plan is a updated schedule that includes the revised production quantities in each cycle, and ensuring the maximization of the total profit in the recovery time window. The number of future cycles allocated to return to the original production schedule from the disrupted cycle, is known as the recovery time Table 1 Summary of different authors' contributions. Banerjee [4] Single-item, perfect production − − − Widyadana and Wee [56] Single-stage, single-item, − − − perfect production Sarkar [44] Single-stage, single-item, − Product Production imperfect production reliability capacity Paul et al. [36] Single-stage, single-item, Single and multiple System Production perfect production production reliability capacity disruptions Cárdenas-Barrón and Sana [7] Single-stage, single-item, − − − perfect production Paul et al. [38] Three-stage, single-item, Production System Production perfect production disruption reliability capacity Kim and Sarkar [20] Multi-stage, multi-item, − − Budget imperfect production constraint Saha et al. [42] Single-stage, single-item, − − − perfect production This model Single-stage, multi-item, Production System Production imperfect production disruption reliability capacity, budget, and space constraint "−" This key word is not available in paper. Journal of Manufacturing Systems 56 (2020) 227-240 window [39] , and is decided by the management of the organization. Fig. 1 (b) presents a disrupted production-inventory system and a recovery plan after the occurrence of production disruption. The production system becomes inoperable for a certain period (T dn ) due to a disruption and after then it operates normally. The recovery plan starts just after the ending point of the disruption period and continues during the recovery time window. It takes M number of production cycles during the recovery time window to returns the production process to its original or normal schedule. Decision variablesX in production quantity of product n for cycle i of the recovery time window, i = 1, 2, .. ., MParametersS tn setup time of machine n for each cycle (year)δ in idle time of machine n for cycle i (year)D n demand of product n per year (units/year)H n holding cost for product n per unit per year ($/unit/year)rreliability of the production process, which is known from historical data of the production systemQ n economic lot size for product n per ideal production cycle with process reliability r (units/cycle)A n setup cost of machine n per setup ($/setup)P n production rate for product n (units/year) in a 100% reliable systemMnumber of cycles in the recovery time window to get to the original production plan (given from the management)T dn disruption period (years)q n pre-disruption production quantity of product n (units)T 0n production time (years) for = q ( ) n of each acceptable unit of product n, this must be greater than 1s n space occupied per unit of product n (meter sq./ unit)Wtotal budget ($)Stotal space(meter sq./ unit)Nmaximum number of machines and products, a positive integer This paper is based on the following assumptions: (1) This model is developed for a multi-product single-machine imperfect production system. In this system, multiple products are produced on the same machine; for example, different types of gears or bearings manufacturing. (2) To fulfill the customer's demand, the production rate is always higher than the demand for each product. Therefore, shortages are not allowed in the model. (3) During the production cycle, the system may go under the disruption due to unexpected factors (i.e. labors strikes, material shortages, or power cut off) and production stops. (4) The recovery cycle starts immediately after a production disruption occurs. The recovery cycle is the time in which management update the current production plan. The management would like to get to the original production plan within the reasonable number of production cycles after a production disruption happens. In this time window, the managers utilize the ideal time to produce products to over come the delay in production process caused by the production disruption during the disruption period [12] . (6) This model is for an imperfect production system. However, all the produced products are well inspected before dispatch. In this system, all defective products are separated and rejected [39] . Hence, no rework is done in the system. (7) The production system has a restriction on the available space for inventory storage. The total available storage space is S and the system should not violate this limit (see [28] ). (8) There is an investment limit for the manufacturer for the production cost. The total available budget for production cost is B and the production cost should not exceed this budget limit. In this Section, the disruption problem for a multi-item productioninventory system is described and formulated. In real life, for multiitem production systems, different types of disruptions may occur, and to manage the system, it is necessary to devise a recovery plan. If a system faces another disruption after the recovery time window of the first disruption, it can still be managed as a single-disruption problem because the second disruption is not affected by the previous disruption and recovery plan. This model considers the multiple products singledisruption problem. Using the lot-for-lot (LFL) policy of Banerjee [4] and the single-item batch production of Sarker and Khan [49] , the proposed model obtains the economic lot size (Q n ). The The total cost for the inventory system (by adding the above-given equations) is given by The mathematical expressions for the setup cost, production cost, holding cost, inspection cost, rejection cost, and interest and depreciation cost are needed to calculate for total cost calculation. These costs are considered in this model for the disruption recovery time window. All the costs involved in the total cost of the production system TC per item are derived as follows: Setup cost (SC) The first step for starting production is devising the production/ manufacturing setup. All the fixed costs, which are experienced every time a single item is produced, are consolidated as the setup cost, and all these fixed costs are mostly associated with the physical activities required within the setup to produce the quantity or lot. These physical activities include setting up the machine, work center, or assembly line. (2) When a manufacturer produces a product or a lot of a given quantity, this requires a certain amount of resources per product. The production cost is the unit production cost for each product multiplied by the total number of units produced for each product. Here, the production costs are calculated as the product of production cost per unit Cp n , production rate P n , and the total production time. Thus, it can be found as follows: All produced products are stored in the warehouse at a cost. To calculate the holding cost, first, an average inventory is calculated. The manufacturing system begins producing products and faces production disruption after some production time T 0n , and from the start to the time when disruption occurs, the produced quantity is q n . The system faces disruption for the time period T dn , the setup time is S tn , the production run time for the first cycle is T 1n , and the system holds the predisruption quantity q n for the entire time of the production cycle (see Fig. 1 ). Thus, the total produced quantity is calculated (See Appendix A) as The unit holding cost multiplied by the total inventory in the recovery time window, which is the area under the curve in the recovery time window, is the total inventory holding cost: All produced products are inspected, and a fixed cost is incurred to inspect all products. The product inspection cost is All produced products are inspected and the defective products are separated from the perfect products. Further, all imperfect products are rejected and thrown out of the manufacturing system at a fixed cost. The process reliability is r and the rate of rejection is (1 − r) . The total rejection cost is the rejection cost per unit rejected multiplied by the number of units rejected after inspection: The total cost of interest and depreciation equation is considered to be a general power function (for instance, see [11] ). It is inversely proportional to the setup cost A n and directly proportional to the process reliability r. Hence, the cost of interest and depreciation is The backorder is the quantity that is not fulfilled at the time of order but is delivered at a later date when the required quantity is available. When the demand quantity is not available at the time of order and the customer will not wait for the stock to be replenished, the lost sales cost is incurred. The lost sales cost is the unit lost sales cost multiplied by the lost sales units. Therefore, the lost sales cost (See Appendix C) is The total cost per item (TC) is calculated by combining setup cost, production cost, inventory holding cost, inspection cost, interest and depreciation cost, backorder cost, and lost sales cost, and is given by The revenue generated by the acceptable items, during the recovery time window, is calculated as selling price per unit multiplied by the demand quantity. The above profit is maximized subject to the following budget, space, demand, capacity, delivery and transportation constraints: In real life production systems, the manufacturer has a limit for the available budget. Thus, the manufacturer never exceeds this budget limit in terms of the costs, i.e., The total production cost is less than or equal to the budget limit, i.e.; ) . Every production setup has a limited storage capacity for the inventory, which usually considered as space constraint. The space constraint ensures that the space occupied by all products is less than or equal to the total space available at the production setup as follows: n N n n 1 in (15) Production capacity constraint When the system functions according to the original production plan without facing any disruption, the production quantity is equal to the desired number. The following constraint ensures the cycle production quantity in an ideal system for all the products:: The following constraints set the production in each cycle of the recovery time window as less than or equal to the production quantity for the ideal system. These two constraints ensures that the production quantity in each cycle of the recovery time window is less than the production quantity in each cycle before the occurrence of the production disruption. This is because of the transportation and delivery requirements. The model is developed considering the recovery window production capacity. The production capacity is limited in that the system can never produce more than the capacity during the disruption recovery period. It gives the recovery time window production capacity constraint as The process reliability is always less than or equal to one and is only equal to 1 when the production system is considered fully reliable. However, in real-life systems, process reliability is always less than 1. Thus, The system has reliability r which affects the overall rate of production of the item and must be higher than the total demand. As the reliability of the system is less than 1 (r < 1), the production rate can be such that the combined rP n is higher than the demand: The following equations represent the ideal time non-negativity conditions: The following two equations ensure the delay time non-negativity: Finally, all the decision variables are non-negative: 1, 2, ..., ; 1, 2. .., . in (27) In literature, various analytical techniques have been used to solve small-and medium-sized problems. However, finding the optimal solution in a complex non-linear constrained problems by using the analytical approaches presented in the literature is very difficult or nearly impossible because of heavy computing overheads. Consequently, there is also a need to reduce the time, and enhance the precision and quality of the solutions for the optimization of complex real-life problems. The proposed multiitem recovery model is a complex constrained mathematical model and it cannot be solved with the analytical procedure because of its complexity. However, this mathematical model is solved by using the two global search metaheurisitc techniques, (1) genetic algorithm (GA) and (2) pattern search (PS). These metaheurisitics are not local search-based and usually these are population-based. These techniques were used because of their ability to effectively deal with a large number of parameters in complex optimization problems. To obtain the solution, we coded the model in MATLAB R2015b and solved using the optimization toolbox. We executed this on an Intel Core i5 with a 3.20 GHz CPU and 8GB RAM. A single disruption problem given in this model can be solved via the help of search algorithms. To validate the model results, two search techniques a genetic algorithm (GA) and pattern search (PS), were implemented in MATLAB R2015b. GA is a natural genetic-based general purpose optimization technique that mostly explores the given search space [15] . Pattern search is also a search technique for use within a given search space and is easy to apply to nonlinear optimization problems with constraints. We compare the results obtained by both search techniques at the end of the next section. Major steps of the proposed solution algorithm for real-time based single-disruption model can be found as follows: Solution algorithm 1 Input all data for the ideal system and get Q n by using the Eq. (1). 2 Assign the values X in0 = Q n 3 Put i = 1 for disruption for all products i = 1, 2, 3 4 Input pre-disruption quantity q n and the disruption period (T dn ). 5 Generate initial population for the recovery cycle from the start of the disruption. 6 Solve the proposed model with two, Genetic algorithm and Pattern search, search 7 Update the X in value as a revised lot size from (6) and note down the revised production plan. 8 Stop. Journal of Manufacturing Systems 56 (2020) 227-240 In this section, we inspect results for the single-disruption case. In this study, we define a disruption situation as the coalition of the predisruption situation and the duration of the disruption. We consider these parameters as uniform random variables in this study. To test the recommended model within the intervals, over 50 disruption test problems are generated by varying the parameters. We solve the proposed model by both the GA and PS search techniques and analyze and compare the obtained results. The ranges of the parameters are given in Table 2 and the input data for single disruption problem is given in Table 3 . The parameters used for the GA and PS to execute the model and obtain the solution are as follows Parameters for GA • All other parameters in the optimization tool box were set as the default. To examine the results, most of the input data from [39] were used: From the above data, X in = Q n can be easily calculated: below the results for each example are given where only the best results out of thirty runs are provided. Example 1. In the first numerical example, a disruption problem with the upper and lower bounds is considered. The results for Example 1 are given in Table 4 for GA and for PS. In the second numerical example, a disruption problem is considered with budget and space constraints, as well as all the upper and lower bounds mentioned in the above mathematical model. The results for Example 2 are given in Table 5 for GA and for PS. Example 3. In the third and last numerical example, a disruption problem is considered with all the constraints as well as the upper and lower bounds already mentioned in the above model and the data given above. The results for Example 3 are given in Table 6 for GA and for PS. To judge the consistency, best results of 24 independent runs have been compared for two different search techniques (GA and PS). All the Table 2 Parameter ranges. test problems have been generated randomly within the ranges given in Table 2 . The obtained average percentage of error or deviation for the two search techniques was 0.00118% and it can be considered as negligible. To calculate the average value of percentage deviation, following equation was used The comparison of results for 24 test runs, obtained from GA and PS, is presented in Table 7 . From the comparison of results, one can see that the results are very much consistent for both, GA and PS, search techniques. In this part, sensitivity analysis for Test instance 1 of example 3 are given. A sensitivity analysis is performed for each of the five key parameters T dn , q n , B n , and L n . The pattern search technique was used for the sensitivity analysis. The total effect of the parameter changes −50%, −25%, +25%, and +50% are calculated, and details are given in Table 8 . For characterizing the impact, the sensitivity analysis is performed different parameters. From the sensitivity analysis table, one can see that the total profit decreases with the increase of the disruption time T dn and unit lost sales cost L n . The total profit increases in the reverse case for disruption time and unit lost sales cost. However, profit is increasing with the increase of the pre-disruption quantity q n and unit backorder cost B n and vice versa. Overall profit is more sensitive for the two parameters, the unit lost sales cost L n and the reliability of the system r, and it changes values very quickly with a smaller change in the parametric value. Figs. 2, 3, and 4 show the changes in the total profit with respect to disruption period, pre-disruption quantity, backorder cost, and lost sales cost. For each analysis, only one parameter is changed and the remaining parameters are kept constant. Authors make these analyses for both the search techniques, GA and PS. Fig. 2 shows the changes in the total profit with the disruption period. The total profit decreases with the increase in the disruption period. The profit decreases linearly when the disruption period is higher than 0.0035 and the authors depict it in Fig. 2. Fig. 3 presents the impact of pre-disruption quantity on the profitability of the system. The total profit increases with the increase in pre-disruption quantity and increment are more when the predisruption quantity is closer to the production quantity X in0 for the ideal cycle. Fig. 4shows the increment in lost sales cost has a strong impact on the profit and causes a linear decrease in overall profit for the production system. The trend is the same for both GA and PS. This production disruption model can be applied to the production systems where the production environment is imperfect and producing multiple items at the same setup. More importantly, in multiple item production systems, the manufacturer must consider the proposed or suitable type of storage space limit and budget limitations over different costs and investments. The managers of any production industry are always anxious about breakdowns (especially random breakdown) due to defective production. Any disruption may occur within a production setup and it can cause defective production or production with low quality. Therefore, the management aims to control the breakdown anyhow. To do this, they generally consider corrective maintenance or preventive maintenance or both, but all these maintenance policies are too costly and it is difficult to maintain and based on literature, only disruptions are available in a single-stage production system. Nowadays, many complex systems are made, where disruptions may occur anytime. The outcomes of this model will benefit industries to maximize profit. Managers can predict the amount of disruption affected products; they can predict how much time the production process will be suspended due to production such that they can commit with other players for selling their products. Finally, by the outcomes, the industry can make proper scheduling of throughput of finished products and the arrival rate of raw material for any multi-product production system. Authors assumed the lot size is fixed throughout the planning horizon. However, the lot size may be split to meet the transportation capacities i.e. truck capacity. We also assumed the demand is known and fixed. While in practical situations the demand may fluctuate in different conditions. Therefore, management should take care of these parameters while applying the proposed approach. The numerical problems are used to develop the methods of making recovery plans after the production disruption occurs which can present managers with some examples to solve the disruption problems in real environment. In a real operation of production systems, there are a number of variables and constants. These numerical are surely not to analyze all possible situations, but they can • demonstrate the effectiveness of the methods, • indicate how to make the recovery plans after the production disruption occurs, • present the advantages of the solution methodology developed for the production system, and • provides a comparison of results and time taken to solve the problem by GA and PS search techniques. In production systems, disruption can occur at any time within the production run-time. Without a recovery plan for such disruptions, the organization may face huge financial and reputation losses. Furthermore, imperfect production lines are very common in real-life and greatly impact a company's profit and loss. The aim of this study was to develop a recovery plan for a disrupted single-stage multi-product imperfect production system. A mathematical model with constraints was developed for single-stage multi-product imperfect production systems. Space and budget constraints were considered, based on real life experiences. For implementation in the real-time recovery planning of a disrupted production system, a few examples were provided. These examples were solved using the pattern search (PS) and genetic algorithm (GA) solution approaches. Similar results were obtained by both solution approaches. In addition to this, sensitivity analyses were performed to show the impact on total profit; these were performed while changing one parameter and keeping all others constant. This model only considers single production disruption or multidisruption without the effect of the previous disruption, which is a major limitation of this model. Immediate extension of this model is possible by considering multiple dependent production disruptions. Further, this model can be extended to different directions related to production or supply chain. The first extension of this model can be for the fuzzy or uncertain environment, parameter like demand can be considered as a triangular or trapezoidal fuzzy number. Several researchers considered different parameters as fuzzy for inventory models [47] , and it would be a new direction if fuzzy parameters are combined with production disruption problem. The fuzzy nature of parameters will have a noticeable impact on the system and profit. Secondly, this model can be extended with the probabilistic proportion for defective products and all imperfect products reworked in the same setup or with separate setup, see for reference [18] . Another possible extension can be, multi-stage production with stochastic lead time demand [20] . Few other future research directions include considering the case of multiple shipments in a supply chain [10] with multi-stage imperfect production and random defective rate [54] . The authors declare no conflicts of interest. The backorder cost is calculated as Coordinating a supply chain system for production, pricing and service strategies with disruptions Rescheduling of parallel machines with stochastic processing and setup times A production inventory model with fuzzy random demand and with flexibility and reliability considerations A joint economic-lot-size model for purchaser and vendor Optimal manufacturing batch size with rework in a singlestage production system-a simple derivation A new approach to solve the multi-product multi-period inventory lot sizing with supplier selection problem A production-inventory model for a two-echelon supply chain when demand is dependent on sales teams? initiatives Multi-item eoq inventory model in a two-layer supply chain while demand varies with promotional effort Optimization of sample size and order size in an inventory model with quality inspection and return of defective items Optimisation of multi-stage jit production-pricing decision: centralised and decentralised models and algorithms An economic production quantity model with flexibility and reliability considerations A disruption recovery model for a single stage production-inventory system A recovery model for a two-echelon serial supply chain with consideration of transportation disruption A recovery mechanism for a two echelon supply chain system under supply disruption Constrained optimization via genetic algorithms Predicting the impacts of epidemic outbreaks on global supply chains: A simulation-based analysis on the coronavirus outbreak (covid-19/sars-cov-2) case Turn bad into good: Using transshipment-before-buyback for disruptions of stochastic demand Impact of random defective rate on lot size focusing work-in-process inventory in manufacturing system An improved way to calculate imperfect items during long-run production in an integrated inventory model with backorders Multi-stage cleaner production process with quality improvement and lead time dependent ordering cost Competition and cooperation in a single-retailer twosupplier supply chain with supply disruption Solving single machine scheduling under disruption with discounted costs by quantum-inspired hybrid heuristics Effect of process unreliability on integrated manufacturing/ production systems A constrained production system involving production exibility and carbon emissions Optimizing a multi-product continuous-review inventory model with uncertain demand, quality improvement, setup cost reduction, and variation control in lead time Coordinating supply-chain management under stochastic fuzzy environment and lead-time reduction A sustainable production-inventory, model for a controllable carbon emissions rate under shortages A hybrid genetic and imperialist competitive algorithm for green vendor managed inventory of multi-item multi-constraint eoq model under shortage A multi-echelon supply chain model for reworkable items in multiple-markets with supply disruption Joint pricing and ordering policy for two echelon imperfect production inventory model with two cycles A multi-echelon production-inventory system with supply disruption A parameter-tuned genetic algorithm to solve multi-product economic production quantity model with defective items, rework, and constrained space A multiproduct single machine economic production quantity model for an imperfect production system under warehouse construction cost Managing sudden transportation disruptions in supply chains under delivery delay and quantity loss Development of a production inventory model with uncertainty and reliability considerations A disruption recovery model in a production-inventory system with demand uncertainty and process reliability Real time disruption management for a two-stage batch production-inventory system with reliability considerations A disruption recovery plan in a three-stage productioninventory system Managing disruption in an imperfect production-inventory system A quantitative model for disruption mitigation in a supply chain Gpso-ls algorithm for a multi-item epq model with production capacity restriction Optimal retailer investments in green operations and preservation technology for deteriorating items A production-inventory model in an imperfect production process An inventory model with reliability in an imperfect production process An economic production quantity model with random defective rate, rework process and backorders for a single stage production system Manufacturing setup cost reduction and quality improvement for the distribution free continuous-review inventory model with a service level constraint Periodic review fuzzy inventory model with variable lead time and fuzzy demand How does an industry reduce waste and consumed energy within a multi-stage smart sustainable biofuel production system? An optimal batch size for a production system operating under periodic delivery policy The resilient enterprise: overcoming vulnerability for competitive advantage Multiproduct single-machine production system with stochastic scrapped production rate, partial backordering and service level constraint Economic production quantity model with repair failure and limited capacity A sustainable development framework for a cleaner multi-item multi-stage textile production system with a process improvement initiative Optimal batch quantity in a cleaner multi-stage lean production system with random defective rate Production rescheduling review: Opportunities for industrial integration and practical applications Optimal deteriorating items production inventory models with random machine breakdown and stochastic repair time The impact of transportation disruptions on supply chain performance Real-time disruption management in a two-stage production and inventory system Resilient adaptive control based on renewal particle swarm optimization to improve production system energy efficiency The work is supported by the National Research Foundation of Korea (NRF) grant funded by the Korea Government (MSIT) (NRF-2020R1F1A1064460). Range of data ) 1 5414 6164 6164 6164 6164 1 2 6042 6803 6892 6892 6892 3148441.09 3 6599 7549 7549 7472 7549 1 5139 6164 6164 6164 6164 2 2 5667 6892 6892 6892 6892 3148065.77 3 6174 7549 7549 7549 7549 1 5664 6164 6164 6164 6164 3 2 6217 6892 6892 6892 6892 3147491.57 3 6749 7549 7549 7549 7549 Test Product type 1 5414 6164 6164 6164 6164 1 2 6042 6892 6892 6892 6892 3148460.78 3 6599 7549 7549 7549 7549 1 5139 6164 6164 6164 6164 2 2 5667 6892 6892 6892 6892 3148065.76 3 6174 7549 7549 7549 7549 1 5664 6164 6164 6164 6164 3 2 6217 6892 6892 6892 6892 3147491.55 3 6749 7549 7549 7549