Optimal Manufacturing-Reconditioning Decisions in a Reverse Logistic System under Periodic Mandatory Carbon Regulation

: Due to environmental concerns, ﬁrms are under increasing pressure to comply with legislations and to take up environmental strategies. This leads researchers and ﬁrms to develop new sustainable supply chains, where a new area has emerged for a manufacturing and reconditioning system. The originality of this work consists in simultaneously considering carbon emissions strategies, carbon tax and mandatory emission in a manufacturing-reconditioning system. The proposed system is composed of two parallel machines, a manufacturing stock, a reconditioning stock and a recovery inventory. In order to make the proposed green manufacturing system more realistic, it is assumed that manufactured (new products) and reconditioned products are distinguishable. The quantity of worn products (used products) depends on the sales in the previous periods, and the repair periods of the machines are stochastic and independent. The aim of this work is to determine the optimal capacities of manufacturing and reconditioning stocks that maximize the total proﬁt, as well as the optimal value of worn products under two carbon emissions’ limitations. An evolutionary algorithm is developed, along with an e ﬃ cient improvement method, to ﬁnd the optimal value of decision variables. Ultimately, numerical results are provided to show the impact of the period of carbon limit and the worn products (returned products) on decision variables. we examine the inﬂuence of the period of mandatory carbon regulation (W) on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, and the percentage of the returned products, p *. This study allows a ﬁrm leader to determine the optimal production and stocking planning when the mandatory carbon period changes. In the third, we analyze the inﬂuence of the limited quantity of carbon emission imposed by the government on optimal values S A *, S Ar *, and p *. This study proposes the optimal stocking management when the government imposes a quantity limit of carbon emissions. In the fourth, we analyze the impact of the carbon emission cost imposed by the government on optimal values S A *, S Ar *, and p *. The used input data


Introduction
Throughout the past few decades, massive carbon emissions have caused serious global environmental damage, such as thick haze and worsening greenhouse gas effects. The Intergovernmental Panel on Climate Change (IPCC), which is the international body for the assessment of climate change, has pointed out in its fifth assessment report that it is necessary to curb the global greenhouse gas (GHG) emissions by 40%-70% from the 2010 level before 2050, and to curb the global GHG emissions to the level of near zero by the end of the 21st century [1]. The IPCC has defined a complete method to standardize the computation of GHG emissions at the national level. Today, most countries are monitoring GHG emissions by using IPCC guidelines to perform annual inventories assessing the quantity of six main gases (CO 2 : carbon dioxide, CH 4 : methane, N 2 O: nitrous oxide, HFCs: hydrofluorocarbons, PFCs: perfluorocarbon, SF 6 : sulfur hexafluoride) [2]. It was found that CO 2 makes up the broad majority of GHG emissions-it represents about three-quarters of total GHG emissions. For instance, in the USA, the emissions from the industry occupies 24% of total emissions. Indeed, GHG emissions from industry primarily come from burning fossil fuels for energy, as well as GHG emissions from certain chemical reactions necessary to produce goods from raw materials [3,4]. To achieve the IPCC targets, drives away the investors. Some European countries are periodically applying the mandatory carbon regulation. For example, the French government applies the mandatory carbon policy on transport and manufacturing activities for a defined period when a high pollution rate is revealed. To the best of our knowledge, there is no work in the literature which deals with a periodic mandatory carbon policy. In this paper, we will combine the carbon tax and mandatory carbon regulation. Indeed, we will consider the carbon tax as a permanent policy to reduce the overall emissions, and the mandatory carbon as a periodic regulation when the emissions exceed a tolerated threshold.
One of the processes to curb carbon emissions in the industrial domain is remanufacturing. Remanufacturing, or refurbing, is considered as a great industrial process, within which worn-out products are restored to new products, providing benefits both environmentally and economically, while at the same time enhancing their image as environmentally responsible, since products are reused instead of being discarded, including curbing carbon emissions [19]. It is supposed that remanufactured products are profitable, at 45%-65% of the price of a new product. Furthermore, they are beneficial under emission regulation [20]. In general, the remanufacturing is a series of industrial processes: disassembly of the worn product on parts, restoring of the reusable parts, replacement of the parts if necessary, and reassembly of the parts [21]. Most of the published papers in the literature have assumed that remanufactured products have the same quality as new ones [22][23][24][25]. However, in practice, the quality of new products is mostly higher than for remanufactured ones. In this case, the remanufactured parts are called "reconditioned parts" and are sold with a lower price than new ones. Over the lasts years, the reconditioning production keeps increasing. In recent years, for example, in the automotive sector, reconditioning firms cover 35% of all production firms in the sector. That in the aerospace sector represents 25%. Concerning the sales, for example, in Asia, the sales of reconditioned smartphones represent 23% of the mobile phones market. In France and according to BACK MARKET and REMADE TECHNOLOGY companies, which produce and sell remanufactured products such as smart phones, laptops and home appliances, the market of remanufactured products posted a growth rate of 7% in 2018, unlike that of new ones, which fell by 6.5%. In the literature, few researchers distinguish between new and reconditioned parts. Gaur et al., [26] have conducted a real case of a battery manufacturer based in India. The authors have proposed a supply chain system, which considers that reconditioned batteries have a lower performance specification and more limited warranty relative to the equivalent new product. They have determined an optimal sales and production plan as well as configuration of the proposed supply chain system. Turki and Rezg [27] have determined the optimal storage, manufacturing, and reconditioning planning, while taking into consideration the difference between new and reconditioned parts. Indeed, they considered that the reconditioned and new parts are distinguishable, and the reconditioned parts are sold at lower price in a market different from that existing for new ones. Therefore, in order to conduce a real-world case of a manufacturing-reconditioning system, we will consider two distinguished markets: the first market for selling new parts and the second for reconditioned ones. Indeed, we will propose a reverse supply chain system that recovers the worn parts from the first market, and then reconditions them to then be sold in the second market. Furthermore, to satisfy markets demands, we will consider two machines: the first for producing new products and the second for reconditioning worn parts. Moshtagh and Taleizadeh [28] have considered a manufacturing-reconditioning system with two machines, one for manufacturing and one for reconditioning. The authors have taken into account inventory costs, production and reconditioning costs, ordering cost, and sales revenue. However, the authors have neglected some characteristics of manufacturing systems, such as machine breakdowns and time to repair [29]. In fact, in real-life manufacturing systems, the time to repair is a manufacturing dead time that cannot be neglected as it causes production delays and raises the risk of stock shortage. Therefore, when the time to repair is stochastic, as in practice, certainly, the control of the manufacturing, reconditioning, and parts' stocks becomes much more complicated. In this paper, in order to make our proposed system closer to the real-life case, we will consider that both machines are subject to random breakdowns and repairs. Furthermore, to control manufacturing and reconditioning processes, we will apply hedging point policy [30] that ensures that the number of parts does not exceed a defined stock threshold. Moreover, concerning the return of the worn parts, most works have assumed that the return of worn products is proportional to the demands, while this assumption is not true and causes a suboptimal in production policies. Therefore, in this work, we assume that the quantity of worn products depends on the sales in the previous periods. The configuration of the proposed system is inspired by the real-life firm, which produces new and reconditioned products such as smart phones or laptops. In addition, due to already existing environmental preoccupations and hard economic concurrence, we assume that the proposed system is operating in a competitive environment. Thus, in order to keep the system competitive, we aim to find the optimal stocking strategy and production planning that maximize the total profit. This consists to determine the optimal capacities of manufacturing and reconditioning stocks and the optimal production control, under carbon tax and periodic mandatory regulation, taking into account stochastic machine breakdowns and repairs. To perform the study and the simulation of the proposed system closer to reality, we adapt the discrete flow model [31] that is adopted in several works. The benefit of the discrete flow model is that it is faithful to describe discrete production systems. In addition, the simulation of the discrete flow model is clear and not painful. However, its optimization takes a huge amount of time to search the solutions. The efficient optimization method chosen was inspired from the evolutionary algorithm [32], and we developed an improvement method based on local search to give better solutions.
Compared with the existing works, this paper contributes mainly in two ways: the first is to examine a manufacturing-reconditioning system simultaneously considering the periodic mandatory carbon emission and carbon tax. The second is to develop a new approach to determine the quantity of new and reconditioned products to produce when the government imposes a limit carbon emission in a random period, the optimal levels of the manufactured-reconditioned stock, and the optimal percentage of returned worn products that maximizes the total profit. The present work is important for implementation of a new circular economy and green deal, which is now on top of the agenda for Europe.
After having stated our purpose, this paper is organized as follows: Section 2 introduces the manufacturing-reconditioning system and the mathematical models. Section 3 presents the developed optimization method inspired from the evolutionary algorithm to find the optimal value of decision variables. Section 4 analyzes the computational results. Finally, in Section 5, we conclude our work.

