A Hybrid Approach to Solve a Model of Closed-Loop Supply Chain

This paper investigates a closed-loop supply chain network, including plants, demand centers, as well as collection centers, and disposal centers. In forward flow, the products are directly sent to demand centers, after being produced by plants, but in the reverse flow, reused products are returned to collection centers and, after investigating, are partly sent to disposal centers and the other part is resent to plants for remanufacturing.The proposed mathematical model is based on mixed-integer programming and helps minimizing the total cost. Total costs include the expenditure of establishing new centers, producing new products, cargo transport in the network, and disposal. The model aims to answer these two questions. (1) What number and in which places the plants, collection centers, and disposal centers will be constructed. (2)What amount of products will be flowing in each segment of the chain, in order to minimize the total cost. Four types of tuned metaheuristic algorithms were used, which are hybrid forms of genetic andfirefly algorithms. Finally an adequate number of instances are generated to analyse the behavior of proposed algorithms. Computational results reveal that iterative sequentialization hybrid provides better solution compared with the other approaches in large size.


Introduction
Nowadays all countries are focusing on sustainable development, which is only achieved by saving and optimal using of the limited and nonrenewable resources in that country.Governments implement a vast range of activities to achieve such purpose, including enacting rules and regulations to encourage using environment-friendly raw materials in industrial complexes, reduce the usage of fossil and oil energy, and recycle the products in both public and private sectors.Sustainable supply chain management includes superseding, reusing, recycling, and defusing materials, as well as improving public wellbeing, and reducing chain costs.
When companies use a forward-reverse system (closedloop), many environmental draw backs will be prevented, for instance, energy waste, high consumption, and transportation costs.Supply chain networks of sustainable development are considered suitable means for designing closed-loop production networks.This method fits the procedures in an economic cycle, which helps economic and environmental development significantly.
Forward activities include improving the design, engineering, production and supply, marketing, sales, and aftersale services of a new product.Reverse activities consist of reverse logistics, usage and disposal places, testing, classifying, refreshing, recycling, marketing, and reselling.Economic approaches, governmental strategies, and customer pressure are 3 main aspects of reverse logistics [1].
There is a key difference between reverse chain and waste management.Waste management means an efficient collection procedure, which does not reuse the products, but in reverse-network, the products are resurrected.
A closed loop network contrives the forward and reverse activities in a unique system.In other words, the integration of forward and reverse supply chains forms closed-loop supply chain.This integration can improve economic and environmental functions in most of the industries.
In recent decades, the pressure of world economy has forced companies to introduce new core competencies.Not only are the companies encouraged by environmental rules to recycle the used products, but also they have learned the economic benefits of this approach; hence, developing reverse logistics and closed-loop networks have become a key focal point for researchers and companies.
The idea of this research arises from a tire manufacturer that aspires to supply the demands of their clients and reduce their costs.Knowing that they felt the need to design a closedloop supply chain, which includes locating and constructing new plants for manufacturing and remanufacturing with two different technologies, collection centers, and disposal centers, therefore, this research tries to develop a mathematical model for the industry in question.Suitable solutions are sought for the same purpose.
The remainder of the paper is organized as follows.In Section 2, we offer a relevant literature on closed-loop supply chain network.Section 3 is dedicated to define and mathematical formulation of the problem.Comprehensive explanations of proposed solutions are discussed in Section 4. Section 5 is dedicated to computational results including data gathering, tuning the parameters and comparisons between solution approaches.Finally in Section 6, conclusions are provided and some future suggestions are given.

