Next Article in Journal
Optimizing Cross-Dock Terminal Location Selection: A Multi-Step Approach Based on CI-DEA–IDOCRIW–MABAC for Enhanced Supply Chain Efficiency—A Case Study
Previous Article in Journal
Bio-Inspired Watermarking Method for Authentication of Fundus Images in Computer-Aided Diagnosis of Retinopathy
Previous Article in Special Issue
Self-Built or Third-Party Blockchain Traceability Strategy in a Dual-Channel Supply Chain Considering Consumers’ Traceability Awareness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Green Closed-Loop Supply Chain Management with Bundling Strategy, Perishable Products, and Quality Deterioration

by
Golnaz Hooshmand Pakdel
1,
Yong He
1,* and
Sina Hooshmand Pakdel
2
1
School of Economics and Management, Southeast University, Nanjing 211189, China
2
Engineering Faculty, Science & Research Branch, Islamic Azad University, Tehran 1477893855, Iran
*
Author to whom correspondence should be addressed.
Submission received: 31 January 2024 / Revised: 25 February 2024 / Accepted: 28 February 2024 / Published: 29 February 2024
(This article belongs to the Special Issue Mathematical Modelling and Optimization of Service Supply Chain)

Abstract

:
This study presents a four-objective mathematical model to improve closed-loop supply chain (CLSC) management. The aim of this research is to reduce the costs of the entire chain, risk, emission of pollutants, and time to deliver the product to the customer in uncertain demand condition. In this paper, the NSGAII algorithm is used to solve the model. In this algorithm, among the answers of each generation, a number of them are selected using the two-run tournament selection method. In the binary selection method, the answers are randomly selected from the population, and then a comparison is made between these two answers, and whichever is better is finally selected. The selection criteria in NSGA-II are, firstly, the rank, and secondly, the crowding distance related to the answer. Also, the performance of the NSGA-II algorithm on the same model and data has been compared with the MOPSO algorithm. In the proposed algorithm, if it encounters an impossible solution, it exits the local mode and solves the problem in global conditions. The results show that the proposed method strikes a better balance between discovery and efficiency criteria and avoids falling into local optima. Therefore, in addition to its effectiveness in discovering optimal answers, the genetic-based method has high speed and subsequently, high convergence and diversity rates compared to the particle swarm method. Also, compared to previous methods in the green closed-loop supply chain, the proposed method is better than the modified genetic algorithm, reducing the costs of the chain by about 2.38%.

1. Introduction

Supply Chain Management (SCM) represents the advancement of developing, implementing, and monitoring supply chain (SC) processes, professionally utilizing information and technology. SCM includes all operations that begin with the f raw material procurement, warehouses, work-in-process inventory, and finished goods, i.e., at consumption and origin points, and ensures organizational productivity while meeting customer demands and fulfilling and satisfying the customers (Alzoubi et al., 2022) [1]. However, in the past few decades, supply chains and their various stages have faced internal and external operational challenges, which may be in the environment, nature, or society. The increasing complexity and uncertainty of the environments in which supply chains operate double the significant role of adaptive planning and control, to guarantee delivery to end consumers with minimum delays and interruptions, avoiding unnecessary costs, and maintain business continuity. To implement compliance-based management principles, real-time coordination of production planning, inventory control, and delivery schedules will be essential, while system control parameters must be dynamically adjusted to minimize costs, maximize revenue, meet target service levels, or pursue any other measurable goals to account for dynamism, instability, and uncertainty (Rolf et al., 2022) [2]. Uncertainty of customer demand and unpredictable disturbances, such as short product life cycles and global sourcing, are challenges or problems that make supply chain management ineffective, unstable, vulnerable, and turbulent (Boskabadi et al., 2022; Chen et al., 2022; Roh et al., 2022) [3,4,5]. The application of supply chain design models has primarily failed to consider carbon emissions throughout cost minimization processes. However, many studies, including Zhu et al. (2015), Large et al. (2016), Abir et al. (2019), and Arani et al. (2020) [6,7,8,9], are among those recently taking into account eco-friendly production, hence considering carbon emissions and optimizing total costs.
The supply chain network has been known as a critically important aspect of supply chain management, affecting the chain’s efficiency and effectiveness for many years. In a general classification, supply chain management is evaluated through two forward and reverse supply chain (FSC and RSC, respectively) approaches. The former includes a set of activities from the raw material to product conversion processes (Jabbarzadeh et al., 2018) [10], while the latter is defined as the returned products’ collection and recovery within SCM. This FSC and RSC combination leads to CLSC formation (Devika et al., 2014) [11], which includes customers, collection centers, recycling, and destruction of used products. This chain focuses on collecting, inspecting, and sending the returned products from the customer to the relevant centers for recycling and destruction. Determining the location of the construction of recycling and destruction centers, along with operational variables such as the flow of returned materials, is one of the decision variables of these types of networks. The issue of integration in CLSC network design involves simultaneously determining the strategic and operational decisions of the two chain types (Farrokh et al., 2018) [12]. The reverse logistics order starts from the customer to first-class warehouses, then to second-class warehouses, and finally to the factory or other destinations. In this reverse logistics process, secondary warehouses play a role in recycling, classification, inspection, and transportation. First-class warehouses take responsibility for concentration, disposal, and transportation (Liu et al., 2018) [13]. The main goal in the CLSC network design problem is to maximize the profit of the entire supply chain by choosing the optimal number and location of facilities, their capacities, and the flow (direct and return) of products between facilities, while at the same time reducing environmental impact (through recycling or destruction of reciprocating products) and maximizing social benefits (Mohtashami et al., 2020) [14].
CLSC design comprises the location, number, and capacity determination of facilities and material flow through the network, all of which affect the flexibility, efficiency, and performance of the chains significantly (Rezae and Kheirkhah, 2017) [15]. Effectively designed CLSCs can help reduce adverse environmental effects caused by human intervention substantially (Garai et al., 2021) [16]. Through CLSC design, the chain of network processes helps reduce waste through reuse or recycling. In line with environmental priorities, economic savings also attract attention to the CLSC.
Considering the role of CLSC in improving the sustainable development process, this paper uses a supply chain model with four objective functions. Thus, the main contributions of the presented research paper include the following:
  • Presenting a mathematical model of a CLSC encompassing the vectors of destruction and recycling of products under cost, time, pollutant, and risk reduction conditions.
  • Considering the parts related to destruction and recycling under conditions of uncertainty, depending on the quality of the products, within a CLSC. Hence, the quality level of the returned products from the customer’s side is evaluated, and based on the evaluation results, decisions are made regarding recycling or destruction.
This paper has been organized as follows: Section 2 discusses the literature review and research investigating CLSC management. Section 3 discusses the mathematical modeling and the developed method. In Section 4, we present the results obtained by implementing the presented method and compare the results with other methods. Ultimately, Section 5 provides the research conclusions.

2. Literature Review