Materials and Methods
In this section, we present and explain the proposed system. In the Table 1 we present parameters and decision variables that are used to formulate the models in this paper: quantity of unmet demand for new products at time t. P Ar (t) quantity of unmet demand for reconditioned products at time t. V A (t) quantity of new products sold at time t. V Ar (t) quantity of reconditioned products sold at time t. R(t) quantity of worn products that are collected and returned from market 1 at time t. p percentage of the returned products (worn products). p* optimal percentage of the returned products. α lifetime of products.
state of the machine M 2 at time t. cv A unit selling price for new products. cv Ar unit selling price for reconditioned products. cu A unit production cost of new products. cu Ar unit production cost of reconditioned products. cs A unit storage cost for new products. cs Ar unit storage cost for reconditioned products. cr unit storage cost for returned product. cp A unit lost sales cost for new products. cp Ar unit lost sales cost for reconditioned products. ct unit carbon emission cost. q l quantity limit of carbon emission. q pm carbon emission quantity of new products. q pr carbon emission quantity of reconditioned products. W length of the mandatory carbon period.
total profit function. F(t)* optimal value of total profit.
The aim of our work is to find the optimal stock capacity for new products, S A *, for reconditioned products, S Ar *, and the percentage of the returned products, p*, that maximizes the total profit, F(t). Figure 1 presents the model system, that is composed of two machines denoted M 1 for manufacturing and M 2 for reconditioning, which are subject to random repairs and failures. The machine M 1 produces new products from raw material, and we supposed that it never starves. This assumption defines the general case of the manufacturing system where the raw materials are usually acquired. The machine produces new products at variable rate u A (t) and they are stored in the stock S A . The variable u A (t) represents the decided amount for producing new parts in the period t. In fact, u A (t) is determined at the period t-∆t according the machine state (θ(t)) and the stock level of new parts s A (t). The production rate u A (t) takes a value between 0 and its maximum U A . The stock S A satisfies the demand of new products d A (t) coming from the Market 1. The sold new products will be used for a period α and then will be collected and returned to the manufacturer, then, after, they are stored in the recovery inventory R, to be reconditioned by the machine M 2 , and finally are stored in S Ar . The production rate of the reconditioning machine (M 2 ) is denoted u Ar (t), which represents the decided amount for producing reconditioned parts in the period t. Indeed, u Ar (t) is determined at the period t-∆t according the machine state (β(t)), the stock level of reconditioned parts s Ar (t) and the stock level of worn products s Ar (t). The production rate u Ar (t) takes a value between 0 and its maximum U Ar . The stock S Ar satisfies the demand of reconditioned products d Ar (t) coming from the Market 2. We supposed that the sold reconditioned products in Market 2 will be destroyed after their end of life. Indeed, in the real case, most reconditioned products are reused one time. In order to study the system under the carbon emission regulations used, we take into consideration the carbon emission emitted by producing the new/reconditioned products. In our model, the carbon emissions represent the quantity of carbon emitted from producing either new or reconditioned products. In other words, we do not consider the carbon footprint, but we directly calculate the carbon quantity generated by the production of new and reconditioned parts.

