An Integrated Inventory-Production-Distribution Model for Crisis Relief Supply Chain Optimization: A Systematic Review and Mixed Integer Programming Formulation

ABSTRACT


INTRODUCTION
Nowadays, most organizations' leaders believe that the optimization of industrial and logistics systems is the perfect tool for their success, due to new requirements of budgetary constraints, competitiveness and global performance.It aims at optimizing supply chain activities from suppliers sourcing until distribution to customers, passing through all the intermediate steps of production, inventory and transportation.To shape the framework of our study, we propose Christopher's insightful definition of supply chain as "a network of companies connected upstream and downstream in various processes and activities that generate value for the end user in the form of goods and services" [1].According to the US Council of Supply Chain Management Professionals, supply chain management (SCM) includes all logistics management tasks from the supplier's supplier to the customer's customer, along with the planning and management of sourcing and procurement processes.Through balancing demand and supply, it consists of the alignment and coordination of the functions of procurement, production, assembly, and storage tracking before delivery to the customer.
The capacity to deploy and support humanitarian relief organizations and the armed forces will always be of utmost importance as long as they are dedicated to serving people worldwide whilst resources are limited.To improve overall responsiveness and force readiness, they need to keep integrating the unique logistic capabilities of all their Services in the most cost-effective way possible by cutting back on handling and storage expenses.Hence, crisis relief supply chains refer to the sustainment of either military or humanitarian operations that share almost the same characteristics.This research seeks to deepen previous and present studies before developing a useful global optimization model of crisis supply chains that supports, facilitates and empowers global and humanitarian organizations or forces' interventions, even far from own countries.Since World War II, military logistics planners proposed valuable sequential and functional optimization models at different stages of the defense supply chain, to support them in making effective and efficient supply chain decisions through analytical tools.However, only few authors like Christopher, Kress and Kiley considered in military and humanitarian relief logistics, the optimization of the whole supply chain as a total process, to help leaders in designing and modeling efficient support strategies in hostile environments [1][2][3].Integrating supply chain echelons and functions is a major turning point in decoupling organizations' outputs and benefits within modern economies and industries.
The drive of this study consists of the review of modeling and solving global supply chain optimization problems to meet customers' demand most efficiently, through the coordination and integration of inventory, production and distribution functions.These models are receiving growing attention from decision-makers of modern society, humanitarian and military organizations, as well as of vehicle and other industries, to build efficient and safe supply chain strategies, how vital and valuable economically, socially and environmentally.It is why the key motivation of this work is to review integrated optimization models of forward/reverse crisis relief supply chains within tactical and operational levels, before developing a Closed-loop Three-Echelon Multi-Period Multi-Commodity Inventory-Production-Distribution Problem (IPDP) formulation, to satisfy efficiently non-stationary deterministic demand of each customer in hostile environment.The review analyses cost minimization approaches through the challenging finding of how to design, deploy, and employ in the most efficient way dedicated supply chains that meet demand in different contexts: the optimal mix of holding required inventories, deploying transportation assets with finished and recycled commodities, and scheduling maintenance.After a deep understanding of existing works, the unique contribution of this work is that this review helps us identify and fill the literature gaps with an inclusive decision model much-awaited by many supply chain scientists and decision-makers, taking into account most of the contingent characteristics of recent crisis relief operations in a hostile dynamic environment.No such optimization formulation has been proposed before, handling simultaneously, besides multi-functions integration and constraints, 1) shipments of finished and recovered repaired commodities, 2) hostility attrition, and 3) both loading and safety constraints of heterogeneous commodities clustering and vehicle-commodity compatibility, to meet efficiently the non-stationary demand of each customer.
Our study is subdivided into four parts.After a detailed description of the related literature review methodology and Distribution in Section 2, plus the state of art of humanitarian and military supply chains, Section 3 will discuss and explore the famous existing works regarding integrated inventoryproduction-distribution models and the research gaps findings valuing the interest of our contribution.Section 4 describes the problem and proposes an integrated mathematical formulation through a generalized mixed integer linear programming model.Potential solution methods will be discussed in section 5, with a focus on decomposition, LINGO solvers and metaheuristic approaches.

MATERIALS AND METHODS
Many problems and models on the integration and coordination of supply chain systems have been analyzed in a huge body of literature.The models were categorized into analytical and mathematical techniques that have been created to combine two or more tasks or functions.Others have integrated and coordinated the supply chain through simulation-based techniques.Stadler and Christoph [4] mentioned two methods to increase supply chain competitiveness: either integration of the layers included in organizations and/or coordination of physical, information and financial flows.

Methodology and classification
It is worth recalling that the literature review includes valuable works developed from the pioneering research [5,6].Within this research, papers are selected according to the following characteristics: mathematical programming, centralized planning models, crisis relief supply chain, one function pure optimization problem, Inventory-Production, Inventory-Distribution, Production-Distribution, and Integrated IPD planning decisions, with an emphasis on the tactical, operational, and potential combinations of these levels with strategic decisions, especially mixed-integer models.For recent researchers and professionals addressing real-world problems, this study can present an outline of the State of the Art in mathematical programming models and solution techniques for supply chain inventory, production, and distribution optimization.

Research process
The proposed literature review methodology applies a "Preferred Reporting Items for Systematic Reviews and Meta-Analyses" PRISMA-like reasoning method [7,8].It follows the progressive steps detailed in Figure 1.
This chronology consists of a search and then an analysis fit for eligibility of the existing body of literature to study, before the description of related works.It helps determine research gaps and the scope of future investigation that needs to be done, as a scientific contribution.First, a total of 127 references spanning 47 years were gathered for this review.In these existing scientific contributions, two groups are identified in Figure 2. The starting block includes books and Ph.D. thesis, general and standard professional planning works with different norms.On the other side, literature reviews and surveys, specific articles and conference proceedings present diverse mathematical programming and planning models, integrated supply chain.
In the studied discipline, the emphasis is to design, deploy and employ a logistics system that responds most efficiently, to the demands of each end-customer in deterministic and uncertain environments.This type pertains to the category of IPDP Optimization Problems with Pick-up/Delivery that pushes us to extend thinking to new trends for finding how to satisfy customers' requirements by optimally using necessary resources despite time, capacities, reliability, compatibility and safety constraints.This section presents a literature review of relevant subjects of the ongoing study.It is structured around three aspects: Supply Chain network, characteristics, Integration Modeling and Optimization methods.

State of art of crisis relief supply chains
Crisis Relief military and humanitarian Supply Chain encompasses the determination of requirements, the building up of inventories and capabilities, and the support of transportation and equipment.So, the general problem of optimization is to design a logistics system that responds efficiently to the needs of each customer, like during military and humanitarian operations in a given theater.