To date, significant research has been allocated to the investigation of CLSC modeling and optimization. In this section, the most important articles from recent years are examined. Jerbiaa et al. (2018) investigated the CLSC network with multi-component recycling options, formulating deterministic problems with integer programming [17]. Their results showed that the solutions for stochastic problems are stable. Also, when using the stochastic model, the profit increased. Mohtashami et al. (2020) introduced a green CLSC design utilizing a queuing system to alleviate environmental impacts and energy consumption [14]. Their paper focused on designing a green supply chain under forward and reverse logistics and utilizing a queuing system to optimize transportation time and fleet network waiting, finally reducing environmental impacts. Santander et al. (2020) investigated a CLSC network for local and distributed plastic recycling for 3D printing [18]. Their research examined the economic and environmental dimensions of distributed plastic recycling from the perspective of logistics, validating their method. Goodarzian et al. (2020) introduced a novel mixed-integer multi-objective linear programming model in a drug supply chain network [19]. Their problem formulation of the Mixed Integer Linear Programming (MILP) model focused on minimizing economic costs and environmental impact while maximizing social outcomes. Zarbakhshnia et al. (2020) presented a multi-product, multi-stage, multi-period, and multi-objective MILP model for a forward and reverse logistics network problem [20]. Nasr et al. (2021) developed a multi-objective, multi-product, multi-period mathematical model within a sustainable CLSC, locating distribution, collection, recycling, and disposal centers while taking into account risk criteria [21]. The main objective of their study was to minimize total costs, along with reducing negative environmental outcomes and promoting social responsibility to provide more job opportunities. A fuzzy inference system was used to model and determine uncertainty in demand and parameters that are dependent on demand. Jian et al. (2021) introduced a green package SC comprising a manufacturer and a retailer and solved it using the Stackelberg game approach [22]. Their results showed that profit-sharing contracts improved the relationship between supply chain members, ensuring sustainable economic and environmental growth. Tavana et al. (2022) presented an inclusive framework for a sustainable CLSC network using the multi-objective MILP (MOMILP) model [23]. This study considered the design of stable CLSC networks with interconnection, location-inventory-routing, time window, supplier selection, order allocation, and transportation with simultaneous pickup and delivery in uncertain conditions. Cheng et al. (2022) examined optimal procedures regarding a CLSC network under economic constraints and greenhouse gas (GHG) emission control [24]. This paper investigated how related parameters, including carbon quotas, consumers’ low carbon preferences, and recovery rates, affected the network status. Kouchaki Tajani et al. (2022) proposed a two-channel network of sustainable CLSC for rice, taking into account energy resources and consumption taxes [25]. Their paper sought to formulate an MILP model optimizing total costs, the number of pollutants, as well as job opportunities throughout the introduced SC network considering cost, supply, and demand uncertainties. Babaeinesami et al. (2022) developed a CLSC network considering the suppliers, assembly centers, retailers, customers, collection centers, refurbishing centers, disassembly centers, and disposal centers [26]. Their paper focused on designing a distribution network according to customer demands for simultaneous total cost and total CO2 emission minimization. Alinezhad et al. (2022) presented a stable CLSC network under uncertainty according to fuzzy theory [27]. Their network was a multi-period multi-product problem formed utilizing a two-objective MILP model with fuzzy demand and return rates for SC profit and customer satisfaction maximization. Bathaee et al. (2023) developed an SLSC network, including used product collection as well as new product distribution [28]. The designed mathematical model included three objective functions: profit maximization, total risk minimization, and product scarcity. Wang et al. (2021) conducted a study on complex manufacturing planning (MP) tasks, aiming to optimize order fulfillment rates and minimize total costs [29]. To address the challenges of large-scale problems, they introduced a novel interactive multi-objective optimization-based MP system. This system utilizes a two-stage multi-objective optimization algorithm (TSMOA). Leung et al. (2020) addressed challenges faced by existing metaheuristic approaches like MOPSOs in solving problems with more than three objectives by introducing HGLSS, a Hybrid Global Leader Selection Strategy [30]. HGLSS incorporates two leader selection mechanisms for exploration and exploitation, enabling each particle to select its global best leader. Bahrampour et al. (2023) presented a novel nonlinear mathematical programming model using the mixed integer approach for sustainable CLSC design problem formulation [31]. Their article evaluated the CLSC model from three aspects of sustainability: social, environmental, as well as economic impacts. Wu et al. (2023) evaluated the choice of recycling channels in CLSCs, taking into account the retailer’s competitive preferences [32]. The authors particularly considered the retailer’s competitive preferences and made conclusions based on the three-channel structure of CLSC. They also examined the manufacturer’s, the retailer’s, and the third-party recycling channels (M, R, and T channels, respectively). Dey and Giri (2023) elaborated on a CLSC with two-channel waste recycling considering corporate social responsibility [33]. In their article, various game theory models were used to optimize design and reduce the economic costs of the model. In the following, the literature review summary is illustrated in Table 1.

3. Mathematical Modeling

Figure 1 highlights the schematic of the problem’s mathematical model.

3.1. Proposed Model Assumptions

The recent suppositions provided below formed the basis of the research model:
  • The problem represents a multi-objective, single-period, and multi-product model.
  • All customer demands must be met.
  • Material flow can only be established between two different levels of the network, and there is no material flow between facilities in one layer.
  • The capacity of the facility is limited.
  • The suppliers’ location, production, distribution, collection, recycling, annihilation centers, and customers are known.
  • A part of returned products is recycled in the reverse supply chain and returned back into the chain, while part of it is excreted and removed from the network.
  • The packaging products’ sales volume is usually greater than that of individual sold goods due to the lower price of packaging products than the total price of individual items combined.
  • The productive–technical risk considered in the proposed network model includes interruptions caused by equipment failure and a shortage of skilled labor.
  • Considering the uncertainties of demand in the real world, customer demand is uncertain in the model as well.
  • The amount of products returned from customers and the rates of recycling and destruction are considered uncertain. Disposal refers to the return of a product that enters the recycling stage, but destruction means returning the product to the destruction center. In fact, a product that cannot be recycled is transferred to the destruction center instead of the disposal center.

3.2. Notations and Formulations