Description of the System
Appl. Sci. 2020, 10, x FOR PEER REVIEW 6 of 20 Market 2. We supposed that the sold reconditioned products in Market 2 will be destroyed after their end of life. Indeed, in the real case, most reconditioned products are reused one time. In order to study the system under the carbon emission regulations used, we take into consideration the carbon emission emitted by producing the new/reconditioned products. In our model, the carbon emissions represent the quantity of carbon emitted from producing either new or reconditioned products. In other words, we do not consider the carbon footprint, but we directly calculate the carbon quantity generated by the production of new and reconditioned parts.

Mathematical Model
As assumed before, the machine M1 is always supplied from raw materials, the state of this machine is either up or down. The machine states are presented in the Equation (1).
The state of machine M2 is given by in the Equation (2).
We supposed that the machines M1 and M2 are subject to a random failure exponentially distributed. When M1 is up, uA(t) is between 0 and its maximum UA and when it is down, uA(t) = 0. It is the same for the machine M2. The production rates of the machines M1 and M2 are given by the Equations (3) and (4).

Mathematical Model
As assumed before, the machine M 1 is always supplied from raw materials, the state of this machine is either up or down. The machine states are presented in the Equation (1).
The state of machine M 2 is given by in the Equation (2).
We supposed that the machines M 1 and M 2 are subject to a random failure exponentially distributed. When M 1 is up, u A (t) is between 0 and its maximum U A and when it is down, u A (t) = 0. It is the same for the machine M 2 . The production rates of the machines M 1 and M 2 are given by the Equations (3) and (4).
The stocks s A (t) and s Ar (t) at time t are equal to the level of the stock in the previous period plus the number of products produced in the previous period, minus the number of products outgoing from the stock. In order to simplify the calculations, it assumed that ∆t = 1. Thus, the stock levels at time t are given by the Equations (5) and (6).
The demands d A (t)/d Ar (t) are satisfied respectively by the stock level in the previous time s A (t − ∆t)/s Ar (t − ∆t) if stock levels are higher or equal to d A (t)/d Ar (t)), otherwise the number of the sold products is equal to the stock levels at t − ∆t. The numbers of new/reconditioned products sold at time t are expressed by by the Equations (7) and (8).
The numbers of unmet demand P A (t) and P Ar (t) equal the demand at time t minus the number of sold products at time t. We assume that unmet demands generate a lost cost. This assumption is considered in practice. Indeed, when the customer is not satisfied, a lost cost equivalent to the loss of the business relevance is considered. The Equations (9) and (10) show the number of lost demands at time t: The number of worn products that are returned at time t is proportional to the quantity of satisfied demand after their lifetime. The number of used products equals a percentage of the number of the products sold at t − α, where α is the lifetime of the products and p represents the percentage of returned products (0 < p < 1). Thus, when t − α < 0, there is no return. Therefore, the quantities of returned products are given by the Equation (11).
The Equation (12) represents the level of the recovery inventory at time t, which equals the number of the recovery inventory at t − ∆t plus the number of worn products that are returned at time t, minus the number of products reconditioned by the machine at t − ∆t.
For the machine M 1 , the production rate at time t (see Equation (13)) depends on the machine state θ(t) and on the stock level of new products at time t. We assumed that the machine is never starved from raw materials. The production rate follows these constraints:

•
Equal to the maximum rate production when the state of the machine is up and the stock level capacity is lower than its maximum level.

•
Equal to difference between the stock level capacity and its maximum level, when the state of the machine is up and the stock level at time t is between zero and its maximum level.

•
Equal to the demand when the state of the machine is up or the stock level at time t has reached its maximum.
• Is null when the state of the machine is down.
For the reconditioned products, the machine M 2 is supplied by the recovery inventory. Then, when it is empty or very low, the machine starves. The production rate for reconditioned items looks like the previous one, but also depends on the recovery inventory R(t). Indeed, the production rate for reconditioned products at time t is given the Equation (14).
In this section, we describe the methodology used to plan the production for manufactured and reconditioned products, in the period where the government sets a limitation for carbon emission on a random period (periodic mandatory carbon emission). First, we present the time horizon discretization for the studied problem (see Figure 2).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 20 • Equal to the demand when the state of the machine is up or the stock level at time t has reached its maximum.