Network, processes and decision levels
To guide our project, a multistage logistics sustainment model is proposed in Figure 3 as a conceptual approach of the complex military sustainment supply chain structure.It includes demand, retail, wholesale and acquisition stages, along with the reverse logistics stage, as the recycling pipeline for intermediate or service repair level as a "value recovery" effort [9].Many types of network configurations have been developed in recent works.From a physical flow perspective, SC structures presented in Figure 4 can be classified as divergent, convergent, complex network (combining divergent and convergent structures), or linear.Complex combined network structure integrates forward divergent flows until customers and reverse convergent flows from customers until suppliers [10].
Kress [2] showed that most supply chain networks are organized like military logistics networks, within a complete k-array tree, rooted in the origin, following the three decision levels in Figure 5.
1) Long-term Strategic decisions: national resources and capabilities, locations; 2) Mid-term Tactical decisions: theater-level deployment and employment, planning of best trade-off between load/capacity; 3) Short-term Operational decisions: combat units' logistics, piloting of flows and scheduling.

Crisis response supply chain topology and characteristics
Before-crisis Supply Chain is like a Commercial SC seeking efficiency, cost minimization and best business practices.But during Crisis response time other factors come into play.In fact, a recent analysis of 24 articles in 2024 dedicated to anti-COVID19 management results has shown that the integration of commercial and disaster relief supply chains has become a necessity [11].
For a better comprehension of the characteristics and challenges of military and humanitarian supply chains, a comparison to business supply chains is presented explicitly in Table 1 [2,12].The goal of a CRSC is to sustain operations through the best trade-off between demand and supply functions, following cyclic sequences of processes and events [2] at the tactical and operational levels.On demand side, operational centers communicate their needs to tactical or strategic logistic sources.On supply side, materials move via the network that connects the source or intermediate nodes to tactical units.To be more precise, in addition to storage and distribution of commodities, the service of production maintenance is integrated in the optimization problem.Demand and supply description.In any crisis relief operation, operational and tactical decisions are made depending on theater and scenario-based factors during wartime.Demand is non-stationary and affected by possible attrition.Deployed combat or humanitarian units consume diverse supplies that vary from fundamental products (food and water) relatively stable and proportional to the number of troops generally constant, to weaponry and ammunition highly scenario-dependent with larger variance.Here, logistics factors manuals prescribe consumptions' baselines according to combat day of supply and combat intensity coefficient depending on enemy attrition.Hence, logisticians strive to switch from the sole push strategy to balancing push/pull resupply, and from resource-based to results-based management as follows: (1) Push strategy: delivery/execution is initiated in response to customers' orders and; (2) Pull strategy: delivery/execution is initiated in anticipation of customers' orders.
Besides a moderately steady base demand forecast by historical data, demand variability and unpredictability with sudden surges remain climax characteristics of humanitarian and military supply chains, posing significant challenges to forecasting the exact nature and scale of demand.Figure 6 displays different emergency demand patterns that can follow either sometimes a base demand with a bell shape or a more frequent severity-based demand surge [13][14][15].The main quantitative examples of these complexity factors are: (1) Sudden surge in demand for the vital supplies of medical products, food and shelters imposing safety stocks like in (a); (2) Bottled water can spike from 10,000 units to 90,000 (800% exponential increase) in a few hours after an unexpected earthquake as in (b); (3) Cultural, Religious and Dietary preferences, as for rice over wheat in some regions; (4) Supply chain disruptions due to damaged roads or blocked ports that can affect the timely delivery of goods for vulnerable populations by 5 days, Guiding principles.In crisis relief management, humanitarian supply chain that is a corollary of logistics sustainment of military operations, needs smart supply chains provided that the following principles, mirrored and explained by United States Army Logistics Doctrine [16,17], are assured.These performance drivers consist of Simplicity, Responsiveness, Flexibility, Attainability, Sustainability, Survivability, Economy (support efficiency), and Integration.Regarding the economy, we distinguish four main categories of costs: inventory holding, maintenance, and transportation, in addition to the disruption costs of shortage and/or backorder penalties.In this context, many national defense organizations issued the "Cost Factors Manual" [18,19].These attributes serve as guidelines for critical thinking and careful planning rather than as a checklist.We want to insist here on the modularity aspect of the sustainment nodes.Modularity implies the grouping of inventory, maintenance and transportation functions in one support facility.Hence, a detailed logistics support topology is elaborated in Figure 7 hereafter, to show the Supply Chain Network, modular nodes and flows' inter-relations.

Evolution of supply chain optimization
According to the savvy scientist Dantzig, the field of mathematical programming began to grow quickly in the 1950s, when the Rand Corporation released the first commercial software based on linear programming [21].
Numerous theoretical concepts transformed linear programming from an appealing mathematical paradigm into strong large-scale methods and tools that revolutionized the process of practical planning.Hence, appeared successively the Compact Inverse before the Upper Bounds and Secondary Constraints in 1955, then the Decomposition Principle in 1960.
Next, Benders discovered the dual form and applied it for the first time in 1962 to solve stochastic programs, and is still widely utilized to solve mixed integer programs.Since 1958, Gomory enriched operations research Integer Programming with the systematic generation of cutting planes, extra necessary conditions, to a current set of inequalities that ensures optimization will result in an integer solution.Its contribution provided another precious approach, different from the earlier work by Dantzig and al on the traveling salesman problem.Gomory theory has been progressively elaborated by many researchers and advanced organizations like IBM, regarding ingenious elimination techniques to resolve 0-1 covering issues.It has been discovered that branch-and-bound is a powerful approach for resolving realistic integer programs.However, the most effective methods up to this point are those that combine branch-andbound with cutting planes [22,23].
The majority of useful planning relations were transformed into a system of linear inequalities, and an objective function took the place of the set of guidelines for choosing ideal strategies [24][25][26].The development of the simplex method, however, is credited with turning the relatively simple economy linear programming problem into a fundamental approach for the NP-hard systems applied planning.For some samples, the optimal solution can quickly be obtained thanks to a basic simplex or interior method and a basic computer.These contributions remain exceptional because they make it possible to address dynamic scheduling and planning issues over time, especially in the face of uncertainty.The present study is interested in methods for continuous and combinatorial optimization as well as dynamic programming.These tools play an important role in solving planning, scheduling, and inventory control problems, and cover most fields of modern optimization: (1) We encountered unconstrained non-linear optimization methods that converge to a local optimum or a saddle point, including specific approaches like conjugategradient ones.However, many scientists found that meta-heuristics were successful for global optimization regarding evolutionary or genetic algorithms, and the differential evolution method.(2) From a constrained linear optimization perspective, most reviewed works considered the revised simplex method for linear programming, intended also for linear network optimization or assignment problems.(3) On another register, authors keen on constrained nonlinear optimization have used the standard theorems of the alternative as necessary conditions for mathematical programming, and sufficient conditions for convex functions.They highlighted the usefulness of Lagrangian multiplier and penalty methods for efficient algorithms.(4) The analysis of reviewed works shows that combinatorial and mixed-integer optimization along with dynamic programming remain by far the most used and preferred methods by supply chain specialists.
It has revolutionized operations research, especially in logistics applied mathematics through the dissemination of Branch-and-Bound techniques and their variants including Branch-and-Price, Branch-and-Cut-and-Price.They provide promising results and impactful solutions for decision-makers.Some complex NP-hard combinatorial optimization problems have also become tractable thanks to successful metaheuristics among which Tabu search and more recent nested partitions methods are very promising techniques.However, due to the challenging real-world volatility of the 1990s, stochastic programming, despite its difficulty, has become an exciting field of research and applied science that has solved some important long-term planning problems in many countries.Hence, great optimization methods have emerged to solve dynamic multi-sector programs, through combining the Nested Decomposition Principle based on staircase structure, the ranking selection and the use of matching processors.Since then, the application of these methods has proven huge results in the building of economic growth strategies [24,26].The breakthrough advancement is the capacity to specify broad goals and figure out the optimal solution for complex practical decision-making problems.
For a good understanding of global supply chain optimization and uncertainty mitigation, the author gives in Table 2, an interesting outline of different modeling methods each with specific type and characteristics, highlighting the conceptual, simulation and mathematical model.In this last one, each decision-maker or researcher must initially define the problem characteristics, network topology (nodes and flows), relevant decision variables to optimize, problem constraints, and the nature of the data set.These inputs help them classify optimization problems and adopt the most suitable Objective Function and Constraints' formulation, mathematical modeling, and powerful solution method.When discussing the logic of optimization models, many academics notice that while constrained programming (CP) permits coming up with possible outputs, linear programming (LP) and Integer Linear Programming (ILP) help reach optimal solutions.However, many authors believe that these programming approaches have specific limitations when it comes to solving large-scale problems with wide solution spaces and accurately representing some of the stochastic and dynamic aspects of optimization problems.Explicitly, records are identified from databases and registers and then screened.After a thorough selection, Retrieved and eligible reports give the studies and respective scientific contents to include in the review and determine gaps' findings.In general, the reviewed literature is structured according to four main specifications, depending on Problem definition; Objective function; Mathematical modeling; Solution methods (LP, MILP, MINLP, Stochastic Mixed Integer Programming (SMIP), and Mixed Integer Goal Programming (MIGP)); and Results (outputs/variables): transportation amount, vehicles and demand satisfaction, inventory, location/allocation, and production quantity.In our case, we'll focus on integrated IPDP with Pickup/Delivery.Table 3 and Figure 8 detail the historical evolution of literature works investigated during our study.Getting informed of past and recent studies by well-known specialized scientists and professionals, the obtained results of this literature review and gaps findings demonstrate the originality of the contribution of the proposed work to the scientific research community.