Literature Review
There have been lots of studies about supply chain network design.In order to focus more precisely, we divided our review in three parts including forward networks, reverse networks, and closed-loop networks.Additionally at the end of this section, we will present 2 tables to provide some comparison between different studies.
Most of studies in supply chain design related to forward network such as [2][3][4][5][6][7][8][9][10] Sabri and Beamon [2] developed a mixed integer linear programming (MILP) model with the aim to minimize cost and maximize service level and flexibility.The mentioned model can be used for simultaneous strategic and operational planning.Syarif et al. [3] considered a multistage logistic problem and formulated it with MILP model.In order to be solved, they utilized spanning treebased genetic algorithm.Jayaraman and Ross [4] focused on planning and execution stages of a PLOT (production, logistics, outbound, and transportation) system.They also presented an MILP model and solved it with a simulated annealing methodology.Miranda and Garrido [5] presented a simultaneous approach to incorporate inventory control decisions in a distribution network design system.In this way, they developed a mixed integer nonlinear programming (MINLP) model and solved it with the Lagrangian relaxation and the subgradient method.Altiparmak et al. [6] developed a MINLP model with some conflicting objectives consisting of cost, service level, and equity of the capacity utilization ratio.For a solution, they used a genetic algorithm approach.Amiri [7] developed a MILP model and provided a heuristic solution procedure based on Lagrangian relaxation for a distribution network.Tsiakis and Papageorgiou [8] determined the optimal configuration of a production and distribution network with respect to operational and financial limitations.In this way they proposed a MILP model with the aim of cost minimization and solved it with an exact method.Pishvaee et al. [9] considered a socially responsible supply chain network design.In addition to dealing with uncertain conditions in a network a robust possibilistic programming was developed.Pishvaee et al. [10] developed a credibility-based fuzzy mathematical programming model with the purpose to minimize total costs and the environmental impacts for a green logistics design under uncertainty.
In this part we focus on reverse logistic networks by review papers [11][12][13][14][15][16][17].Govindan et al. [11] classified 382 papers by problem classes, solution, and modeling approaches.Mentioned paper is generally divided into 4 parts including reverse logistics, closed-loop, sustainable and green but more focused on the first two parts.Jayaraman et al. [12] considered a reverse distribution network and developed an MILP model for it.Apart from that, they introduced a heuristic solution methodology.Listes ¸and Dekker [13] developed a stochastic approach for product recovery network design and applied it to a real case study on recycling sand from demolition waste in The Netherlands.Aras et al. [14] developed a mixed-integer nonlinear facility location-allocation model under pick-up strategy with capacitated vehicles.The solution approach was based on a Tabu search method.Schweiger and Sahamie [15] considered a combined continuous and discrete facility location problem containing home-recycling techniques, as well as selling and disposing technologies.In addition a hybrid Tabu search algorithm is used to solve it.Eskandarpour et al. [16] proposed a biobjective MILP model for postsales reverse logistics network operated by a 3PL.They solved their model with a novel multistart variable neighborhood search with three new encoding-decoding mechanisms.Roghanian and Pazhoheshfar [17] considered multiproduct, multistage reverse logistics network under stochastic environment and proposed a probabilistic mix integer linear model for it.They also presented a priority based genetic algorithm to solve their mathematical model.
Papers [18][19][20][21][22][23][24][25][26][27][28] belong to closed loop supply chain.Fleischmann et al. [18] presented a MILP facility location model and used it to analyze the impact of product return flows on logistics networks.Ko and Evans [19] developed a MINLP model for a dynamic integrated forward/reverse logistics network and solved it by a genetic algorithm.Pishvaee et al. [20] tried to minimize the total costs and maximize responsibility.They applied mixed-integer programming with two objectives, to find a set of nondominated solutions.Pishvaee and Torabi [21] used stochastic approach in closed-loop supply chains with uncertainty conditions.They aimed to prevent suboptimizing the proposed model.This study combined the design decisions in both forward and reverse networks and also used fuzzy approach to solve the model.Amin and Zhang [22] considered a supply chain network with the aim to maximize profit, minimize the rate of defects, and maximize importance of suppliers.Their model included 2 phases.In the first phase, the selection factor of suppliers is determined in a reverse flow.The second phase introduced a mixed-inte ger programming model.This was the first attempt to simultaneously select the suppliers and allocate the orders in a closed-loop network.Pishvaee and Razmi [23] applied a multiobjective fuzzy programming model to design an environment-friendly supply chain.Mentioned model had several parameters of uncertainty and intended to minimize the costs and environmental effects.
2013 and 2014 are remarkable years with regard to the researches on closed-loop systems.Numerous articles, such as [24][25][26][27][28], have been published in these years.Amin and Zhang [24] used a mixed-integer linear model for minimizing the total costs and at the same time considering environmental factors.They applied Weighted Sum Method and constraint Method to develop the model.They also tried to probe the effect of demand uncertainty and return uncertainty in the network.The mentioned article used stochastic programming technique for this purpose.Ramezani et al. [25] devised a multiobjective stochastic model to design a forward-reverse logistic chain in uncertainty environment.This model incorporates 3 levels in forward flow and 2 levels in reverse flow and attempted to maximize the interest, responsiveness to clients, and quality in the logistic network.Soleimani et al. [26] developed a mix integer linear model for a multiechelon, multiproduct, and multiperiod closed loop network.In their model plants, warehouses and distributors can stratify customer demand.For a solution they utilized CPLEX and genetic algorithm.Faccio et al. [27] introduced a linear programming model for a closed-loop problem, with the aim to minimize the total costs of this chain.Parametric study and robust analysis of this model have been contrasted with the classic forward-flow approach and social responsibility.Devika et al. [28] developed a mixed integer programming model in closed-loop supply chain based on triple bottom line approach.They also presented three novel hybrid metaheuristic methods based on adapted imperialist competitive algorithms and variable neighborhood search to solve it.
These mentioned researches indicate that metaheuristic algorithms are becoming the most efficient approaches to solve complex problems in supply chain network design.
Table 1 is dedicated to coding system and Table 2 shows the mentioned models in different studies.