The research model is introduced using the following related notations:
Sets:
T The entire sales period
s Suppliers set ( s = 1 , , N )
i Products set ( i = 1 , , N )
f Production centers ( f = 1 ,   2 )
μ Customers ( μ = 1 , , N )
d c Distribution-collection centers ( d c = 1 , , N )
d i Annihilation centers ( d i = 1 ,   2 )
r Recycling centers ( r = 1 ,   2 )
k Staff ( k = 1 , , N )
l Equipment set ( l = 1 , , N )
m Vehicle types ( m = 1 , , N )
t Decision cycle, pricing, and packaging phases ( t = 1 , , N )
Parameters:
n i , t Number of type i products produced in period t
c i Inventory cost of type i product per decision period
S i , t Product i inventory in period t
q i Product i logistics cost
f s Purchasing costs for each raw material unit from supplier s
t r m s f Raw material transportation costs from supplier s to production center f utilizing type m   vehicle depending on distance
t r m f d c Product transportation costs from the production center f to the distribution-collection center d c by type m vehicle depending on distance
t r m d c μ Product transportation costs from the distribution–collection center d c to customer μ by type m vehicle depending on distance
t r m μ d c Transportation costs for the returned products from customer μ to the distribution–collection center d c by type m vehicle depending on distance
t r m d c r Transportation costs for the returned products from the distribution–collection center d c to the recycling center r by type m vehicle depending on distance
t r m d c d i Transportation costs for the returned products from the distribution–collection center d c to the annihilation center d i by type m vehicle depending on distance
t r m r f Transportation costs for the recycled products from the recycling center r to the production center f by type m vehicle depending on distance
t r m r d i Transportation costs for the returned products from the recycling center r to the annihilation center d i by type m vehicle depending on distance
C h Maintenance costs of type i products in the distribution–collection center d c
C N Equipment repair and maintenance costs
C t r Cost of staff training
C r i Recycling costs of type i products at the recycling center r
C d i Costs of annihilating type i products in the annihilation center d i
C d Delay costs per time
f R Inspection costs for each returned product unit in the distribution–collection center d c
f H Purchasing costs for each unit of type i   product from customer μ
C r e n t m Renting costs for a car type m
C b m Purchasing costs for a car type m
T i f 1 d Person–hour needed to manufacture each type i product unit at the first production center
T i f 2 d Person–hour needed to manufacture each type i product unit at the second production center
T w 1 Working hours during the order period at the first production center
T w 2 Working hours during the order period at the second production center
k a 1 Number of existing workers in the first production center
k a 2 Number of existing workers in the second production center
k e 1 Number of hired workers in the first production center
k e 2 Number of hired workers in the second production center
e s f Transport-related  C O 2  emission rates from supplier s to the production center of f per product unit
e f d c Transport-related  C O 2  emission rates from the production center f to the distribution–collection center d c per product unit
e d c μ Transport-related  C O 2  emission rates from the distribution–collection center d c to customer μ per product unit
e μ d c Transport-related  C O 2  emission rates from customer μ to the distribution–collection center d c   per product unit
e d c r Transport-related  C O 2  emission rates from the distribution–collection center d c to the recycling center r per product unit
e r f Transport-related  C O 2  emission rates from the recycling center r to the production center f per product unit
e r d i Transport-related  C O 2  emission rates from the recycling center r to the annihilation center d i per product unit
e d c d i Transport-related  C O 2  emission rates from the distribution–collection center d c to the annihilation center d i per product unit
e i C O 2 emission rates for producing each type i product unit
e r i C O 2 emission rates for the recycling process of each type i returned product unit
e d i i C O 2 emission rates for the annihilation process of each type i returned product unit
e d c C O 2 emission rates for the returned product inspection process in the distribution–collection center d c
d s f Supplier s to the production center f distance
d f d c Production center f to the distribution–collection center d c distance
d d c μ Distribution–collection center d c to client μ distance
d μ d c Client μ to the distribution–collection center d c distance
d d c r Distribution–collection center d c to the recycling center r distance
d r f Recycling center r to the production center f distance
d r d i Recycling center r to the annihilation center d i distance
d d c d i Distribution–collection center d c to the annihilation center d i distance
C a s Maximum capacity for supplier s
C a f Maximum capacity for the production center f
C a d c Maximum capacity for the distribution–collection center d c
C a r Maximum capacity for the recycling center r
C a d i Maximum capacity for the annihilation center d i
G Maximum allowable C O 2 emissions from production
M R Maximum allowable production–technical risk
t k Maximum acceptable delivery time for customer μ
d μ Customer demand
d f Raw materials’ demand of the production center f from supplier s
R l Device L failure risk
T e i C O 2 emission rate of total production at the production center
t t Total time from supplying raw materials from supplier s to sending final product to customer μ
t s f Length of time for supplying raw materials from supplier s to the production center f
t f d c Length of time for finished products to be sent from the production center f to the distribution–collection center d c
t d c μ Length of time for sending the finished product from the distribution–collection center d c to customer μ
α Product return rate
β Product return rate from the distribution–collection center d c to the annihilation center d i
γ Product return rate from the distribution–collection center d c to the recycling center r
λ Recycled product sending rate from the recycling center r to the production center f
τ Recycled product sending rate from the recycling center r to the annihilation center d i
q s μ Customer μ expected quality
q s d c d i Quality level of the product that leads to transfer from the distribution–collection center d c to the annihilation center d i
q s d c r Quality level of the product that leads to transfer from the distribution–collection center d c to the recycling center r
q s r d i Quality level of the product level that leads to transfer from the recycling center r to the annihilation center d i

Model’s Decision Variables

The proposed mathematical model has two types of decision variables. The first category includes non-zero continuous decision variables to determine the material and product flow between different facilities and the number of employees in each production facility, while also determining the time period of product delivery between various facilities in the direct supply chain. The second category consists of zero and one variables, which are used to select facilities, automobiles, and equipment for production centers.
Continuous Decision Variables:
p i , t Sale price of product i in cycle t
Y s f Raw material amounts sent from supplier s to the production center f
Y i f d c Type i product amounts sent from the production center f to the distribution–collection center d c
Y i d c μ Type i product amounts sent from the distribution–collection center d c to customer μ
Y i μ d c Type i product amounts sent from customer μ to the distribution–collection center d c
Y i d c r Type i product amounts delivered from the distribution–collection center d c to the recycling center r
Y i r f Recycled product amounts delivered from the recycling center r to the production center f
Y i r d i Product amounts delivered from the recycling center r to the annihilation center d i
Y i d c d i Returned product amounts delivered from the distribution–collection center d c to the annihilation center d i
t s f Time to send raw materials from supplier s to the production center f
t f d c Time to send product from the production center f to the distribution–collection center d c
t d c μ Time to send product from the distribution–collection center d c to customer μ
r μ i Returned rate of product i with quality level q s from customer μ
x μ d c . q s Returned product amounts considering quality level q s from customer μ to the distribution–collection center d c
Binary Variables:
X s If supplier s is chosen, it is equal to 1 ; otherwise, it equals 0 .
X f If the production center f is chosen, it is equal to 1 ; otherwise, it equals 0 .
X d c If the distribution–collection center d c is chosen, it is equal to 1 ; otherwise, it equals   0 .
X μ If customer μ is chosen, it is equal to 1 ; otherwise, it equals 0 .
X r If the recycling center r is chosen, it is equal to 1 ; otherwise, it equals 0 .
X d i If the annihilation center d i is chosen, it is equal to 1 ; otherwise, it equals 0 .
X v If the vehicle type m is sent, it is equal to 1 ; otherwise, it equals 0 .
X r m If the vehicle type m is rented, it is equal to 1 ; otherwise, it equals 0 .
X i m If the vehicle type m is purchased, it is equal to 1 ; otherwise, it equals 0 .
X l If equipment l is used, it is equal to 1 ; otherwise, it equals 0 .

3.3. Objective Functions