Integration: Climax of supply chain revolution
The integration of supply chain echelons and functions is a major turning point in decoupling organizations' outputs and benefits within modern economies and industries.When optimizing a situation for one issue but neglecting the other, some of the costs of the disregarded issue or function might be very high.Decisions that significantly affect one another must be coordinated to guarantee the supply chain operates efficiently and agile [27].Shapiro, among other well-known authors, contributes widely in elevating the value of logistics science in strategic planning and management from an unbiased perspective [28].
In our supply chain model, decisions aim at improving the analysis approach of sustainment supply chains from a sequential functional dimension to more integrated and global management perspective that is still at the beginning stages in scientific literature, especially in military, humanitarian and other crisis management operations.We investigate four types of Supply Chain Integration: (1) Space-time or Geographical Integration; (2) Functional Integration; (3) Inter-organizations Integration: From a domain-based to Inter-domains structure or Vertical Integration [29]; (4) Decisional Integration combining strategic, tactical and operational decisions in one inclusive model.

Space-time integration
It refers to a dynamic geographical integration from local to multi-echelon or worldwide supply chain network and flows, where multi-period planning means finding the optimal solution for the model in a finite or infinite planning horizon, in place of a single period one.Given the scale and complexity of such optimization problems, the interesting reviews of Bhatnagar, Chen and al. address the issues of supply chain coordination and integration under two types [9,10].The first concerns alignment within the same function at several stages of the supply chain, a Multi-Plant Coordination Problem.

Functional integration
This transition from one function-dominated to a flowdominated supply chain, promoting coordination between functions in one or more echelons, is the second type proposed by the same authors, as a General Coordination Problem [9,10].The integration of Inventory, Production and Distribution decisions in one global optimization problem like in our study is a relevant illustration of this more inclusive integrated model.Table 4 illustrates most dedicated integrated supply chain optimization models, studied in literature.Scientists started with network design or one function pure optimization problem, and have progressively investigated Inventory-Distribution, Inventory-Production or Production-Distribution planning.Others investigated sustainable resilient supply chain in emerging economies [6].However, only a few works discussed Integrated Inventory-Production-Distribution planning decisions concentrating on the tactical and operational levels, and/or incorporating potential integration with strategic location decisions.
Reviewed works imply that there are two main reasons why integrated optimization of these large-scale supply chain problems is challenging.Throughout the supply chain, various facilities and/or choices may have opposing objectives.A lot of studies show that sequential analysis yields only to locally optimal decisions that could play as a structural constraint to the global supply chain.But, modern successful organizations appeal for global optimization that allows service level satisfaction while minimizing total sustainment cost through integrating and aligning simultaneously inventory, production and distribution processes and decisions with aggregate demand [30].This is the reason behind the study of research works focusing on modeling and resolution of integrated Inventory Production Distribution problems (IPDP).

Closed-loop supply chain optimization
The research project surveys particular planning problems for supply chain optimization including forward, reverse, decoupled, integrated, or coupled types.This section's goals are to organize the relevant literature review and categorize optimization models, upon four general characteristics e.g.problem description, supply chain structure, goals or outputs, and solution methods.This taxonomy paves the way and motivates our detailed description and contribution to an inclusive integration inventory-production-distribution decision model.Depending on planning and scheduling objectives, optimization aims at finding the best order of tasks' fulfillment and resource allocation to maximize outputs or minimize total cost during the shortest time.This approach fits more likely for production processes as presented by Hillier and Lieberman [31].The majority of literature works consider sequential or decoupled optimization regarding inventory, production, distribution and MILP facilities location models, evaluating from one period, single commodity, simple uncapacitated facility location model to more elaborated models.
However, Jayaram et al. [32] and Fleischmann et al. [33] have surveyed specific integrated forward/reverse supply chain problems that include capacitated multi-period, multistage and multi-commodity models through either MILP or MINLP.Hence, to solve integrated supply chain problems, optimization models and algorithms moved progressively from single-function pure inventory problems, production, distribution, network design [34] or vehicle routing problems [35], to more integrated approaches, such as location-routing [36], location-allocation [37,38], inventory routing, production routing or inventory production distribution (IPDP) problems that prove to be more complicated as in [39,40].

