Next Article in Journal
Exact Solutions to the Maxmin Problem max‖Ax‖ Subject to ‖Bx‖≤1
Previous Article in Journal
The Bourguignon Laplacian and Harmonic Symmetric Bilinear Forms
Previous Article in Special Issue
Cooperative Co-Evolution Algorithm with an MRF-Based Decomposition Strategy for Stochastic Flexible Job Shop Scheduling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sustainable Closed-Loop Supply Chain Design Problem: A Hybrid Genetic Algorithm Approach

1
Department of Business Administration, Chosun University, 309 Pilmun-daero, Dong-gu, Gwangju 61452, Korea
2
Fuzzy Logic Systems Institute, Tokyo University of Science, 1-3 Kagurazaka Shinjiku-ku, Tokyo 162-8601, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(1), 84; https://doi.org/10.3390/math8010084
Submission received: 4 November 2019 / Revised: 19 December 2019 / Accepted: 2 January 2020 / Published: 4 January 2020

Abstract

:
In this paper, we propose a solution to the sustainable closed-loop supply chain (SCLSC) design problem. Three factors (economic, environmental, and social) are considered for the problem and the three following requirements are addressed while satisfying associated constraint conditions: (i) minimizing the total cost; (ii) minimizing the total amount of CO2 emission during production and transportation of products; (iii) maximizing the social influence. Further, to ensure the efficient distribution of products through the SCLSC network, three types of distribution channels (normal delivery, direct delivery, and direct shipment) are considered, enabling a reformulation of the problem as a multi-objective optimization problem that can be solved using Pareto optimal solutions. A mathematical formulation is proposed for the problem, and it is solved using a hybrid genetic algorithm (pro-HGA) approach. The performance of the pro-HGA approach is compared with those of other conventional approaches at varying scales, and the performances of the SCLSC design problems with and without three types of distribution channels are also compared. Finally, we prove that the pro-HGA approach outperforms its competitors, and that the SCLSC design problem with three types of distribution channels is more efficient than that with a single distribution channel.

1. Introduction