•
Is null when the state of the machine is down.
For the reconditioned products, the machine M2 is supplied by the recovery inventory. Then, when it is empty or very low, the machine starves. The production rate for reconditioned items looks like the previous one, but also depends on the recovery inventory R(t). Indeed, the production rate for reconditioned products at time t is given the Equation (14).

t S and s t S s t or if t and s t S ands t d t if t and s t
In this section, we describe the methodology used to plan the production for manufactured and reconditioned products, in the period where the government sets a limitation for carbon emission on a random period (periodic mandatory carbon emission). First, we present the time horizon discretization for the studied problem (see Figure 2). The horizon [0, T] is the finite horizon discretized into periods of time Δt, where Δt is the simulation time step. The occurrence of a period under mandatory carbon emission in the system and the length of the period are generated by truncated normal distribution.

•
In order to determine the production quantity for manufacturing and reconditioning products under carbon regulation, we have provided a method to determine the quantity that respects the limited carbon emission and answers to the demand according to the determined production quantity under hedging point policy. We recall qpm and qpr present the carbon emissions of manufacturing and reconditioning respectively, and ql is the limit quantity of carbon emission to not exceed. We have developed two algorithms (please see Algorithms 1 and 2 that are given in Appendix A and B) that build a function system capable to calculate the proportion of the production rate of the nearest manufactured and reconditioned products under the constraint to not exceed the limit of carbon emissions. The steps are as follows: We calculate the proportion of the production for new and reconditioned products: The horizon [0, T] is the finite horizon discretized into periods of time ∆t, where ∆t is the simulation time step. The occurrence of a period under mandatory carbon emission in the system and the length of the period are generated by truncated normal distribution.

•
In order to determine the production quantity for manufacturing and reconditioning products under carbon regulation, we have provided a method to determine the quantity that respects the limited carbon emission and answers to the demand according to the determined production quantity under hedging point policy. We recall q pm and q pr present the carbon emissions of manufacturing and reconditioning respectively, and q l is the limit quantity of carbon emission to not exceed. We have developed two algorithms (please see Algorithms 1 and 2 that are given in Appendices A and B) that build a function system capable to calculate the proportion of the production rate of the nearest manufactured and reconditioned products under the constraint to Appl. Sci. 2020, 10, 3534 9 of 20 not exceed the limit of carbon emissions. The steps are as follows: We calculate the proportion of the production for new and reconditioned products: • We search, with two counters i and j running from 1 to the upper bound (u A (t) and u Ar (t)), to respect the constraint of not exceeding the limit of carbon emission: i.q pm + j.q pr < = q l (17) • We calculate the proportion: Then, we developed a function that gives the nearest values corresponding to the sum value P m2co and P R2CO closest to P SC and P AC . The total profit function is determined by the difference between the total revenue and the total costs. The total revenue includes the total cost from the new products T t=0 V A (t) · cv A , and the total cost for reconditioned products T t=0 V Ar (t) · cv Ar . The total cost includes: the total production costs for new products T T=0 u A (t) · cu A and the total production costs for reconditioned products T T=0 u Ar (t) · cu Ar , the total storage cost for new products T T=0 s A (t) · cs A , the total storage cost for reconditioned products T T=0 s Ar (t) · cs Ar , and the total storage cost for returned products T T=0 s r (t) · cs r . The total return cost for worn products T T=0 R(t) · cr, the total lost sales for new products T T=0 P A (t) · cp A , and the total lost sales for reconditioned products T T=0 P Ar (t) · cp Ar . The total amount of carbon emitted from producing new products over the horizon [0, T] is T T=0 u A (t) · ct · q pm and from reconditioned products is T T=0 u Ar (t) · ct · q pr . We try to maximize the total profit function given by:

Optimization Method
The mathematical model presented above is clearly based on simulation. Machines states are simulated, random failures are also simulated according to an exponential distribution law, and time flow is also simulated. From all of this, values of stock levels, number of new and reconditioned products, lost demands, return products quantities, production rates of both machines, limit of carbon emission constraint, production plan, and total profit function are calculated, according to Equations (5) to (20). It is obviously not an integer linear mathematical model, but we have developed an optimization method to cope with the non-linearity and random variables challenges, and managed to obtain a total profit function maximization.
The total profit function was programmed in the language C++ with the free software Dev-C++. The optimization method was also developed in this language. In order to facilitate the method handling, we have put the decision variables in a vector. Each value of the vector can be between a minimum and a maximum value. The optimization process is illustrated in Figure 3.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 20 an optimization method to cope with the non-linearity and random variables challenges, and managed to obtain a total profit function maximization. The total profit function was programmed in the language C++ with the free software Dev-C++. The optimization method was also developed in this language. In order to facilitate the method handling, we have put the decision variables in a vector. Each value of the vector can be between a minimum and a maximum value. The optimization process is illustrated in Figure 3. The model is stochastic by construction. Therefore, we calculate the average of three simulations to obtain a more precise value of the objective function, in order to be able to have the best surface possible to perform optimization. The simulation time horizon is fixed to 10 7 time units. The optimization method is inspired from evolutionary algorithms and returns the optimal values of the vector which maximizes the total profit function. The quality of the results obtained by this method depends on the number of tests. The larger the number of tests, the better the results. However, the simulation time can also increase too much. Thus, we developed an improvement method based on local search that can give better solutions by testing the function the least amount possible. The steps are the following: Initialization: random individuals are created, with each vector of values respecting the boundaries of each of its parameters.
Evaluation: In addition to the pre-determined constraints on variables, the following constraint is written: SA > SAr. This equation imposes that the storage of new products should be higher than that of reconditioned ones in order to ensure the return of worn products. Consequently, unadapted individuals are rejected.
Selection: We retain the best adapted individuals, those who give satisfying total profit function results.
Mutation: Among the best adapted individuals, we try some neighborhood tests, thanks to a deviation of some values of the vector describing an individual, to see if its characteristics are better adapted.
Evaluation: If a new generated individual meets the constraints, it is evaluated with its total profit value.
Replacement: If the new total profit function is better than the previous, it is retained. In the end, the final solution is represented by the following vector (SA, SAr, p). The second part of the proposed optimization method is composed of local improvement procedures, which are very useful to perform a better optimization [32]. We developed three main procedures. The first one changes the value of a parameter, chosen at random, in one sense or in the other, here again chosen at random (positive or negative). It can go up until the appropriate bound, but the variation is also at random (it may not "hit" the bound). The set of parameters is then inside each variable bounds, by construction, but the constraints may be violated anyway. Then, after a constraint checking, its objective function is evaluated. If it improves the current solution, it is kept, and if not, the counter of tested solutions is incremented. We allow to test up to 100 valid neighborhoods.
The second neighborhood iteratively repeats the procedure described below, until there is no improvement between two successive iterations. Positive and negative alterations are tested on each (S A * ,S Ar * ,p * ) 1   The model is stochastic by construction. Therefore, we calculate the average of three simulations to obtain a more precise value of the objective function, in order to be able to have the best surface possible to perform optimization. The simulation time horizon is fixed to 10 7 time units. The optimization method is inspired from evolutionary algorithms and returns the optimal values of the vector which maximizes the total profit function. The quality of the results obtained by this method depends on the number of tests. The larger the number of tests, the better the results. However, the simulation time can also increase too much. Thus, we developed an improvement method based on local search that can give better solutions by testing the function the least amount possible. The steps are the following: Initialization: random individuals are created, with each vector of values respecting the boundaries of each of its parameters.
Evaluation: In addition to the pre-determined constraints on variables, the following constraint is written: S A > S Ar . This equation imposes that the storage of new products should be higher than that of reconditioned ones in order to ensure the return of worn products. Consequently, unadapted individuals are rejected.
Selection: We retain the best adapted individuals, those who give satisfying total profit function results.
Mutation: Among the best adapted individuals, we try some neighborhood tests, thanks to a deviation of some values of the vector describing an individual, to see if its characteristics are better adapted.
Evaluation: If a new generated individual meets the constraints, it is evaluated with its total profit value.
Replacement: If the new total profit function is better than the previous, it is retained. In the end, the final solution is represented by the following vector (S A , S Ar , p). The second part of the proposed optimization method is composed of local improvement procedures, which are very useful to perform a better optimization [32]. We developed three main procedures. The first one changes the value of a parameter, chosen at random, in one sense or in the other, here again chosen at random (positive or negative). It can go up until the appropriate bound, but the variation is also at random (it may not "hit" the bound). The set of parameters is then inside each variable bounds, by construction, but the constraints may be violated anyway. Then, after a constraint checking, its objective function is evaluated. If it improves the current solution, it is kept, and if not, the counter of tested solutions is incremented. We allow to test up to 100 valid neighborhoods.
The second neighborhood iteratively repeats the procedure described below, until there is no improvement between two successive iterations. Positive and negative alterations are tested on each of the chromosomes of an individual. If any best improves the objective function, it is kept. Since this neighborhood is iterative, it is only launched once.
The third neighborhood we have implemented is also iterative. For each couple of values of the decision variables vector, we test an alteration of each possible direction's combinations. Since the vector we have to optimize is composed of 3 values (S A *, S Ar *, p*), at each time this neighborhood method is called, it evaluates (c 2 3 × 4 = 12) combinations of values around a current vector. A Pseudo-code algorithm of this iterative neighborhood is given in Appendix C (Algorithm 3). The optimization method we have developed combines both the necessity to consider a stochastic function, with the uncertainty and weakness associated, and the need of efficiency with a relatively long computing, time-consuming objective function. As well as the objective function, this method is also stochastic, and its results cannot be considered as optimal. Nevertheless, they allow a quick convergence towards the optimal region by testing an infinitesimal part of the solution space, thanks to the power of the neighborhoods developed. Moreover, in a stochastic solution space, all but convex, and where quasi-optimal solutions have values in a relatively narrow range, the method we have proposed, designed, and developed, optimizes correctly. This allows us first to work on this exhaustive model and to exploit its richness to come to interesting scientific results. This also proves the relevance of this method to treat optimization problems with high computing, time-consuming objective functions.