The study addresses a multi-objective, single-period, and multi-product model within a closed-loop supply chain network. The network design integrates environmental factors and imposes restrictions on material flow across different network levels. Facilities have limited capacities, and a dual-purpose location is implemented to minimize network costs for product distribution and collection. The locations of suppliers, production centers, distribution and collection centers, recycling centers, annihilation centers, and customers are known. The network consists of two production centers, two annihilation centers, two recycling centers, and various suppliers and distribution and collection centers. The reverse supply chain partially recycles returned products while removing the remaining ones from the network. Meeting specific inventory targets within a fixed sales period is crucial, particularly for perishable products. Packaging products have higher sales volumes due to their lower price compared to individual items. The total costs encompass several components, including sales, logistics, production, raw material procurement, transportation, delay, inventory, staff training, equipment maintenance, inspection, and recycling and annihilation costs. The model incorporates productive–technical risks associated with equipment failure and labor shortages. Uncertainty is considered in customer demand, product returns, recycling rates, and destruction. Greenhouse gas emissions are accounted for in production, transportation, and inspection processes throughout the supply chain network. The proposed model of the present study has four objective functions, all of which aim for minimization. Objective function (1) presented below, calculates and minimizes all network costs.
M i n   Z 1 = t N i N p i , t n i , t + i N S i , 0 q i + t N i N S i , t n i , t c i + s S f F f s Y s f X s X f + i N f F c i n i X f + i N μ N d c D C f H Y i μ d c X d c + i N μ N d c D C f R Y i μ d c X d c + i N r R d c D C C r i Y i d c r X r + i N d i i C d i Y i d c d i + Y i r d i X d i + l L C N X l + μ N i N d c D C C d max 0 , t t t k Y i d c µ + m M s S f F t r m s f d s f Y s f X s X f + m M i N d c D C f F t r m f d c d f d c Y i f d c X d c X f + m M i N μ N d c D C t r m d c μ d d c μ Y i d c μ X d c X μ + m M i N μ N d c D C t r m μ d c d μ d c Y i μ d c X d c X μ + m M i N r R d c D C t r m d c r d d c r Y i d c r X d c X r + m M r R i N f F t r m r f d r f Y i r f X r X f + m M i N r R d i D I t r m r d i d r d i Y i r d i X d i X r + m M i N r R d c D C t r m d c d i d d c d i Y i d c d i X d i X d c + i N d c D C f F C h Y i d f c X d c + m M C r e n t m X r m + m M C b m X i m + k e 1 K k e 2 K C t r ( k e 1 + k e 2 )
The first statement represents the cost of sales and logistics of product i based on inventory levels, costs, and raw material purchasing costs from supplier s . The second statement represents the production costs of each product unit. The third and fourth statements represent the purchasing costs of each product unit returned from customers and the inspection costs of each returned product unit in the distribution–collection center, respectively. The fifth and sixth statements represent the recycling costs of each returned product unit and the annihilating costs of each returned and unrecyclable product unit, respectively. The seventh and eighth statements represent the cost of equipment repairs and maintenance, and the cost of delaying product delivery per unit of time, respectively. The ninth to sixteenth statements, respectively, represent the raw material transportation costs from suppliers to production centers, products transported from production to distribution–collection centers, products from distribution–collection centers to customers’ locations, returned products transported from customers to distribution–collection centers, products from distribution–collection centers to recycling centers, products from recycling centers to production centers, products from recycling centers to annihilation centers, and products from distribution–collection centers to annihilation centers. The seventeenth statement represents maintenance costs in distribution–collection centers, and the eighteenth and nineteenth statements indicate the cost of renting and purchasing cars for transporting raw materials and products between facilities, respectively. Lastly, the twentieth statement is equivalent to the cost of training staff.
M i n   Z 2 = i N n i e i + i N r R Y i d c r e r i + i N d i D I r R ( Y i d c d i + Y i r d i ) e d i i + i N d c D C Y i μ d c e d c + m M s S f F Y s f e s f X v + m M i N d c D C f F Y i f d c e f d c X v + m M i N d c D C c C Y i d c µ e d c μ X v + m M i N μ N d c D C Y i μ d c e μ d c X v + m M i N d c D C r R Y i d c r e d c r X v + m M i N r R f F Y i r f e r f X v + m M i N r R d i D I Y i r d i e r d i X v + m M i N d c D C d i D I Y i d c d i e d c d i X v
Objective function (2) aims to minimize greenhouse gas emissions from different supply chain network processes. The first and second statements indicate the amount of emissions caused by the production process of products in production centers and the amount of emissions caused by the process of recycling, respectively. The third and fourth statements indicate the amount of emissions caused by the annihilation process and the amount of emissions caused by the inspection process of returned products in distribution–collection centers, respectively. The fifth to twelfth terms also represent the amount of emissions from the transportation process: from suppliers to production centers, from production to distribution–collection centers, from distribution–collection centers to customers, from customers to distribution–collection centers for transporting returned products, from distribution–collection centers to recycling centers, from recycling centers to production centers, from recycling centers to annihilation centers, and from distribution–collection centers to annihilation centers.
M i n   Z 3 = l L R l X l
Objective function (3) aims to minimize the risk of device failure. This function is aimed at minimizing production–technical risk, which is the most critical risk in the green supply chain. As mentioned in the model assumptions, two factors of interruption, caused by equipment failure and a shortage of skilled labor, are considered as the main causes of production–technical risk. The repair and maintenance of network equipment seem necessary and inevitable to minimize the interruption caused by equipment failure. Therefore, the costs for repairing and maintaining equipment in the network to prevent equipment failure or reduce downtime are considered and modeled in objective function (1). Also, staff training is considered an effective solution to solve the problem of the shortage of skilled labor, so costs for training employees to compensate for this shortfall and increase the skills of employees are modeled in objective function (1):
M i n   Z 4 = s S f F t s f X f + f F d c D C t f d c X d c + d c D C μ N t d c μ X μ
The fourth objective function aims to minimize product delivery time. The first to last statements indicate time taken for raw materials to be sent from the supplier to the production center, the time taken for products to be sent from production to distribution–collection centers, and the time taken for products to be sent from distribution–collection centers to customers, respectively.

3.4. Constraints

In the following, the constraints of the whole system is described in detail.
s S X s Y s f C a s
The amount of raw materials sent from suppliers should be smaller than their capacity.
f F i N d c D C μ N ( Y i f d c + Y i µ d c ) X d c C a d c
The amount of products entering distribution–collection centers should be smaller than their capacity.
μ N i N d c D C Y i d c μ X d c μ N d μ
All customer requests must be satisfied.
i N r R d c D C Y i d c r X r C a r
The total amount of goods sent from distribution–collection to recycling centers should be smaller than the capacity of recycling centers.
i N r R f F Y i r f X r C a f
The total amount of products delivered from recycling centers to production centers must be smaller than the production center’s capacity.
i N d c D C r R d i D I Y i d c d i X d c X d i + Y i r d i X d i X r C a d i
The total amount of goods sent to annihilation centers should be smaller than the capacity of annihilation centers.
f F s S X s Y s f d f
The amount of raw materials sent from the suppliers to the production centers should be large enough to meet the demands of the production center.
i N d c D C μ N Y i μ d c = μ N i N d c D C α Y i d c μ
The total amount of returned products from customers to distribution–collection centers should be α % of the total products sent to customers.
i N d c D C f F Y i f d c = i N N i
The total number of type i products delivered from the production centers to distribution centers must equal the total production of type i products.
i N d c D C d i D I Y i d c d i = i N μ N d c D C β Y i μ d c
The total product delivered from the distribution–collection centers to annihilation centers can be represented as β % of the total returned products sent from customers to distribution–collection centers.
i N d c D C r R Y i d c r = i N μ N d c D C γ Y i μ d c
The total number of products sent from distribution–collection centers to recycling centers can be indicated as γ % of the total returned products sent from customers to distribution–collection centers.
d i D I i N r R Y i r d i = i N d c D C r R τ Y i d c r
The total number of products sent from recycling centers to annihilation centers should be τ % of the total products sent from distribution–collection centers to recycling centers.
i N d c D C μ N f F Y i f d c + Y i μ d c i N d c D C μ N r R d i D I Y i d c μ + Y i d c r + Y i d c d i
The sum of products that enter distribution–collection centers must exceed the output sum from distribution–collection centers.
i N r R f F Y i r f = i N d c D C r R λ Y i d c r
The total product sent from recycling centers to production centers should be λ % of the total product amounts sent from distribution–collection centers to recycling centers.
i N d c D C r R Y i d c r = f F i N r R d i D I Y i r f + Y i r d i
The total product amounts entering recycling centers must equal the total output of recycling centers.
i N d c D C μ N Y i d c μ = i N d c D C f F Y i f d c  
The sum of type i products delivered from distribution–collection centers to customers must equal the type i product total number sent from production centers to distribution–collection centers.
m M X v X r m + X i m
The total number of vehicles on the journey must be less than or equal to the total purchased and rental vehicles in the organization.
l L R l M R
Production–technical risk must be smaller than the maximum production–technical risk allowed.
i N T e i G
The total C O 2 emission rate of production in production facilities should be less than the maximum allowable rate of C O 2 emissions from production.
t t = t s f + t f d c + t d c µ
The product delivery total time to the customer equals the total time for product transfer from the supplier to the production centers, the product transfer time from production centers to distribution–collection centers, and the product transfer time from distribution–collection centers to the customers.
β + γ + λ + τ = 1
The above constraint ensures that the value of 1 is obtained for the sum of the returned products’ coefficients.
d c D C μ N x μ d c . q s = r μ i  
The above statement ensures the collection of all returned products from customer centers throughout the return process.
q s d c d i + q s d c r + q s r d i = q s μ
The above statement ensures that the total product quality in the whole chain (except the product transfer stage from recycling centers to manufacturing centers) equals the original product quality.
l L i N n i X l = d c D C f F i N Y i f d c
The total production of type i equipment in each production center equals the number of products delivered from the desired production center to distribution–collection centers.
f F i N T i f 1 d Y i f d c k a 1 k k e 1 k k a 1 + k e 1 T w 1
Person–hours needed to manufacture each product type i unit in the first production center should be supplied by the existing labor force and by hiring new labor in case of shortages.
f F i N T i f 2 d Y i f d c k a 2 k k e 2 k k a 2 + k e 2 T w 2
Person–hours needed to manufacture each product type i unit in the second production center should be supplied by the existing labor force and by hiring new labor in case of shortages.
Y s f , Y i f d c , Y i d c µ , Y i µ d c , Y i d c r , Y i r f , Y i r d i , Y i d c d i , t s f , t f d c , t d c µ , r μ i , x μ d c . q s > 0
X s , X f , X d c , X μ , X r , X d i , X v , X r m , X i m , X l 0,1
These two constraints ensure that the mentioned parameters and variables are positive and between zero and one, respectively.
t s f . m i n t s f t s f . m a x
This constraint ensures that the time spent on raw material sending from supplier s to the production center f must be within the time frame set by the production center f .
t f d c . m i n t f d c t f d c . m a x
This constraint ensures that the time it takes to send final products from the production center f to the distribution–collection center d c must be within the time frame set by the distribution–collection center d c .
t d c µ . m i n t d c µ t d c µ . m a x
This constraint ensures that the time it takes to send final products from the distribution–collection center d c to customer μ must be within the time frame set by customer μ .
X s , X f , X d c , X r , X d i , X v , X r m , X i m , X l 1
s S , f F , d c D C , r R , d i D I , m M , l L
This constraint ensures that at least one facility is used during the product transfer to the customer.
q i + c i p i , t p i m a x
This constraint ensures that pricing must be higher than the cost and below the maximum price reserved by consumers.
n i , t 0
This constraint ensures that the volume of product sales in each decision cycle must be non-negative.
n i , t S i
This constraint is an indirect limiting condition to ensure that the sales volumes of the two types of products must be less than the total inventory.
m a x p 1 , t , p 2 , t < p 3 , t < p 1 , t + p 2 , t
According to the bundling strategy, this constraint represents a package pricing that must be less than the sum of the separated pricing of two products and higher than the highest price of the two separate products sold. Otherwise, advertising through the bundling strategy is meaningless.
n 3 , t < min S i , t
This constraint ensures that the sales volume of the bundling product must be less than the minimum inventory of two products.
n 1 , t + n 3 , t S i ,   n 2 , t + n 3 , t S i
This constraint is also a limited requirement to ensure that the sales volumes of two product types must be less than the total inventory.