Pure inventory optimization model
The optimization of pure inventory problems, or the role of inventory management, serves to identify the best mix of anticipation, cycle, safety, pipeline or/and decoupling stocks.In 2014, Graves Stephen [24] and Chen [41] investigated successively many types of deterministic and stochasticdemand inventory models: economic order quantity (EOQ); single-period-like newsvendor model; and multi-period with periodic review, base stock policy and continuous review, (R,Q) policy.
If not, stock-outs result in lost sales very quickly.Table 5 summarizes the main inventory models depending on the following characteristics: (1) Demand.Constant, deterministic, stochastic.
(4) Capacity.Limited or uncapacitated order/inventory limits.(5) Service.Satisfy all demand, allowed shortages, or meet all demand points but with a certain level.Being the cornerstone of inventory management, stocking control helps ensure at the lowest possible cost, that the proper quantity of stock is available for sustaining the organization's targeted fill rate in the market.In order to ensure that customers are served within lead time throughout the supply chain, decision-makers must establish and oversee specific service levels and logistics strategies.
The evaluation of results has shown that decoupled approaches are rapidly deceiving when inventory, production and distribution costs belong to the same scale.So, these functions must be considered simultaneously within the same model.

Multi-functional optimization
If one pure inventory optimization problem is difficult to solve, it is evident that large-scale integrated closed-loop multi-echelon multi-function planning models, such as Inventory Routing, Production Routing, or Inventory Production Distribution (IPDP) problems, probe to be more complicated.We are interested in integrated optimization approaches that consider simultaneously two or more functions to minimize the total cost of the dynamic global supply chain.
Table 6 summarizes reviewed works regarding respectively single and multi-echelon supply chain optimization.No need to mention that multi-echelon Inventory Production Distribution Problems probe to be more complex to resolve as they imply many more variables and constraints to deal with.Chandra and Fisher [42] figure among the prior investigations valuing the integration of production and distribution decisions in comparison with decoupled optimization of production planning and transportation problems.The findings of computational analysis highlighted that the coupled model offers potential cost savings of 5 to 20%.Fumero and Vercellis [43] also developed a synchronized production-inventory-distribution schedule where they used Lagrangian relaxation to solve the integrated MILP problem.In 1996, Pirkul and Jayaraman [53] modeled a cost minimization problem to simulate an integrated multicommodity, tri-echelon, plant capacity, and warehouse location system.They employed a heuristic approach to solve the issue and Lagrangian relaxation (LR) to determine the lower bound.After three years, Özdamar and Yazgac [59] resolved an integrated production-inventory-routing problem through a hierarchical method.Following collective aggregated decision-making, they improved the planning model to a daily schedule.
From then, many scientists proposed some integrated optimization models for production and distribution planning, to coordinate optimally relevant interrelated supply chain decisions, like capacity management, inventory allocation, and vehicle routing.Forma et al. [60] defined and formulated an integrated SCM problem, including procurement, inventory, production and transportation decisions.It is now used as a foundation for contrasting different functional decomposition techniques.Additionally, it has modeled the potential savings from combining tours linked to the transportation of finished goods and raw materials together, as well as the potential to use both owned and third-party transporters.This strategy proved to be more economic and environmental.Qu et al. [61] developed a decomposition method for an integrated inventory-transportation model with stochastic demand for multi-commodity.They adopted isolated computation for inventory and routing decisions before synchronizing them suitably.
However, many researchers like Federgruen and Tzur [62] introduced effective time-partitioning heuristics to solve multi-commodity production/distribution Supply Chain problems, and decomposed the planning horizon into smaller ones, with possible extension to geographic application.They applied it to a multi-location dynamic lot-sizing problem, elaborating then a heuristic for the one warehouse multiretailer model of a two-echelon distribution network delivering distinct items.Computational results have shown that the partitioning heuristics were efficient, leading to a close-to-optimal solution within 1.5% of optimality, even with up to 150 periods of planning horizon and large numbers of retailers and commodities.An interesting optimization model was developed by Berman and Wang [63] to integrate simultaneously inventory and distribution functions with decisions regarding the supply chain network and location of cross-docking sites.Lately, Toptal and Çetinkaya [64] aligned the supply chain demonstrating through analytical and numerical results, interesting cost enhancement.But, they discovered that this performance comes from the variation in cost criteria.In the framework of integrated supply chain optimization, alignment of inventory and distribution system was conducted in 2007 by Archetti et al [44].They performed a branch and cut algorithm for a MILP and compared the results of two relaxations of the problem.Lei et al. [65] proposed a smart effective two-phase approach to solve an integrated production-inventory-distribution routing problem through cost minimization while satisfying customers' demand.After solving the original problem when routing costs were estimated only by direct shipment costs, they resolved the integrated supply chain problem through the consolidation of the inventory, production and direct shipment decisions, already obtained in the first phase with routing costs and constraints.
An in-depth study is given in part V dedicated to discussing potential solution methods.From then, reverse logistics captivated substantial consideration from researchers who published many literature reviews, surveys and models.Amin and Zhang [66] proposed in 2012, a closed-loop supply chain network configuration based on return-recovery recycling through a MILP model.The model aimed at maximizing profit from the optimal quantity of remanufactured products resent to the secondary market.Fleishman et al. [67] reviewed and categorized quantitative mathematical models for reverse logistics distinguishing three key functions optimization: distribution planning, inventory and production planning.Jamsa [68] and Rubio et al. [69] provided a literature review of interesting works and issues, along with case studies and opportunities of research, regarding reverse logistics planning and network design published since the 1990's.

Forward/reverse supply chain optimization
It is only during the three last decades that many integrated optimization models were meticulously investigated and developed either for reverse or closed-loop supply chain planning and/or network design as illustrated in the practical case presented by Chekoubi et al [39].Jayaraman [54,70].were among the top contributors in this field.In 1999, they provided a MILP decision model for the design of reverse logistics networks with a pull system determined by customers' demand for returned and repaired products [32,71].Fleischmann et al. [33] showed that an integrated approach optimizes simultaneously the forward and reverse flows providing considerable cost savings.In the wake of this work, Salema et al. [72] expanded Fleischmann's model to multiproduct network design under demand uncertainty paving the way for very promising outcomes.

Crisis relief supply chain optimization
Some key papers proposed valuable inventory, location, production and distribution models applying useful mathematical optimization techniques and algorithms.To boot, the premier practical reference book [12] expands the academic understanding of crisis response supply chain through sharing applications lessons learned by professionals in daily logistics and worst case scenarios during recent natural and manmade disasters alleviation operations.In 2022, Zhu et al. validated a collaborative fuzzy-based optimization method that ensures both integrative emergency suppliers' evaluation and a sufficient reserve material safety inventory to avoid stockouts and minimize human and material damages, through uncertainty mitigation by a fuzzy interval control ante-, during and post-phases [13].Habib et al. also offered an exhaustive systematic review of most relevant scientific contributions in the two last decades [15].
If deterministic models are more utilized, some authors like Salmeron and Apte elaborated a two-level stochastic locationdistribution-inventory optimization model but used the exact solution method, branch-and-bound algorithm, for scenariodependent equipment and budget allocation prepositioning before a natural crisis starts so that to satisfy variable demand and minimize the expected number of casualties [73].Noyan et al. [56] and Ozdamar [74] developed disaster relief preparedness and response models addressing successively strategic mitigation decisions of warehouse and inventory prepositions to reduce response time and both tactical fleet size decisions and operational last-mile distribution decisions and casualties' evacuation.Throughout Scenario-based 2-stages Stochastic Programming models, the author maximized the anticipated accessibility score using an adapted branch-andcut decomposition algorithm, applied successively in the recent Turkey earthquake.
However, the advanced contributions in the studies [57,58] greatly inspired us on the design and formulation of our integrated optimization model.They elaborated multi-echelon disaster relief supply chain optimization models and algorithms, through progressive 2-stage deterministic-tostochastic programming approach.The total cost minimization of facility location, inventory holding, distribution and shortages subject to some realistic constraints was expressed as a MILPdeterministic equivalent model, tackled by a heuristic approach based on Lagrangian Relaxation and by vehicle symmetry breaking constraints that value the impact of branch-and-bound technique.
Most recently, Guide Jr. and Van Wassenhove [75] gave an exhaustive summary of the development of closed-loop supply chain design and planning research throughout the following five successive stages: From remanufacturing as a technical problem to estimating the reverse logistics steps, before aligning the reverse supply chain, and closing the loop.After a thorough selection, 88 retrieved and eligible reports give the studies and respective scientific contents to include in the review and determine gaps' findings.