Problem Definition
The problem in this paper has been taken from a closed-loop supply chain, including plants, collection centers, demand centers, and disposal centers.Plants are able to manufacture new products and remanufacture the returned products.In forward flow, the products are dispatched to demand centers, after being produced by factories.In reverse flow, the returned products are returned to collection centers.Collection centers are in charge of the following:  (iii) sending renewable returned products to plants (for remanufacturing), (iv) sending nonrenewable returned products to disposal centers.
The model tries to answer this question.In which places and what number of plants, collection centers, and disposal centers will be established, and what amount of products must go there to minimize the total cost.Figure 1 provides a pictorial representation of the closed-loop supply chain.In addition to defining the problem precisely we consider the following assumption and limitations.
(i) All the returned products from demand centers are kept in collection centers.The demand centers have fixed locations.
(ii) The potential locations, capacity of plants, collection centers, and disposal centers are already determined.
(iii) The needs of every demand center can be met with the help of several plants.When founding a new plant, diverse manufacturing and remanufacturing technologies are available.The needs of demand centers will be met completely.
(iv) There is no cargo-flow between the same level facilities.

Proposed Mathematical Formulation. The mathematical model of this problem becomes
Min: TC The objective is minimizing the total costs.There are 7 types of cost.The first 3 are the constant expenses of establishing new plants, new collection centers, and disposal centers.The 4th part includes the transportation and manufacturing costs of producing new products.The 5th expenditure will be recycling and transporting the returned products.The 6th segment is associated with the total costs of recycling and transporting returned products from collection centers to factories.Finally, the 7th part of target function represents the transportation and disposal costs.Constraint (2) guarantees that the total number of manufactured products for each demand center is equal to or greater than their demand rate.Constraint (3) shows the capacity of plants.Constraint (4) shows that the flow rate in forward method is higher than reverse flow.Constraint (5) represents the disposal amount.Constraint (6) presents the limited capacity of collection centers.Equation (7) indicates that the amount of returned products from customer centers is equal to the amount of returned products to plants and the amount of products sent to disposal centers from every collection center for every product.Constraint (8) illustrates the limited capacity of disposal centers.Equation ( 9) refers to the returned products.Constraint (10) demonstrates that every plant, if established, enjoys only one type of technology.In other words no combination of technologies will be accepted.The two final constraints in the model focus on the decision variables.Three variables are associated with establishing and not establishing the plants, collection centers, and disposal centers, respectively, in their potential locations.These variables are binary.Other decision variables in the last constraint are defined as nonnegative.