Numerical Results
This section conducts experimental computations to illustrate the proposed model and explore different results to investigate the impact of interesting parameters on the optimal values of the decision variables. In the beginning of this section, we introduce numerical data, then we present four subsections. In the first, we investigate the impact of the returned products percentage, p, on the optimal stock capacity for new products, S A *, and reconditioned products, S Ar *. This study analyzes the impact of the quantity of returned worn products on the manufacturing, reconditioning, and stocking decisions. In the second, we examine the influence of the period of mandatory carbon regulation (W) on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, and the percentage of the returned products, p*. This study allows a firm leader to determine the optimal production and stocking planning when the mandatory carbon period changes. In the third, we analyze the influence of the limited quantity of carbon emission imposed by the government on optimal values S A *, S Ar *, and p*. This study proposes the optimal stocking management when the government imposes a quantity limit of carbon emissions. In the fourth, we analyze the impact of the carbon emission cost imposed by the government on optimal values S A *, S Ar *, and p*.
The used input data are: • The time simulation, T = 10 7 periods • The life cycle of new product, α = 100 periods The carbon emission quantity for producing new products, q pm = 100 carbon units • The carbon emission quantity for producing reconditioned products, q pr = 10 carbon units • The quantity limit of carbon emission, q l = 1000 carbon units • The unit carbon emission cost, ct = 0.01 monetary unit (for example, dollars or euros) • The maximum production rate, U A = 2500 products/period • The maximum production rate, U Ar = 2300 products/period The demand d A (t) is generated by truncated normal distribution, in which the average = 15 and the standard deviation = 7. For d Ar (t), the average = 7 and the standard deviation = 4. The generation of time to repair and time between failures are exponentially distributed. The occurrence of the period of carbon emission limit and length of the period denoted (W) are generated by truncated normal distribution. For length of the period W, the lower truncated is 40 and the upper truncated is 60, with the average = 50 and standard deviation = 10, and for the occurrence of the period, the upper truncated is 4100 and the lower truncated is 3900, and the standard deviation = 100.

Impact of the p on S A *, S Ar *, and F(T)
In this subsection, the we investigate the impact of the returned products percentage, p, on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, and the total profit, F(T). Thus, we vary the value of p, and by using the optimization algorithm, we determine the corresponding S A *, S Ar *, and F(T). The simulations' results are illustrated in Table 2. In order to improve the visualization of the obtained results in the analyzed part, we add two Figures, Figures 4 and 5, to the presented Table. • The unit storage cost for returned product cr = 0.0003 monetary unit.
The demand dA(t) is generated by truncated normal distribution, in which the average = 15 and the standard deviation = 7. For dAr(t), the average = 7 and the standard deviation = 4. The generation of time to repair and time between failures are exponentially distributed. The occurrence of the period of carbon emission limit and length of the period denoted (W) are generated by truncated normal distribution. For length of the period W, the lower truncated is 40 and the upper truncated is 60, with the average = 50 and standard deviation = 10, and for the occurrence of the period, the upper truncated is 4100 and the lower truncated is 3900, and the standard deviation = 100.

Impact of the p on SA*, SAr*, and F(T)
In this subsection, the we investigate the impact of the returned products percentage, p, on the optimal stock capacity for new products, SA*, reconditioned products, SAr*, and the total profit, F(T). Thus, we vary the value of p, and by using the optimization algorithm, we determine the corresponding SA*, SAr*, and F(T). The simulations' results are illustrated in Table 2. In order to improve the visualization of the obtained results in the analyzed part, we add two Figures, Figures 4  and 5, to the presented Table.   As shown in Table 2, the total profit is, at the first time, proportional to the percentage of returned products. Then, as the percentage of returned products increases, the profit again decreases and leads to losses. Indeed, it is observed that we have an optimum when p is equal to 50%. This is explained by the fact that when p increases, the number of worn products returned, R(t), increases and replenishes the recovery inventory. Thus, when p is below 50%, the value of profit function increases, which is normal, and the quantity of returned products increases and replenishes the recovery inventory, then the reconditioned products are satisfied. Nevertheless, when p it is above 50%, many of the returned products will be stored in the recovery inventory, which leads to overstocking and generates a high storage cost, and thus explains the losses. The optimal manufacturing stock, SA*, remains constant, which means that it does not depend on p (see Figure 4). In fact, the production of new parts is independent from the production of reconditioned parts. The optimal reconditioning stock, SAr*, increases when p increases (see Figure 5), and this is explained by the need of the worn parts for supplying the demand of reconditioned parts in Market 2. If more worn parts are returned to the production system, more reconditioned parts will be available in the stock, SAr, to satisfy the demand, dAr(t). This study allows a firm leader to manage, in an optimal way, the manufacturing, reconditioning, and stocking when the quantity of returned worn products varies. Furthermore, when the leader has the option to set the quantity of the returned worn products, he can improve the system management by proposing the optimal quantity of the returned worn products that maximize the profit. Therefore, in the next studies, we add p to the decision variables.

Impact of the W on SA*, SAr*, p*, and F(T)
In this subsection, we investigate the impact of the length of the mandatory carbon period, W, on the optimal stock capacity for new products, SA*, reconditioned products, SAr*, returned products percentage, p, and total profit, F(T). Thus, we vary the value of W, and by using the optimization algorithm, we determine the corresponding SA*, SAr*, p*, and F(T). The simulations' results are illustrated in Table 3. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 6 to the presented Table.  Table 2, the total profit is, at the first time, proportional to the percentage of returned products. Then, as the percentage of returned products increases, the profit again decreases and leads to losses. Indeed, it is observed that we have an optimum when p is equal to 50%. This is explained by the fact that when p increases, the number of worn products returned, R(t), increases and replenishes the recovery inventory. Thus, when p is below 50%, the value of profit function increases, which is normal, and the quantity of returned products increases and replenishes the recovery inventory, then the reconditioned products are satisfied. Nevertheless, when p it is above 50%, many of the returned products will be stored in the recovery inventory, which leads to overstocking and generates a high storage cost, and thus explains the losses. The optimal manufacturing stock, S A *, remains constant, which means that it does not depend on p (see Figure 4). In fact, the production of new parts is independent from the production of reconditioned parts. The optimal reconditioning stock, S Ar *, increases when p increases (see Figure 5), and this is explained by the need of the worn parts for supplying the demand of reconditioned parts in Market 2. If more worn parts are returned to the production system, more reconditioned parts will be available in the stock, S Ar , to satisfy the demand, d Ar (t). This study allows a firm leader to manage, in an optimal way, the manufacturing, reconditioning, and stocking when the quantity of returned worn products varies. Furthermore, when the leader has the option to set the quantity of the returned worn products, he can improve the system management by proposing the optimal quantity of the returned worn products that maximize the profit. Therefore, in the next studies, we add p to the decision variables.