Research gap findings
So as outlined, the great part of the reviewed literature on supply chain optimization planning models has targeted forward and reverse logistics separately and hasn't fully considered the integrated closed-loop configuration.However, few studies have investigated the integrated forward/reverse logistics functions problems of multi-echelon multicommodity CRSC optimization in a multiple time horizon.Our review investigated some of the most important methods and tools, previously used for situation-based modeling and solving supply chain management problems.So, it inspires us to describe, simplify then formulate a real-world problem for possible analysis, before designing efficient algorithms to obtain high-quality solutions for impactful decision-making.
Further significantly to our knowledge, none of the previous research works has tried to find efficient strategies for a multiperiod multi-echelon multi-commodity supply chain with non-stationary demand through an integrated Inventory-Production-Distribution planning model with pickup/delivery.Based on Table 6, the proposed model presents the characteristics detailed in Table 7.
To motivate our future contribution and model development, this literature review inspired us to situate the proposed integrated model among three classical optimization subproblems: Multi-Bin Packing Problem (MBPP) as in the work [76], Inventory Distribution Problem and Resource-Constrained Production Scheduling Problem (RCSP).In that Multi-Stage Multi-Period IPD Problem as described in Figure 9, the presence of non-stationary demand, forward and reverse flows and hostility attrition, along with commodity clustering and commodity-vehicle compatibility constraints, make this optimization model unique and valuable for military and humanitarian logistics planners, as well as electronic and vehicles industries.It seeks to enhance both supply chain efficiency and safety with utmost customers' satisfaction.As far as we know, this problem is the first to address diverse new real-world features of military and humanitarian supply chains has not been dealt with before in literature works.
Hence, Figure 9 shows a simplified schematic representation of the multi-period model's structure: (1) 1 Source of Support Node with consolidation, 2 Intermediate Support Nodes with consolidation and 6 Forward Organizational Support Node with End-User Demand, considered Customer Demand Node.(2) Both Warehouses are modular nodes with two or more functions integration.

MATHEMATICAL MODEL: CIPDND
This study considers an integrated logistics sustainment model that coordinates the in-theater core functions of Supply, Maintenance and Transportation upstream and downstream, taking into consideration most of relevant parameters and constraints linked to the required sustainment performance level of fighting units.So, we will try to answer this question: Given a dynamic three-echelon multi-period military supply chain network,  a set different support edges locations:

Problem description
A three-echelon multi-period logistics sustainment network model is adopted where each node is composed of inventory, maintenance and transportation facilities.This tactical supply chain falls within the category, of Inventory-Production-Distribution on nodes for it aims at satisfying the requirements of a set of customers that can be represented by nodes on a network.It is a dynamic multi-stage multi-period IPDP Problem as detailed in Figures 7 and 9.The situation we look at is a single Source of Supply that delivers logistic support to warehouses, then retailers through the Military/Humanitarian Logistics Sustainment Network.We consider four types of flows: transportation trucks and helicopters, new commodities, and returned commodities repaired by the maintenance workforce.We also consider delivery disruption with backorders (or shortages) cost penalties.The concerned closed-loop divergent supply chain network is an integrated three-echelon network with both forward and reverse flows through one Theater Source Support Node (SSF), two Intermediate Modular Support Nodes (ISF) and three Forward Organic Support Nodes (FOSF) for each ISN.Due to safety concerns, lateral supply between both the ISFs (Warehouses) and the FOSF (Retailers) is not allowed.An ISF is to be assigned at least one Retailer, but each Retailer must be assigned exactly one ISF.

Coding of the problem: Model's features
Capacitated Flows.Flows in the logistics network comprise two types.Commodities: new commodities such as: ammunition, food, water, spares and POL and recycled commodities referring to repaired returned commodities.And Transportation vehicles that carry the supplies: Trucks and Helicopters of heterogeneous classes and two modes.
Distribution.A weight per pallet (W) is associated to each pallet and a payload (PC) and bulk (BC) capacity is associated with transportation assets.Arrows of the bottom of figure show some potential transportation relationships between commodities on one hand and between commodities and means of transportation on the other.For instance, some products are never transported together because of technical and tactical danger (e.g mixing ammunition and oxygen bottle in one 5 Tons Expansible Van Truck).
Production.A service labor capacity is associated with battlefield commodity failure, implying maintenance workforce skill at each location.Maintenance service encompasses three-level maintenance system comprising organic maintenance level, intermediate maintenance level and theater level operations.This multi-level system provides timely repairs and necessary evacuations, ensuing in rapid servicing of equipment and returned products, and quick return of items to customers' units in an operational status, as these organizations are developing new multifunctional maintenance facilities to meet 21st century requirements [77].
Limited resources.They include only most critical and restricted resources in the planning model.They consist of vehicles and technicians' workforce that refer respectively to transportation and production.
Demand.Deterministic, non-stationary (advancing) but must be periodically revised according to environmenttechnical and tactical attrition.It consists of multi-commodity: finished and returned repaired ones.
Balancing Push and Pull Resupply.The dual behavior of demand motivates us to reflect on moving progressively from push strategy to the best trade-off between push and pull resupply for customers' requirements satisfaction.At the forward flow, commodities with valid usage rates require pushed automatic resupply, suitable for getting stocks of common-user items.It is better to use pull strategy requisitions for variable usage rate needs of commodities, drawn through a divergent network.In reverse logistics, flows of defective returned products are picked up and shipped following push principles through a semi-convergent network.A proportional quantity of demand from each customer node in each period satisfied is reversed items with a known disposal fraction.
Contemporary logistic inventiveness seeks to reduce costs, develop pull resupply for maximum efficiency, and grow speed, forecasting and visibility of resupply, even for products like returned parts and main items [16,17].
Multi-period Planning Horizon.A CRSC is normally stratified into time-periods (days, hours) needing different levels of supply as shown in Figure 8.The time-dependent dynamics of sustainment flows are considered in our proposed model with the use of two modes of transportation means: trucks and helicopters of heterogeneous classes.The schedule or sequences of production and distribution runs are assigned to any period.
Facility Capacity.Capacitated facilities except for the Source Support Node which can supply all requirements.Only most critical or limiting resources in the optimization problem will be considered like transportation resources (vehicles), production resources (workforce), limited inventory capacity and safety stock, be it known that facilities to be opened are determined.
Relevant Costs.Inventory holding costs, variable production cost, variable delivery vehicle cost, and imperfect customer service penalty cost: when demand is back-ordered or isn't filled on time.There is a linear relationship between production quantity and resource usage.
And road freight cost is estimated as a function of distance and load.The Cost per driving hour is defined as the hourly cost of operating ground mobility or transportation vehicle that covers fuel costs, consumables like washers and bearings, and the repair of major systems and subsystems, such as engines at Intermediate-level repairable.So, the variable transportation cost is calculated with its hourly cost rate, cruising speed and traveled distance or time.In short, we distinguish three main categories of logistics costs: Planning Decisions to include in the model.allocation quantities, distribution decisions in terms of fleet size (number of vehicles), transportation and demand satisfaction quantities, location decisions for the two Intermediate Support Nodes, production quantities and schedule for the repair of returned commodities, inventory holding quantities Modelling Approach with MILP Programming, supporting both forward supply, recovery and resupply tasks, adapted to military/humanitarian operations' sustainment and different kinds of industries like in vehicle and electronic industries surveyed by Üster et al. [78].

