Next Article in Journal
Developing a Sustainable Logistics Service Quality Scale for Logistics Service Providers in Egypt
Previous Article in Journal
The Effects of Different Supply Chain Integration Strategies on Disruption Recovery: A System Dynamics Study on the Cheese Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Designing a Sustainable Green Closed-Loop Supply Chain under Uncertainty and Various Capacity Levels

Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA 02115, USA
*
Author to whom correspondence should be addressed.
Submission received: 17 February 2021 / Revised: 12 March 2021 / Accepted: 16 March 2021 / Published: 7 April 2021

Abstract

:
The ever-increasing concerns of the growth in the volume of waste tires and new strict government legislations to reduce the environmental impact of the end-of-life (EOL) tires have increased interest among companies to design a sustainable and efficient closed-loop supply-chain (CLSC) network. In the real world, the CLSC network design is subject to a variety of uncertainties, such as random and fuzzy (epistemic) uncertainties. Designing a reliable and environmentally cautious CLSC with consideration of risks and the uncertainty of the parameters in the network is necessary for a successful supply-chain network. This study proposes a sustainable and environmentally cautious closed-loop supply-chain network for the tire industry, by considering several recovery options, including retreading, recycling, and energy recovery. This study aims to design and develop a robust multi-objective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network under hybrid uncertainty. There are two types of uncertainties associated with the parameters in the network. There is an uncertainty associated with the demand, which is expressed in some future scenarios according to the probability of their occurrences, as well as fuzzy-based uncertainty associated with return rates, retreading rates, recycling rates, procurement, and production costs, which are expressed with possibilistic distributions. In order to deal with this hybrid uncertainty, a robust fuzzy stochastic programming approach has been proposed, and the proposed mixed integer programming model is applied to a case study in the tire industry to validate the model. The result indicates the applicability of the proposed model and its efficiency to control the hybrid uncertainties and the risk level in the network.

1. Introduction

In today’s business world, companies need to pay close attention to their supply-chain management systems and logistics network, in order to effectively compete in a highly competitive market. A great supply-chain system is a necessity for any manufacturing organization. Most manufacturers put their focus primarily on the forward supply chain; however, the key to success is to have a holistic view of the entire forward and reverse logistic network. Reverse logistics is a network where selected products transfer up the stream from the end users for repair, remanufacturing, or recycling [1]. Environmental concerns, government regulations, and the potential economic benefits are some of the motives for manufacturers to move toward this direction [2,3]. Besides environmental regulations, consumer’s awareness and their responsibility against the environment are some of the motives. The tire industry is one of the industry sectors moving toward this direction. Worldwide, more than a billion tires are manufactured every year, and a significant number of them are thrown away after being used. According to the US Environmental Protection Agency (EPA), in the United States alone, 280 million tires are being scraped every year. The average passenger tire weighs about 22 pounds, and the average commercial semi-truck tire weighs about 110 pounds, and these tires contain oil and other combustible carbon compounds, which create a fire hazard. Figure 1. Tire fires are difficult to extinguish and can break the tire down into hazardous compounds, including gases, heavy metals, and oil. Tire piles can also create a breeding ground for pests such as rodents and mosquitos that contribute to the spread of diseases [4,5]. Although the forward tire supply chain is well developed, there is no systemic and comprehensive reverse logistics system in place. Most collection, retreading, and recycling centers are small individual shops with a limited capacity. Government regulations and requirements vary from state to state. This lack of a comprehensive plan can be seen by a large volume of tires still ending up in the landfills.
There are many factors that are considered by the decision maker (DM) when designing a closed-loop supply-chain (CLSC) network, such as costs, environmental impacts, social dimensions, uncertainty in the parameters, governmental regulations, disruption and operational risks, reliability of the network, etc. Uncertainty with the parameters is probably the biggest challenge facing the supply-chain network (SCN) design. Random and fuzzy uncertainties are generally the two types of uncertainty in the design of a supply-chain network [6]. Possibilistic, stochastic, and robust programming are typically the methods utilized to control these uncertainties. Stochastic programming can be used when there are sufficient data available to estimate the probability distribution of the uncertain parameters that is used for modeling the uncertain parameters under different scenarios in the network [7]. For instance, several scenarios for demand may be related to economic conditions, occurrence of natural disasters, changes in customer preferences, development of new technologies, etc. However, the value of specific parameters in each scenario could be imprecise and can be determined by possibilistic distribution.
Possibilistic programming is used when there are no sufficient historical data available, but the possibility distribution of uncertain parameters can be estimated by using the expert’s knowledge in the field. Robust programming is a method that allows the DMs to adjust the risk level of output results with respect to the uncertainty of parameters [8]. Combining the stochastic, possibilistic, and robust programming can supplement the weaknesses of each individual method and provide more reliable results for the DMs in the supply-chain network.
In this study, robust fuzzy stochastic programming was used to cope with hybrid uncertainties in the model. Moreover, in this paper, we provide a holistic and realistic modeling approach for the coordination of a forward and reverse tire management system that extends the life span of re-manufacturable tires, recovers the value of scrap tires, and eliminates the disposal in the landfills. In that respect, we examine different recovery options simultaneously, such as retreading, recycling, and energy recovery. Our model aims to maximize the total CLSC profit, minimize the total greenhouse gas (GHG) generated along the CLSC and maximize the social responsibilities which is presented in the form of the total number of jobs created by expanding the network in the forward direction and establishing a new reverse logistics network.
As the body of literature about CLSC network design shows, most models ignore the fact that products can be remanufactured multiple times and remain in the supply-chain network to satisfy the customer demand, rather than being recycled prematurely. We refer to this as a multi-cycle supply-chain network. This can impact model assumptions and outputs with respect to overall market demand, the production rate of manufacturing centers, the return rate of end-of-life (EOL) products, and the capacity of the facilities in the network. Also, a genuinely green supply-chain network must consider and plan for the return of EOL remanufactured products for final processing. This has also been neglected in most previous studies.
This study aims to design a closed-loop supply-chain network considering the economic, environmental, and social factors with the following contributions:
  • Proposing a new sustainable, robust, multi-objective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network by considering the real-world situation where used tires are retreaded multiple times, and the retread tires are being returned to the collection center for final processing.
  • Presenting a novel Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) model that controls scenario-based uncertainty (disruption risks) and fuzzy-based uncertainty (operational risks) in the network simultaneously and provide a more reliable output results for the managers and the DMs.
  • Evaluating the possibility and feasibility of expanding the existing forward logistic network at different capacity levels to cope with market demand under various scenarios, and also the feasibility of establishing a new reverse logistic network at different capacity levels to cope with the return rate of EOL tires.
  • Validation of the proposed model by sensitivity analysis through a case study.
The paper is structured as follows. In Section 2, the relevant literature on closed-loop supply-chains network design and tire industry is reviewed; some works concerning uncertainty in the supply chain are also analyzed. In Section 3, the main problem characteristics, assumptions, and mathematical formulation are presented. In Section 4, a proposed Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) approach and augmented epsilon constraint method are described in detail. In Section 5, the proposed model is applied to an example problem inspired by a case in the U.S. The computational results and the sensitivity analysis are also presented in this section. Finally, in Section 6, some conclusions and future work are discussed.

2. Literature Review

Designing a robust and flexible supply-chain network (SCN) that includes optimizing the number, capacity, and location of network facilities has received a lot of attention these days [9]. A closed-loop supply-chain network (CLSC) considers both forward and reverse supply-chain networks in an integrated manner. There have been many papers published in this subject area and especially green closed-loop supply-chain network. Govindan et al. (2014) and Devika et al. (2015) provided a comprehensive survey of the CLSC network [10,11]. Moreover, a comprehensive review on environmentally conscious manufacturing and product recovery was provided by Ilgin and Gupta (2010) [9]. Jindal and Sangwan (2016) presented the design and optimization of a multi-objective closed-loop supply chain, considering the economic and environmental factors under uncertainty using the interactive ε-constraint method [12]. Sarkar et al. (2017) studied the effect of transportation cost and carbon emission cost on remanufacturing. They proposed a multi-echelon closed-loop supply-chain model with third-party logistics (3PL) [13]. In another study, they developed a multi-attribute closed-loop supply-chain model for self-healing polymers-based returnable transport packaging. Their proposed multi-objective supply-chain model aimed to maximize profit while minimizing carbon emission of the system [14].
Dulman and Gupta (2018) studied the financial impact of sensors on closed-loop supply-chain systems [15]. Zhou and Gupta (2019) provided a pricing strategy for new and remanufactured high-tech products. Their model aimed to maximize the profit associated with customer demand based on the manufacturer, retailer, and joint supply chain [16]. In another study, they proposed partial least square method to explore the factors that affect value depreciation rate and price differentiation between new and remanufactured iPhone and iPad [17]. Aldoukhi and Gupta (2019) proposed a new model for designing a closed-loop supply-chain network by considering downward product substitution policy under four carbon-emission regulation policies. They used robust optimization to deal with uncertainty related to product demand and number of returned products [18]. In another study, they proposed a multi-objective model to design a closed-loop supply-chain (CLSC) network, with the aim of minimizing the total cost, minimizing the carbon emission, and maximizing the service level of the retailers. They used physical programing approach to model the problem [19]. Fadhel and Gupta (2019) evaluated the food waste valorization alternatives from a sustainability point of view. They estimated energy utilization and GHG emission reduction for each potential food waste processing technique [20]. Gupta et al. (2019) studied the optimal control of production and maintenance operations in smart custom manufacturing systems, using a heuristic algorithm [21]. De and Giri (2020) studied a closed-loop supply chain (CLSC), focusing on managing, scheduling, and routing problems, to achieve economic and environmental sustainability [22]. Mogale et al. (2018) designed a dual market facility network by considering manufacturing the common components centrally and locally [23]. Oliveira and Machado (2021) provided a literature review on the application of optimization methods in a closed-loop supply chain [24]. Zahedi et al. (2021) proposed a closed-loop supply-chain network by considering multi-task sales agencies and multi-mode transportation, and they solved their model by using metaheuristic and hybrid algorithm [25].
Guarnieri et al. (2020) proposed a theoretical model to assess reverse logistics using multi-criteria decision-making for e-waste [26]. Vieira et al. (2020) studied the barriers in implementing reverse logistics for e-waste and investigated the widely applied multi-criteria decision aid (MCDA) methods used to prioritize these barriers [27].
There have been many researches in the area of sustainability [28]. Mogale et al. (2019) designed a sustainable food-grain supply-chain distribution system with the aim of minimizing the cost and carbon dioxide emission simultaneously [29]. Ullah et al. (2020) investigated the optimal remanufacturing strategy and reusable packaging capacity under stochastic demand and return rate [30]. Ghadge et al. (2020) studied obstacles in sustainability implementation in cold-food supply-chain networks by considering a case of UK artisan cheese producers [31]. Gupta et al. (2020), presented a compilation of six recent papers on a variety of topics to demonstrate the pioneer research activity within responsible and sustainable manufacturing [32]. Sazvar et al. (2021) designed a sustainable CLSC for pharmaceutical by considering manufacturer’s brand and waste management [33].
Inconsistent, unreliable, and uncertain data are one of the biggest challenges facing the SCN [34,35,36]. In the real world, there are many parameters in SCN that are uncertain, such as customer demands, costs of opening new facilities at a specific location, manufacturing cost, cost of raw material, transportation cost, etc. [37]. According to Mula et al. (2006/2007), the two main types of uncertainty are randomness and epistemic [38,39]. Stochastic programming and possibilistic programming are the most common methods to deal with randomness and epistemic (fuzzy) uncertainty, respectively [40,41]. A flexible robust optimization approach for scenario-based stochastic programming models was proposed by Mulvey et al. (1995) [42]. Their approach was later developed by Yu and Li, 2000 [7]. Pishvaee et al. (2018) also developed a new approach, called robust possibilistic programming [43]. Mousazadeh et al. used the robust possibilistic approach for a health-service network design [44]. Tozanli and Gupta et al. (2017) provided a literature review paper for environmentally concerned logistics operation (ECLO) in fuzzy environment. They utilized over 800 papers between 1994 and 2017 and performed a detailed analysis of ECLO with an emphasis on fuzzy application [45]. Vahdat et al. (2017) provided a two-stage stochastic programming modeling to design a multi-period, multistage, and single-commodity CLSC network under uncertainty [46]. Parakash et al. (2018) proposed a CLSC network under risks and demand uncertainty [47]. Yolmeh and Saif (2020) provided a supply-chain network design with assembly and disassembly line balancing under uncertainty [48]. Wang et al. (2021) studied various incentive mechanism in green supply chain under demand uncertainty [49]. Tordecilla et al. (2021) provided a review of simulation optimization methods for designing a resilient supply chain under uncertainty [50].
There has been some research on the tire CLSC network. In 1997, Ferrer published a paper and estimated the optimal number of times each tire can be retreaded [51]. Sasikumar et al. (2010) provided an optimization model for a tire remanufacturing case [52]. According to Chopra and Meindl (2015), between 2009 and 2011, only 3 percent of total sales by US firms within the tire industry were retread tires [53]. One of the challenges facing the tire industry is the coordination of the entire forward and reverse logistics network and managing the flow of the new and used tires [54]. Amin et al. (2017) examined a case study for CLSC network for tire industry in Toronto, Canada. They developed a mixed integer linear programming model, under uncertainty, with a goal of maximizing the profit [55]. Subulan et al. (2015) examined tire remanufacturing case study in Turkey using a fuzzy mixed integer programming approach [56]. Derakhshan et al. (2017) presented a technique for recycling the tire waste, as well as reuse of the waste [57]. O’Brien and North (2017) investigated the emission of pollutants from tire waste. They studied the types of gases that are produced from car tire waste [58]. Fathollahi Fard et al. (2018) proposed a tri-level programming model for tire CLSC to identify the decision variable, using real datasets [59]. Pedram et al. (2018) used stochastic programing to design a single objective tire CLSC network [60]. Lokesh et al. (2018a) studied a tire-retreading supply-chain network under carbon tax policy. In another study, Lokesh et al. (2018b) used a fuzzy goal programming approach for Brownfield tire retreading case study under carbon tax policy [61,62]. SahebJamnia et al. (2018) also modeled a robust closed-loop supply chain of the tire industry and used hybrid metaheuristic algorithms to solve the model [63]. Abdolazimi et al. (2019) proposed a multi-objective closed-loop supply chain by integrating on time delivery and the cost [64]. Mehrjerdi at al. (2020) studied sustainability and resiliency of a closed-loop supply chain, using multiple sourcing strategy [65]. Lokesh et al. (2020) also designed a sustainable tire supply chain with a fuzzy goal planning approach [66].