4. Numerical Example and Results

4.1. Numerical Example

As a result of the NP-hard nature of the model and its computational complexity, it is not possible to solve it with exact methods. Hence, we will use a meta-heuristic algorithm to solve the model. Among multi-objective meta-heuristic algorithms, the NSGA-II algorithm has been selected in this research due to its significant advantages. The nature of this algorithm’s random search in the problem space is considered a parallel search, because each of the random chromosomes generated by the algorithm is considered a new starting point for searching a part of the problem state space, and the search is conducted across all of them simultaneously. Also, due to the extensive dispersion of the points that are searched, it obtains favorable results for problems that have large search spaces. It is also considered a targeted random search and will reach different answers through different paths. In addition, it does not face any restrictions in the search and selection of random answers. Finally, due to competition among answers and selection of the best among the population, it will reach the global optimal point with a high probability. The proposed NSGA-II algorithm starts with a population initialization process that randomly generates populations with the values of each gene bounded by given values. Then, it uses a selection process to form a pool of parents to generate offspring using a crossover process. The genetic algorithm used in this research uses a simple crossover process that randomly selects two parents as parent1 and parent2 and then generates a random variable (r) from 0 to 1. The weight of this random variable is defined as R. We considered R as 0.8 and used c h i l d = p a r e n t 1 + r R ( p a r e n t 2 p a r e n t 1 ) to create each gene from a child. This is shown in Figure 2.
In the proposed NSGA-II algorithm, in order to improve population diversity and convergence speed, a combination of disorder mapping and conflict-based learning methods is used to generate the initial population. Also, an approach based on a penalty function is used for solutions that do not meet the time limit. This problem shows the differences between the proposed method and the classical NSGA-II and MOPSO methods.
The problem at hand is a multi-objective, single-period, multi-product green CLSC model. The settings for the NSGAII algorithm are presented in Table 2. It should be noted that the implementation of the proposed method was carried out using MATLAB version 2023 software, and the coding was performed on a system with a core i5 processor and 8 GB of RAM.
Table 3 shows the initial values of some main parameters. These data were collected from a food packaging company in Tehran, Iran.
During the process of solving the model, some solutions were infeasible, which is due to the placement of the algorithm in local solutions. In order to overcome this issue, the proposed algorithm shifts from local mode and solves this problem under global conditions. The difference between feasible and infeasible solutions lies in establishing the necessary restrictions and conditions to solve the problem that is considered in the model.
Now, results related to all four objective functions have been obtained by setting the population number to 10, the operators’ probability to 0.8, and the mutation operators’ probability to 0.3. Nevertheless, as the developed model represents several different levels with 10 customers, five distributors, and two manufacturers (Table 3), the flow of products is determined considering the amount of customer demand. Hence, Table 4 illustrates the amount of demand for the 10 target customers separately.
According to the above table, the first customer demands 176 pieces, the second demands 329 pieces, and so on until the end. It is necessary to note that the maximum time needed for product delivery to each customer without causing dissatisfaction is set at 120 h. Additionally, given the completely random determination of initial model values in problem-solving by the genetic algorithm (GA), unreasonable or infeasible solutions may arise, prompting the algorithm to promptly address the problem and find a practical answer. The above-mentioned process should be considered by the algorithm during 300 iterations. Table 5 presents the model implementation results with the genetic multi-objective algorithm for all four specific objectives.
As shown in Table 5, the desired results are presented for the proposed model with a numerical example. The point that should be mentioned, regarding the presentation of results related to risk, is that in this research, supply chain risk assessment is related to two interruption factors: equipment failure and the lack of skilled labor, both considered the main factors contributing to production–technical risk. Any failure leads to an interruption in production of the product and increases the costs of the entire chain. Therefore, equipment maintenance costs and staff training costs are presented in the form of a cost model to compensate for this deficiency and increase the skills of employees. Therefore, by reducing risk, in addition to minimizing the delay in the production and shipping of the product, the costs of the entire chain are also reduced.
Table 6 and Table 7 present the product transfer amounts from manufacturers to distributors and from distributors to customers to meet their expectations throughout the considered interval. The presented values highlight the desired figures to achieve the target functions’ optimal amounts.
Based on Table 6, the first producer only transferred 1231 units of products to the fifth distributor, and the second producer only transferred 1853 units to the first distributor, 2 units to the third distributor, and 272 units to the fourth distributor. However, there have been no product deliveries to the second distributor. The issue that needs consideration is that the total customer demand, as indicated in Table 6, is 3358 units. As mentioned before, all demands should be answered. Therefore, the total amount of the provided products in Table 6 should equate to 3358 units as well. In addition, the maximum production capacities for the first and second producers were 1577 and 3071 units, adopting 1231 and 2127 units for each, respectively. In other words, the product amounts sent from manufacturers to distributors showed lower values than the manufacturers’ maximum capacities. Table 7 presents the amounts of product transfers from each distributor to customers.
According to this table, the first customer takes all the demand of 176 units from the fourth distributor, and the ninth customer receives 9 and 225 product units from the first and second distributors, respectively, to satisfy its request of 234 units. As shown, each row’s sum equals the total number of customer demands met.

4.2. Different Sample Size Solution

Due to the complexity of the model and its limitations in the real world, this section aims to use different numerical multiples across different dimensions to make applications in the real world easier. So, in this section, we will examine the model’s solution for various examples in small, medium, and large dimensions. Table 8, Table 9 and Table 10 highlight the findings for all four objective functions for all three problem dimensions.
Figure 3, Figure 4, Figure 5 and Figure 6 are presented to compare the objective functions of cost, time, pollutant emission, and risk level, respectively, to better understand the results.

4.3. Sensitivity Analysis