Assumptions: Simplified supply chain problem
To limit the complexity of our optimization model, the present study integrated the most impacting functions of inventory, production and distribution in terms of efficiency and customers' satisfaction, which is supposed to be a largescale problem.From then on, the following assumptions are considered: • Demand satisfaction rate (level) α must be more than the lowest possible limit β as a threshold determined depending on the techno-operational requirements of each period.• Transportation routes have no capacity constraints.One transportation mission is allowed for each vehicle between the support facility and a DP which means one arc in one time period.And travel time between each supplier and respective customer is the same in both forward and reverse directions in each period.• While the location of the Source Supplier remains the same throughout the planning horizon, those of Warehouses and Retailers change each period.But they are known as they are predetermined for each period.Consequently, the unit transportation cost for shipment of commodity k from a Source Supplier to a Warehouse, and between Warehouse and Retailer (forward and reverse) are given accordingly to each period.(Cjnkt Unit transportation cost for shipment of commodity k from intermediate support facility j to demand location n ($/Pallet), in time-period t).
• Loads are standardized (Pallet) when possible, to maximize lift capabilities.Then, we will try to find feasible support plans that improve the efficiency, responsiveness and robustness of the model.

Problem formulation
The optimization problem is to determine the most efficient Closed-Loop Inventory-Production-Distribution Plan to meet Non-stationary Demand (CIPDPND) within a global supply chain.It is a MILPproblem that minimizes the total cost of an integrated forward/reverse supply chain while satisfying customers' demands.We provide an optimization model where we define the objective function, decision variables, and constraints through a MILP formulation.Nomenclature defines the used Sets, Indices and Parameters.Objective function.The objective function minimizes the sum of the inventory holding costs at source and intermediate support facilities, the forward transportation costs from source to intermediate support facilities and from these facilities to the demand locations, and the shortages penalty costs from each demand location to its intermediate support facility, for all commodities; in addition to the reverse recovery costs from each demand location to its intermediate support facility and the production costs at this facility, for the commodity k', over the planning horizon of T periods.Shortage penalty costs are considered linear in the number of shortages from each demand location in each No need to mention that we used a combination of many extensions to the basic IPDP model before arriving at the proposed formulation.

Demand constraints Eqs. (2)-(4):
( ) The demand of each customer must be satisfied more than the desired level.Deterministic demand is aggregated into a set of demand locations towards higher-stage support facilities ISFs.Shortage quantity for each commodity from each demand location has to remain very low so that the determined demand satisfaction level α is assured to be more than the minimum tolerated rate β.
The supply of a commodity at each facility is either held in inventory or routed to a demand point to satisfy demand.It refers to the supply of a commodity in a period with its demand or usage taking into account the shortages that behave like negative inventory.
A terminal constraint on the shortages ( 7) is introduced at the end of the planning horizon, showing that all demand is finally met over the T-period planning horizon thanks to unlimited inventory holding capacity of SSF.In the same way, initial backorder is dropped from the formulation as it can be included in initial demand dnjk1.

Transportation constraints Eqs. (9)-(18):
The set of constraints Eqs. ( 9)-( 12), (17) guarantees the feasibility of forward delivery and reverse recovery flows compared to the available transportation capacity of SSF and ISFs, or vehicle capacity.They notify as well, that the number of transportation vehicles in charge of commodities' distribution belongs to the regional support facility.Besides, commodity-vehicle compatibility constraints are also considered.In fact, 5 Tons Cargo Trucks (Vs1; Vj1) transport all commodities (∀k≠2) except for fuel commodity (k=2) transported exclusively on 6 Tons Tank Trucks (Vj2).
Constraints (13)(14)(15)(16) prescribe that only required number of transportation trucks necessary for commodities' distribution forward to ISFs should be at SSF, and those between each ISFs to its associated demand locations should be at respective ISF.However, constraints (17)(18) assure the feasibility of reverse recovery flows where only commodities of class k=3, main item equipment (e.g.Ambulance or Combat System respectively in Humanitarian or Military Supply chain), can be recovered from demand locations to their associated ISF for repair.The number of 5 Tons cargo trucks used for forward shipment of commodities from ISF to a demand location can recover all defective commodities of class k=3 from this location to the same ISF on the same arc.

Production constraints Eqs. (19)-(22):
Maintenance resource constraints are related to the labor workforce resource consumption due to the quantity of recovered commodities repaired at ISF, that remain lower than ISF maintenance production capacity or t-period time duration.In the same way, the existence of production in a period assures of a production setup.
Only required number of maintainers necessary for repair of recovered commodities of class k=3 should be at each ISF on one hand, and only main item commodities of class k=3 are repaired at ISF, on the other.Quantity of class k3 commodities repaired at ISFj cannot exceed quantity of commodities recovered from its 3 associated demand locations.

Model enhancement
To strengthen the distribution part, along with technooperational risk mitigation, we could add hereafter two interesting inequalities regarding commodities clustering compatibility (25)(26) and transportation suitability constraints (27)(28), grasped from the optimization model by Minic et al. [57].

Commodity-commodity clustering compatibility constraints onboard the vehicle:
-From SSF s to ISF j: -From ISF j to Customer n: Constraints (25)(26) are derived based on real world considerations of safety and risk assessment.They ensure that commodities k and k' cannot be assigned to the same vehicle (in equation left side) unless they are compatible according to the clustering compatibility parameter γvkk, based on safety and regulatory concerns.Example: Chemicals, Flammable commodities and fresh food must be transported separately.