In general, closed-loop supply chain (CLSC) design problems are multi-stage networks endowed with various facilities of forward logistics (FL) and reverse logistics (RL) at each stage. In FL, finished products are provided to customers via manufacturers, distribution centers (DCs), and retailers. On the other hand, in RL, some of the products returned by customers are recovered after performing classification and recovery processes at the collection center and the recovery center, respectively, and the rest are disposed at the disposal center.
Among the major trends in CLSC design problems, the principle of sustainable development has become popular in conventional literatures [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. According to this principle, economic, environmental, and social factors are usually considered for effectively constructing a sustainable CLSC (SCLSC) design problem.
The economic factors considered include the maximization of the total profit or the minimization of the total cost resulting from construction and operation processes of the CLSC design problem [4,5,6,7,8,16,17,20,22]. Wang and Hsu [5] developed a CLSC design problem to minimize the total cost, which is the aggregate of the total production cost, the total transportation cost, and the total fixed cost. In the CLSC design problem, the supplier, the manufacturer, the DC, and the customer are considered in FL and the dismantler is considered in RL. Following Wang and Hsu [5]. The CLSC design problem by Devika et al. [1] also seeks to minimize the total cost which is the aggregate of the total manufacturing cost, the total transportation cost, the total handling cost, the total purchasing cost, the total opening cost, and the total sale revenue. On the other hand, the CLSC design problem by Chen et al. [6] seeks to maximize the total profit, which is the aggregate of the total revenue (=the total sale revenue of new product) and the total cost (=total processing cost + the total transportation cost + the total fixed cost).
The environmental factors usually considered include the minimization of the total amount or the cost of CO2 or carbon emitted during the production and transportation of materials and products at each stage of the CLSC design problem [3,9,10,18,19,21]. The CLSC design problem by Paksoyet et al. [9] seeks to minimize the total cost of CO2 emissions during the transportation of materials or products at each stage by raw material suppliers, plants, DCs, warehouses, and customers in FL and collection centers, repairing centers, dismantlers, decomposition centers, and final disposals in RL. Talaei et al. [10] showed a simple CLSC design problem to minimize the total amount of CO2 emissions during the establishment of facilities, production and transportation of products, and disposal of waste products at each stage by the manufacturers and markets in FL and disposal centers, collection centers, and inspection centers in RL. Following Talaei et al. [10], Özceylan et al. [3] also showed a CLSC design problem for the automotive industry in Turkey. They minimized the total amount of CO2 emissions during the transportation of finished products and returned products at each stage. The CLSC design problem by Fahimnia et al. [19] seeks to minimize total cost of carbon emitted during the production and transportation of products at each stage by manufacturers, warehouses, and end-users in FL and collection centers, recycling centers and disposal centers in RL.
The relevant social factors usually considered are metrics like the number of job opportunities created by introducing new technology, the number of lost days due to damage to work, and general unemployment [1,2,3,23]. Devika et al. [1] constructed a CLSC design problem to optimize various social influences (e.g., the number of newly created job opportunities, their job types, and damage to work caused either during the establishment of facilities or during the manufacture and handling of products). Özceylan et al. [3] also considered some social influences (e.g., the number of new job opportunities by introducing new technology and the number of lost days due to damage to work) in their CLSC design problem.
As mentioned above, economic, environmental, and social factors have been individually considered in various studies. However, few papers have considered all three factors simultaneously for the SCLSC design problem [1,3,11].
The SCLSC design problem presented by Devika et al. [1] represents a multi-objective optimization problem by considering the three factors to be objective functions. For the economic factor, the problem focused on the minimization of the total cost (=the sum of the costs of manufacture, transportation, handling, purchase and opening − the sum of sale revenues). For the environmental and social factors, the metrics considered were the minimization of the total cost (=the sum of environmental damages − the sum of environmental advantages) and the maximization of the total number (=the sum of the number of the new job opportunity − the sum of the number of lost days due to damage to work), respectively. To solve the multi-objective optimization problem, a meta-heuristic approach using an imperialist competitive algorithm and a variable neighborhood search was used. By performing numerical experiments based on various problem scales, it was proven that the meta-heuristic approach outperforms certain conventional hybrid approaches. However, the search speed of the meta-heuristic approach became slower than those of the conventional hybrid approaches as the problem scales were increased. On the other hand, the SCLSC design problem by Özceylan et al. [3] was represented by a single-objective optimization problem, by using the economic factor as the only objective function, and the remaining two factors (environmental and social factors) as constraints. The maximization of total profit (=the sum of total revenues − the sum of total costs) served as the economic factor. The total amount of CO2 emissions as well as the number of job opportunities created and the number of lost days due to damage to work were considered as the environmental and social factors, respectively.
One of the motivations behind the constructions of CLSC design problems is to distribute materials and products efficiently via the CLSC network. Consideration of various distribution channels reinforces the efficiency of production and distribution in a CLSC design problem. In general, three types of distribution channels (normal delivery: NRD, direct delivery: DRD, and direct shipment: DRS) are considered in such a problem [6]. The NRD is generally recognized to be the most basic distribution channel and is used to send materials or products from one stage to an adjoining one. For example, products are sent from the manufacturer to the DC via the NRD in a CLSC design problem that consists of a supplier, a manufacturer, a DC, a retailer, a customer, a collection center and a recovery center. The distribution channel in which a certain quantity of the products can be sent from DC to the customer directly, without routing it through the retailer, is called a DRD. The channel that allows direct shipment from the manufacturer to the customer, without any involvement of the DC and the retailer is called a DRS.
CLSC design problems with the three types of distribution channels mentioned above have been considered in conventional literature [4,8,13,14,15]. The CLSC design problem by Cardoso et al. [14] consists of a plant, a warehouse, a retailer and a market. Certain products manufactured at the plant are sent to market through the warehouse and the retailer via the NRD, and some others are sent from the plant to the market, without being routed through the warehouse or the retailer, via DRS. The rest of the products at the warehouse are sent to the market via the DRD, without involving the retailer. The CLSC design problem by Son et al. [8] consists of an external supplier, a pre-manufacturing processor, a parts manufacturer, an end-product manufacturer, and a customer in FL, and a collector, a recycler, and a disposal unit in RL. Certain materials and products are sent from one stage to the next via the NRD, some others are sent from the external supplier to the parts manufacturer via the DRD without involving the pre-manufacturing processors. Varsei and Polyakovskiy [15] designed a CLSC design problem for the wine industry in Australia. The CLSC design problem consisted of a supplier, a winery, a bottling plant, a DC and a demand point at each stage. Certain raw materials (e.g., grapes and bottles) and wine products were sent from one stage to the next via the NRD. Some others were sent from the supplier to the bottling plant via the DRD, without routing them through the winery.
As mentioned above, existing works have considered various distribution channels in their CLSC design problems. Unfortunately, most of them, except for Soleimani and Kannan [13], do not use the three distribution channels simultaneously.
Table 1 presents a brief review of conventional works in the literature on SCLSC design problems with regard to three factors of sustainability (economic, environmental, and social factors) and three distribution channels (NRD, DRD, and DRS).
As is evident from Table 1, although there are several works that have considered either all three factors for sustainability or all three types of distributions channels, no work has considered both sets simultaneously.
In this paper, we propose a SCLSC design problem considering all three factors for sustainability and all three types of distribution channels simultaneously. The conceptual structure of the SCLSC design problem is described in Section 2.1 As the three factors for sustainability are used as three separate objective functions, the proposed SCLSC design problem is a multi-objective optimization problem and it can be represented by the mathematical formulation discussed in Section 2.2 To effectively solve the multi-objective optimization problem, a hybrid genetic algorithm (pro-HGA) approach is proposed in Section 2.3 that combines the genetic algorithm (GA) with a cuckoo search (CS). Via the numerical experiments presented in Section 3, the SCLSC design problem is explored on five different scales. The performance of the pro-HGA approach is compared with those of other conventional approaches, and the performances of the SCLSC design problems with and without three types of distribution channel are also compared. Finally, some conclusions are summarized and directions for future study are mentioned in Section 4.

2. Materials and Methods

2.1. Conceptual Structure of SCLSC Design Problem

The conceptual structure of the SCLSC design problem is depicted in Figure 1. The manufacturer produces products and sends a certain proportion of them (α1%) to the DC using the NRD and the remaining products (α2% = 1 − α1%) are sent to customer directly via the DRS. At the DC, a certain proportion (β1%) of the incoming products are sent to the retailer via the NRD and the rest (β2% = 1 − β1%) are sent to customer directly via the DRD. The retailer sends products to the customer via the NRD. A proportion (γ1%) of the products returned by the customer are transported to the recovery center via the collection center and the rest (γ2% = 1 − γ1%) are sent to the disposal center to be disposed. At the recovery center, a certain proportion (δ1%) of the returned products are recycled into materials and then returned to the manufacturer, and the rest (δ2% = 1 − δ1%) are recovered and then are resold to the customer.
Figure 1 also depicts the use of the three factors of sustainability at each stage. The total cost (=the sum of the total fixed costs incurred by the manufacturer, the DC, the retailer, the collection center, the recovery center, and the disposal center + the total transportation costs for all distribution channels + the total handling costs for the manufacturer, the DC, the retailer, the collection center, the recovery center and the disposal center) is considered as the economic factor. The total amount of CO2 emissions during the manufacture of products at the manufacturer, during the recovery of returned products at the recovery center, and during transportation of products between each stage is considered to constitute the environmental factor. The number of job opportunities created, general unemployment, and the number of lost days due to damage to work at the manufacturer are together used as the relevant social factors.

2.2. Mathematical Formulation

The following assumptions are made to implement the proposed SCLSC design problem.
  • Products of a single type are considered.
  • The numbers of facilities considered at each stage at the manufacturer, the DC, the retailer, the customer, the collection center, the recovery center, and the disposal center are fixed and known beforehand.
  • Only one facility is opened at each stage at the manufacturer, the DC, the retailer, the collection center, and the recovery center. However, all facilities available at the site of the customer and the disposal center are always considered open.
  • The costs to operate the facilities considered at each stage at the manufacturer, the DC, the retailer, the collection center, and the recovery center are all constant, different from each other in value, and known beforehand.
  • The unit handling costs of the facilities considered at each stage at the manufacturer, the DC, the retailer, the collection center, and the recovery center are different from each other in value and known beforehand.
  • The unit transportation costs among the manufacturer, the DC, the retailer, the customer, the collection center, the recovery center, and the disposal center are different from each other in value and known beforehand.
  • The unit amount of CO2 emitted during transportation at each stage, and those emitted during manufacture and recovery at the manufacturer and the recovery center, respectively, are different from each other in value and known beforehand.
  • The proposed SCLSC design problem is considered to be in a steady-state situation.
Indices, parameters, and decision variables are defined as follows:
Index Set
m: index of manufacturer; M: set of manufacturer, m M
d: index of DC; D: set of DC, d D
r: index of retailer; R: set of retailer, r R
c: index of customer; C: set of customer, c C
l: index of collection center; L: set of collection center, l L
e: index of recovery center; E: set of recovery center, e E
p: index of disposal center; P: set of disposal center, p P
Parameter
fm: fixed cost at manufacturer m
fd: fixed cost at DC d
fr: fixed cost at retailer r
fl: fixed cost at collection center l
fe: fixed cost at recovery center e
hm: unit handling cost at manufacturer m
hd: unit handling cost at DC d
hr: unit handling cost at retailer r
hl: unit handling cost at collection center l
he: unit handling cost at recovery center e
tsmd: unit transportation cost from manufacturer m to DC d
tsmc: unit transportation cost from manufacturer m to customer c via DRS
tsdr: unit transportation cost from DC d to retailer r
tsdc: unit transportation cost from DC d to customer c via DRD
tsrc: unit transportation cost from retailer r to customer c
tscl: unit transportation cost from customer c to collection center l
tscp: unit transportation cost from customer c to disposal center p
tsle: unit transportation cost from collection center l to recovery center e
tsec: unit transportation cost from recovery center e to customer c
tsem: unit transportation cost from recovery center e to manufacturer m
dmd: distance between manufacturer m and DC d (unit: kilometer)
dmc: distance between manufacturer m and customer c via DRS (unit: kilometer)
ddr: distance between DC d and retailer r (unit: kilometer)
ddc: distance between DC d and customer c via DRD (unit: kilometer)
drc: distance between retailer r and customer c (unit: kilometer)
dcl: distance between customer c and collection center l (unit: kilometer)
dcp: distance between customer c and disposal center p (unit: kilometer)
dle: distance between collection center l and recovery center e (unit: kilometer)
dec: distance between recovery center e and customer c (unit: kilometer)
dem: distance between recovery center e and manufacturer m (unit: kilometer)
wc: weight allocated to the created job opportunity
wl: weight allocated to the number of lost days due to damage to work
wu: weight allocated to unemployment
sem: number of the job opportunities created due to the use of technology t at manufacturer m
slm: number of lost days due to damage to work due to use of technology t at manufacturer m
sum: number of the unemployed workers due to the use of technology t at manufacturer m
CO2V: amount of CO2 emitted from a vehicle per kilometer
CO2m: unit amount of CO2 emitted during production process at manufacturer m
CO2e: unit amount of CO2 emitted during recovery process at recovery center e
vc: capacity that can be shipped in a vehicle
Decision Variables
cm: capacity of manufacturer m
cd: capacity of DC d
cr: capacity of retailer r
cc: capacity of customer c
cl: capacity of collection center l
ce: capacity of recovery center e
cp: capacity of disposal center p
qmd: quantity transported from manufacturer m to DC d
qmc: quantity transported from manufacturer m to customer c via DRS
qdr: quantity transported from DC d to retailer r
qdc: quantity transported from DC d to customer c via DRD
qrc: quantity transported from retailer r to customer c
qcl: quantity transported from customer c to collection center l
qcp: quantity transported from customer c to disposal center p
qle: quantity transported from collection center l to recovery center e
qec: quantity transported from recovery center e to customer c
qem: quantity transported from recovery center e to manufacturer m
xm: takes the value 1 if manufacturer m is opened and 0 otherwise
xd: takes the value 1 if DC d is opened and 0 otherwise
xr: takes the value 1 if retailer r is opened and 0 otherwise
xl: takes the value 1 if collection center l is opened and 0 otherwise
xe: takes the value 1 if recovery center e is opened and 0 otherwise
tm: takes the value 1 if technology t is used at manufacturer m and 0 otherwise
The first objective function (F1), corresponding to the economic factor, is defined as follows to minimize the total cost. The total cost is the aggregate of the total fixed cost, the total handling cost, and the total transportation cost.
M i n .   F 1   =   T o t a l   f i x e d   c o s t   ( T F C )   +   T o t a l   h a n d l i n g   c o s t   ( T H C )   +   T o t a l   t r a n s p o r t a t i o n   c o s t   ( T T C )
T F C   =   ( m f m · x m ) + ( d f d · x d ) + ( r f   r · x r ) + ( l f   l · x l ) + ( e f e · x e )
T H C = ( m h m · c m · x m ) + ( d h d · c d · x d ) + ( r h r · c r · x r ) + ( l h l · c l · x l ) + ( e h e · c e · x e )
T T C = ( m d t s m d · ( α 1 % · c m ) · x m · x d ) + ( m c t s m c · ( α 2 % · c m ) · x m ) + ( d r t s d r · ( β 1 % · c d ) · x d · x r ) + ( d c t s d c · ( β 2 % · c d ) · x d ) + ( c l t s c l · ( γ 1 % · c c ) · x l ) + ( c p t s c p · ( γ 2 % · c c ) ) + ( l e t s l e · c l · x l · x e ) + ( e m t s e m · ( δ 1 % · c e ) · x e · x m ) + ( e c t s e c · ( δ 2 % · c e ) · x e )
The term TFC in Equation (2) is calculated as the sum of the fixed costs when the facilities at all manufacturers, DCs, retailers, collection centers and recovery centers are opened. The term THC in Equation (3) is the sum of the handling costs incurred to handle products at the facilities opened at manufacturers, DCs, retailers, collection centers, and recovery centers. The term TTC in Equation (4) is the sum of the transportation costs of products from one stage to another.
The second objective function (F2) corresponding to the environmental factor is defined as follows to minimize the total amount of CO2 emissions during production, recovery and transportation processes.
M i n .   F 2   =   T o t a l   a m o u n t   o f   C O 2   e m i t t e d   d u r i n g   p r o d u c t i o n   a n d   r e c o v e r y   p r o c e s s e s   ( T C P )   +   T o t a l   a m o u n t   o f   C O 2   e m i t t e d   d u r i n g   t r a n s p o r t a t i o n   p r o c e s s e s   ( T C T )
T C P   =   ( m c m ·   x m · C O 2 m )   + ( e c e ·   x e · C O 2 e )
T C T = C O 2 V · [ ( m d d m d · x m · x d · c m · α 1 % v c ) + ( m c d m c · x m · c m · α 2 % v c ) + ( d r d d r · x d · x r · c d · β 1 % v c ) + ( d c d d · x d · c d · β 2 % v c ) + ( r c d r c · x r · c r v c ) + ( c l d c l · x l · c l · γ 1 % v c ) + ( c p d c p · c c · γ 2 % v c ) + ( l e d l e · x l · x e · c l v c ) + ( e m d e m · x e · x m · c e · δ 1 % v c ) + ( e c d e c · x e · c e · δ 2 % v c ) ]
The term TCP in Equation (6) denotes the total amount of CO2 emitted during production and recovery processes at manufacturers and recovery centers. The term TCT denotes the total amount of CO2 emitted during transportation between each pair of stages. In Equation (7), the parameters values of vc and CO2V are randomly generated according to the methods described in [3,9].
The third objective function (F3) corresponding to the social factor is defined as follows to maximize the social influence, which consists of (i) the number of job opportunities created, (ii) the number of lost days due to damage to work, and (iii) the number of unemployed workers at the manufacturers. The three components are already mentioned and studied by many conventional literatures [24,25,26,27]. Equation (8) defines F3.
M a x .   F 3 = ( w e i g h t   a l l o c a t e d   t o   c r e a t e d   j o b   o p p o r t u n i t y · n u m b e r   o f   c r e a t e d   j o b o p p o r t u n i t y   a t   m a n u f a c t u r e r ) ( w e i g h t   a l l o c a t e d   t o   l o s t   d a y   c a u s e d   b y   w o r k s d a m a g e · n u m b e r   o f   l o s t   d a y   c a u s e d   b y   w o r k s   d a m a g e   a t   m a n u f a c t u r e r ) ( w e i g h t a l l o c a t e d   t o   u n e m p l o y m e n t · n u m b e r   o f   u n e m p l o y m e n t   a t   m a n u f a c t u r e r ) = ( w c · m s e m · x m · t m ) ( w l · m s l m · x m · t m ) ( w u · m s u m · x m )
The three aforementioned objective functions should be optimized with respect to the following constraints.
( α 1 % · m d q m d · x m · x d ) ( d c d · x d ) 0  
( α 1 % · m d q m d · x m · x d ) ( d c d · x d ) 0  
( β 1 % · d r q d r · x d · x r ) ( r c r · x r ) 0
( β 2 % · d c q d c · x d ) ( c c c ) 0
( r c q r c · x r ) ( c c c ) 0
( γ 1 % · c l q c l · x l ) ( l c l · x l ) 0
( γ 2 % · c p q c p ) ( p c p ) 0
( l e q l e · x l · x e ) ( e c e · x e ) 0
( δ 1 % · e m q e m · x e · x m ) ( m c m · x m ) 0
( δ 1 % · e c q e c · x e ) ( c c c ) 0
m x m = 1
d x d = 1
r x r = 1
l x l = 1
e x e = 1
x m = { 0 , 1 } ,     m M
x d = { 0 , 1 } ,     d D
x r = { 0 , 1 } ,     r R
x l = { 0 , 1 } ,     l L
x e = { 0 , 1 } ,     e E
c m , c d , c r , c c , c l , c e , c p 0 ,   m M , d D , r R , c C , l L , e E , p P
Equation (9) indicates that the amount of transported goods from a manufacturer to a DC is less than or equal to the capacity of the DC. In a similar vein, inequalities (10) to (18) represent the quantity limitation for transportation between each pair of stages. Equations (19)–(23) impose the restriction that only one facility may be opened at each stage at each manufacturer, DC, retailer, collection center, and recovery center. Relations (24)–(28) represent the fact that each decision variable should take the values of 0 or 1. Inequality (29) represents non-negativity of the variables concerned.

2.3. Proposed HGA Approach

Most complicated network problems are known to be NP-complete [16,28,29], and the SCLSC design problem proposed in this paper is no exception. Existing works in the literature [4,29,30,31] have shown that meta-heuristics-based approaches such as GA, CS, and Tabu search (TS) have been adapted to solve complicated network problems efficiently. However, there exist many situations in which most conventional single-based meta-heuristics approaches do not perform particularly well [28,31,32]. To address this gap, various hybrid approaches have been developed that use GA or other similar approaches [29,31,32,33,34,35,36,37].
Zhang et al. [35] and Xinyu and Liang [36] developed the hybrid GA (HGA) approaches that combine GA with TS to reinforce the hybrid search ability by utilizing GA for global searches and TS for local searches. However, their HGA approaches suffer from weaknesses in applications to complicated network problems. First, the search schemes used in their HGA approaches require considerable computation time, as all individuals in a current population are adapted to produce the TS list and this process is time-consuming. Second, the search quality for locating global optimal solutions via GA sometimes deteriorates, because the Tabu list in TS is produced based on the initial population. Lin et al. [33] developed an HGA approach that utilizes GA for global searches and the iterative hill climbing method for local searches. The approach was adapted to solve complicated network problems of several sizes. Although the authors proved the superior efficiency of the approach compared to that of a conventional GA, its search speed is significantly slower than that of GA, especially for problems of large size.
Therefore, this paper proposes an improved HGA approach, called the pro-HGA approach, that reinforces search speed and search quality to solve the SCLSC design problem efficiently. The pro-HGA approach combines GA for global searches with CS for local searches. First, a new population is produced by using three operators (crossover, mutation, and selection) of GA. Then, CS is applied to the new population produced by GA to locate more respective solutions. To enhance the search speed and search quality of the pro-HGA approach, the respective solutions are continuously produced, and they are used to carry out continuous improvement during GA and CS search processes. For this strategy, the Lévy flight scheme [38] of CS is adapted to all solutions of the population produced during the GA search process. The implementation procedure of the pro-HGA approach is depicted in detail as follows.
procedure: pro-HGA approach
input: problem data, parameters
output: Pareto optimal solutions
begin
t ← 0  //t: generation number
 initialize parent population PP(t) by encoding routine;
 calculate each objective function Fi, i = 1, 2, 3 of PP(t) by decoding routine;
 create Pareto optimal solutions E(P) by non-dominated routine;
while (t < max generation)
  produce offspring OP(t) from PP(t) by adapting 2X crossover operator [30] and random mutation operator [30];
  calculate Fi of OP(t) by decoding routine;
  find current E(P) by non-dominated routine;
  keep best solution set GLbest using current E(P);
  for each solution xi of OP(t) do
   generate a new solution xnew from xi by adapting Lévy flight scheme [38];
   randomly select another solution xi in OP(t);
   if (F(xnew) > F(xi)) then CP(t) ← xnew//CP(t): CS population
  end for
  worst solutions with fraction rate (fr) are abandoned;
  randomly regenerate new solutions xr_new as many as fr;
  CP(t) ← xr_new
  calculate Fi of CP(t) by decoding routine;
  find current E(P) by non-dominated routine;
  keep best solution set CLbest using current E(P);
  if (F(GLbest) > F(CLbest))
   then update E(P) using GLbest by non-dominated routine
   else update E(P) using CLbest by non-dominated routine
  end if
   reproduce PP(t + 1) using OP(t) and CP(t) by adapting elitist selection scheme [30];
   tt + 1;
end
output Pareto optimal solutions E(P);
end;

3. Results and Discussions

Via numerical experiments, the mathematical formulation of the SCLSC design problem suggested in Section 2.2 is implemented at five different scales as presented in Table 2.
To evaluate the performance of the pro-HGA approach comparatively, two conventional approaches (GA: conventional GA by Gen and Cheng [28], HGA: conventional HGA by Kanagaraj et al. [39]) are chosen. All approaches were programmed using MATLAB Version 2014b and run under the same computational environment (IBM compatible PC 1.3 GHz processor-Intel core I5-1600 CPU, 4GB RAM, and OS-X EI).
The parameter settings to implement the GA, HGA, and pro-HGA approaches are as follows: Total number of generations is 1000, population size is 20, crossover rate is 0.5, and mutation rate is 0.3. The number of host nests is 20, α = 1 and p a = 0.25 for the search of CS used in the HGA and pro-HGA approaches. These parameter values were obtained by fine-tuning the procedure of each approach. The parameter values to implement the mathematical formulation in Section 2.2 are defined as follows: CO2V = 3, CO2m = 0.5, CO2e = 0.3, vc = 10, wc = 0.15, wl = 0.15, wu = 0.15, sem = U [90, 100], slm = U [20, 30], sum = U [15, 20]. To eliminate the randomness in the implementation of each approach, each method is evaluated over 30 independent runs. Various measures of performance are used to compare the performances of the approaches, as recorded in Table 3.
The values of |Sj|, RNDS(Sj), and DIR(Sj) presented in Table 3 are the averages obtained over 30 independent runs. The three sub-problems generated from the original problem which had been suggested in Section 2.2 are used to aid the analyses of computation results as follows [31]. The experimental results are summarized in Table 4, Table 5, Table 6, Table 7 and Table 8. The boldfaces in Table 4, Table 5, Table 6, Table 7 and Table 8 indicate the best solutions for each measure.
  • Problem 1: min. F1 and min. F2
  • Problem 2: min. F1 and max. F3
  • Problem 3: min. F2 and max. F3
In Table 4, the number of Pareto optimal solutions of Problem 1 obtained via the GA approach is greater than that obtained via the HGA and pro-HGA approaches for the same problem in terms of the values |Sj|. However, in Problem 2, the number of the Pareto optimal solutions obtained via the pro-HGA approach is the same as that obtained via the HGA approach, and they are both greater than the number of solutions obtained via the GA approach. In Problem 3, the pro-HGA approach locates two Pareto optimal solutions, but the GA and HGA approaches fail to locate any. The values of |Sj| also influence the rate of Pareto optimal solutions in S*. In other words, in terms of RNDS(Sj), the rate corresponding to the GA approach is higher than those corresponding to the HGA and pro-HGA approaches in Problem 1. In contrast, the rates corresponding to the HGA and pro-HGA approaches are higher than that corresponding to the GA approach in Problems 2 and 3. In terms of DIR(Sj), the pro-HGA approach exhibits significantly better results than the GA and HGA approaches in Problems 1 and 3, which can be inferred from the fact that the average distance between the Pareto optimal solutions and elements of S* is inversely proportional to the quality of the approach.
According to the data presented in Table 5, the pro-HGA approach performs significantly better in terms of |Sj| compared to the GA and HGA approaches in Problems 1, 2, and 3. Similar results are observed in terms of RNDS(Sj) and DIR(Sj). The computation results based on |Sj|, RNDS(Sj), and DIR(Sj) indicate that the search scheme used in the pro-HGA approach is more efficient than those used in the GA and HGA approaches.
According to the data presented in Table 6, the pro-HGA approach outperforms the GA and HGA approaches in terms of |Sj| and RNDS(Sj) in Problems 1 and 3, but not for Problem 2. On the other hand, in terms of DIR(Sj), the performance of the GA approach is slightly better than that of the HGA and pro-HGA approaches in Problem 1. However, for Problems 2 and 3, the performance of the pro-HGA approach is superior to those of the others.
In terms of the average CPU times required, as recorded in Table 4, Table 5 and Table 6, no significant differences among the approaches are observed, even though the computation speeds of the GA approach were slightly quicker than those of the HGA and pro-HGA approaches.
In the data corresponding to Problems 1, 2, and 3 presented in Table 7, the pro-HGA approach performs the best in terms of |Sj|, RNDS(Sj), and DIR(Sj) compared to the GA and HGA approaches. A similar observation also follows from the results tabulated in Table 8 in all the Problems, the performances of the pro-HGA approach are significantly superior to those of the GA and HGA approaches in terms of |Sj|, RNDS(Sj), and DIR(Sj). The search speeds of the GA approach were observed to be slightly quicker than those of the HGA and pro-HGA approaches, but no significant difference was noticed.
Figure 2, Figure 3 and Figure 4 depict the Pareto optimal solutions of the GA, HGA, and pro-HGA approaches in Problems 1, 2, and 3 at Scale 5.
From Figure 2, it is evident that five Pareto optimal solutions obtained via the pro-HGA approach coincide with the S*, and that the GA and HGA approaches do not show any advantages compared to it in this regard. Similar results are observed regarding the number of Pareto optimal solutions in terms of |Sj| for Problem 1 as presented in Table 8. Similar situations are also depicted in Figure 3 and Figure 4. The pro-HGA approach exhibits two Pareto optimal solutions in Figure 3 and three solutions in Figure 4 that coincide with those of S*, but none of the solutions obtained via the GA and HGA approaches coincide in position with elements of S*.
Based on the experimental results presented in Table 4, Table 5, Table 6, Table 7 and Table 8 and Figure 2, Figure 3 and Figure 4, we can reach the following conclusion.
  • The pro-HGA approach has performed better than the GA and HGA approaches in terms of the metrics |Sj|, RNDS(Sj), and DIR(Sj), as the problem sizes were increased from scales 1 to scale 5. However, the former does not exhibit any advantage in terms of the CPU times required compared to the latter pair. With respect to the Pareto optimal solutions obtained via each approach in comparison with the S*, the performance of the pro-HGA approach is superior to those of the GA and HGA approaches, which proves that the search scheme used in the pro-HGA approach is more efficient than those used in the GA and HGA approaches.
As recorded in Table 1, the SCLSC design problem proposed in this paper utilizes various distribution channels (i.e., NRD, DRD, and DRS) and this feature distinguishes it from methods presented in other works in the literature. Therefore, we can compare the performances of the following SCLSC design problems.
  • . SCLSC-S: the SCLSC design problem with a single distribution channel, i.e., NRD
  • . SCLSC-V: the SCLSC design problem with various distribution channels, i.e., NRD, DRD, and DRS
Table 9 and Table 10 present the computation results for SCLSC-S and SCLSC-V at scales 4 and 5. The boldfaces in Table 9 and Table 10 indicate the best solutions for each measure.
From the data presented in case of Problem 1 in Table 9, the performances of SCLSC-V are observed to be significantly better than those of the SCLSC-S in terms of the metrics |Sj|, RNDS(Sj), and DIR(Sj). In Problem 2, SCLSC-V again exhibits slightly better performances than SCLSC-S in terms of |Sj| and RNDS(Sj), but not for DIR(Sj). However, in Problem 3, SCLSC-V no longer exhibits any comparative merits in terms of |Sj| and RNDS(Sj). Similar results are gleaned by analyzing the results presented in Table 10. In Problems 1 and 3, SCLSC-V significantly outperforms SCLSC-S in terms of |Sj|, RNDS(Sj), and DIR(Sj), even though the latter performs slightly better in terms of |Sj| and RNDS(Sj) in Problem 2. In Problems 1, 2 and 3 of Table 9 and Table 10, the search speeds of SCLSC-S are slightly quicker than those of SCLSC-V as is apparent from the average CPU times required.
Figure 5, Figure 6 and Figure 7 depict the Pareto optimal solutions for SCLSC-S and SCLSC-V with respect to Problems 1, 2, and 3 at Scale 5.
From Figure 5, it is evident that four Pareto optimal solutions out of five coincide with S* for SCLSC-V, while the same is true for two out of seven solutions for SCLSC-S. These results are the same as the distributions of Pareto optimal solutions in terms of |Sj| for Problem 1 as presented in Table 10. From Figure 6, we can observe that all Pareto optimal solutions for SCLSC-V and SCLSC-S coincide with S*. However, in Figure 7, all Pareto optimal solutions for SCLSC-V coincide with S*, but none of the solutions for the SCLSC-S do.
Based on the experimental results presented in Table 9 and Table 10 and Figure 5, Figure 6 and Figure 7, the following conclusion can be reached.
  • SCLSC-V outperforms SCLSC-S with respect to most performance measures, although the former has slightly lower search speeds than the latter. This indicates that the SCLSC design problem with various distribution channels, i.e., NRD, DRD, and DRS is more efficient in locating optimal solutions than that with a single distribution channel, i.e., NRD.

4. Conclusions

In this paper, we have proposed an SCLSC design problem. Three factors (economic, environmental, and social factors) and three distribution channels (NRD, DRD, and DRS) have been considered for the construction of the problem. The economic factor is concerned with minimizing the total cost, which is the aggregate of the total fixed cost, the total handling cost, and the total transportation cost. The environmental factor also is concerned with minimizing the total amount of CO2 emissions during the transportation and manufacturing processes. On the other hand, the social factor is concerned with maximizing the social influence of the process, which consists of the number of job opportunities created, the number of unemployed workers and the number of lost days due to damage to work. As these three factors have been considered in the proposed SCLSC design problem simultaneously, the problem has been formulated as a multi-objective optimization problem. A pro-HGA approach that combines GA with CS has been proposed to solve the SCLSC design problem.
The SCLSC design problem has been presented in five different scales with the help of numerical experiments and the performance of the pro-HGA approach has been compared with those of two conventional approaches (GA and HGA). Various measures of performance, such as Pareto optimal solutions, have been used for the comparison. The experimental results have conclusively shown that the pro-HGA approach is more efficient than the GA and HGA approaches with respect to various measures of performance, although the former shows no advantage over the latter with respect to the average CPU time required to complete the process.
SCLSC design problems with and without various distribution channels have also been compared in the paper. SCLSC-S (the SCLSC design problem with single distribution channel, i.e., NRD) and SCLSC-V (the SCLSC design problem with various distribution channels, i.e., NRD. DRD, and DRS) have been compared at two separate scales. Experimental results demonstrated that SCLSC-V is more efficient in locating optimal solutions than SCLSC-S with respect to most performance measures, except the average CPU time required.
With regard to the direction of our future studies, we wish to present SCLSC design problems of a larger scale using real field data to compare the performances of SCLSC-V and SCLSC-S. We also wish to consider other meta-heuristic approaches, such as particle swam optimization, ant colony optimization, etc. to construct a more efficient pro-HGA approach.

Author Contributions

Investigation, Y.Y.; Methodology, Y.Y. and M.G.; Software, A.C.; Validation, Y.Y.; Writing—original draft, Y.Y. and M.G.; Review and editing, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Devika, K.; Jafarian, A.; Nourbakhsh, V. Designing a sustatinable closed-loop supply chain network based on triple bottom line approach: A comparison of metaheuristics hybridixation techniques. Eur. J. Oper. Res. 2014, 235, 594–615. [Google Scholar] [CrossRef]
  2. Eskandarpour, M.; Dejax, P.; Miemczyk, J.; Péton, O. Sustainable supply chain network design: An optimization-oriented review. Omega 2015, 54, 11–31. [Google Scholar] [CrossRef]
  3. Özceylan, E.; Demirel, N.; Çetinkaya, C.; Demirel, E. A closed-loop supply chain network design for automotive industry in Turkey. Comput. Ind. Eng. 2017, 113, 729–745. [Google Scholar] [CrossRef]
  4. Min, H.; Ko, C.S.; Ko, H.J. The spatial and temporal consolidation of returned products in a closed-loop supply chain network. Comput. Ind. Eng. 2006, 51, 309–320. [Google Scholar] [CrossRef]
  5. Wang, H.F.; Hsu, H.W. A closed-loop logistic model with a spanning tree based genetic algorithm. Comput. Oper. Res. 2010, 37, 376–389. [Google Scholar] [CrossRef]
  6. Chen, Y.T.; Chan, F.T.S.; Chung, S.H. An integrated closed-loop supply chain model with location and allocation problem and product recycling decisions. Int. J. Prod. Res. 2015, 53, 3120–3140. [Google Scholar] [CrossRef]
  7. Jabbarzadeh, A.; Haughton, M.; Khosrojerdi, A. Closed-loop supply chain network design under distription risks: A robust approach with real world application. Comput. Ind. Eng. 2018, 116, 178–191. [Google Scholar] [CrossRef]
  8. Son, D.; Kim, S.; Park, H.; Jeong, B. Closed-loop supply chain planning model of rare metals. Sustainability 2018, 10, 1061. [Google Scholar] [CrossRef] [Green Version]
  9. Paksoy, T.; Bektaş, T.; Özceylan, E. Operational and environmental performance measures in a multi-product closed-loop supply chain. Transp. Res. Part E 2011, 47, 532–546. [Google Scholar] [CrossRef]
  10. Talaei, M.; Moghaddam, B.F.; Pishvaee, M.S.; Bozorgi-Amiri, A.; Gholamnejad, S.A. A robust fuzzy optimization model for carbon-efficient closed-loop supply chain network design problem: A numerical illustration in electronics industry. J. Clean. Prod. 2016, 113, 662–673. [Google Scholar] [CrossRef]
  11. Zhalechian, M.; Tavakkoli-Moghaddam, R.; Zahiri, B.; Mohammadi, M. Sustainable design of a closed-loop location-routing-inventory supply chain network under mixed uncertainty. Transp. Res. Part E 2016, 89, 182–214. [Google Scholar] [CrossRef]
  12. Sahebjamnia, N.; Fathollahi-Fard, A.M.; Hajiaghaei-Keshteli, M. Sustainable tire closed-loop supply chain network design: Hybrid metaheuristic algorithms for large-scale networks. J. Clean. Prod. 2018, 196, 273–296. [Google Scholar] [CrossRef]
  13. Soleimani, H.; Kannan, G. A hybrid particle swarm optimization and genetic algorithm for closed-loop supply chain network design in large-scale network. Appl. Math. Model. 2015, 39, 3990–4012. [Google Scholar] [CrossRef]
  14. Cardoso, S.R.; Barbosa-Póvoa, A.P.; Relvas, S. Integrating financial risk measures into the design and planning of closed-loop supply chains. Comput. Chem. Eng. 2016, 85, 105–123. [Google Scholar] [CrossRef]
  15. Varsei, M.; Polyakovskiy, S. Sustainable supply chain network design: A case of the wine industry in Australia. Omega 2017, 66, 236–247. [Google Scholar] [CrossRef]
  16. Savaskan, R.C.; Bhattacharya, S.; Van Wassenhove, L.V. Closed-loop supply chain models with product remanufacturing. Manag. Sci. 2004, 50, 239–252. [Google Scholar] [CrossRef] [Green Version]
  17. Lee, Y.J.; Baker, T.; Jayaraman, V. Redesigning an integrated forward-reverse logistics system for a third party service provider: An empirical study. Int. J. Prod. Res. 2012, 50, 5615–5634. [Google Scholar] [CrossRef]
  18. Özkir, V.; Başlıgıl, H. Modelling product-recovery processes in closed-loop supply-chain network design. Int. J. Prod. Res. 2012, 50, 2218–2233. [Google Scholar] [CrossRef]
  19. Fahimnia, B.; Sarkis, J.; Dehghanian, F.; Banihashemi, N.; Rahman, S. The impact of carbon pricing on a closed-loop supply chain: An Australian case study. J. Clean. Prod. 2013, 59, 210–225. [Google Scholar] [CrossRef]
  20. Faccio, M.; Persona, A.; Sgarbossa, F.; Zanin, G. Sustainable SC through the complete reprocessing of end-of-life products by manufacturers: A traditional versus social responsibility company perspective. Eur. J. Oper. Res. 2014, 233, 359–373. [Google Scholar] [CrossRef]
  21. Jindal, A.; Sangwan, K.S. Multi-objective fuzzy mathematical modelling of closed-loop supply chain considering economic and environmental factors. Ann. Oper. Res. 2017, 257, 95–120. [Google Scholar] [CrossRef]
  22. Jerbia, R.; Boujelben, M.K.; Sehli, M.A.; Jemai, Z. A stochastic closed-loop supply chain network design problem with multiple recovery options. Comput. Ind. Eng. 2018, 118, 23–32. [Google Scholar] [CrossRef]
  23. Lie, 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] [Green Version]
  24. Pishvaee, M.; Razmi, J.; Torabi, S. Robust possibilistic programming for socially responsible supply chain network design: A new approach. Fuzzy Sets Syst. 2012, 206, 1–20. [Google Scholar] [CrossRef]
  25. Pishvaee, M.; Razmi, J.; Torabi, S. An accelerated Benders decomposition algorithm for sustainable supply chain network design under uncertainty: A case study of medical needle and syringe supply chain. Transp. Res. Part E 2014, E67, 14–28. [Google Scholar] [CrossRef]
  26. Pishvaee, M.; Zanjirani, F.R.; Dullaert, W. A memetic algorithm for bi- objective integrated forward/reverse logistics network design. Comput. Oper. Res. 2010, 37, 1100–1112. [Google Scholar] [CrossRef]
  27. Saffari, M.M.; Shakouri, G.H.; Razmi, J. A new multi objective optimization model for designing a green supply chain network under uncertainty. Int. J. Ind. Eng. Comput. 2015, 6, 15–32. [Google Scholar] [CrossRef] [Green Version]
  28. Gen, M.; Cheng, R. Genetic Algorithms and Engineering Optimization; John-Wiley & Sons: New York, NY, USA, 2000. [Google Scholar]
  29. Gen, M.; Lin, L.; Yun, Y.S.; Inoue, H. Recent advances in hybrid priority-based genetic algorithms for logistics and SCM network design. Comput. Ind. Eng. 2018, 115, 394–412. [Google Scholar] [CrossRef]
  30. Gen, M.; Cheng, R. Genetic Algorithms and Engineering Design; John Wiley & Son: New York, NY, USA, 1997. [Google Scholar]
  31. Yun, Y.S.; Chung, H.S.; Moon, C.U. Hybrid genetic algorithm approach for precedence-constrained sequencing problem. Comput. Ind. Eng. 2013, 65, 137–147. [Google Scholar] [CrossRef]
  32. Yun, Y.S. Hybrid genetic algorithm with adaptive local search scheme. Comput. Ind. Eng. 2006, 51, 128–141. [Google Scholar] [CrossRef]
  33. Lin, L.; Gen, M.; Wang, X. Integrated multistage logistics network design by using hybrid evolutionary algorithm. Comput. Ind. Eng. 2009, 56, 854–873. [Google Scholar] [CrossRef]
  34. Yun, Y.S.; Gen, M.; Hwang, R.K. Adaptive genetic algorithm to multi-stage reverse logistics network design for product resale. Inf. Inte. Interdiscip. J. 2012, 15, 6117–6138. [Google Scholar]
  35. Zhang, Q.; Manier, H.; Manier, M.A. A genetic algorithm with Tabu search procedure for flexible job shop scheduling with transportation constraints and bounded processing times. Comput. Oper. Res. 2012, 39, 1713–1723. [Google Scholar] [CrossRef]
  36. Xinyu, L.; Liang, G. An effective hybrid genetic algorithm and Tabu search for flexible job shop scheduling problem. Int. J. Prod. Econ. 2016, 174, 93–110. [Google Scholar]
  37. Gen, M.; Zhang, W.; Lin, L.; Yun, Y.S. Recent advance in hybrid evolutionary algorithms for multiobjective manufacturing scheduling. Comput. Ind. Eng. 2017, 112, 616–633. [Google Scholar] [CrossRef]
  38. Yang, X.S.; Deb, S. Cuckoo search via Levy flights. In Proceedings of the World Congress on Nature and Biologically Inspired Computing (NaBIC 2009), Coimbatore, India, 9–11 December 2009; pp. 210–214. [Google Scholar]
  39. Kanagaraj, G.; Ponnambalam, S.G.; Jawahar, N. A hybrid cuckoo search and genetic algorithm for reliability-redundancy allocation problems. Comput. Chem. Eng. 2013, 66, 1115–1125. [Google Scholar] [CrossRef]
  40. Ishibushi, H.; Yoshida, T.; Murata, T. Balance between genetic search and local search in memetic algorithms for multiobjective permutation flowshop scheduling. IEEE Trans. Evol. Comput. 2003, 7, 204–223. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Conceptual structure of the proposed SCLSC design problem.