Solutions
For the optimization problems with large scale and high complexity, classic solutions will not be reliable.The solutions must search the feasible space in a more sophisticated way.As mentioned in literature review nowadays many researchers devise metaheuristic algorithms for complicated problems of supply chain.
Before applying metaheuristic we must show that the closed-loop problem is NP hard.In this way we mentioned the fact that the closed-loop problem is a capacitated location-allocation problem and also can be considered as a multiple-choice Knapsack problem, for it is known to be NPhard [3,29,30].So metaheuristic algorithm can be one of the best candidate solvers for such problems.The problem in question becomes more complicated by increasing the number of facilities and product technologies.This might make the accurate and heuristic solutions inapplicable; hence, metaheuristic algorithms become the focal point.Among the metaheuristic solving techniques, hybridization methods play an important role in optimizing the case.

Genetic Algorithms.
The main idea of genetic algorithms arises from nature.These are effective random searching techniques in mega-scale target areas.They search in different directions to find the optimum answer.
Genetic algorithms usually generate several random solutions to the problem initially.This set of solutions is called "The initial population, " and each member of the set is titled "Chromosome." In order to find better chromosomes, the operators of genetic algorithm cross the chromosomes over each other.This cross-over leads to mutation, so that new chromosomes are generated.The best members are separated.This is how the second generation is formed.In other words, the operators of selecting, crossing over, and mutating every generation produce a new generation of chromosomes.The new generation is different from the previous one.The whole process is repeated over and over for the next generations, until the ultimate solution is obtained.The steps of this algorithm are explained below in detail, respectively: (1) creating and apprising the initial population that is generated randomly, (2) selecting the parents with one of the selection methods (roulette wheel, tournament, and random), (3) applying crossover operator on the selected chromosomes as parents, (4) selecting the chromosomes for mutating, (5) applying the mutation operator, (6) creating a new population from the best members of all the 3 main populations: crossed-over, mutated, and initial population, (7)   the simplest method for presenting the solution structure is using a string in problems.Matrix view is sometimes applied as well.
With modern calculation tools, such as Matlab Software, stronger and more efficient structures have become available for presenting the solutions.Data structure is one of them.In data structure, every answer is associated with a structure.The structure will also have substructures based on the number of variables in the problem.For the current problem, every answer is initially presented with one structure.There are 7 key variables, so every answer has 7 substructures, which are also divided into substructures according to the problem scale (number of indexes for manufacturers, collection centers, demand centers, disposal centers, etc.).The scheme of one answer is shown in Figure 2.

Selection.
For every member, a pair of parents is selected.The aim is to choose the best elements.Even the weakest element has the chance to be selected.This prevents stopping on local solutions.Numerous selection operators have been introduced so far, for example, operators roulette wheel, tournament, and random.This paper has coded all the 3 methods, but the vast range of tests proved that method roulette wheel is the best selection technique for our model.

Crossover.
Genetic algorithms show the probability for crossover operator, which is between 0 and 1.This presents the probability of creating offspring.This is the likelihood of chromosomes being crossed over again.The crossover of every pair of chromosomes will create an offspring, which is a new member of the next generation.This operator tries to find suitable potential answers in the next generation.
There are various cross-over methods, such as methods partial-mapped crossover (PMX), order crossover (OX), multipoint crossover, and cycle crossover (CX) [29].The current study uses multipoint crossover technique, which is a developed version of single-point crossover.Two structures are primarily considered.Then several locations are selected in the structure randomly.The elements of these locations are substituted with each other.A scheme of this crossover technique comes in Figure 3.This paper creates the initial solutions with the technique that prevents generating any type of unfeasible answer.During cross-over and mutation steps, if the changes in elements lead to such answers, again a correction method heals the answers towards the feasible space; hence, penalty function will not be necessary; therefore, the fitness function in this study will be exactly the same as target function in (1).The suggested correction method is divided into 3 techniques according to the constraints of this model.The first one is "less than or equal to constraints." Assume constraint ∑    ≤   is applied and, for certain reasons, variable   that is generated in chromosome structure has violated feasibility conditions (∑    ≻   ).In such case, we find the maximum number in the structure of   and suppose it is equal to the minimum number in the structure.This is repeated over and over until we get back to feasibility.
The second case is "greater than or equal to constraints." Assume constraint ∑    ≥   exists and, for certain reasons, the generated   created by the constraint is not feasible (∑    ≺   ).In such case, we find the minimum number in the structure of   and make it equal to the maximum number in the structure.This is repeated until we get back to feasibility.
The 3rd case is "equality constraint." Assume the existence of constraint ∑    =   .The generated variable   does not fit equality conditions (∑    ̸ =   ).In such case again, we create a matrix of random numbers between 0 and 1 and called it    .The new   is obtained from the following formula: It is noteworthy that all cases 1, 2, and 3 can be generalized to several summations.
4.1.5.Mutation.Genetic algorithms have small and constant change probability, which is between 0 and 1.According to this probability, child chromosomes change randomly and come across mutation.There are various methods of mutating, such as methods inversion, insertion, and uniform.
In the suggested chromosomes' structure, every variable sits in a fixed location; therefore, using mutation operators, that substitute the location of one variable with another, leads to losing the answer structure; hence, methods such as uniform can be proper choices for mutating.
The proposed mutation operator in the current study selects one of the 7 variables randomly and causes some change in it.Figure 4 shows a scheme of mutating method.