Commodity-class of vehicle suitability compatibility:
-From SSF s to ISF j: -From ISF j to Customer n: By incorporating vehicle suitability constraints, these constraints are typically derived within the realistic feasibility of transportation.They ensure through the binary commodityvehicle suitability parameter hkv that only suitable vehicles are considered for transporting specific commodities.E.g. perishable commodities require Refrigerated trucks, Tanker truck for fuel distribution and Cargo truck for spare parts.
Our real-world multi-period multi-echelon multi-product closed-loop supply chain system is now described, simplified and modeled.We will discuss hereafter what are the most suitable potential solution methods for this integrated problem.Both existing proven exact and heuristic methods, and designing efficient hybrid or new algorithms could be investigated to obtain high-quality solutions for impactful decision-making.

DISCUSSION OF POTENTIAL SOLUTION METHODS AND FINDINGS
Once the problem is identified and its key components defined, we examine technical feasibility and modeling approach to solve the problem which means adapted programming methods of resolution in reasonable time.This problem is a mixed-integer linear model with too many continuous and binary variables, in addition to too many constraints.The size of the optimization problem generates a challenge for applying and maintaining such a decision model.For moderate-scale problems containing small number of binary decision variables, this planning model can be effectively resolved by commercial optimization software.However extended models intensify the magnitude and complexity of the problem.This MILP has a deterministic time-hard aspect which is supposed to be a large-scale problem as it considers different new real-world aspects of military and humanitarian supply chains.
For this problem, several MILP formulations use different variables.So, to facilitate its resolution and find out the optimal or near-optimal solution, we will oversee good reformulation changing the variable space of a problem before evaluating its quality or potential depending on: • the obtained relaxation value: If it allows an improvement to the previous formulation, whenever approaching the integer optimum, the formulation is more likely to perform well in a branching algorithm.
• the efficient possibilities of branch techniques: avoiding the presence of too many symmetric solutions, with the same cost and/or quasi-similar structures.
• the formulation structure: when a program can be almost decoupled into independent sub-programs.Decomposition techniques make it possible to use this structure to construct particular algorithmic principles.
Most performing solutions approaches consist generally of Danzig-Wolfe decomposition and column generation to deal with too many variables, Bender's decomposition for problem with too many constraints [46], or Lagrangian relaxation [79].
Considering the MILP (P) defined with the following Vectors and Matrixes, in practical problems, the coefficients matrix M has often the block-diagonal structure of Figure 10 where the gray areas represent the non-zero coefficients, sometimes after simultaneous permutations of rows and columns.While the bottom block represents the coupling constraints and concerns all the variables of the problem, the other diagonal blocks are smaller and concern only a limited batch of variables where they appear exclusively.Many authors solved IPDP without using this matrix composition, coordinating simultaneously the optimization of inventory, production and delivery, and routing in a single model.But, lots of times, this difficult approach couldn't be resolved or blocked many times, pushing analysts to the decomposition methods that use these multiple blocks to solve it by dividing the resolution into several programs and subprograms.The master problem obtained from the original one, which gives the resolution of the initial problem, generally abides with coupling constraints and ensures a pilot role in the decomposition technique.Besides, satellite subproblems help finalize solving the whole problem and are most often connected with the isolated blocks of the coefficients' matrix M. The entire master program and satellite programs are then integrated into an algorithmic process converging to a lower bound or to the optimal solution of the problem.The Lagrangian duality gives a motivating setting for relaxation and decomposition.Relaxation often offers very good relaxation value, habitually better than continuous relaxation.Furthermore, a dedicated study is frequently required for the implementation of decomposition techniques that pave the way for great success on very large industrial problems.
Many scientists as Adulyasak et al. [46] reviewed existing formulations, decomposition methods and solution algorithms for inventory, production and distribution problems, discussing not only Lagrangian Relaxation, Lower bound, and Branch and Cut techniques but also exact methods, heuristics, stochastic, fuzzy optimization.Highlighting modern research findings, Ioannis [80] offered some interesting methods and resources to model and solve supply chain management problems, which propose adapted efficient and effective algorithms.In this section, this article summarizes and discusses the useful formulations, reformulations, and solution techniques for our integrated supply chain cost minimization problem IPDP, along with algorithmic and computational issues.

Approaches to compute lower bounds
In general, the IPDP is a complex combinatorial optimization problem with too many variables and constraints.Normally, the lower bound of the basic IPDP formulation attained through the resolution of the LP relaxation is not of good quality.Nananukul [47] demonstrated in his doctoral thesis that the LP relaxation is not suitable for offering relaxed solutions in exact algorithms (e.g.branch-and-bound) or in finding further near-solution methods' quality.Therefore, they investigated substitutes for these relaxation methods to handle complexity and get better lower bounds.It concerns Lagrangian Relaxation techniques illustrated by Fisher [25], Fumero and Vercellis [43], then Solyali and Süral [52], or Column Generation approaches by Nananukul [47] and Lübbecke and Desrosiers [81], as detailed in Appendix.

Decoupled versus integrated solution method
As the IPDP embraces the structure of the LSP and Distribution or VRP problems, it is useful to give a summary of the different formulations and promising reformulations for these two problems.Pochet and Wolsey [82] highlighted the weakness of the basic LSP formulation with poor quality lower bounds.Many reformulations have been used to strengthen it, including the shortest path [83] and facility location [84] ones.Hence, both of these formulations possess the property of integrality in the case of single-level incapacitated LSP, and feasible mixed-integer solutions are attained through solving the relaxed Linear Programming (LP) formulation.Ruokokoski et al. confirmed the efficiency of multiple Lot sizing reformulations in the IPDPs.In the Distribution or VRP, many formulations are investigated to find solutions to the problems with diverse features [48].While a classical simple formulation is appropriate for a homogeneous fleet and reduced quantity of vehicles, a developed formulation should be employed to manage a heterogeneous fleet with different numbers of transporters, consumption or costs, including vehicle indexation and multi-commodity flow consideration.This mathematical method is also normally formulated as the master problem in a column generation technique.Like in the distribution problems, some models were included in the IPDP to handle specific issues, such as the use of vehicle indexation in formulation considering a heterogeneous fleet by Lei et al. [49], and the path-based approach in the column generation method [47,50,51].
Regarding the integrated approach, we found out that the contributions of Chandra and Fisher [42] and Lei et al. [65] offered inspiring integrated solution approaches for our problem.Lei et al. [65] considered the integrated Production, Inventory, and Distribution Routing cost minimization Problem (PIDRP), aligning production, inventory, and delivery operations to satisfy customer demand most efficiently.It also includes heterogeneous vehicles with noninstantaneous traveling times and several client demand locations with different inventory capacities.Due to its combinatorial nature, it was hard for them to optimally solve such an integrated problem especially when vehicle routing is included; and used a two-phase solution approach.Stage I solved a mixed-integer program addressing all the constraints in the initial model but restricting distribution to direct shipments between facilities and demand points.The optimal solution they get in stage I is most of the time feasible to the original model.Stage II deepened the opportunity of a delivery consolidation problem, modeled as a capacitated transportation problem with supplementary constraints and addressed using a heuristic method.
This two-phase approach is unique as it gives the huge opportunity of coordinating simultaneously the production, inventory, and transportation functions throughout the whole planning horizon, without aggregating the demand or relaxing the constraints on transportation capacities.It is a virtuous paradigm shift from the other classical decoupled approach that optimizes separately even successively the production lot sizes before the distribution problem.This enabled them to rapidly pinpoint a high-quality suboptimal solution to the original complex problem and correct this suboptimality through the consolidation effort of stage II.The structural and performance assessment of this proposed two-phase approach shows its effectiveness and encourages our interest in this integrating method to enrich our developed model and facilitate its application to real-world supply chain networks.
In their coordinated inventory, production and distribution planning, Chandra and Fisher compared two management approaches to this system.In the first method, the production scheduling and vehicle routing problems are addressed distinctly, and in the next one, the coupled method, they are integrated within a single model.Computational results highlighted that the coupled model offers potential savings of 5 to 20% [42].