Tire Industry and Recovery Options

According to Modern Tire Dealer (MTD), in 2019, 267.8 million new passenger tires, 37.9 new light truck tires, and 25.4 million new medium and heavy truck tires were shipped within the United States. EOL tires due to the large quantities being manufactured and the durability of the material being used, pose a significant challenge for the government and the supply chain. Disposing scrap tires in landfills consumes valuable space and can also result in the penetration of pollutants and metal contents into the underground water. With recent research and advancement in the composition of tire, tires have become even more durable. When the tire treads are worn out, the majority of the time, the tire casing is still in good condition and can be reused many times. Disposing of in the landfill or incineration for energy recovery only shortens the life span of the tire casing. There are few recovery options available for EOL tires, some of which are as follows:
Reuse or repurposing the tire is the most environmentally friendly option that extends the life cycle of the tires. Tires with sufficient thread depth can be sold in the used tire market. Whole tires can also be used for a number of applications, including erosion control, playground equipment, crash barriers on the side of the highway, and boat bumpers at the docks.
Retread or recap is a process in which the worn thread of the tire is removed from the tire casing and replaced with new treads. Study shows 25 gallons of oil needed to produce a new truck tire compared to only 7 gallons of oil to retread a used one. That results in a significant reduction in greenhouse gas emission, preserving raw material, and minimizing the amount of waste disposed of in the landfills.
Waste tires are typically shredded before they are recycled. Recycled tires can be used as crumb rubber in asphalt pavement or as an aggregate in Portland cement concrete. Tires are also cut up and used in garden beds to hold in the water and to prevent weeds from growing [67].
A large percentage of scrap tires are used as fuel. Tire-derived fuel (TDF) is a fuel derived from scrap tires. Scrap tires, either shredded or as a whole, is used to supplement coal in power plants, cement kilns, paper mills, etc. Research has shown that atmospheric contamination dramatically increases when tire rubber is used as fuel [68,69,70,71,72].
Disposal in the landfill is the most harmful option due to the large volume and high percentage of void space inside the tire. Tires can trap gas that causes them to become buoyant and damage the landfill covers.
Designing an economically optimized and environmentally friendly CLSC network is a prerequisite for tire manufacturers not only to earn a profit but also to decrease waste, preserve natural resources, and save landfills and the environment, with the ultimate goal of sustainable development. This has motivated this study to take a closer look at the challenges facing the tire industry and come up with a solution that creates a win-win situation for companies in the tire industry network and the environment. Based on the reviewed papers reported in Table 1, in this study, we aimed to identify the uncertain parameters in the logistics network and find the solution for dealing with these uncertainties in order to reach a robust modeling in the logistics network. We introduce a robust multi-objective, multi-product, multi-echelon, multi-cycle, multi-capacity, green closed-loop supply-chain network under hybrid uncertainty with the aim of reducing costs, reducing environmental impacts, and increasing the social responsibilities and the reliability of the network.

3. Problem Description and Formulation

As shown in Figure 2, we are considering a multi-product and multi-echelon CLSC network consisting of suppliers, production centers, distribution centers, collection centers, retreading centers, recycling, and energy-recovery centers. In the forward flow, various raw materials are procured from the available choice of suppliers. After tires are manufactured at the production facilities, they are distributed through the distribution centers and the retailers to the end users. The location of retailers is fixed, and there might be unsatisfied demands due to uncertainty in demand. In the reverse flow, used tires are transferred to collection centers, where the initial inspection is performed, and based on the condition of the tires, they are either routed to retreading center or the recycling centers.
The main objective in the CLSC network design problem is to maximize the total profit of the supply chain by choosing the optimal number and location of facilities, their capacities, and the flow of products between the facilities, while reducing the environmental impacts and maximizing the social benefits. In real life, the product demand could be a function of economic condition; therefore, in this study, the future demand was expressed under a few scenarios. The return rates, retreading rates, recycling rates, procurement, and production costs are fuzzy; therefore, the possibilistic approach was utilized to deal with these parameters.

3.1. Problem Assumption and Notation

The main characteristics and assumptions are as follows:
  • The potential location of new manufacturing centers, distribution centers, retreading centers, and recycling centers is fixed and predefined.
  • Production, distribution, retreading, and recycling capacities are limited, and all are subject to capacity level U.
  • Capacity of the energy recovery centers is unlimited, but it has the most negative environmental impact.
  • The return, recycling, and retreading rates of new and retread tires are fuzzy.
  • Raw material procurement costs and production costs are fuzzy, and their fuzzy expected value are utilized.
  • Customer demand is stochastic and is expressed in different scenarios.
  • Ground transportation is the only transportation method used in the network.
  • There is cost associated with the collection of used tires. This provides an incentive for the end users to return the EOL tires.
  • Tires can be retreaded two times, retread-1 and retread-2, and each tire will be identified how many times it has been retreaded.
  • The value of the tires that were retreaded once (retread-1) is greater than the value of the tires that were retreaded two times (retread-2), and both are less than the value of the new tires.
  • The demand for new tires must be satisfied, and there is a penalty for unsatisfied demand. There is no obligation to meet the demand for retread tires.
Table 2 shows the list of indices, parameters and variables related to the network under study.

3.2. Problem Formulation