Finishing Conditions.
It is obvious that the designed algorithm, which searches through the answer space, has to stop somewhere and somehow; otherwise, the search procedure will grow too long and sometimes does not come to any result at all.
Several methods have been introduced and devised so far about the finishing criteria of metaheuristic algorithms, for instance, gaining an acceptable level in the answer and passing through a certain number of generations in the algorithm and a particular number of loops without improvement in answer.Some researchers use a time limit for ending the algorithm.They restrict the run time or the maximum time without gaining an improved answer.Among all introduced criteria, time-based techniques seem to be unreliable because various computers have different processing speeds; besides, every year new models with higher speed and vaster features come to market.This paper regards the number of generations as finishing criterion for the algorithm.Every algorithm stops after cycling through a certain number of generations.

Firefly Algorithm.
Firefly algorithm was introduced firstly in 2007 by someone named Yang [31].This is one of the newest optimization algorithms, whose idea arises from nature, and is based upon swarm intelligence.In FF algorithm, a few artificial fireflies are initially generated randomly, with a structure similar to the chromosomes in genetic algorithm.
After generating and distributing the fireflies, every firefly radiates some light, whose intensity represents the optimality level of its location.In other words, the firefly with higher fitness is more attractive and grabs more attention.The attractiveness of fireflies is based on their brightness and light strength; hence, the closer they are to the light source, the more attractive they are.Additionally, if a firefly is weak but close to the light source, again it can enjoy suitable attractiveness.The light intensity of every firefly is continuously compared with the light strength of other fireflies, and the difference of light level stimulates the firefly to move.This means if the firefly has higher fitness than the contrasted firefly, it moves to the new location, which is an average between the location of higherfitness firefly and the contrasted one.
If the problem objective is minimization, lighter firefly is attracted towards the dimmer firefly.In this algorithm, even the weakest firefly moves randomly inside the problem domain with the aim of increasing the chance of finding global optimum.
The decrease in light level is caused and determined by light absorption factor and controlled within a certain scale.The information exchange among fireflies is done via light radiation towards each other.Finally, this whole group procedure causes the fireflies to tend towards the more optimal point.

4.2.1.
The Procedure of FF.The steps of this algorithm are shown below.
(1) Generate the initial population of fireflies.
(2) Calculate light intensity  for every firefly (we can always consider objective function of a firefly as light intensity of it [31]): (3) For every firefly (as "") and for every other firefly (as ""), if the received light intensity of firefly "" is less than the received light of firefly "" (  ≺   ), we produce new solution according to following formula: is the location of firefly ,   is the location of firefly ,  is the mutation coefficient (it is between 0 and 1 [31]). 0 is attraction coefficient base value (usually considered 1 [31]),  is a random vector (usually considered sign [rand−1/2]), and  is a light absorption coefficient (it typically varies from 0.01 to 100 [31]).
In addition  is the Cartesian distance to light source and calculated by the following equation [31]: (4) Evaluate the new fireflies.
(5) Determine the best detected answer.
(6) Loop to step (3), if the ending conditions are not true.
Again, the structure of each firefly is similar to the chromosome in Figure 2. Cost function is similarly calculated as well.