Sensitivity analysis is a topic that can provide proper insight into solving problems. In other words, sensitivity analysis determines how much the dependent variable will change if the value of an independent variable changes in a specific and defined situation, assuming that other variables are constant. Here, the values presented in Table 3 are changed to find out how an increase or decrease in each parameter affects the final cost of the model. Based on sensitivity analyses performed on the main model parameters, the supplier, collection–distribution center, and producer numbers have the greatest impact on the target functions. It should be noted that the time required for product delivery to customers has an inverse relationship with the increase in these parameters. Hence, increased supplier, collection–distribution center, and producer numbers would lead to a decrease in the time required to transfer products to customers. Additionally, parameters related to destruction and recycling centers only have an effect on pollutant levels and chain costs, with no effects on objective functions related to time and risk. Thus, with an increase in destruction centers, the amount of pollutants increases more, while an increase in recycling centers results in a higher chain costs. The increased number of customers and produced products does not affect risk. Since the impact of other problem parameters, such as the costs of producing each product unit, purchasing raw materials, training, recycling, destruction, and inspection, as presented in Table 3, only affect chain costs and have no effects on other objective functions, they were not investigated in the sensitivity analysis section.

4.4. Comparison

We also use the MOPSO algorithm and compare its performance with the NSGA-II algorithm in solving large-scale problems to evaluate the developed methodology. Table 11 presents the comparison results.
As revealed by comparisons, the genetic optimization method (NSGA-II) has higher diversity in discovering solutions due to its ability to identify local and global optimal solutions and its continuity. Compared to MOPSO, a better balance is established between discovery and efficiency criteria to prevent the risk of falling into local optima. Therefore, in addition to its high power in discovering optimal answers, the genetics-based method shows higher speed and subsequently a higher convergence rate and diversity than MOPSO.
In the following analysis, we use the MGA, as employed by Gholizadeh and Fazlollahtabar (2020), to maximize the single objective model associated with the CLSC total costs [34]. The objective function presented here seeks to maximize the total profit of the closed-loop green SC. Gross profit equals the difference between income and expenses. Revenue sources are products sold to customers in both Tier 1 and Tier 2. The company’s total costs also comprise operational and transportation costs. Thus, each period’s operating costs within the future flow are equal to raw material purchasing costs, production of grade 1 products, assembly of grade 1 and 2 products, and the distribution centers’ operating costs. A reverse supply chain also involves paying for the purchase of used products from customers, separation costs of returned products, quality checking of separated parts in the parts separation center, and waste disposal costs in the reverse flow. Hence, Table 12 presents the modeling results of the modified GA and the proposed model.
Based on Table 12, among the 10 examples of the solved problem, the algorithm introduced in this research has favorable performance in determining chain costs and enhacing profits, with the highest difference equal to 4.72% and the average difference equal to 2.38% compared to MGA.
The Non-dominated Sorting Genetic Algorithm II (NSGA-II) outperforms the modified genetic algorithm due to several key factors inherent in its design and operational mechanisms. NSGA-II’s enhanced performance, which led to a 2.38% reduction in costs, can be attributed to its superior handling of multiple objectives, preservation of diversity among solutions, and efficient sorting and selection process.

5. Discussion and Conclusions

In this paper, the multi-objective genetic optimization method was used to find the optimal answer to the four-objective closed-loop supply chain nonlinear programming problem. The study involved comparing and analyzing the optimal values obtained from each method across different dimensions. The proposed algorithm was performed in several phases. In the first phase, the problem was modeled and initialized, encoding the supply chain network in question with a string of real numbers. Each solution (chromosome) in the genetic algorithm is equivalent to the components of suppliers, distribution centers, producers, repair centers, collection, renovation, destruction, and recycling, and the values of each sub-component of these parameters are randomly initialized within the intervals determined for each. Genetic algorithms usually use a higher quality population to speed up convergence. In the second phase, a criterion must be defined to evaluate the members of the population and enable the recognition of the better organisms in the population. This work, i.e., determining the suitability of an entity, is called the evaluation of that entity. The fitness function is equivalent to the planning problem discussed in Section 3, which deals with profit maximization and environmental impact minimization. The determined values for each of the components of suppliers, distribution centers, producers, and other investigated centers are placed in the objective function to calculate the objective functions and solve the problem according to the considered limitations. In the next step, Pareto solutions are extracted depending on the genetic algorithm model, representing the optimal values of regulatory parameters, including the number of disposals, renovation, repair, recycling, distributor, and supplier centers. The experimental results showed that the genetic algorithm had a favorable performance in finding solutions in all three dimensions, small, medium, and large, in finding the optimal solution, maximizing profit, and reducing the effects on the environment, risk and product transfer time. It was also determined by sensitivity analysis that the parameters related to the number of suppliers, collection–distribution centers, and manufacturers have the greatest impact on the performance of the proposed model’s objectives. It was also found that the product delivery time to the customer will have the opposite effect with the increase in these parameters, and with the increase in the number of suppliers, collection–distribution centers, and manufacturers, the product transfer time to the customer will decrease. The results of the comparison of the proposed method with the MOPSO algorithm showed that the genetic optimization method has more diversity in discovering solutions due to finding local and global optimal solutions and its continuous nature, and a better balance is established between the discovery and efficiency criteria to avoid failing to local optima. It was also found that compared to the MOPSO method, it has a high speed and subsequently higher convergence and diversity in achieving optimal solutions. In other words, the advantage of the proposed method compared to other methods is in automatic subset creation and finding global optima and local optima to consciously maximize the profit from each stage of the supply chain. In addition, the proposed method has wider applicability, and we intend to focus on the design of supply chain systems with other objectives in the uncertain and ambiguous periods in the future. Also, in the future, the authors will attempt to perform a quantitative assessment of the environmental benefits resulting from pollutant reduction strategies, possibly including life cycle assessments (LCAs) and economic analyses of cost-saving measures, potentially incorporating a break-even analysis to assess financial feasibility.

Author Contributions

Conceptualization, G.H.P. and Y.H.; Methodology, G.H.P.; Software, G.H.P.; Validation, G.H.P.; Formal analysis, G.H.P.; Investigation, G.H.P.; Resources, G.H.P. and Y.H.; Data curation, G.H.P.; Writing—original draft, G.H.P.; Writing—review & editing, Y.H. and S.H.P.; Visualization, G.H.P.; Supervision, Y.H.; Project administration, Y.H.; Funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