In the following, we introduce three objective functions z 1 , z 2 , and z 3 , under scenario ξ , representing the total profit of the network, the total environmental impacts, and the social benefits, respectively.
max Z 1 = [ d c p ( ρ p 0 Q 1 d , c . p , ξ D C + ρ p 1 Q 1 d , c , p , ξ D C + ρ p 2 Q 2 d , c , p , ξ D C ) + j p ρ p r Q j , p , ξ J ]             ( [ f F b u fc f , u F . x f , u F + d D b u fc d , u D . x d , u D + i u fc i , u I . x i . u I           + j u fc j , u J . x j , u J + k u fc k , u K . x k , u K ]           + [ s fc s S . x s S + m s f mc m , s , ξ S ˜ . Y m , s , f , ξ ] + [ f p pc p , f F ˜ . V f , p , ξ ]           + [ c i p cc p , c ( Q c , i , p , ξ C I + Q 1 c , i , p , ξ C I + Q 2 c , i , p , ξ C I ) ]           + [ k d p ( rtc 1 p , k k Q 1 k , d , p , ξ K D + rtc 2 p , k k Q 2 k , d , p , ξ K D ) ] + [ j p rc pj J Q j , p , ξ J ]           + [ w p ic p , w W Q w , p , ξ W ]           + [ m s f tr m , s , f S F . Y m , s , f , ξ + f d p tr p , f , d F D . Q f , d , p , ξ F D           + d c p tr p , d , c D C . ( Q d , c , p , ξ D C + Q 1 d , c , p , ξ D C + Q 2 d , c , p , ξ D C )           + c i p tr p , c , i C I ( Q c , i , p , ξ C I + Q 1 c , i , p , ξ C I + Q 2 c , i , p , ξ C I )           + i k p tr p , i , k I K ( Q i , k , p , ξ I K + Q 1 i , k , p , ξ I K ) + i j p tr p , i , j I J ( Q 2 i , j , p , ξ I J )           + k d p tr p , k , d K D ( Q 1 k , d , p , ξ K D + Q 2 k , d , p , ξ K D )           + k j p tr p , k , j K J ( Q 1 k , j , p , ξ K J + Q 2 k , j , p , ξ K J )           + j w p tr p , j , w J W . Q j , W , p , ξ J W ] ) ξ   Ξ
Equation (1) maximizes the total profit of the supply chain. The components of the Equation (1) include the following: the overall cost of opening new centers, cost of procurement of raw material, overall manufacturing cost, collection cost, overall retreading cost, overall recycling cost, overall energy recovery cost, total transportation cost, revenue from sales of new tires, revenue from sales of tires that are retreaded once, revenue from sales of tires that are retreaded twice, and revenue from sales of recycled material.
min Z 2 = [ f p e f , p F . V f , p , ξ ] + [ k d p e k , p K ( Q 1 k , d , p , ξ K D + Q 2 k , d , p , ξ K D ) ]          + [ j p e j , p J Q j , p , ξ J ] + [ w p e w , p W Q w , p , ξ W ]          + [ m s f p e m , s , f S F . Y m , s , f , ξ + f d p e p , f , d F D . Q f , d , p , ξ F D          + d c p e p , d , c D C   ( Q d , c , p , ξ D C + Q 1 d , c , p , ξ D C + Q 2 d , c , p , ξ D C )          + c i p e p , c , i C I   ( Q c , i , p , ξ C I + Q 1 c , i , p , ξ C I + Q 2 c , i , p , ξ C I )          + i k p e p , i , k I K ( Q i , k , p , ξ I K + Q 1 i , k , p , ξ I K )          + i j p e p , i , j I J ( Q 2 i , j , p , ξ I J )          + k d p e p , k , d K D ( Q 1 k , d , p , ξ K D + Q 2 k , d , p , ξ K D )          + k j p e p , k , j K J ( Q 1 k , j , p , ξ K J + Q 2 k , j , p , ξ K J )          + j W p e p , j , k J W . Q j , W , p , ξ J W ] ξ   Ξ
Equation (2) minimizes the overall environmental impact, which is expressed as the amount of greenhouse gases (GHGs) generated at production centers, retreading centers, recycling centers, and energy recovery centers, as well as during the transportation of raw material and finished tires between centers. Among greenhouse gases, this research is focusing on carbon dioxide (CO2) emission.
max Z 3 = [ f F b u LBR f , u F . x f , u F + d D b u LBR d , u D . x d , u D + i u LBR i , u I . x i . u I     + j u LBR j , u J . x j , u J + k u LBR k , u K . x k , u K ]   ξ   Ξ
Equation (3) maximizes the number of job opportunities created as a result of expanding the forward logistics and also opening new facilities in reverse logistics.
Constraints:
f   Y m , s , f , ξ     cap m , s S x s S ;     m M ,   s S ,   ξ Ξ
p   V f , p , ξ     u   cap u F . g f , u F ;     f F a ,   ξ Ξ
p   V f , p , ξ     u   cap u F . x f , u F ;     f F b ,   ξ Ξ  
u   x f , u F     1 ;     f F b ,   ξ Ξ  
p   r m F V f , p , ξ     s   Y m , s , f , ξ ;     m M ,   f F ,   ξ Ξ
V f , p , ξ = d Q f , d , p , ξ F D ;     f F ,   p P ,   ξ Ξ
f p Q f , d , p , ξ F D + k p Q 1 k , d , p , ξ K D + Q 2 k , d , p , ξ K D     u   cap u D . g d , u D ;     d D a ,   ξ Ξ
f p Q f , d , p , ξ F D + k p Q 1 k , d , p , ξ K D + Q 2 k , d , p , ξ K D     u   cap d D . x d , u D ;     d D b ,   ξ Ξ
u   x d , u D     1 ;     d D b
f Q f , d , p , ξ F D = c Q d , c , p , ξ D C ;     d D ,   p P ,   ξ Ξ
k Q 1 k , d , p , ξ K D = c Q 1 d , c , p , ξ D C ;     d D ,   p P ,   ξ Ξ
k Q 2 k , d , p , ξ K D = c Q 2 d , c , p , ξ D C ;     d D ,   p P ,   ξ Ξ
d Q d , c . p , ξ D C = dem c , p , ξ C ˜ ;     c C ,   p P ,   ξ Ξ
d Q 1 d , c , p , ξ D C     dem 1 c , p , ξ C ˜ ;     c C ,   p P ,   ξ Ξ
d Q 2 d , c , p , ξ D C     dem 2 c , p , ξ C ;     c C ,   p P ,   ξ Ξ
i Q c , i , p , ξ C I = α c ˜ d Q d , c , p , ξ D C ;     c C ,   p P ,   ξ Ξ
i Q 1 c , i , p , ξ C I = α 1 c ˜ d Q 1 d , c , p , ξ D C ;     c C ,   p P ,   ξ Ξ
i Q 2 c , i , p , ξ C I = α 2 c ˜ d Q 2 d , c , p , ξ D C ;     c C ,   p P ,   ξ Ξ
c p Q c , i , p , ξ C I + c p Q 1 c , i , p , ξ C I + c p Q 2 c , i , p , ξ C I     u   cap i I . x i , u I ;     i I ,   ξ Ξ
c Q c , i , p , ξ C I = k Q i , k , p , ξ I K ;     i I ,   p P ,   ξ Ξ
c Q 1 c , i , p , ξ C I = k Q 1 i , k , p , ξ I K ;     i I ,   p P ,   ξ Ξ
c Q 2 c , i , p , ξ C I = j Q 2 i , j , p , ξ I J ;     i I ,   p P ,   ξ Ξ
i p Q i , k , p , ξ I K + Q 1 i , k , p , ξ I K     u   cap u K . x k , u K ;     k K ,   ξ Ξ
d Q 1 k , d , p , ξ K D = θ 1 ˜ i Q i , k , p , ξ I K ;     k K ,   p P ,   ξ Ξ
d Q 2 k , d , p , ξ K D = θ 2 ˜ i Q 1 i , k , p , ξ I K ;     k K ,   p P ,   ξ Ξ
j Q 1 k , j , p , ξ K J   =   i Q i , k , p , ξ I K d Q 1 k , d , p , ξ K D ;     k K ,   p P ,   ξ Ξ
j Q 2 k , j , p , ξ K J   =   i Q 1 i , k , p , ξ I K d Q 2 k , d , p , ξ K D ;     k K ,   p P ,   ξ Ξ
k p Q 1 k , j , p , ξ K J + k p Q 2 k , j , p , ξ K J + i p Q 2 i , j , p , ξ I J     u   cap j J . x j , u J ;     j J ,   ξ Ξ
Q j , p , ξ J = β ˜ k Q 1 k , j , p , ξ K J + k Q 2 k , j , p , ξ K J + i Q 2 i , j , p , ξ I J ;     j J ,   p P ,   ξ Ξ
w Q j , w , p , ξ J W   =   k Q 1 k , j , p , ξ K J + k Q 2 k , j , p , ξ K J + i Q 2 i , j , p , ξ I J Q j , p , ξ J ;     j J ,   p P ,   ξ Ξ
Q w , p , ξ W = j Q j , w , p , ξ J W ;     w W ,   p P ,   ξ Ξ
x * , * . , Y m , s , f , ξ     0 , 1 V f , p , ξ ,   Q * , * . .     0
Constraint (4) ensures that the total quantity of raw material sent to tire manufacturing centers is less than or equal to the maximum capacity of each supplier for that raw material. Constraints (5) and (6) ensure that the total quantity of tires produced at each manufacturing center is less than or equal to the maximum capacity of each existing and potential manufacturing center. Constraint (7) ensures that only one production facility at capacity level u can be open at location f. Constraint (8) guarantees that the raw material needed at production centers are satisfied. Constraint (9) is an equilibrium constraint indicating the total quantity of tires transferred from each manufacturing center f to distribution center d is equal to the quantity of tires manufactured at that center. Constraints (10) and (11) ensure that the maximum capacity of the existing and potential distribution centers are not violated. Constraint (12) guarantees that only one distribution center at capacity level u can be open at any potential location d. Constraints (13)–(15) are equilibrium constraints to ensure the total quantity of new and retread (retread-1 and retread-2) tires transferred out of each distribution center is equal to the total quantity of tires enter that center. Constraint (16) ensures that the market demand for new tires is satisfied. Constraints (17) and (18) indicate that the demand for retread-1 and retread-2 tires does not have to be entirely satisfied. It should be noted that the demand parameters are stochastic-scenario-based. Constraints (19)–(21) calculate the total quantity of tires returned to each collection center. It should be noted α c ˜ , α 1 c ˜ , α 2 c ˜ are fuzzy parameters that will be discussed in the next section. Constraint (22) controls the maximum capacity of the collection centers. Constraints (23)–(25) are equilibrium constraints which create flow balance at each collection center. Constraint (26) ensures the flow of products into and out of each retreading center do not exceed the capacity of that center. Constraints (27) and (28) calculate the quantity of the tires that are being successfully retreaded at each retreading center. Moreover, ( θ 1 ˜ ) and ( θ 2 ˜ ) are fuzzy parameters. Constraints (29) and (30) calculate the total quantity of new and retread tires that cannot be retreaded at retreading centers and are being transferred to recycling centers. Constraint (31) controls the capacity of the recycling centers. Constraint (32) calculates the total quantity of products transfers to each recycling center. It should be noted that β ˜ is a fuzzy parameter. Constraint (33) creates equilibrium at recycling centers. Constraint (34) calculates the total number of tires being incinerated at each energy recovery center. Ultimately, constraint (35) enforces the binary and non-negative constraints on the corresponding decision variables.
The proposed tire CLSC model is a Mixed Integer Linear Programming (MILP) problem under hybrid uncertainty. Fuzzy uncertainty related to return rate, retreading rate, recycling rate, procurement cost, production cost, and scenario-based uncertainty related to demand parameters. In the following, we discuss a method that can simultaneously deal with these hybrid uncertainties.

4. The Proposed Approach

Robust Fuzzy/Possibilistic Stochastics Programming (RFSP) is a novel approach in mathematical programming that deals with hybrid uncertain parameters [6]. It is an expansion of possibilistic programing [43] in which the possible scenarios of the disturbance in the system are additionally being considered and addressed using stochastic scenario-based programming method [42]. An example of a disturbance in the system would be a change in customer behavior that creates a noticeable change in demand in comparison to the initial forecast. These changes can create a new scenario in the system with respect to demand parameters.
Since the proposed tire CLSC problem faces both scenario-based uncertainty related to demand parameters and fuzzy uncertainty associated with return rates, retreading rates, recycling rates, procurement, and production costs, Robust Fuzzy/Possibilistic Stochastics Programming method is proposed in this research to deal with hybrid uncertain parameters.
In the following, we first describe the basics of robust stochastic programming and possibilistic programming. Then, the proposed RFSP approach is explained; and at the end, this approach is used to control the uncertainties in the proposed tire CLSC problem.

4.1. Stochastic Programming

In general, the scenario-based stochastic programming approach can be described as follows:
max ξ Ξ π ξ . z ξ z ξ = c ξ . x ξ + d ξ y      ξ Ξ A ξ x ξ + K ξ y b ξ   ξ Ξ R y = q y Y ,   x ξ 0
where z ξ is the objective function under scenario ξ Ξ , x ξ are the scenario dependent variables, y are the variables independent of the scenario (i.e., binary variables), c ξ , d ξ , A ξ , K ξ , b ξ are the parameters under scenario ξ   Ξ , (R, q) are crisp parameters and π ξ is the probability of occurrence of scenario ξ   Ξ [7,73].
Since scenario-based programming models are not sensitive to deviations of the value of the objective function, in order to prevent excessive deviation of the of the objective function under different scenarios, the variance of the objective function is also considered. Therefore, the resulting robust stochastic programming model would be as follows:
max ξ Ξ π ξ . z ξ     λ . ξ Ξ π ξ 2 θ ξ z ξ ξ Ξ π ξ . z ξ ω ξ Ξ π ξ . Δ ξ z ξ = c ξ . x ξ + d ξ y     ξ Ξ θ ξ z ξ ξ S π ξ . z ξ     ξ Ξ A ξ x ξ + K ξ y b ξ + Δ ξ   ξ Ξ R y = q y Y ,   x ξ 0 θ ξ 0
where λ denotes the importance or weighing scale to measure the tradeoff between optimality and the cost under different scenarios. θξ denotes the absolute deviation from the mean under each scenario. Δξ represent the deviation for violation of the control constraint, and ω is the robustness coefficient or feasibility robustness multiplier which its value can be equal to the penalty cost of the shortage. It should be noted that the above model is known as the Mulvey Robust stochastic model [42], although the above model is the linear version of the model based on the research of Yu and Li [7].
In the above model, the following linearization was utilized:
min ξ Ξ π ξ z ξ ξ Ξ π ξ . z ξ Linearization min ξ Ξ π ξ 2 θ ξ z ξ ξ S π ξ . z ξ θ ξ z ξ ξ Ξ π ξ . z ξ     ξ Ξ