Hybridization of Algorithms.
With the help of hybridizing process, two or more algorithms can be combined to create new methods, which inherit part of the characteristics of each creator algorithm and form a new algorithm with new different characteristics [32].The goal of algorithms' hybridization is improving their characteristics or creating new characters that are not available in each algorithm alone.That is why in the recent years more and more hybrid algorithms are being born and applied.Among different hybrids, combining a search method within a genetic algorithm can upgrade the search performance and also it is a chance to hybrid optimization to attain their target [33].
The range of hybridization methods is enormous, for example, sequentialization and parallel methods.In current study, four hybridization methods have been devised, and with the help of Matlab Software, all of them were coded; besides, a summary about the functionality of each of the 4 hybridization methods is explained with the related chart.

Algorithm PHGF.
The first suggested algorithm is a simple parallelization technique.In this technique, the calculations of genetic and firefly algorithms are performed in a parallel way with similar cycles, and the final output is chosen from the best answers of both algorithms.Figure 5 shows an illustration of simple parallelization hybridization.
The steps of PHGF are shown below.

Algorithm ISHGF.
In this technique, the output of genetic algorithm is used as the input of firefly algorithm, and then the output of firefly algorithm will be used as the input of genetic algorithm.This will be done over and over until finishing conditions come true.Such new algorithm is obtained from hybridization with iterative sequentialization method and is called ISHGF.
In ISHGF, according to the sizes of GA and FF population members, 3 scenarios might occur.The first scenario is   =   (  is the number of individuals in the population and   is the total number of fireflies).In such case, the answers are sent continually from GA to FF algorithm and from FF to GA, until gaining the desired answer.The second scenario is   >   .In such case, the best members of GA population sit in the place of FF members, and the remaining group of members whose number is calculated by subtracting   from   are thrown away.Of course when the algorithm has not reached the desired answer, GA receives all the members of FF and compensates its lacking members by generating random answers.The 3rd scenario is   <   .Hereon, all GA members go into FF algorithm, and new random answers are generated as many as   −   .Of course in this scenario again, if the ending conditions do not arise, the algorithm returns to its beginning, so GA selects only as many as   members among the best answers of   and throws the rest away.Figure 7 shows an illustration of iterative sequentialization hybridization.
The steps of ISHGF are shown below (  >   ).
(  (d) Rank the solution, update the best and also generate (  −   ) members randomly and add them to   (Added members +   →   ).
End while (3) Select best of final generation as ISHGF solution.

Algorithm WHGF.
The fourth suggested algorithm is simple workload division.This technique is very similar to simple parallelization, meaning that every algorithm works in parallel with and independent of other algorithms, but in the end, the collected answers appear in the final answer according to the determined share for every algorithm.For example, firefly algorithm has 60% and genetic algorithm has 40% share in the ultimate answer.Obviously, the final solution is the best answer among all the answers collected from both participant algorithms in the hybrid process.Figure 8 shows an illustration of simple workload division hybridization.The steps of WHGF are shown below.
(1) Randomly generate an initial population (  ) (2) While (not stop condition) (a) For ℎ = 1:   (  is the number of individuals in the population) (a.

Data Gathering.
In order to investigate the performance of suggested algorithms, numerous tests are required.First of all, certain numbers of the parameters required in the problem must be determined.The performance of suggested solutions is to be appraised in various scales: small, medium, and large.Knowing that, direct input of numbers into the parameters, and changing the numbers in every test will make the solution too complicated; hence, uniform statistical distribution was applied here.Table 3 shows the values of required parameters for solving this problem.

Tuning the Parameters.
For improving the performance of suggested hybrid algorithms, a particular process called parameter tuning is required.Parameter tuning keeps the required parameters of the problem in the best possible range; so that every parameter performs in the best possible way.Numerous techniques have been used so far for setting the parameters of metaheuristic algorithms, including Trial And Error, Full-Factorial, Taguchi, Response Surface, and Neural Networks.This paper aims to set the parameters of all 4 suggested algorithms with Taguchi Method, which was innovated by Taguchi in 1986 [34].The first step is finding and introducing the influential parameters in algorithm.Since the suggested solutions are all based on genetic and firefly algorithms, the influential factors in performance of their hybrid also include the parameters of those two.
The parameters of genetic algorithm include mutation rate, crossover rate, number of iterations, and number of chromosomes in the population.The selected parameters in firefly algorithm also include the number of artificial fireflies, light absorption coefficient, mutation coefficient, and the maximum number of iterations.The parameter values are defined in 3 levels: small, medium and large.These values are shown with indexes 1, 2, and 3, respectively.Since the mathematical model of this problem is not similar to other studies, the value of parameters cannot be compared with other researches; therefore, these levels have been determined with trial-and-error, within their permitted range.Taguchi solution approach looks at the number of factors and the levels and determines the required number of tests accordingly.Since this paper has selected 8 factors and 3 levels, 27 problems are to be solved.Taguchi technique in this problem has been performed by Minitab Software, on a PC with 3 GB RAM.Table 4 shows the considered parameters and their levels.
For tuning the parameters, the obtained values should become normalized, so as to make them independent of range.Numerous methods are available for making the parameters dimensionless, including methods fuzzy normalization, and linear normalization.This paper uses related percentage deviation (RPD) technique, with the following formula: Table 5 shows the results of 27 times problem solving, with different levels of factors, and in RPD form.The ultimate aim of parameter tuning is adjusting each of the algorithms in their best possible position.This is why a sample problem is initially selected and solved with various levels of parameters, and then the results are entered in a column.The best answer is the minimum value (because this is a minimizing problem).According to RPD formula, the best answer turns into zero, and the other answers are evaluated according to their deviation from the best; therefore, since every algorithm is evaluated separately, every column has 1 zero.The obtained numbers in every algorithm were entered into Minitab separately.The following graphs show the results of suggested levels for parameters, in Taguchi Method.Minitab Software introduces two graphs for every algorithm: one associated with mean and the other one with SNR.Every Mean diagram consists of 8 graphs, indicating The other diagram is SNR, which illustrates the optimum level of every parameter in another way.The higher SNR level is better.Here is a caveat to keep in mind.If the parameters are tuned properly, the results of both Means and SNR diagrams will be equal; otherwise, if the interpretation of every parameter in average diagram is different from its optimum level in SNR diagram, the tests must be repeated.
Note that we have one type of population in the whole EHGF algorithm; therefore, the number of fireflies does not define.
Table 6 shows the tuned levels of parameters in every algorithm, according to Figures 9,10,11,and 12.After entering the tuned value of every parameter, the algorithm will be used for solving 40 problems in small,     Interval plot 95% CI for the Mean Small size medium, and large scale.Table 7 demonstrates the resulted values.In addition we consider our base algorithms GA and FF to investigate effect of hybridization on the test problems.

Comparisons.
The values obtained from all 6 solutions are contrasted according to statistical comparison techniques, with the aim to find the best solution.For each problem, the proposed algorithms are run 3 times, and then the average value is considered.In order to determine the significant difference between different values obtained by algorithm, ANOVA Statistical Test was devised with SPSS Software.Note that in order to compare more precisely each size considered separately.Tables 8, 9, and 10 and Figures 13, 14, and 15 show the comparisons between algorithms.As shown in Figure 13 and also from ANOVA Table 8 (sig > 0.05) there is no significant difference between algorithms in solving small size problems.
Although it seems that there is a difference between algorithms in Table 9 (Sig < 0.05) because of overlapping between plots in Figure 14, we cannot certainly report which one is better in solving medium size problems.
Value obtained for sig lower than 0.05 is indicating that there is a significant difference between algorithms.In addition as shown in Figure 15, it is obvious that ISHGF achieves lower cost in large size than the other algorithms.

Conclusion
This paper introduces a new model for reducing the costs of closed-loop supply chain.The model includes several plants that are able to manufacture and remanufacture various types of products, using diverse technologies.In forward flow, the demanded products of customer centers were sent to them directly, without any go-between.In reverse flow, the collection centers are in charge of collecting and evaluating the products, so as to decide about sending the products to plants or disposal centers.The proposed mathematical model in this problem grew more and more complicated by increasing the number of facilities, products, and technologies, so that classic methods became unreliable for solving the problem; hence, this paper After tuning the optimum parameters of those two algorithms, using Taguchi Approach in Minitab Software, the parameters were used for solving 40 problems with small, medium, and large scale.Finally, ANOVA Statistical Test was devised with SPSS Software, and this shows that iterative sequentialization technique gave the best answers in large size.
Closed-loop supply chain has considerable economic, environmental, and social effects; therefore, it has become a focal point in recent years.Such model can consider social and environmental factors, which are influential in the network.It can be formulated in multiobjective form.
Increasing the number of levels and participants of the supply chain design can help improve its compatibility with the needs of that industry and society.
Aside from that, considering uncertainty in parameters, and investigating the flexibility of model with regard to gradual changes can also help make the model significantly more robust.
Using other metaheuristic methods and artificial neural networks can also be considered among the future solution ideas.

Figure 2 :
Figure 2: An illustration of chromosomal structure.

Figure 4 :
Figure 4: An illustration of mutation operations.

Figure 8 :
Figure 8: An illustration of simple workload division hybridization.

Figure 9 :
Figure 9: Main effects plots of Means (a) and SNR (b) for PHGF.

Figure 10 :
Figure 10: Main effects plots of Means (a) and SNR (b) for EHGF.

Figure 12 :
Figure 12: Main effect plots for Mean (a) and SNR (b) of WHGF.

Figure 13 :
Figure 13: Comparison between algorithms in small size.

Figure 14 :
Figure 14: Comparison between algorithms in medium size.

Figure 15 :
Figure 15: Comparison between algorithms in large size.

Table 1 :
Coding of literature review.
: quantity of returned product  shipped from collection-inspection center  to plant ,

Table 2 :
Review of previously published literature.  : distance between plant  and customer center , ℎ  : distance between collection-inspection center  to plant ,   : distance between collection-inspection center  to disposal center ,   : distance between collection-inspection center  to customer center ,   : demand of customer center  for product ,   : amount of returned product  from customer center ,   : minimum disposal rate for product ,   : maximum manufacturing capacity of plant  for product  with technology , d Facilities to be located in the network.: quantity of returned product  shipped from collection-inspection center  to disposal center ,   : 1 if the plant center  is to be established with technology , 0 otherwise,   : 1 if the collection-inspection center  is to be established, 0 otherwise,   : 1 if the disposal center  is to be established, 0 otherwise.3.1.3.Parameters  : unit cost of manufacturing product  in plant  with technology ,   : unit cost of transportation product  from plants to customer centers (Per km),   : unit cost of transportation product  from customer centers to collection-inspection centers (Per km),   : unit cost of transportation product  from collection-inspection centers to plants (Per km),   : unit cost of transportation product  from collection-inspection centers to disposal centers (Per km),   : fixed cost of establishing plant  with technology ,    : fixed cost of establishing collection-inspection center ,    : fixed cost of establishing disposal center ,   : unit saving cost of product ,   : unit disposal cost of product  at center ,  : maximum capacity of collection-inspection center  for product ,    : maximum capacity of disposal center  for product .
Mutate   percent of population (  ) to produce   offspring and then evaluate them.
(b) For  = 1:   (  is the total number of fireflies) (b.1)For  = 1:   If   <   Generate new solution and update it (  ) according to formulation in (15).(1) Randomly generate an initial population () (2) While (not stop condition) (d) Merge population generated via mutation and population from firefly algorithm in a unique set (  ∪   ) and then sort them.(e) Truncate best members of merged set and generate new population of size .(  ∪   ) →  End while (3) Select best of final generation () as EHGF solution.
1) Recombine   percent of population to produce    offspring and then evaluate them.(a.2) Mutate   percent of population to produce Save   members of best solutions obtain by GA in final For  = 1:   (  is the total number of firefly) (b.1)For  = 1:   If   <   Generate new solution and update it (  ) according to formulation in (15) Save   members of best solutions obtain by FF in final ff (7) Merge % best of GA and 1 − % best of FF (% (final ga ) ∪ 1 − % (final ff )) → final ℎ (8) Rank the member of final ℎ and select best of them as WHGF solution.

Table 3 :
Values of model parameters.

Table 4 :
Levels of algorithm parameters.

Table 5 :
Values of proposed algorithms in parameter tuning.

Table 7 :
Result of objective function in different sizes.

Table 8 :
ANOVA of small size.

Table 9 :
ANOVA for medium size.

Table 10 :
ANOVA for large size.