Impact of the W on S A *, S Ar *, p*, and F(T)
In this subsection, we investigate the impact of the length of the mandatory carbon period, W, on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, returned products percentage, p, and total profit, F(T). Thus, we vary the value of W, and by using the optimization algorithm, we determine the corresponding S A *, S Ar *, p*, and F(T). The simulations' results are illustrated in Table 3. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 6 to the presented Table.  Table 3. Impact of length of the mandatory carbon period W on S A *, S Ar *, p*, and F(T).   As it is shown in Table 3, the larger the mandatory carbon period is, the more the total profit decreases, which is explained by lower production in a large period. Also, the optimal of the length of the period of carbon is given (the bold line). Indeed, when the mandatory carbon period is low, the limitation for producing either for new or reconditioned products is low, thus the production is less limited, and the number of satisfying products increases, and then the lost sales costs decreases (see Figure 6). Consequently, to hold the produced products, the optimal storage capacities for both types of products increase. Thus, the demands will be satisfied more, and of course, the profit increases. This study allows a company manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when a mandatory carbon period is imposed. Indeed, the government imposes a mandatory carbon period when the pollution reaches a level that is considered harmful to humans. In this period, the manager has to respect a fixed limit of carbon emissions. Consequently, the production is limited, and the manager has to find the optimal planning of production, stocking, and returning of worn products.

Impact of the ql on SA*, SAr*, p*, and F(T)
In this subsection, the we investigate the impact of the limited quantity of carbon emission, ql, on the optimal stock capacity for new products, SA*, reconditioned products, SAr*, returned products percentage, p, and the total profit, F(T). Thus, we vary the value of ql, and by using the optimization algorithm, we determine the corresponding SA*, SAr*, p*, and F(T). The simulations' results are As it is shown in Table 3, the larger the mandatory carbon period is, the more the total profit decreases, which is explained by lower production in a large period. Also, the optimal of the length of the period of carbon is given (the bold line). Indeed, when the mandatory carbon period is low, the limitation for producing either for new or reconditioned products is low, thus the production is less limited, and the number of satisfying products increases, and then the lost sales costs decreases (see Figure 6). Consequently, to hold the produced products, the optimal storage capacities for both types of products increase. Thus, the demands will be satisfied more, and of course, the profit increases. This study allows a company manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when a mandatory carbon period is imposed. Indeed, the government imposes a mandatory carbon period when the pollution reaches a level that is considered harmful to humans. In this period, the manager has to respect a fixed limit of carbon emissions. Consequently, the production is limited, and the manager has to find the optimal planning of production, stocking, and returning of worn products.
4.3. Impact of the q l on S A *, S Ar *, p*, and F(T) In this subsection, the we investigate the impact of the limited quantity of carbon emission, q l , on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, returned products percentage, p, and the total profit, F(T). Thus, we vary the value of q l , and by using the optimization algorithm, we determine the corresponding S A *, S Ar *, p*, and F(T). The simulations' results are illustrated in Table 4. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 7 to the presented Table.   Table 4. Impact of the limited quantity of carbon emission, q l , on S A *, S Ar *, p*, and F(T).  Table 4. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 7 to the presented Table.  As it observed, when the limited quantity of carbon increases, the profit slightly increases, but both SA* and SAr* decrease (see Figure 7). In fact, when the limited quantity of carbon emission is high, the production either for the new or reconditioned products are less limited, and then both types of products are more available, thus, both demands dA(t) and dAr(t) are better satisfied. Moreover, when both types of products are available, the stocking decreases as the demands are satisfied. However, when the limited quantity of carbon emission is low, the production either for the new or reconditioned products are very limited. Thus, both types of products are less available and then it shall increase the storing to satisfy the demands. This study allows a firm manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when the government imposes a limited quantity of carbon emission.