4.2. Possibilistic Programming

If some of the data in an optimization problem are imprecise due to insufficient data or unavailability of the information, but their approximate values can be expressed in fuzzy numbers, fuzzy mathematical programming can be used [74]. If fuzzy data in the optimization problem are expressed as Trapezoidal Fuzzy Number (TFN) and the possibility of occurrence of the possibilistic parameters are measured by Pos. and Nec., then possibilistic programming can be used, where Pos. is the maximum level of occurrence of parameters and Nec. is the minimum possibility level.
The general form of possibilistic programming can be expressed as follows:
max   Z = E c ˜ x s . t . Me ( Ax F ˜ ) 𝕡 x X
where Me(.) is a fuzzy measure; ˜ = ( f ( 1 ) , f ( 2 ) , f ( 3 ) ,   f ( 4 ) ) and c ˜ = ( c ( 1 ) , c ( 2 ) , c ( 3 ) ,   c ( 4 ) ) are the fuzzy parameters in the form of TFN, Figure 3; 𝕡 is the minimum confidence level; and E c ˜ x is the fuzzy expected value of the objective function, which can be expressed as follows:
E c ˜ x = c 1 + 2 c 2 + 2 c 3   +   c 4 6 x
In order to control fuzzy constraints, one of the most widely used fuzzy measures is the credibility measure (Cr), which is expressed as the average of the optimistic (Pos) and pessimistic (Nec) measures:
Cr . = 1 2 ( Pos . + Nec . )
And the following relationship can exist [70]:
Cr ( Ax F ˜ ) 𝕡     Ax     f 3   +   2 α 1 f 4 f 3   :   𝕡 0.5 Cr ( Ax F ˜ ) 𝕡     Ax     f 2     2 α 1 f 2 f 1   :   𝕡 0.5
One of the disadvantages of classical programming is that the 𝕡 value is determined subjectively by the decision-maker (DM) and is usually set inactively or ultimately interactively and may not be optimal. This issue was first addressed in the research conducted by Pishvaee (2012) [43].
In their proposed model, taking into account the two concepts of feasibility robustness and optimality robustness, the optimal levels of 𝕡 value are obtained proactively and objectively. This approach has been used in several studies [74,75,76,77,78,79,80].

4.3. Robust Fuzzy/Possibilistic Stochastic Programming