Inspiring optimization solvers
Once we have built our closed-loop multi-period multiechelon inventory, production and distribution mathematical model aiming at minimizing the integrated supply chain cost, we will thrive to solve our complex MIP problem through a synergic prescriptive analytics approach.It combines the abovementioned simplifying reformulation technique with powerful commercial optimization software taking advantage of recent and developing computer technology improvements.The recent tremendous improvement of commercial optimization solvers, along with seeking almost the best compatibility between the appropriate solver and the analyzed model can significantly upsurge our capacity to compute satisfactory solutions to various problems with solvers like CPLEX, GUROBI, and IBM Optimization Solutions and Library (OSL) solver.
However, some recent software like General Algebraic Modeling System (GAMS) and LINGO, annually extended and enhanced, offer high-level model development environments that support the analysis and solution of an extensive range of optimization problems with a significant worldwide user base.
In our study, after simplified reformulation, we will use LINGO thanks to its flexibility, modeling language, and especially its treatment capacity, emphasized by both Hillier and Lieberman [85], and Sithole et al. [86], through the following points: • LINGO is very much useful to solve hard optimization Problems, by using qualitative branching and relaxing methodologies, performing better than meta-heuristic algorithms in certain linear, nonlinear and MILP programming problems.• It can be used to verify and compare the results with the traditional and meta-heuristic optimization methods.Hence, the last version of LINGO demonstrated lately to be a comprehensive modeling language and optimizer with high performance in building (including data codification) and solving large-scale Linear, Nonlinear (convex, non-convex / Global), Quadratic, Quadratically Constrained, Second Order Cone, Semi-Definite, Stochastic, and Integer optimization models faster, easier and more efficiently [86,75].

CONCLUSION AND FUTURE RESEARCH
The optimization of the proposed closed-loop multi-period multi-echelon supply chain consists of minimizing the total cost of inventory and back-orders, production and distribution functions from support facilities to all non-stationary demand points over T periods within a determined Satisfaction Performance Level (SPL).The general structure of the considered network can be applied to different kinds of humanitarian and military organizations, as well as to vehicle and electronic industries.Getting informed of 72 past and recent studies by well-known specialized scientists and professionals, the obtained results of this literature review and gaps findings demonstrate the originality of the contribution of the proposed work to the scientific research community.
These integrating models are garnering increasing global attention as contemporary society, institutions, and governments recognize more acutely not only the critical importance of this knowledge for survival but also for its sustainable economic value.The problem formulation of the model led to a challenging MINLP problem, with too many variables and constraints.Some possible solutions are proposed to illustrate such an integrated approach, which can be improved in future research through risk pooling and uncertainty mitigation to face attritions of real-world environments.
Having reviewed and classified related literature works based on problem characteristics, supply chain structure, objective, modeling approach, decision variables and solution methods, the purpose of this study consists also of describing and proposing an optimal push-pull integrated inventoryproduction-distribution model with a discussion of potential solution methods (reformulation, branching and relaxation) and promising LINGO algorithms.It is a prelude to their future application and computation for results validation and outcomes comparison with the standard push supply chain model.This prescriptive analysis introduces a model that empowers decision-makers to simultaneously optimize product and customer allocations, fleet size and inventory holding within a multi-echelon, multi-period, multicommodity forward/reverse supply chain system.The problem formulation of the model led to a challenging MILP problem, with too many variables and constraints.Some possible solution methods are proposed to illustrate such an integrated approach.
A comparative study of LINGO integrated and decoupled approaches' outputs of the model would be an interesting avenue of reflection, that can be extended with predictive analytics, especially demand forecasting like the Mamdanifuzzy logic decision model just recently proposed by three authors of the present article in [87] and last-mile logistics [88].Indeed, we point out stimulating research perspectives that can enhance the proposed supply chain planning model, including sustainability-based risk pooling [89] and fuzzy optimization to satisfy and adapt scenario-based demand for uncertainty mitigation in a real-world environment.By switching from consumption-based deterministic to fuzzy demand considering simultaneously the three relevant indicators of Consumption Severity, Market Sensitivity, and Commodity Importance, the integration of Artificial Intelligence would provide dynamic visibility and efficiency into volatile supply chain operations.

Figure 2 .
Figure 2. Distribution of literature works per type

Table 7 .
Research gap: Defining our problem : Single or Multi ; Det/Sto : Deterministic or Stochastic Demand; Distribution DS/Rout : Direct Shipment ou VRP; P/DC/C : Inventory at Plant, Distribution Center or Customer.

Figure 9 .
Figure 9. Multi-period supply chain network

( 1 )
Subject to the following constraints:• Satisfaction of each customer demand of commodities more than a given sustainment performance level; • No exceeding the capacity of storage, distribution and production of each source and intermediate support facilities.
k, v, n, j, t

Table 1 .
Commercial versus crisis relief supply chains

Table 2 .
Methods of supply chain modeling

Table 4 .
Main integrated supply chain planning models Distribution modes (DS: Direct Shipment; Rout: Routing)

Table 6 .
Reviewed works regarding single-echelon and multi-echelon integrated supply chain optimization • No flow transshipment is authorized among DPs for safety concerns.It is a common practice at the humanitarian relief and military tactical logistics level, as risks' prevention from lives and equipment losses in a hostile environment.• Initial Stock inventory Ifk0 is known.The Supplier Inventory capacity is Unlimited, and by the end of the planning horizon, all demands are fulfilled.• Maintenance resources exist only in ISF or modular Warehouses and production cost estimates are proportional to the necessary time-to-repair regarding the Intermediate level maintenance costs.Stability at each ISF node is supposed to be more than repair time or respective time period.• We assume that disruption (shortage) penalty cost is linear according to the number of shortages in each period.
Xsjkevt =1 if eth vehicle of class v of SSF s is used to ship commodity of class k to ISF j in period t for e=1,…, Lvs, v∈Vs, j∈J, k∈Kv, t∈T; 0 otherwise.Xjnkevt =1 if eth vehicle of class v of ISF j is used to ship commodity of class k to customer n in period t for e =1,.
, Lvj, v ∈ Vj, j ∈ J, n ∈ N, k ∈ Kv, t ∈ T; 0 otherwise.Yjot Number of maintainers o needed for maintenance of commodity k at intermediate support facility j 4.2.2Objective function and constraints