Impact of the ct on SA*, SAr*, p*, and F(T)
In this subsection, we investigate the impact of the cost of carbon emission, ct, on the optimal stock capacity for new products, SA*, reconditioned products, SAr*, returned products percentage, p, and the total profit, F(T). Thus, we vary the value of ct, and by using the optimization algorithm, we determine the corresponding SA*, SAr*, p*, and F(T). To study the influence of the cost of carbon As it observed, when the limited quantity of carbon increases, the profit slightly increases, but both S A * and S Ar * decrease (see Figure 7). In fact, when the limited quantity of carbon emission is high, the production either for the new or reconditioned products are less limited, and then both types of products are more available, thus, both demands d A (t) and d Ar (t) are better satisfied. Moreover, when both types of products are available, the stocking decreases as the demands are satisfied. However, when the limited quantity of carbon emission is low, the production either for the new or reconditioned products are very limited. Thus, both types of products are less available and then it shall increase the storing to satisfy the demands. This study allows a firm manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when the government imposes a limited quantity of carbon emission.
4.4. Impact of the ct on S A *, S Ar *, p*, and F(T) In this subsection, we investigate the impact of the cost of carbon emission, ct, on the optimal stock capacity for new products, S A *, reconditioned products, S Ar *, returned products percentage, p, and the total profit, F(T). Thus, we vary the value of ct, and by using the optimization algorithm, we determine the corresponding S A *, S Ar *, p*, and F(T). To study the influence of the cost of carbon emission, we changed the parameters' values for cs A = 0.00002, cs Ar = 0.00002, and cr = 0.0000001 monetary units. The simulations' results are illustrated in Table 5. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 8 to the presented Table.  Table 5. In order to improve the visualization of the obtained results in the analyzed part, we add the Figure 8 to the presented Table.  The first observation in Table 5 is that the total profit obviously decreases when the cost of carbon increases. Indeed, when the cost of carbon emission increases, the system produces less new parts in order to decrease the cost of emissions, and then, optimal stock, SA*, decreases (see Figure 8). On the other hand, in order to decrease the cost of emissions, the system produces more reconditioned parts instead of new ones. That explain the increase of optimal stock, SAr*, when the cost of carbon emission increases (see Figure 8). This study allows a firm manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when the government increases the carbon emission cost.

Conclusions
This paper has proposed a new design of a manufacturing and reconditioning system, taking into account two carbon regulations: carbon tax and periodic mandatory emission. New and reconditioned products are distinguished and sold in different markets. Our work provides, in the framework of this system, a new and quick strategy to quantify the production of parts under the periods of carbon limitation. To make this green manufacturing more realistic, we have considered the failure and repair of the machines, and the amount of returned products (worn products) depends on the previous sales. The demands are stochastic. To simulate the proposed system and to faithfully  Table 5 is that the total profit obviously decreases when the cost of carbon increases. Indeed, when the cost of carbon emission increases, the system produces less new parts in order to decrease the cost of emissions, and then, optimal stock, S A *, decreases (see Figure 8). On the other hand, in order to decrease the cost of emissions, the system produces more reconditioned parts instead of new ones. That explain the increase of optimal stock, S Ar *, when the cost of carbon emission increases (see Figure 8). This study allows a firm manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when the government increases the carbon emission cost.

Conclusions
This paper has proposed a new design of a manufacturing and reconditioning system, taking into account two carbon regulations: carbon tax and periodic mandatory emission. New and reconditioned products are distinguished and sold in different markets. Our work provides, in the framework of this system, a new and quick strategy to quantify the production of parts under the periods of carbon limitation. To make this green manufacturing more realistic, we have considered the failure and repair of the machines, and the amount of returned products (worn products) depends on the previous sales. The demands are stochastic. To simulate the proposed system and to faithfully describe the system behavior, we have developed a mathematical model based on a discrete flow model. Then, an optimization method based on an evolutionary algorithm was developed, to determine the optimal quantity of the returned worn products, and the stock capacity for new and reconditioned products. Four numerical results were obtained by using the optimization algorithm, which help a firm manager to make decisions on the production and storing of new and reconditioned parts. The results have provided optimal strategies of production, storing, and recovery of worn parts in function of the mandatory carbon period length, emission limit, and carbon cost.
In future research, we will consider the mentioned limitations and the carbon emissions exceeded by transportation from stores to markets.

Discussion
This work concerned the consideration of a periodic mandatory carbon policy that is combined with carbon tax. A new approach was developed to determine the quantity of new and reconditioned products to produce when the government imposes a limited carbon emission in a random period, the optimal levels of the manufactured-reconditioned stock, and the optimal percentage of returned worn products that maximizes the total profit. Furthermore, the proposed approach combined different options that make the study more realistic, such as: distinction between new and reconditioned parts, the dependence of worn products on the sales in the previous periods, and the variability of the machines repair periods. An efficient optimization method based on evolutionary algorithms was developed that returns the optimal values of the vector, which maximizes the total profit function. The studies provided in this work allow a company manager to propose optimal decisions on the manufacturing, reconditioning, stocking, and quantity of returned worn products when a mandatory carbon period is imposed or when the government increases the carbon emission cost. This work is not exempt from limitations. First, it is considered that the reconditioning process is assured by one machine, but in practice, usually it consists of several machines that are subject to random failures. Second, only one type of product was considered, whereas in practice, almost all companies are producing multiple products. Therefore, in our future research, we will extend the model by assuming that the manufacturing process as well as the reconditioning process consists of several activities carried out at several workplaces. Also, we will consider a multi-product system of production and reconditioning.
Author Contributions: S.T. is the responsible for conducting this research. S.T. proposed the main paper contribution. S.S. developed the mathematical model. C.S. developed the optimization method and algorithms. N.S. supervised this work and contributed to revise this paper. All authors have read and agreed to the published version of the manuscript. Production plan (u A (t), u Ar (t)) 01: k = 0; 02: If (u A (t)! 0 and u Ar (t)! = 0) 03: P SC = u A (t)/(u A (t) + u Ar (t)) 04: P AC = u Ar (t)/(u A (t) + u Ar (t)) 05: For (i = 1 to u A (t)) 06: For (j = 1 to u Ar (t)) 07: If (i·q pm + j·q pr < = q L ) 08: P m2co = i/(i + j) 09 P R2CO = j/(i + j) 10: h