The current work received support from the National Natural Science Foundation of China (No. 72171047).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alzoubi, H.M.; Ghazal, T.M.; Sahawneh, N.; Al-kassem, A.H. Fuzzy assisted human resource management for supply chain management issues. Ann. Oper. Res. 2022, 326, 125004644. [Google Scholar]
  2. Rolf, B.; Jackson, I.; Müller, M.; Lang, S.; Reggelin, T.; Ivanov, D. A review on reinforcement learning algorithms and applications in supply chain management. Int. J. Prod. Res. 2022, 61, 7151–7179. [Google Scholar] [CrossRef]
  3. Boskabadi, A.; Mirmozaffari, M.; Yazdani, R.; Farahani, A. Design of a distribution network in a multi-product, multi-period green supply chain system under demand uncertainty. Sustain. Oper. Comput. 2022, 3, 226–237. [Google Scholar] [CrossRef]
  4. Chen, J.; Wang, H.; Fu, Y. A multi-stage supply chain disruption mitigation strategy considering product life cycle during COVID-19. Environ. Sci. Pollut. Res. 2022, 1–15. [Google Scholar] [CrossRef]
  5. Roh, T.; Noh, J.; Oh, Y.; Park, K.S. Structural relationships of a firm’s green strategies for environmental performance: The roles of green supply chain management and green marketing innovation. J. Clean. Prod. 2022, 356, 131877. [Google Scholar] [CrossRef]
  6. Abir, A.S.; Bhuiyan, I.A.; Arani, M.; Billal, M.M. Multi-Objective Optimization for Sustainable Closed-Loop Supply Chain Network Under Demand Uncertainty: A Genetic Algorithm. arXiv 2020, arXiv:2009.06047. [Google Scholar]
  7. Arani, M.; Liu, X.; Abdolmaleki, S. Scenario-based simulation approach for an integrated inventory blood supply chain system. In Proceedings of the 2020 Winter Simulation Conference (WSC), Orlando, FL, USA, 14–18 December 2020; pp. 1348–1359. [Google Scholar]
  8. Large, R.O.; Thomsen, C.G. Drivers of Green Supply Management Performance: Evidence from Germany. J. Purch. Supply Manag. 2016, 17, 176–184. [Google Scholar] [CrossRef]
  9. Zhu, Q.; Sarkis, J.; Lai, K.H. Green Supply Chain Management Practices and Sustainability Performance. J. Purch. Supply Manag. 2015, 19, 106–117. [Google Scholar] [CrossRef]
  10. Jabbarzadeh, A.; Haughton, M.; Khosrojerdi, A. Closed-loop Supply Chain Network Design under Disruption Risks: A Robust Approach with Real World Application. Comput. Ind. Eng. 2018, 116, 178–191. [Google Scholar] [CrossRef]
  11. Devika, K.; Jafarian, A.; Nourbakhsh, V. Designing a sustainable closed-loop supply chain network based on triple bottom line approach: A comparison of metaheuristics hybridization techniques. Eur. J. Oper. Res. 2014, 235, 594–615. [Google Scholar] [CrossRef]
  12. Farrokh, M.; Azar, A.; Jandaghi, G.; Ahmadi, E. A novel robust fuzzy stochastic programming for closed loop supply chain network design under hybrid uncertainty. Fuzzy Sets Syst. 2018, 341, 69–91. [Google Scholar] [CrossRef]
  13. Liu, M.; Liu, R.; Zhu, Z.; Chu, C.; Man, X. A bi-objective green closed loop supply chain design problem with uncertain demand. Sustainability 2018, 10, 967. [Google Scholar] [CrossRef]
  14. Mohtashami, Z.; Aghsami, A.; Jolai, F. A green closed loop supply chain design using queuing system for reducing environmental impact and energy consumption. J. Clean. Prod. 2020, 242, 118452. [Google Scholar] [CrossRef]
  15. Rezae, S.; Kheirkhah, A. A comprehensive approach in designing a sustainable closed-loop supply chain network using cross-docking operations. In Computational and Mathematical Organization Theory; Springer Science + Business Media: New York, NY, USA, 2017. [Google Scholar]
  16. Garai, A.; Chowdhury, S.; Sarkar, B.; Roy, T.K. Cost-effective subsidy policy for growers and biofuels-plants in closed-loop supply chain of herbs and herbal medicines: An interactive bi-objective optimization in T-environment. Appl. Soft Comput. 2021, 100, 106949. [Google Scholar] [CrossRef]
  17. Jerbiaa, R.; Kchaou Boujelbenb, M.; Amine Sehlia, M.; Jemaia, Z. A stochastic closed-loop supply chain network design problem with multiple recovery options. Comput. Ind. Eng. 2018, 118, 23–32. [Google Scholar] [CrossRef]
  18. Santander, P.; Sanchez, F.A.C.; Boudaoud, H.; Camargo, M. Closed loop supply chain network for local and distributed plastic recycling for 3D printing: A MILP-based optimization approach. Resources. Conserv. Recycl. 2020, 154, 104531. [Google Scholar] [CrossRef]
  19. Goodarzian, F.; Hosseini-Nasab, H.; Fakhrzad, M.B. A Multi-objective Sustainable Medicine Supply Chain Network Design Using a Novel Hybrid Multi-Objective Metaheuristic Algorithm. Int. J. Eng. 2020, 33, 1986–1995. [Google Scholar] [CrossRef]
  20. Zarbakhshnia, N.; Kannan, D.; Mavi, R.K. Annals of Operations Research A novel sustainable multi-objective optimization model for forward and reverse logistics system under demand uncertainty. Ann. Oper. Res. 2020, 295, 843–880. [Google Scholar] [CrossRef]
  21. Nasr, A.K.; Tavana, M.; Alavi, B.; Mina, H. A novel fuzzy multi-objective circular supplier selection and order allocation model for sustainable closed-loop supply chains. J. Clean. Prod. 2021, 287, 124994. [Google Scholar] [CrossRef]
  22. Jian, J.; Li, B.; Zhang, N.; Su, J. Decision-making and coordination of green closed-loop supply chain with fairness concern. J. Clean. Prod. 2021, 298, 126779. [Google Scholar] [CrossRef]
  23. Tavana, M.; Kian, H.; Nasr, A.K.; Govindan, K.; Mina, H. A comprehensive framework for sustainable closed-loop supply chain network design. J. Clean. Prod. 2022, 332, 129777. [Google Scholar] [CrossRef]
  24. Cheng, P.; Ji, G.; Zhang, G.; Shi, Y. A closed-loop supply chain network considering consumer’s low carbon preference and carbon tax under the cap-and-trade regulation. Sustain. Prod. Consum. 2022, 29, 614–635. [Google Scholar] [CrossRef]
  25. Kouchaki Tajani, E.; Ghane Kanafi, A.; Daneshmand-Mehr, M.; Hosseinzadeh, A.A. A Robust Green Multi-Channel Sustainable Supply Chain based on RFID technology with considering pricing strategy and subsidizing policies. J. Ind. Eng. Int. 2022, 18, 43–78. [Google Scholar]
  26. Babaeinesami, A.; Tohidi, H.; Ghasemi, P.; Goodarzian, F.; Tirkolaee, E.B. A closed-loop supply chain configuration considering environmental impacts: A self-adaptive NSGA-II algorithm. Appl. Intell. 2022, 52, 13478–13496. [Google Scholar] [CrossRef]
  27. Alinezhad, M.; Mahdavi, I.; Hematian, M.; Tirkolaee, E.B. A fuzzy multi-objective optimization model for sustainable closed-loop supply chain network design in food industries. Environ. Dev. Sustain. 2022, 1–28. [Google Scholar] [CrossRef]
  28. Bathaee, M.; Nozari, H.; Szmelter-Jarosz, A. Designing a new location-allocation and routing model with simultaneous pick-up and delivery in a closed-loop supply chain network under uncertainty. Logistics 2023, 7, 3. [Google Scholar] [CrossRef]
  29. Wang, Z.; Zhen, H.L.; Deng, J.; Zhang, Q.; Li, X.; Yuan, M.; Zeng, J. Multiobjective optimization-aided decision-making system for large-scale manufacturing planning. IEEE Trans. Cybern. 2021, 52, 8326–8339. [Google Scholar] [CrossRef]
  30. Leung, M.F.; Coello, C.A.C.; Cheung, C.C.; Ng, S.C.; Lui, A.K.F. A hybrid leader selection strategy for many-objective particle swarm optimization. IEEE Access 2020, 8, 189527–189545. [Google Scholar] [CrossRef]
  31. Bahrampour, P.; Najafi, S.E.; Edalatpanah, A. Designing a Scenario-Based Fuzzy Model for Sustainable Closed-Loop Supply Chain Network considering Statistical Reliability: A New Hybrid Metaheuristic Algorithm. Complexity 2023, 2023, 1337928. [Google Scholar] [CrossRef]
  32. Wu, Z.; Qian, X.; Huang, M.; Ching, W.K.; Wang, X.; Gu, J. Recycling channel choice in closed-loop supply chains considering retailer competitive preference. Enterp. Inf. Syst. 2023, 17, 1923065. [Google Scholar] [CrossRef]
  33. Dey, S.K.; Giri, B.C. Corporate social responsibility in a closed-loop supply chain with dual-channel waste recycling. Int. J. Syst. Sci. Oper. Logist. 2023, 10, 2005844. [Google Scholar] [CrossRef]
  34. Gholizadeh, H.; Fazlollahtabar, H. Robust optimization and modified genetic algorithm for a closed loop green supply chain under uncertainty: Case study in melting industry. Comput. Ind. Eng. 2020, 147, 106653. [Google Scholar] [CrossRef]