Figure 1. Conceptual structure of the proposed SCLSC design problem.
Mathematics 08 00084 g001
Figure 2. Pareto optimal solutions for each approach in Problem 1 at Scale 5.
Figure 2. Pareto optimal solutions for each approach in Problem 1 at Scale 5.
Mathematics 08 00084 g002
Figure 3. Pareto optimal solutions for each approach in Problem 2 at Scale 5.
Figure 3. Pareto optimal solutions for each approach in Problem 2 at Scale 5.
Mathematics 08 00084 g003
Figure 4. Pareto optimal solutions for each approach in Problem 3 at Scale 5.
Figure 4. Pareto optimal solutions for each approach in Problem 3 at Scale 5.
Mathematics 08 00084 g004
Figure 5. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problem 1 at Scale 5.
Figure 5. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problem 1 at Scale 5.
Mathematics 08 00084 g005
Figure 6. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problems at Scale 5.
Figure 6. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problems at Scale 5.
Mathematics 08 00084 g006
Figure 7. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problems 3 at Scale 5.
Figure 7. Pareto optimal solutions for SCLSC-S and SCLSC-V in Problems 3 at Scale 5.
Mathematics 08 00084 g007
Table 1. A brief review of conventional works in the literature on SCLSC design problems with regard to three factors for sustainability and three distribution channels.
Table 1. A brief review of conventional works in the literature on SCLSC design problems with regard to three factors for sustainability and three distribution channels.
Three Factors for SustainabilityThree Types of Distribution Channels
Eco.Env.Soc.NRDDRDDRS
Savaskan et al. [16]
Min et al. [4]
Wang and Hsu [5]
Paksoy et al. [9]
Lee et al. [17]
O¨zkir & Başlıgıl [18]
Fahimnia et al. [19]
Devika et al. [1]
Faccio et al. [20]
Chen et al. [6]
Soleimani & Kannan [13]
Cardoso et al. [14]
Talaei et al. [10]
Zhalechain et al. [11]
JinDal & Sangwan [21]
Özceylan et al. [3]
Sahebjamnia et al. [12]
Jabbarzadeh et al. [7]
Jerbia et al. [22]
Lie et al. [23]
Son et al. [8]
This study
Table 2. The SCLSC design problem at five different scales.
Table 2. The SCLSC design problem at five different scales.
ScaleNumber of ManufacturersNumber of DCsNumber of RetailersNumber of CustomersNumber of Collection CentersNumber of Recovery CentersNumber of Disposal Centers
15881851
2253030130251
3506060160501
410080801801001
515010010011001501
Table 3. Measures of performance.
Table 3. Measures of performance.
MeasureDescription
|Sj|Number of Pareto optimal solutions in reference solution set (S*) [40]
RNDS(Sj)Rates of Pareto optimal solutions in the S* [40]
DIR(Sj)Average distance between Pareto optimal solutions and the set S* [40]
CPU time (sec.)Average CPU time required for each run
Table 4. Computation results for each approach in Scale 1.
Table 4. Computation results for each approach in Scale 1.
Scale 1
Problem 1Problem 2Problem 3
MeasureGAHGApro-HGAGAHGApro-HGAGAHGApro-HGA
|Sj|211022002
RNDS(Sj)0.500.250.250.000.500.500.000.001.00
DIR(Sj)574171541524601096496450
CPU time13.714.014.113.714.014.113.714.014.1
Table 5. Computation results for each approach in Scale 2.
Table 5. Computation results for each approach in Scale 2.
Scale 2
Problem 1Problem 2Problem 3
MeasureGAHGApro-HGAGAHGApro-HGAGAHGApro-HGA
|Sj|104003003
RNDS(Sj)0.200.000.800.000.001.000.000.001.00
DIR(Sj)8472560828437012865850
CPU time13.113.414.613.113.414.613.113.413.4
Table 6. Computation results for each approach in Scale 3.
Table 6. Computation results for each approach in Scale 3.
Scale 3
Problem 1Problem 2Problem 3
MeasureGAHGApro-HGAGAHGApro-HGAGAHGApro-HGA
|Sj|025122003
RNDS(Sj)0.000.290.710.200.400.400.000.001.00
DIR(Sj)2413943736022922465234400
CPU time13.213.514.513.213.514.513.213.514.5
Table 7. Computation results for each approach in Scale 4.
Table 7. Computation results for each approach in Scale 4.
Scale 4
Problem 1Problem 2Problem 3
MeasureGAHGApro-HGAGAHGApro-HGAGAHGApro-HGA
|Sj|007003002
RNDS(Sj)0.000.001.000.000.001.000.000.001.00
DIR(Sj)359406046923506985730
CPU time41.342.743.641.342.743.641.342.743.6
Table 8. Computation results for each approach in Scale 5.
Table 8. Computation results for each approach in Scale 5.
Scale 5
Problem 1Problem 2Problem 3
MeasureGAHGApro-HGAGAHGApro-HGAGAHGApro-HGA
|Sj|005002003
RNDS(Sj)0.000.001.000.000.001.000.000.001.00
DIR(Sj)569557070748306706840
CPU time41.842.144.041.842.144.041.842.144.0
Table 9. Computation results for SCLSC-S and SCLSC-V at Scale 4.
Table 9. Computation results for SCLSC-S and SCLSC-V at Scale 4.
Scale 4
Problem 1Problem 2Problem 3
MeasureSCLSC-SSCLSC-VSCLSC-SSCLSC-VSCLSC-SSCLSC-V
|Sj|181232
RNDS(Sj)0.110.890.330.660.60.4
DIR(Sj)12,6700940112400
CPU time32.843.632.843.632.843.6
Table 10. Computation results for SCLSC-S and SCLSC-V at Scale 5.
Table 10. Computation results for SCLSC-S and SCLSC-V at Scale 5.
Scale 5
Problem 1Problem 2Problem 3
MeasureSCLSC-SSCLSC-VSCLSC-SSCLSC-VSCLSC-SSCLSC-V
|Sj|243203
RNDS(Sj)0.330.660.600.400.001.00
DIR(Sj)88997870013,6540
CPU time33.044.033.044.033.044.0

Share and Cite

MDPI and ACS Style

Yun, Y.; Chuluunsukh, A.; Gen, M. Sustainable Closed-Loop Supply Chain Design Problem: A Hybrid Genetic Algorithm Approach. Mathematics 2020, 8, 84. https://doi.org/10.3390/math8010084

AMA Style

Yun Y, Chuluunsukh A, Gen M. Sustainable Closed-Loop Supply Chain Design Problem: A Hybrid Genetic Algorithm Approach. Mathematics. 2020; 8(1):84. https://doi.org/10.3390/math8010084

Chicago/Turabian Style

Yun, YoungSu, Anudari Chuluunsukh, and Mitsuo Gen. 2020. "Sustainable Closed-Loop Supply Chain Design Problem: A Hybrid Genetic Algorithm Approach" Mathematics 8, no. 1: 84. https://doi.org/10.3390/math8010084

APA Style

Yun, Y., Chuluunsukh, A., & Gen, M. (2020). Sustainable Closed-Loop Supply Chain Design Problem: A Hybrid Genetic Algorithm Approach. Mathematics, 8(1), 84. https://doi.org/10.3390/math8010084

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