The RFSP approach is an extension of the combination of scenario-based stochastic programming and possibilistic programming [43] in which random scenarios and fuzzy data are controlled simultaneously [6]. The RFSP approach can be summarized as follows:
{ max ξ Ξ π ξ . z ξ     λ . ξ Ξ π ξ ( 2 θ ξ     ( z ξ ξ Ξ π ξ . z ξ ) ) ω ξ Ξ π ξ . Δ ξ     Ψ . ϑ ( 𝕡 ) z ξ = E ( c ξ . x ξ + d ξ y ) = ( c ( 1 ) + 2 c ( 2 ) + 2 c ( 3 ) + c ( 4 ) ) 6 x ξ + d ξ   y      ξ Ξ A ξ x ξ + K ξ y b ( 2 ) ξ ( 2 𝕡 1 ) ( b ( 2 ) ξ b ( 1 ) ξ ) + Δ ξ       ξ Ξ R y = q θ ξ z ξ s S π ξ . z ξ   ξ Ξ y Y ,   x ξ 0 θ ξ 0 𝕡 0.5
where Ψ is the coefficient for estimating fuzzy numbers, ϑ 𝕡 is a non-negative descending function that calculates the violation of fuzzy constraints with confidence level 𝕡 . Moreover, ϑ 𝕡 can be expressed as the difference between the substituted crisp numbers, ( b 2 ξ 2 𝕡 1 b 2 ξ b 1 ξ ) , and the worst-case value of the imprecise parameters (TFNs) can be expressed as follows:
ϑ 𝕡 = b 2 ξ 2 𝕡 1 b 2 ξ b 1 ξ b 1 ξ = b 2 ξ b 1 ξ 2 𝕡 1 b 2 ξ b 1 ξ = 1 2 𝕡 + 1 b 2 ξ b 1 ξ = 2 1 𝕡 b 2 ξ b 1 ξ
Δ ξ calculates shortage under each scenario. It should be noted that λ denotes the weighing scale to measure the tradeoff between optimality and the cost under different scenarios, θ ξ denotes the absolute deviation from the mean under each scenario. Δ ξ represents the deviation for violation of the control constraint, and ω is the model robustness coefficient or feasibility robustness multiplier.

4.4. Coping with Multi-Objective Functions

In order to solve multi-objective decision-making (MODM) optimization problems, several approaches have been proposed, such as Weighted Sum Method (WSM), Epsilon Constraint (EC), Augmented Epsilon Constraint (AEC), Goal Programing (GP), Lexicographic (Lex), and so on. The general form of a MODM problem is as follows:
Min   f 1 x , f 2 x , ,   f n x x X Min   f 1 x , f 2 x , ,   f n x x X

4.5. Epsilon Constraint

Let us assume the first objective function is the main objective, and the remaining objectives will have an upper bound of epsilon and will be presented as constraints of the model. In this case, the EC method is being used, and the resulting single objective model can be described as follows:
Min   f 1 x f i x     e i   i = 2 , 3 , , n x X
In the above model, the first objective is considered as the main objective, and the remaining objectives are considered as constraints and are bound by maximum e i . By changing the e i , different results can be obtained, which may be weakly efficient or not efficient.
This issue can be resolved with minor modification, which is known as the AEC method [81]. To better implement the AEC method, the appropriate range of epsilons e i can be obtained by using the Lex method [82]. In the AEC method, the appropriate range of e i is first determined and then the Pareto front obtained for different values of e i .

4.6. Appropriate Interval for e i with Lex Method

In order to find the appropriate interval for e i related to the objective function f i (i = 2, …, n), the following optimization problem is first being solved:
PayOff jj = Min   f j x j = 1 , 2 , ,   n x X
where PayOff jj = f j x j , * is the optimum value of f j and x j , * is the vector of decision variables that optimizes f j .
Then, with the solution that optimizes the objective function f j , the value of the remaining objective functions is obtained as follows:
PayOff ij = Min   f i x f j x   =   PayOff jj x X j = 1 , 2 , , n ;   j i
where P a y O f f j j = f i x i , j , * is the optimum value of f i and x i , j , * is the vector of decision variables that optimizes f i . Using the Lex method, the following payoff matrix is obtained:
PayOff   =   payOff ij
After calculating the payoff table for f i i = 1, …, n, the following can be defined:
Min f i   =   Min j payOff ij   =   payOff ii Max f i   =   Max j payOff ij R f i   =   Max f i   Min f i
Using the above definition, the appropriate rage for each e i can be obtained by using the Lex method:
e i     Min f i , Max f i
The value R f i is used to normalize the objectives in the AEC objective functions.

4.7. Augmented Epsilon Constraint

The AEC model can be expressed as follows:
Min   f 1 x     i = 2 n ϕ i s i f i x   +   s i = e i   i = 2 , 3 , , n x X s i 0
where s i are non-negative variables, ϕ i is a parameter for normalizing the first objective with respect to f i and ϕ i = R f 1 R f i .
In the proposed AEC method in this study, first e i     Min f i , Max f i is obtained using Lex method, then Equation (45) is solved, which provides an optimum solution.

5. Computational Experiment

To check the validity, practicality, and performance of the model, this model was applied for the design and optimization of a network in Northeast USA, as illustrated in Figure 4. The relevant data presented in Table 3, Table 4, Table 5 and Table 6 were generated from the literature, field data, interviewing experts in the field, and some organizations involved in this field, including the EPA, US Tire Manufacturers Association, Tire Retread & Repair Information Bureau (TRIB), Modern Tire Dealer (MTD), and Institute of Scrap Recycling Industries (ISRI).
In this research, there are ten types of medium/heavy truck tires considered. There are five groups of raw materials (synthetic polymers, natural rubber, steel, fillers, and curing compounds) required to manufacture a tire and five suppliers to purchase the raw material from, which are all located out of the state. There are also two existing production facilities and five distribution centers at medium-capacity level, two energy recovery centers, and forty commercial retail locations. We are considering the feasibility of expanding the network by adding two additional manufacturing facilities, five additional distribution centers, ten collection centers, four retreading centers, and five recycling centers. We believe that no tire should end up in the landfills, and by providing appropriate incentives, we can encourage customers to return the EOL tires for collection and processing. These facilities are selected such that they are in the proximity of the facilities they are linked to in the network flow and also to provide maximum coverage to the entire region. The goal is to optimize the number, location, and capacity of the facilities, increase the efficiency and profitability of the network and reduce the amount of greenhouse gases generated in the network while maximizing the total number of jobs that can be created. The demand for truck tires depends on the performance of the logistics and transportation industry in the country, which is directly correlated to the condition of the economy. We assume three scenarios for future demand, corresponding to demand at the time of poor economy (recession), during a good steady economy, and when the economy is growing. Ground transportation is the only mean of transferring the tires and raw material between the facilities. Assuming a 22-foot full truckload on each transfer, the cost of transferring one unit of tire per kilometer is $0.02 and $0.002 for one unit of raw material. The amount of Co2 generated per unit of tire transferred between facilities, using ground transportation is 0.002 kilogram per kilometer and is 0.0002 kg per kilometer for one unit of raw material being transferred.

5.1. Computational Results and Sensitivity Analysis

In order to evaluate the performance of the proposed model, first, the bi-objective version of the model is solved by considering the maximization of the profit and minimization of the environmental impact for different values of λ . The model is solved using GAMS software and a PC with the following configuration: Intel® Core™ i7-1165G7 Processor (2.80 GHz) and 16 GB DDR4 3200MHz.
Figure 5 and Figure 6 illustrate the changes in objective functions with respect to changes to the risk factor ( λ ). It can be observed that by increasing λ , meaning, by penalizing the deviation of the objective function from the mean value (reducing the risk taken by DM), both total profit and scenario variability are decreased. The optimality robustness of the model can be controlled by adjusting the risk factor ( λ ) to DM’s preferences. Reducing the risk factor to the lowest ( λ = 0) results in the highest profit and highest variability.
Table 7 shows the optimum configuration of the network, using the mean approach ( λ = 0). In this configuration, two additional production centers will be open at high capacity and five additional distribution centers will be open, namely two at low, two at medium, and one at high capacity, to cope with the demand. In the reverse logistics, six collection centers will be open at low capacity; three retreading centers will be open, one at medium and two at high capacity; and two recycling centers open at low capacity. The network flow for this configuration is illustrated in Figure 7. The penalty cost for unsatisfied demand or shortage ( ω ) is another tool at DM’s disposal to manage shortage in the system. An increase in the value of ω results in a decrease in mean profit and an increase in feasibility robustness. It can be seen in the optimum configuration of the network that the appropriate capacity levels of the open centers were determined in order to eliminate the shortage. The optimum confidence level of the fuzzy constraints ( 𝕡 ) was also determined based on the objective functions and the constraints which resulted in zero violation of fuzzy constraints ( ϑ 𝕡 ).
Further sensitivity analysis is performed to evaluate the impact of changes to the problem parameters on the optimal solution and the model behavior. In that respect, we evaluate the impact of increasing the demand by 20% without making any changes to the number of potential facilities. As illustrated in Table 8, an increase in demand has resulted in an increase in total profit and GHG production which is the result of an increase in production. The capacity of one of the distribution centers has also increased from low to high level. It can be noted that there is some unsatisfied demand due to the capacity constraint at production centers.
We further evaluate the model by increasing the demand by 20% and also considering opening an additional production facility in order to eliminate the bottle neck at production centers. As it is shown in Table 9, the total profit further increased, the new production facility is open at a high-capacity level, and the capacity of the other two production facilities decreased to medium, in order to adjust with the overall demand, which resulted in zero unsatisfied demand. There was also an increase in the capacity of distribution centers, collection centers, recycling centers, and retreading centers. This is reflected in the total number of the job opportunities created along the network.
In this section, the model is evaluated by optimizing all three objective functions, which are (1) maximizing the total profit, (2) minimizing the environmental impacts, and (3) maximizing the total number of jobs created along the network. As illustrated in Table 10, the total number of job opportunities has been increased by opening more centers at a high-capacity level. This has resulted in a reduction in the total profit due to the opening cost of the new facilities. The total amount of GHG remains unchanged, as it is directly correlated with the demand and the quantity of products being processed in each center. The overall unsatisfied demand remains at zero. The optimum confidence level of fuzzy constraints is obtained, which results in zero violation of fuzzy constraints ( ϑ 𝕡 ).

5.2. The Managerial Implication

There are many decision factors facing the DMs in designing a complex closed-loop supply-chain network. Consideration of disruption risk in addition to operational risk is one of the biggest challenges facing the DMs and the managers in supply-chain network. An example of disruption risk could be changes to economic condition, introduction of new technologies, changes in customer habit and preferences. Operational risks can be related to uncertainty in demand, procurement, production and transportation cost. RFSP is a powerful tool that provides flexibility to the DM to control these risks in the design of CLSC network simultaneously. By adjusting the risk factors (λ) to the DM’s preferences and their degree of conservatism, the optimality robustness of the network can be controlled. Feasibility robustness and the acceptable degree of violation of the constraints of the network can also be controlled by adjusting the value of Ψ and ω.
It is evident that this robust optimization method is capable of controlling scenario variability, possibilistic variability, and unsatisfied demand by adjusting the weighting scale based on the DM’s preferences.

6. Conclusions

Recovery of EOL tires has been a challenge due to the large volume being produced and the material composition of the tires that makes them difficult to decompose. With all recent advancements in recycling options for EOL tires, a significant volume of them still ends up in the landfill. This indicates the importance of a need for a systemic, sustainable CLSC network for the tire industry. In this study, we developed a mixed-integer linear programming model for a CLSC network for the tire industry by considering several recovery options and taking into account the hybrid uncertainty associated with some parameters in the network, with the aim of optimizing the total profit, while minimizing the environmental impact and maximizing the social responsibility of the network. A robust fuzzy stochastic programming method was proposed to cope with the scenario-based uncertainty associated with the demand and fuzzy-based uncertainty associated with return rate, retreading rate, recycling rate, procurement cost, and production cost simultaneously.
To illustrate the performance and practicality of the proposed model, this model was applied for the design and development of a tire CLSC network in the Northeast USA. The result indicates the efficiency and applicability of the model in real-world problems and the feasibility of designing a profitable and environmentally friendly CLSC network for the tire industry. This model and proposed solution can be further extended and applied to the supply-chain network design of other industries.
One of the limitations of exact solution methods is their inability to solve large-scale problems; therefore, for a future study, a larger-scale problem using metaheuristics methods, with a focus on the impact of low-quality, low-cost imported tires and governmental intervention, regulations, and incentives on the performance and profitability of the network, can be considered.

Author Contributions

Investigation, M.T.; Methodology, M.T.; Supervision, S.M.G.; Validation, M.T.; Writing—original draft, M.T.; Writing—review & editing, S.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the reviewers for the insightful comments and feedback to enhance the quality of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Brito, M.P.; Dekker, R. A Framework for Reverse Logistics. In Reverse Logistics; Springer: Berlin, Heidelberg, 2004; pp. 3–27. [Google Scholar]
  2. Daugherty, P.J.; Autry, C.W.; Ellinger, A.E. Reverse logistics: The relationship between resource commitment and program performance. J. Bus. Logist. 2001, 22, 107–123. [Google Scholar] [CrossRef]
  3. Linton, J.; Klassen, R.; Jayaraman, V. Sustainable supply chains: An introduction. J. Oper. Manag. 2007, 25, 1075–1082. [Google Scholar] [CrossRef]
  4. Selbes, M.; Yilmaz, O.; Khan, A.A.; Karanfil, T. Leaching of DOC, D.N., and inorganic constituents from scrap tires. Chemosphere 2015, 139, 617–623. [Google Scholar] [CrossRef] [PubMed]
  5. Lin, C.; Huang, C.L.; Shern, C.-C. Recycling waste tire powder for the recovery of oil spills. Resour. Conserv. Recycl. 2008, 52, 1162–1166. [Google Scholar] [CrossRef]
  6. Farrokh, M.; Azar, A.; Jandaghi, G.; Ahmadi, E. A novel robust fuzzy stochastic programming for closed loop supply chain network design under hybrid uncertainty. Fuzzy Sets Syst. 2018, 341, 69–91. [Google Scholar] [CrossRef]
  7. Yu, C.S.; Li, H.L. A robust optimization model for stochastic logistic problems. Int. J. Prod. Econ. 2000, 64, 385–397. [Google Scholar] [CrossRef]
  8. Pishvaee, M.S.; Torabi, S.A. A possibilistic programming approach for closed loop supply chain network design under uncertainty. Fuzzy Sets Syst. 2010, 161, 2668–2683. [Google Scholar] [CrossRef]
  9. Ilgin, M.A.; Gupta, S.M. Environmentally conscious manufacturing and product recovery (ECMPRO): A review of the state of the art. J. Environ. Manag. 2010, 91, 563–591. [Google Scholar] [CrossRef]
  10. Govindan, K.; Jafarian, A.; Khodaverdi, R.; Devika, K. Two echelon multiple vehicle location–routing problem with time windows for optimization of sustainable supply chain network of perishable food. Int. J. Prod. Econ. 2014, 15, 9–28. [Google Scholar] [CrossRef]
  11. Devika, K.; Jafarian, A.; Nourbakhsh, V. Designing a sustainable closed loop supply chain network based on bottom line approach: A comparative metaheuristic hybridization technique. Eur. J. Oper. Res. 2014, 253, 594–615. [Google Scholar] [CrossRef]
  12. Jindal, A.; Sangwan, K.S. Multi-objective fuzzy mathematical modelling of closed-loop supply chain considering economical and environmental factors. Ann. Oper. Res. 2017, 257, 95–120. [Google Scholar] [CrossRef]
  13. Sarkar, B.; Tayyab, M.; Kim, N.; Habib, M.S. Optimal production delivery policies for supplier and manufacturer in a constrained closed-loop supply chain for returnable transport packaging through metaheuristic approach. Comput. Ind. Eng. 2019, 135, 987–1003. [Google Scholar] [CrossRef]
  14. Sarkar, B.; Ullah, M.; Kim, N. Environmental and economic assessment of closed-loop supply chain with remanufacturing and returnable transport items. Comput. Ind. Eng. 2017, 111, 148–163. [Google Scholar] [CrossRef]
  15. Dulman, M.; Gupta, S.M. Evaluation of Maintenance and EOL Operation Performance of Sensor-Embedded Laptops. Logistics 2018, 1, 3. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, L.; Gupta, S.M. A Pricing and Acquisition Strategy for New and Remanufactured High-Technology Products. Logistics 2019, 3, 8. [Google Scholar] [CrossRef] [Green Version]
  17. Zhou, L.; Gupta, S.M. Value depreciation factors for new and remanufactured high-technology products: A case study on iPhones and iPads. Int. J. Prod. Res. 2020, 58, 7218–7249. [Google Scholar] [CrossRef]
  18. Aldoukhi̇, M.; Gupta, S. A robust closed loop supply chain network design under different carbon emission policies. Pamukkale University. J. Eng. Sci. 2019, 25, 1020–1032. [Google Scholar]
  19. Aldoukhi, M.; Gupta, S.M. Use of Maximal Covering Location Problem to Design a Closed Loop Supply Chain Network Under Product Substitution. Appl. Manag. Sci. 2020, 20, 71–96. [Google Scholar]
  20. Fadhel, A.W.; Gupta, S.M. Carbon Emissions and Energy Balance in the Design of a Sustainable Food Waste Network. Int. J. Ind. Eng. Oper. Manag. 2020, 2, 63–77. [Google Scholar]
  21. Gupta, E.V.; Mogale, D.G.; Tiwari, M.K. Optimal Control of Production and Maintenance Operations in Smart Custom Manufacturing Systems with Multiple Machines. IFAC PapersOnLine 2019, 52, 241–246. [Google Scholar] [CrossRef]
  22. De, M.; Giri, B.C. Modeling a closed-loop supply chain with a heterogeneous fleet under carbon emission reduction policy. Transp. Res. Part E Logist. Transp. Rev. 2020, 133, 101813. [Google Scholar] [CrossRef]
  23. Mogale, D.G.; Lahoti, G.; Jha, S.B.; Shukla, M.; Kamath, N.; Tiwari, M.K. Dual Market Facility Network Design under Bounded Rationality. Algorithms 2018, 11, 54. [Google Scholar] [CrossRef] [Green Version]
  24. Oliveira, L.S.; Machado, R.L. Application of optimization methods in the closed-loop supply chain: A literature review. J. Comb. Optim. 2021, 41, 357–400. [Google Scholar] [CrossRef]
  25. Zahedi, A.; Salehi-Amiri, A.; Hajiaghaei-Keshteli, M.; Diabat, A. Designing a closed-loop supply chain network considering multi-task sales agencies and multi-mode transportation. Soft Comput. 2021, 1–33. [Google Scholar]
  26. Guarnieri, P.; Camara, E.; Silva, L.; Vieira, B.O. How to Assess Reverse Logistics of e-Waste Considering a Multicriteria Perspective: A Model Proposition. Logistics 2020, 4, 25. [Google Scholar] [CrossRef]
  27. Vieira, B.O.; Guarnieri, P.; Nofal, R.; Nofal, B. Multi-Criteria Methods Applied in the Studies of Barriers Identified in the Implementation of Reverse Logistics of E-Waste: A Research Agenda. Logistics 2020, 4, 11. [Google Scholar] [CrossRef]
  28. Mogale, D.G.; Kumar, S.K.; Tiwari, M.K. Green food supply chain design considering risk and post-harvest losses: A case study. Ann. Oper. Res. 2020, 295, 257–284. [Google Scholar] [CrossRef]
  29. Mogale, D.G.; Cheikhrouhou, N.; Tiwari, M.K. Modelling of sustainable food grain supply chain distribution system: A bi-objective approach. Int. J. Prod. Res 2019. [Google Scholar] [CrossRef] [Green Version]
  30. Ullah, M.; Asghar, I.; Zahid, M.; Omair, M.; AlArjani, A.; Sarkar, B. Ramification of remanufacturing in a sustainable three-echelon closed-loop supply chain management for returnable products. J. Clean. Prod. 2021, 290, 125609. [Google Scholar] [CrossRef]
  31. Ghadge, A.; Kara, M.; Mogale, D.G.; Choudhary, S.; Dani, S. Sustainability implementation challenges in food supply chains: A case of U.K. artisan cheese producers. Prod. Plan. Control. 2020. [Google Scholar] [CrossRef]
  32. Gupta, S.M.; Güngör, A.; Govindan, K.; Özceylan, E.; Can Berk, K.; Piplani, R. Responsible & sustainable manufacturing. Int. J. Prod. Res. 2020, 58, 7181–7182. [Google Scholar]
  33. Sazvar, Z.; Zokaee, M.; Tavakkoli-Moghaddam, R.; Al-sadat Salari, S.; Nayeri, S. Designing a sustainable closed-loop pharmaceutical supply chain in a competitive market considering demand uncertainty, manufacturer’s brand and waste management. Ann. Oper. Res. 2021, 1–32. [Google Scholar] [CrossRef]
  34. Meade, L.; Sarkis, J. A conceptual model for selecting and evaluating third-party reverse logistics providers. Supply Chain Manag. Int. J. 2002, 7, 283–295. [Google Scholar] [CrossRef]
  35. Ahmadi, A.; Amin, S. An integrated chance-constrained stochastic model for a mobile phone closed-loop supply chain network with supplier selection. J. Clean. Prod. 2019, 226, 988–1003. [Google Scholar] [CrossRef]
  36. Fazli-Khlaf, M.; Chaharsoohi, S.; Pishvaee, S. A new robust possibilistic programming model for reliable supply chain network design: A case study of lead- acid battery supply chain. RAIRO Oper. Res. 2019, 53, 1489–1512. [Google Scholar] [CrossRef] [Green Version]
  37. Sabouhi, F.; Pishvaee, M.S.; Jabalameli, M.S. Resilient Supply Chain Design under Operational and Disruption Risks Considering Quantity Discount: A Case Study of Pharmaceutical Supply Chain. Comput. Ind. Eng. 2018, 126, 657–672. [Google Scholar] [CrossRef]
  38. Mula, J.; Poler, R.; Garcia, J.P. MRP with flexible constraints: A fuzzy mathematical programming approach. Fuzzy Sets Syst. 2006, 157, 74–97. [Google Scholar] [CrossRef]
  39. Mula, J.; Poler, R.; Garcia-Sabater, J.P. Material Requirement Planning with fuzzy constraints and fuzzy coefficients. Fuzzy Sets Syst. 2007, 158, 783–793. [Google Scholar] [CrossRef]
  40. Pishvaee, M.S.; Razmi, J. Environmental supply chain network design using multi-objective fuzzy mathematical programming. Appl. Math. Model. 2012, 36, 3433–3446. [Google Scholar] [CrossRef]
  41. Pishvaee, M.S.; Torabi, S.A.; Razmi, J. Credibility-based fuzzy mathematical programming model for green logistics design under uncertainty. Comput. Ind. Eng. 2012, 62, 624–632. [Google Scholar] [CrossRef]
  42. Mulvey, J.M.; Vanderbei, R.J.; Zenios, S.A. Robust optimization of large-scale systems. Oper. Res. 1995, 43, 264–328. [Google Scholar] [CrossRef] [Green Version]
  43. Pishvaee, M.S.; Razmi, J.; Torabi, S.A. Robust possibilistic programming for socially responsible supply chain network design: A new approach. Fuzzy Sets Syst. 2012, 206, 1–20. [Google Scholar] [CrossRef]
  44. Mousazadeh, M.; Torabi, S.A.; Pishvaee, M.S.; Abolhassani, F. Health service network design: A robust possibilistic approach. Int. Trans. Oper. Res. 2018, 25, 337–373. [Google Scholar] [CrossRef]
  45. Tozanli, O.; Duman, G.M.; Kongar, E.; Gupta, S.M. Environmentally Concerned Logistics Operations in Fuzzy Environment: A Literature Survey. Logistics 2017, 1, 4. [Google Scholar] [CrossRef] [Green Version]
  46. Vahdat, V.; Vahdatzad, M.A. Accelerated Benders’ Decomposition for Integrated Forward/Reverse Logistics Network Design under Uncertainty. Logistics 2017, 1, 11. [Google Scholar] [CrossRef] [Green Version]
  47. Prakash, S.; Kumar, S.; Soni, G.; Jain, V.; Rathore, A. Closed-loop supply chain network design and modeling under risks and demand uncertainty: An integrated robust optimiza-tion approach. Ann. Oper. Res. 2020, 290, 837–864. [Google Scholar] [CrossRef]
  48. Yolmeh, A.; Saif, U. Closed-loop supply chain network design integrated with assembly and disassembly line balancing under uncertainty: An enhanced decomposition approach. Int. J. Prod. Res. 2020, 1–18. [Google Scholar] [CrossRef]
  49. Wang, W.; Zhang, Y.; Zhang, W.; Gao, G.; Zhang, H. Incentive mechanisms in a green supply chain under demand uncertainty. J. Clean. Prod. 2018, 279, 123636. [Google Scholar] [CrossRef]
  50. Tordecilla, R.D.; Juan, A.A.; Montoya-Torres, J.R.; Quintero-Araujo, C.L.; Panadero, J. Simulation-optimization methods for designing and assessing resilient supply chain networks under uncertainty scenarios: A review. Simul. Model. Pract. Theory 2021, 106, 102166. [Google Scholar] [CrossRef]
  51. Ferrer, G. The economics of tire remanufacturing. Resour. Conserv. Recycl. 1997, 19, 221–255. [Google Scholar] [CrossRef]
  52. Sasikumar, P.; Kannan, G.; Noorul Haq, A. A multi-echelon reverse logistics network design for product recovery—A case of truck tire remanufacturing. Int. J. Adv. Manuf. Technol. 2010, 49, 1223–1234. [Google Scholar] [CrossRef]
  53. Chopra, S.; Meindl, P. Supply Chain Management: Strategy, Planning and Operation; Pearson: Boston, MA, USA, 2013; Volume 232. [Google Scholar]
  54. Yadollahinia, M.; Teimoury, E.; Paydar, M. Tire forward and reverse supply chain design considering customer relationship management. Resour. Conserv. Recycl. 2018, 138, 215–228. [Google Scholar] [CrossRef]
  55. Amin, S.H.; Zhang, G.; Akhtar, P. Effects of uncertainty on a tire closed-loop supply chain network. Exp. Syst. Appl. 2017, 73, 82–91. [Google Scholar] [CrossRef] [Green Version]
  56. Subulan, K.; Tas, A.S.; Baykasoğlu, A. Designing an environmentally conscious tire closed-loop supply chain network with multiple recovery options using interactive fuzzy goal programming. Appl. Math. Model. 2015, 39, 2661–2702. [Google Scholar] [CrossRef]
  57. Derakhshan, Z.; Ghaneian, M.; Mahvi, A.; Gea, O.C.; Faramarzian, M.; Dehghani, M.; Ferrante, M. A new recycling technique for the waste tires reuse. Environ. Res. 2017, 158, 462–469. [Google Scholar] [CrossRef] [PubMed]
  58. O’Brien, B.O.; North, B. A Review of Waste Tyre Gasification. J. Environ. Chem. Eng. 2017, 5, 10. [Google Scholar] [CrossRef]
  59. Fathollahi-Fard, A.M.; Hajiaghaei-Keshteli, M. Hybrid optimizers to solve a tri-level programming model for a tire closed-loop supply chain network design. Appl. Soft Comput. 2018, 70, 701–722. [Google Scholar] [CrossRef]
  60. Pedram, A.; Yusoff, N.; Udoncy, O.; Mahat, A.; Pedram, P.; Babolola, A. Integrated forward and reverse supply chain: A tire case study. Waste Manag. 2017, 60, 460–470. [Google Scholar] [CrossRef]
  61. Lokesh, K.S.; Jain, P.K.; Sharma, A.K. Tactical supply chain planning for Tyre remanufacturing considering carbon tax policy. Int. J. Adv. Manuf. Technol. 2018, 97, 1505–1528. [Google Scholar]
  62. Lokesh, K.S.; Jain, P.K.; Sharma, A.K. A fuzzy goal programme with carbon tax for Brownfield tyre remanufacturing supply chain planning. J. Clean. Prod. 2018, 198, 737–753. [Google Scholar]
  63. Sahebjamnia, N.; Fathollahi-Fard, A.; 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]
  64. Abdolazimi, O.; Esfandarani, M.; Salehi, M.; Shishebori, D. Robust design of a multi-objective closed-loop supply chain by integrating on-time delivery, cost, and environmental aspects, case study of a Tire Factory. J. Clean. Prod. 2020, 264, 121566. [Google Scholar] [CrossRef]
  65. Mehrjerdi, Y.Z.; Shafiee, M. A Resilient and sustainable closed-loop supply chain using multiple sourcing and information sharing strategies. J. Clean. Prod. 2020, 289, 125141. [Google Scholar] [CrossRef]
  66. Lokesh, K.S.; Jain, P.K.; Sharma, A.K. A fuzzy goal programme–based sustainable Greenfield supply network design for tyre retreading industry. Int. J. Adv. Manuf. Technol. 2020, 108, 2855–2880. [Google Scholar]
  67. Gupta, V.K.; Nayak, A.; Agarwal, S.; Chaudhary, M.; Tyagi, I. Removal of Ni (II) ions from water using scrap tire. J. Mol. Liq. 2014, 190, 215–222. [Google Scholar] [CrossRef]
  68. Derakhshan, Z.; Ehrampoush, M.H.; Faramarzian, M.; Ghaneian, M.T.; Mahvi, A.H. Waste tire chunks as a novel packing media in a fixed-bed sequence batch reactor: Volumetric removal modeling. Desalin. Water Treat. 2017, 64, 40–47. [Google Scholar] [CrossRef]
  69. Battista, M.; Gobetti, A.; Agnelli, S.; Ramorino, G. Post-consumer tires as a valuable resource: Review of different types of material recovery. Environ. Technol. Rev. 2021, 10, 1–25. [Google Scholar] [CrossRef]
  70. Conti, G.O.; Heibati, B.; Kloog, I.; Fiore, M.; Ferrante, M. A review of AirQ Models and their applications for forecasting the air pollution health outcomes. Environ. Sci. Pollut. Res. Int. 2017, 24, 6426–6445. [Google Scholar] [CrossRef]
  71. Sciacca, S.; Conti, G.O. Mutagens and carcinogens in drinking water. Med. J. Nutr. 2009, 2, 157–162. [Google Scholar]
  72. Dehghani, M.; Keshtgar, L.; Javaheri, M.R.; Derakhshan, Z.; Conti, G.O.; Zuccarello, P.; Ferrante, M. The effects of air pollutants on the mortality rate of lung cancer and leukemia. Mol. Med. 2017, 15, 3390–3397. [Google Scholar] [CrossRef] [Green Version]
  73. Salehi, F.; Mahootchi, M.; Husseini, S.M.M. Developing a robust stochastic model for designing a blood supply chain network in a crisis: A possible earthquake in Tehran. Ann. Oper. Res. 2019, 283, 679–703. [Google Scholar] [CrossRef]
  74. Kacprzyk, J.; Orlovski, S.A. Optimization Models Using Fuzzy Sets and Possibility Theory; Springer Science & Business Media: Berlin, Germany, 2013; Volume 4. [Google Scholar]
  75. Xu, J.; Zhou, X. Approximation based fuzzy multi-objective models with expected objectives and chance constraints: Application to earth-rock work allocation. Inf. Sci. 2013, 238, 75–95. [Google Scholar] [CrossRef]
  76. Zahiri, B.; Jula, P.; Tavakkoli-Moghaddam, R. Design of a pharmaceutical supply chain network under uncertainty considering perishability and substitutability of products. Inf. Sci. 2018, 423, 257–283. [Google Scholar] [CrossRef]
  77. Savadkoohi, E.; Mousazadeh, M.; Torabi, S.A. A possibilistic location-inventory model for multi-period perishable pharmaceutical supply chain network design. Chem. Eng. Res. Des. 2018, 138, 490–505. [Google Scholar] [CrossRef]
  78. Heydari, J.; Sabbaghnia, A.; Razmi, J. A dynamic bi-objective model for after disaster blood supply chain network design; a robust possibilistic programming ap-proach. J. Ind. Syst. Eng. 2018, 11, 16–28. [Google Scholar]
  79. Sadghiani, N.S.; Torabi, S.; Sahebjamnia, N. Retail supply chain network design under operational and disruption risks. Transp. Res. Part E Logist. Transp. Rev. 2015, 75, 95–114. [Google Scholar] [CrossRef]
  80. Zhang, P.; Zhang, W.-G. Multiperiod mean absolute deviation fuzzy portfolio selection model with risk control and cardinality constraints. Fuzzy Sets Syst. 2014, 255, 74–91. [Google Scholar] [CrossRef]
  81. Mavrotas, G. Effective implementation of the ε-constraint method in multi-objective mathematical programming problems. Appl. Math. Comput. 2009, 213, 455–465. [Google Scholar] [CrossRef]
  82. Aghaei, J.; Amjady, N.; Shayanfar, H.A. Multi-objective electricity market clearing considering dynamic security by lexicographic optimization and augmented epsilon constraint method. Appl. Soft Comput. 2011, 11, 3846–3858. [Google Scholar] [CrossRef]
Figure 1. Tire composition.
Figure 1. Tire composition.
Logistics 05 00020 g001
Figure 2. Schematic view of the proposed CLSC.
Figure 2. Schematic view of the proposed CLSC.
Logistics 05 00020 g002
Figure 3. Function of Trapezoidal Fuzzy Numbers.
Figure 3. Function of Trapezoidal Fuzzy Numbers.
Logistics 05 00020 g003
Figure 4. Location of the existing and potential facilities.
Figure 4. Location of the existing and potential facilities.
Logistics 05 00020 g004
Figure 5. Value of objective functions for different value of λ.
Figure 5. Value of objective functions for different value of λ.
Logistics 05 00020 g005
Figure 6. Mean value and scenario variability for different value of λ.
Figure 6. Mean value and scenario variability for different value of λ.
Logistics 05 00020 g006
Figure 7. Optimum network configuration and product flow.
Figure 7. Optimum network configuration and product flow.
Logistics 05 00020 g007
Table 1. Research comparison.
Table 1. Research comparison.
Network TypeSupplier SelectionMulti Raw MaterialMulti ProductMulti Capacity LevelNetwork Expansion
ForwardReverseCLSC
Sasikumar et al. (2010) [52] x
Pishvaee et al. (2012) [43]x
Subulan et al. (2015) [56] xxxxx
Jindal et al. (2016) [12] xx
Amin et al. (2017) [55] xx x
Pedram et al. (2018) [60] x x
Fazli-Khalaf et al. (2018) [36] xxx x
Fathollahi-Fard et al. (2018) [59] xx x
Sabouhi et al. (2018) [37]x xxxx
Lokesh et al. (2018a) [61] x xx
Lokesh et al. (2018b) [62] x x
Sahebjamnia et al. (2018) [63] xx x
Yadollahinia et al. (2018) [54] x x
Ahmadi and Amin (2019) [35] xxxx
Mogale et al. (2019) [29]x
De et al. (2020) [22] x x
Abdolazimi et al. (2020) [64] xxxx
Mehrjerdi et al. (2020) [65] xxxx
Lokesh et al. (2020) [66] x xxx
This Paper xxxxxx
ObjectiveGreenUncertainty TypeRobust OptimizationMulti Cycle NetworkBuyback IncentivePenalty CostCase Study
SingleMultiDeterministicStochasticFuzzy
Sasikumar et al. (2010) [52]x x Tire
Pishvaee et al. (2012) [43] x xx xMedical
Subulan et al. (2015) [56] x x x Tire
Jindal et al. (2016) [12] xx x
Amin et al. (2017) [55]x x Tire
Pedram et al. (2018) [60]x x xTire
Fazli-Khalaf et al. (2018) [36]x xx xBattery
Fathollahi-Fard et al. (2018) [59] xxx
Sabouhi et al. (2018) [37]x x Pharma
Lokesh et al. (2018a) [61] xxx Tire
Lokesh et al. (2018b) [62] xx x Tire
Sahebjamnia et al. (2018) [63] xxx Tire
Yadollahinia et al. (2018) [54] x x x Tire
Ahmadi and Amin (2019) [35] xx x Phone
Mogale et al. (2019) [29] xxx Grain
De et al. (2020) [22] xxx
Abdolazimi et al. (2020) [64] xx x x Tire
Mehrjerdi et al. (2020) [65] xx x Tire
Lokesh et al. (2020) [66] xx x Tire
This Paper xx xxxxxxTire
Table 2. The list of model indices, parameters and variables.
Table 2. The list of model indices, parameters and variables.
ΞSet of Possible Scenarios
S Set of raw material suppliers
M Set of raw materials
P Set of types of tires
F a Set of existing production centers
F b Set of potential production centers
D a Set of existing distribution centers
D b Set of potential distribution centers
C Set of customer delivery points
I Set of collection centers
J Set of recycling centers
K Set of retreading centers
W Set of energy recovery centers
U Set of capacity levels
Parameters: Tilde symbol (∼) is to distinguish the uncertain parameters.
f c f , u F Fixed opening cost of production center at location f at capacity level u
f c d , u D Fixed opening cost of distribution center at location d at capacity level u
f c i , u I Fixed opening cost of collection center at location i at capacity level u
f c j , u J Fixed opening cost of recycling center at location j at capacity level u
f c k , u K Fixed opening cost of retreading center at location k at capacity level u
f c s S Fixed ordering cost of raw material from supplier s
m c m , s S ˜ Unit cost of raw material m from supplier s (fuzzy)
p c p , f F ˜ Production cost of a unit of tire p at production center f (fuzzy)
c c p c Buyback price of a unit of EOL tire p from customer c
r t c 1 p , k k Retreading cost of a unit of tire p for first time at retreading center k
r t c 2 p , k k Retreading cost of a unit of tire p for second time at retreading center k
r c p , j J Recyclincost of a unit of tire p at recycling center j
i c p , w W Incineration cost of a unit of tire p at incineration center w
t r m , s , f S F Transportation cost for unit of raw material m from center s to center f
t r p , f , d F D Transportation cost for a unit of tire p from center f to d
t r p , d , c D C Transportation cost for a unit of tire p from center d to c
t r p , c , i C I Transportation cost for a unit of tire p from center c to i
t r p , i , j I J Transportation cost for a unit of tire p from center i to j
t r p , i , k I K Transportation cost for a unit of tire p from center i to k
t r p , k , d K D Transportation cost for a unit of tire p from center k to d
t r p , k , J K J Transportation cost for a unit of tire p from center k to j
t r p , j , W J W Transportation cost for a unit of tire p from center j to w
c a p m , s S Maximum capacity of the supplier s for material m
c a p u F Maximum c a p m , s S capacity of production center f at capacity level u
g f , u F Binary parameter defining capacity level of existing center f F a
c a p u D Maximum cap F a city of distribution center d at capacity level u
g d , u D Binary parameter defining capacity level of existing center d D a
c a p u I Maximum capacity of collection center i at capacity level u
c a p u K Maximum capacity of retreading center k at capacity level u
c a p u J Maximum capacity of recycling center j at capacity level u
e f , p F GHG generated at center f for producing a unit of tire p
e j , p J GHG generated at center j for recycling a unit of tire p
e k , p K GHG generated at center k for retreading a unit of tire p
e w , p W GHG generated at center w for incineration of a unit of tire p
e m , s , f S F GHG generated for transporting a unit of material m from center s to f
e f , d , p F D GHG generated for transporting a unit of tire p from center f to d
e d , c , p D C GHG generated for transporting a unit of tire p from center d to c
e c , i , p C I GHG generated for transporting a unit of tire p from center c to i
e i , k , p I K GHG generated for transporting a unit of tire p from center i to k
e k , j , p K J GHG generated for transporting a unit of tire p from center k to j
e k , d , p K D GHG generated for transporting a unit of tire p from center k to d
e i , j , p I J GHG generated for transporting a unit of tire p from center i to j
e j , W , p J W GHG generated for transporting a unit of tire p from center j to w
r m , p F Quantity of material m needed to produce one unit of tire p
d e m c , p , ξ C ˜ Market demand for new tires p (scenario-based)
d e m 1 c , p , ξ C ˜ Market demand for retread-1 tires p (scenario-based)
d e m 2 c , p , ξ C ˜ Market demand for retread-2 tires p (scenario-based)
α c ˜ Return rate of EOL new tire from center c (fuzzy)
α 1 c ˜ Return rat α c ˜ of EOL retread-1 tire from center c (fuzzy)
α 2 c ˜ Return rate of EOL retread-2 tire from center c (fuzzy)
θ 1 ˜ Retread rate of EOL new tires (fuzzy)
θ 2 ˜ Retread rate of EOL retread-1 tires (fuzzy)
β ˜ Recycling rate (fuzzy)
ρ p 0 Unit price of the new tire p
ρ p 1 Unit price of the retread-1 tire p
ρ p 2 Unit price of the retread-2 tire p
ρ p r Value of the recycled material for unit tire p
L B R f , u F Number of employees at new manuf. center f F b at capacity u
L B R d , u D Number of employees at new dist. center d D b at capacity u
L B R k , u K Number of employees at new retreading center k at capacity u
L B R i , u I Number of employees at new collection center i at capacity u
L B R j , u J Number of employees at new recycling center j at capacity u
Variables
Y m , s , f , ξ Quantity of material m shipped form supplier s to f under scenario ξ
V f , p , ξ Total production quantity of tire p at the center f under scenario ξ
Q f , d , p , ξ F D Quantity of new tire p shipped from center f to d under scenario ξ
Q d , c , p , ξ D C Quantity of new tire p shipped from center d to c under scenario ξ
Q 1 d , c , p , ξ D C Quantity of retread-1 tire p shipped from center d to c under scenario ξ
Q 2 d , c , p , ξ D C Quantity of retread-2 tire p shipped from center d to c under scenario ξ
Q c , i , p , ξ C I Quantity of new tire p shipped from center c to i under scenario ξ
Q 1 c , i , p , ξ C I Quantity of retread-1 tire p shipped from center c to i under scenario ξ
Q 2 c , i , p , ξ C I Quantity of retread-2 tire p shipped from center c to i under scenario ξ
Q i , k , p , ξ I K Quantity of new tire p shipped from center i to k under scenario ξ
Q 1 i , k , p , ξ I K Quantity of retread-1 tire p shipped from center i to k under scenario ξ
Q 2 i , j , p , ξ I J Quantity of retread-2 tire p shipped from center i to j under scenario ξ
Q 1 k , d , p , ξ K D Quantity of retread-1 tire p shipped from center k to d under scenario ξ
Q 2 k , d , p , ξ K D Quantity of retread-2 tire p shipped from center k to d under scenario ξ
Q 1 k , j , p , ξ K J Quantity of retread-1 tire p shipped from center k to j under scenario ξ
Q 2 k , j , p , ξ K J Quantity of retread-2 tire p shipped from center k to j under scenario ξ
Q j , p , ξ J Total quantity of recycled tire p at recycling center j under scenario ξ
Q j , W , p , ξ J W Quantity of tire p shipped from center j to w under scenario ξ
Q w , p , ξ W Total quantity d D b f tire p incinerated at center w under scenario ξ
Binary variables:
x f , u F Indicating if center f at capacity level u opened at location f F b or not
x d , u D Indicating if center d at capacity level u opened at location d D b or not
x i , u I Indicating if center i at capacity level u opened or not
x k , u K Indicating if center k at capacity level u opened or not
x j , u J Indicating if center j at capacity level u opened or not
x s S Indicating if raw material is purchased from supplier s or not
Table 3. Values for scenario-based parameters (demands).
Table 3. Values for scenario-based parameters (demands).
Scenarios ( ξ )RecessionGoodGrowing
Scenario Probability0.200.500.30
dem c , p , ξ C ˜ U (5000–5100)U (5100–5300)U (5300–5500)
dem 1 c , p , ξ C ˜ U (2400–2700)U (2700–2800)U (2800–2900)
dem 2 c , p , ξ C ˜ U (500–600)U (600–700)U (700–800)
Table 4. Values for fuzzy parameters.
Table 4. Values for fuzzy parameters.
Parameters c 1 c 2 c 3 c 4
α c   ˜ U (0.7–0.75)U (0.75–0.8)U (0.80–0.85)U (0.85–0.9)
α 1 c   ˜ U (0.65–0.7)U (0.7–0.75)U (0.75–0.8)U (0.8–0.85)
α 2 c   ˜ U (0.55–0.6)U (0.6–0.65)U (0.65–0.7)U (0.7–0.75)
θ 1 ˜ U (0.6–0.65)U (0.65–0.7)U (0.7–0.75)U (0.75–0.8)
θ 2 ˜ U (0.5–0.55)U (0.55–0.6)U (0.6–0.65)U (0.65–0.7)
β ˜ U (0.88–0.91)U (0.91–0.93)U (093–0.96)U (0.96–0.99)
pc p , f F   ˜ U ($100–$110)U ($110–$120)U ($120–$130)U ($130–$140)
mc m = 1 , s S ˜ U ($9, $9.5)U ($9.5, $10)U ($10, $10.5)U ($10.5, $11)
mc m = 2 , s S ˜ U ($43, $44)U ($44, $45)U ($45, $46)U ($46, $47)
mc m = 3 , s S ˜ U ($9, $9.5)U ($9.5–$10)U ($10–$10.5)U ($10.5–$11)
mc m = 4 , s S ˜ U ($9, $9.5)U ($9.5–$10)U ($10–$10.5)U ($10.5k–$11)
mc m = 5 , s S ˜ U ($10, $10.5)U ($10.5, $11)U ($11, $11.5)U ($11.5, $12)
Table 5. Parameters at multi-capacity levels.
Table 5. Parameters at multi-capacity levels.
Capacity Level (U)
LowMediumHigh
fc f , u F U ($3M–$4M)U ($4M–$6M)U ($6M–$10M)
fc d u D U ($0.2M–$0.3M)U ($0.3M–$0.5M)U ($0.5M–$0.7M)
fc iu I U ($0.2M–$0.3M)U ($0.3M–$0.4M)U ($0.4M–$0.5M)
fc j , u J U ($0.3M, $0.4M)U ($0.4M, $0.5M)U ($0.5M–$0.7M)
fc k , u K U ($0.8M–$1.5M)U ($1.5M–$2.5M)U ($2.5M–$4M)
cap u F U(200K–400K)U (400K-600K)U (600K–1000K)
cap u D U (100K–200K)U (200K–300K)U (300K–400K)
cap u I U (200K–400K)U (400K–600K)U (600K–800K)
cap u K U (100K–400K)U (400K–600K)U (600K–800K)
cap u J U (100K–400K)U (400K–600K)U (600K–800K)
L B R f , u F U (100–150)U (200–300)U (300–500)
L B R d , u D U (30–60)U (60–100)U (100–150)
L B R k , u K U (30–50)U (50–70)U (70–100)
L B R i , u I U (10–30)U (30–50)U (50–70)
L B R j , u J U (10–30)U (30–50)U (50–70)
Table 6. Other parameters.
Table 6. Other parameters.
r m , p F m1, U (4–6); m2, U (14–17); m3, (4–6); m4, U (9–12); m5, U (10–12)
fc s S U ($1000–$1200)
cc p , c U ($30–$50)
rtc 1 p , k k U ($40–$50)
rtc 2 p , k k U ($45–$55)
rc p , j J U ($5–$10)
ic p , w W U ($3–$8)
e f F U (9–10) kg
e w W U (10–12) kg
e k K U (6–7) kg
e j J U (1–2) kg
ρ p 0 U ($450–$550)
ρ p 1   U ($200–$300)
ρ p 2 U ($150–$200)
ρ p r U ($25–$30)
Table 7. Summary of optimum results, using mean approach.
Table 7. Summary of optimum results, using mean approach.
Value of Coefficients λ = 0 ,   ω = 0 ,   Ψ = 0
Objective 1 (profit, million dollar)1004
Objective 2 (total GHG =kg)324,893
Objective 3 (total job opportunities)1600
Production CentersProduction Centers 3 and 4 open at high capacity
Distribution CentersDistribution Centers 6 and 7 open at low capacity
Distribution Centers 8 and 9 open at medium ca-pacity
Distribution Center 10 open at high capacity
Collection CentersCollection Centers 1, 2, 4, 6, 8, and 9 open at low capacity
Retreading CentersRetreading Center 2 open at medium capacity
Retreading Centers 1 and 3 open at high capacity
Recycling CentersRecycling Centers 2 and 3 open at low capacity
Suppliers being utilizedSuppliers 1, 3, 4, and 5
θ ξ (scenario variability)164.926
Δ ξ (shortage)0
ϑ 𝕡 0
𝕡 (confidence level of fuzzy constraints)𝕡1 = 1/𝕡2 = 0.5/𝕡3 = 0.5/𝕡4 = 1/𝕡5 = 1/𝕡6 = 1
Table 8. Summary of optimum result (demand increased by 20%).
Table 8. Summary of optimum result (demand increased by 20%).
Value of Coefficients λ = 0 ,   ω = 0 ,   Ψ = 0
Objective 1 (profit, million dollar)1115
Objective 2 (total GHG = kg)348,253
Objective 3 (total job opportunities)1715
Production CentersProduction Centers 3 and 4 open at high capacity
Distribution CentersDistribution Centers 6 open at low capacityDistribution Centers 8 and 9 open at medium capacityDistribution Centers 7 and10 open at high capacity
Collection CentersCollection Centers 1, 2, 4, 6, 8, and 9 open at low capacity
Retreading CentersRetreading Center 2 open at medium capacityRetreading Centers1 and 3 open at high capacity
Recycling CentersRecycling Centers 2 and 3 open at low capacity
Suppliers being utilizedSuppliers 1, 3, 4, and 5
θ ξ (scenario variability)184.695
Δ ξ (shortage)107,151
ϑ 𝕡 0
𝕡 (confidence level of fuzzy constraints)𝕡1 = 1/𝕡2 = 1/𝕡3 = 1/𝕡4 = 1/𝕡5 = 1/𝕡6 = 1
Table 9. Summary of optimum result (demand increase by 20% and considering opening an additional production center).
Table 9. Summary of optimum result (demand increase by 20% and considering opening an additional production center).
Value of Coefficients λ = 0 ,   ω = 0 ,   Ψ = 0
Objective 1 (profit, million dollar)1160
Objective 2 (total GHG, kg)358,321
Objective 3 (total job opportunities)1980
Production CentersProduction Centers 3 and 4 open at medium capacityProduction Centers 5 open at high capacity
Distribution CentersDistribution Centers 6, 7, 8, 9, and 10 open at high capacity
Collection CentersCollection Centers 1, 2, 4, 8, and 9 open at low capacityCollection Center 6 open at medium capacity
Retreading CentersRetreading Center 2 open at low capacityRetreading Center 2 open at medium capacityRetreading Centers 1 and 3 open at high capacity
Recycling CentersRecycling Centers 1, 2, 3, and 5 open at low capacity
Suppliers being utilizedSuppliers 1, 2, 3, 4, and 5
θ ξ (scenario variability)190.680
Δ ξ (shortage)0
ϑ 𝕡 0
𝕡 (confidence level of fuzzy constraints)𝕡1 = 1/𝕡2 = 0.5/𝕡3 = 0.5/𝕡4 = 1/𝕡5 = 1/𝕡6 = 1
Table 10. Summary of optimum result (maximizing the social benefits).
Table 10. Summary of optimum result (maximizing the social benefits).
Value of Coefficients λ = 0 ,   ω = 0 ,   Ψ = 0
Objective 1 (profit, million dollar)998
Objective 2 (total GHG =kg)320,563
Objective 3 (total job opportunities)2452
Production CentersProduction Centers 3 and 4 open at high capacity
Distribution CentersDistribution Centers 6, 7, 8, 9, and 10 open at high capacity
Collection CentersCollection Centers 1, 2, 4, 6, 8, and 9 open at high capacity
Retreading CentersRetreading Centers 1, 2, 3, and 4 open at high capacity
Recycling CentersRecycling Centers 1, 2, 3, 4, and 5 open at high capacity
Suppliers being utilizedSuppliers 1, 3, 4, and 5
θ ξ (scenario variability)168.635
Δ ξ (shortage)0
ϑ 𝕡 0
𝕡 (confidence level of fuzzy constraints)𝕡1 = 1/𝕡2 = 1/𝕡3 = 1/𝕡4 = 0.5/𝕡5 = 0.5/𝕡6 = 0.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tehrani, M.; Gupta, S.M. Designing a Sustainable Green Closed-Loop Supply Chain under Uncertainty and Various Capacity Levels. Logistics 2021, 5, 20. https://0-doi-org.brum.beds.ac.uk/10.3390/logistics5020020

AMA Style

Tehrani M, Gupta SM. Designing a Sustainable Green Closed-Loop Supply Chain under Uncertainty and Various Capacity Levels. Logistics. 2021; 5(2):20. https://0-doi-org.brum.beds.ac.uk/10.3390/logistics5020020

Chicago/Turabian Style

Tehrani, Mohsen, and Surendra M. Gupta. 2021. "Designing a Sustainable Green Closed-Loop Supply Chain under Uncertainty and Various Capacity Levels" Logistics 5, no. 2: 20. https://0-doi-org.brum.beds.ac.uk/10.3390/logistics5020020

Article Metrics

Back to TopTop