Figure 1. The schematic of the developed model.
Figure 1. The schematic of the developed model.
Mathematics 12 00737 g001
Figure 2. Crossover and mutation operators.
Figure 2. Crossover and mutation operators.
Mathematics 12 00737 g002
Figure 3. Comparison of costs in each dimension.
Figure 3. Comparison of costs in each dimension.
Mathematics 12 00737 g003
Figure 4. Comparison of times in each dimension.
Figure 4. Comparison of times in each dimension.
Mathematics 12 00737 g004
Figure 5. Comparison of pollution emissions in each dimension.
Figure 5. Comparison of pollution emissions in each dimension.
Mathematics 12 00737 g005
Figure 6. Comparison of risks in each dimension.
Figure 6. Comparison of risks in each dimension.
Mathematics 12 00737 g006
Table 1. Literature review summary.
Table 1. Literature review summary.
ReferencesObjectivesReturn System TypeConditionProduct Type
CostTimeEnvironmentRiskSocialRecyclingDestructionCertainUncertaintyPerishableImperishable
Jerbiaa et al. (2018) [17]
Mohtashami et al. (2020) [14]
Santander et al. (2020) [18]
Goodarzian et al. (2020) [19]
Zarbakhshnia et al. (2020) [20]
Nasr et al. (2021) [21]
Jian et al. (2021) [22]
Tavana et al. (2022) [23]
Cheng et al. (2022) [24]
Kouchaki Tajani et al. (2022) [25]
Babaeinesami et al. (2022) [26]
Alinezhad et al. (2022) [27]
Bathaee et al. (2023) [28]
Bahrampour et al. (2023) [31]
Wu et al. (2023) [32]
Dey and Giri (2023) [33]
Present Study
Table 2. Initial parameters of NSGA-II.
Table 2. Initial parameters of NSGA-II.
ParameterValue
Population size60
No. of generations100
Selection methodRandom
Non-dominant choiceTournament DCD
Crossover methodLaplace
Probability of crossover95%
Mutation methodPower
Probability of mutation0.5%
Table 3. Initial model parameters.
Table 3. Initial model parameters.
ParameterSymbolValue
Supplier s 4
Production center f 2
Distribution-collection center D c 5
Destruction center D i 2
Recycling center r 2
Customer C 10
Staff k 5
Set of products p 3
Set of equipment l 4
All vehicle types m 3
Number of products produced N p 1000
Production costs per product unit C p 50,000
The costs of each raw material unit purchasing from the suppliers f s 10,000
The employee training costs C t r 1,200,000
Product recycling costs C r p 20,000
Product destruction costs C d i 13,000
Inspection fee per returned product unit F c h p 30,000
Table 4. Customer demand values.
Table 4. Customer demand values.
Customers12345678910
Demand176329427729102356449224234332
Table 5. Problem-solving results with the genetic algorithm (NSGAII).
Table 5. Problem-solving results with the genetic algorithm (NSGAII).
The Objective FunctionCost Objective Function ($)Pollutant Emission Objective Function (PPM)Risk Objective Function (%)Time Objective Function (h)Execution Time (s)
Value10,981,18511,7440.2711714.5963
Table 6. Transferring the products (units) from manufacturers to distributors.
Table 6. Transferring the products (units) from manufacturers to distributors.
Distributors
The first producer00001231
The second producer1853022720
Table 7. Product transfer (units) from distributors to customers.
Table 7. Product transfer (units) from distributors to customers.
Customers
First Distributor004276910227449090
Second Distributor0000000000
Third Distributor0200000000
Fourth Distributor1760009600000
Fifth Distributor03270386790227225332
Table 8. Objective function results for the problem in small dimensions.
Table 8. Objective function results for the problem in small dimensions.
Samples12345678910
Costumer No.1244778899
Distributer No.2223344455
Cost103,813104,190104,938109,889110,440120,444133,912145,530145,838147,781
Pollution (ppm)35038439945045963288593510021027
Time (h)781013192539446875
Risk (%)0.0580.1190.130.2980.1510.2240.3240.1330.2370.276
Table 9. Objective function results for the problem in medium dimensions.
Table 9. Objective function results for the problem in medium dimensions.
Samples12345678910
Costumer No.12141415151616181818
Distributer No.6677889121314
Cost1,086,8331,180,5321,209,7451,222,1651,232,8361,295,8911,302,7851,339,5291,344,9991,377,401
Pollution (ppm)1178120813211342157416221700176519322120
Time (h)50515262646973778084
Risk (%)0.0590.3410.3670.2280.3090.230.3990.2020.2170.389
Table 10. Objective function results for the problem in large dimensions.
Table 10. Objective function results for the problem in large dimensions.
Samples12345678910
Costumer No.30455050606565707070
Distributer No.20202530323841455060
Cost1,520,3982,026,5272,142,3702,165,4512,255,8482,257,9462,342,9302,411,2602,412,8172,595,378
Pollution (ppm)13,52613,65814,00214,23015,20016,21016,89017,53218,05019,850
Time (h)112114123129135139145160175191
Risk (%)0.1770.3460.4430.3030.3910.4080.4790.5870.5990.631
Table 11. Comparing the performance of NSGA-II and MOPSO in model solving.
Table 11. Comparing the performance of NSGA-II and MOPSO in model solving.
Distribution No.Costumers No.NSGA-IIMOPSO
CostPollution (ppm)Risk (%)Time (h)CostPollution (ppm)Risk (%)Time (h)
203015,203,98013,5260.17711215,878,00915,3210.201122
204520,265,27013,6580.34611416,508,71416,8070.359126
255021,423,70114,0020.44312321,046,48517,2430.481131
305021,654,51414,2300.30312921,245,25917,2960.329135
326022,558,48715,2000.39113522,857,77217,6330.405141
386522,579,46016,2100.40813923,429,83818,8280.415154
416523,429,30116,8900.47914523,968,08319,3090.502161
457024,112,60717,5320.58716024,099,88619,6550.592167
507024,128,17218,0500.59917525,479,45219,7380.621181
Table 12. Comparison of the developed algorithm (NSGA-II) and MGA, proposed by Gholizadeh and Fazlollahtabar (2020) [34].
Table 12. Comparison of the developed algorithm (NSGA-II) and MGA, proposed by Gholizadeh and Fazlollahtabar (2020) [34].
Problem No.12345678910
Modified GA11,005,41310,363,21615,344,75820,483,95616,050,51028,996,98920,668,67922,790,12323,998,90814,999,575
Present Study11,152,82110,852,52415,817,90820,854,21516,642,75829,528,21021,217,93123,269,60424,752,04215,241,798
Diff. (%)1.344.723.081.813.691.832.662.103.141.61
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pakdel, G.H.; He, Y.; Pakdel, S.H. Multi-Objective Green Closed-Loop Supply Chain Management with Bundling Strategy, Perishable Products, and Quality Deterioration. Mathematics 2024, 12, 737. https://0-doi-org.brum.beds.ac.uk/10.3390/math12050737

AMA Style

Pakdel GH, He Y, Pakdel SH. Multi-Objective Green Closed-Loop Supply Chain Management with Bundling Strategy, Perishable Products, and Quality Deterioration. Mathematics. 2024; 12(5):737. https://0-doi-org.brum.beds.ac.uk/10.3390/math12050737

Chicago/Turabian Style

Pakdel, Golnaz Hooshmand, Yong He, and Sina Hooshmand Pakdel. 2024. "Multi-Objective Green Closed-Loop Supply Chain Management with Bundling Strategy, Perishable Products, and Quality Deterioration" Mathematics 12, no. 5: 737. https://0-doi-org.brum.beds.ac.uk/10.3390/math12050737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop