Fuzzy multi-objective optimization model to design a sustainable closed-loop manufacturing system

View article
PeerJ Computer Science

Introduction

It is an eminent fact that major environmental concerns have arisen because of global megatrends such as the exponential growth of population and pollution due to the disposal of untreated and unavoidable manufacturing waste, carbon dioxide CO2 emissions, energy consumption, and depletion of natural resources. Manufacturing companies represent the backbone of contemporary industrialized civilization and have long been regarded as the fulcrum of the global economy. It has been indicated by Karuppiah, Sankaranarayanan & Lo (2024) that advancement in industrialization is adversely impacting the environment and human life on earth but it is yet uncertain how to evaluate their long-term viability. Sheng et al. (2023) identified that the manufacturing sector is a major contributor to huge energy consumption, waste generation, and greenhouse emissions, therefore researchers are required to propose/develop such kind of industrial model that may ensure enhanced resource efficiency, minimize waste, and reduce carbon emission. The waste generation and carbon emission rate in Fig. 1 (Kaza et al. (2018)) and Fig. 2 (IEA (2023)) indicate the scale of environmental degradation caused by these factors. Although new terms like environmentally conscious manufacturing, green manufacturing, and reverse manufacturing are thought to be connected to sustainability, there is a wide range of views and interpretations, due to which it has become difficult to find a definitive definition for sustainability in manufacturing firms (Sartal et al., 2020).

Projected waste generation (million tons per year) by region.

Figure 1: Projected waste generation (million tons per year) by region.

Carbon emission from industrial and transportation sector.

Figure 2: Carbon emission from industrial and transportation sector.

It is incumbent that industrial communities develop such modern enterprises which may ensure that enacting a sustainable system is an essential requirement (Anwar et al., 2023). Most definitions of sustainability include terms like green, ‘clean, keep, retain, stability, ecological balance, natural resources, and the environment. It can be difficult to define the term “sustainability,” usually it is understood as a way of nonlethal manufacturing products, and natural resources preservation, and frugal. In the meantime, there is still a disagreement over what constitutes and defines sustainability; promoting sustainable development has emerged as a significant worldwide goal. The phrase sustainable development (SD) is widely used to describe contemporary ideas of sustainability. SD has been regarded as a form of development that can gratify the needs of human beings devoid of jeopardizing the capability of future generations to gratify their desires. It is crucial to consider SD principles in every activity or choice made to achieve responsible and sustainable development. One of the most crucial topics for advancing the overall goal of SD is sustainable manufacturing. Designing for sustainability and sustainable development becomes the design of sustainable manufacturing firms. To accomplish sustainability, it is significant to appreciate the models and interrelated needs sustainably.

Over the last century, sustainability has been attracting both researchers and the industrial sector from almost all over the world (Bello-Pintado, Machuca & Danese, 2023). To start along a more sustainable road, several objectives and targets have been established. To encourage the concept of sustainability by confronting the three primary global concerns, the United Nations has defined seventeen sustainable development targets (SDGs) (see Table 1). The goals for sustainable development are interconnected, which means that achieving one requires addressing problems that are frequently linked with another. Just as the unadorned stressful situation on earth is caused by numerous interconnected global megatrends, so are the goals for sustainable development.

Table 1:
Seventeen sustainable development goals (UN, 2015).
1. No poverty 2. Zero hunger 3. Good health and well-being 4. Quality education
5. Gender equality 6. Clean water and sanitation 7. Affordable and clean energy 8. Decent work and economic growth
9. Industry, innovation and infrastructure 10. Reduced inequalities 11. Sustainable cities and communities 12. Responsible consumption and production
13. Climate action 14. Life below water 15. Life on land 16. Peace, justice and strong institutions
17. Partnership for the goals
DOI: 10.7717/peerj-cs.2591/table-1

One of the important goals of SDGs is SDG9 which focuses on industry, innovation, and infrastructure. Here, it is important to mention that an important element of SD is sustainable industrialization. In developing countries, industrialization has produced several environmental impacts and sustainable industrialization is being promoted at a rapid pace; as it has been indicated by Chen et al. (2024) that due to it there is an increase in carbon emission in many countries. Regarding this, action has been taken on the defined principles to encourage fuel (energy) saving and the reduction in the emissions of CO2 in many manufacturing industries. The expanded efforts (cost, energy, and emissions) of CO2 at all levels and elements, from the process, product, and the linked system through the complete product life cycle, must be made to make a manufacturing firm more sustainable. It is necessary to define a sustainable manufacturing system to minimize environmental effects. Utilizing mathematical optimization approaches, the issues with sustainable industrialization are resolved. When solving problems with sustainable manufacturing, many goal functions must be taken into account at once. For instance, in the sugar sector, boosting sugar output while decreasing waste is a problem. Because of this, no one solution can be referred to as the best and a variety of trade-off options must be taken into account. These issues serve as illustrations of the class of optimization problems known as multi-objective optimization problems.

To handle the range of optimization approaches used in real-world circumstances, multi-objective optimization is crucial. For addressing several objectives and offering solutions that are more in line with system requirements, there are numerous specialized strategies. The construction of the non-inferior (Pareto optimum) solution is a general objective in multi-objective optimization due to the multi-criteria of such issues; a non-inferior solution is one in which an improvement in one objective necessitates the deterioration of another. Typically, several trade-off solutions (referred to as the Pareto set) are discovered using the stated concept of optimality. Some data used for analysis and other purposes in the actual world are questionable.

However, decision-makers must take this unpredictability into account while designing their networks. The creation of a multi-objective optimization helps to address uncertainty resulting from the impreciseness of data. Fuzzy multi-objective optimization is a method that relies only on the formulation of fuzzy information in terms of membership functions to address the optimality of the fuzziness model utilizing currently available multi-optimization tools and methodologies (Paksoy, Pehlivan & Özceylan, 2012).

Supply chain operations now inevitably influence several company sectors and international trade. Customers play a fundamental part in these processes, and markets work to meet their needs by taking lead time minimization, product quality improvement, quick delivery, and their prospective expectations into account. The lack of precious resources, the state of the economy, social and environmental effects, and the rise in customer requests aimed at the finest products, lead policymakers towards deliberation and reverse logistics networks (Soon et al., 2022; Kazancoglu et al., 2022). Reverse logistics advantages in this area let organizations best the most of their businesses through material recycling, return, and reuse in multi-level systems that help improve the environment (Tavana et al., 2022).

The modern closed-loop manufacturing system will unavoidably have an impact on several industrial sectors as well as global trade. The primary goal of the research is to minimize the cost, total energy consumption, and overall CO2 emissions of a suggested closed-loop manufacturing system (CLMS) by choosing the overall quantity of existing services and significantly reducing the overall cost within multiple levels. Each place in this design seems to have a limited scope, and other places for each place are defined. The supply rates to every type of customer are stated for each period. Each warehouse has a certain capacity, and after each period, the inventory cost is taken into consideration for the items that are still supplied. Since the multi-objective issue has been taken into account, the goal of the research is to use optimization approaches to validate the model and produce the best outcomes. The closed-loop supply chain (CLSC) is thought to be a tactic for developing a sustainable supply network, as Roudbari, Fatemi Ghomi & Eicker (2024) stated that it can play two major roles; satisfying customers’ demands and effectively recapturing the value of the returned product. The following are the primary goals of the research in the manufacturing system:

  • 1)

    Constructing a closed-loop, sustainable supply network for the manufacturing phase in the industrial sector.

  • 2)

    Adopting a linear programming methodology to reduce the overall expenses of the manufacturing logistics chain.

  • 3)

    To supply materials and transportation of returned goods for subsequent use, take into consideration forward and reverse logistics.

The methodology initially evaluates the various providers under sustainable standards. The network of closed-loop supply chains is assessed to get the best possible amount from the potential providers chosen based on sustainable standards to satisfy consumer requests. Supplier (SU) then ships the quantity qSU.PLRaw through plant (PL) to retailer (RE), and then transport it to customer (CU). The product as waste material is collected by collection center (CC) by CU and transported to disassembly center (DC) to refurbishing center (RC) back on the road to PL using various modes of transportation. The next section of this study presents a comprehensive literature review along with the current study’s motivation, aim, and objectives. “Closed-Loop Manufacturing System Model Formulation” develops the thorough aspects of dilemma formulation, taking into consideration objective function and problem constraint. “Solution Methodology” presents the solution methodology followed by a case study, results, comparison of techniques, and conclusion.

Literature review

There have been several studies in supply chain systems, closed-loop supply chains, logistics system planning, and remanufacturing (Lo et al., 2021; Nunes, Causer & Ciolkosz, 2020; De & Giri, 2020). However, there are relatively few chronological writings on manufacturing supply chains. Initially, Lim et al. (2021) offers industrial sector approaches and created supply chain topics to draw scholars by providing an example of a specific supply chain network. Marić & Opazo-Basáez (2019) supports sustainable industrialization and ensures that providing an industrial supply chain for recycling products, results in more quality in the manufacturing industry. The planned measure was then given environmental and economic investigation in the same study. A detailed analysis of industrial practices and approaches for a design intended for sustainable actions is presented by Shekarian et al. (2022). They provided a new and enhanced classification of the work indicating more routes so it may be used to analyze how different firms reach sustainability.

The implementation of systems for closed-loop supply chain processes and reverse logistic (RL) has been studied in depth by a significantly as well as an increasing number of studies. RL is a mechanism in which products and other products are moved from their ultimate place to further destinations to be recycled, disposed of, or used in other operations (Fathollahi Fard & Hajaghaei-Keshteli, 2018; Cheraghalipour, Paydar & Hajiaghaei-Keshteli, 2017). Since RL promotes better growth and worth including both company owners and consumers, it has been adopted by the majority of emerging nations in various industries and enterprises (Hsu, Tan & Mohamad Zailani, 2016).

Keeping in view these advantages, previous research utilized reverse or closed-loop supply chain systems in their networks (Ritchie et al., 2000; Govindan, Soleimani & Kannan, 2015), and discussed closed-loop supply chain systems reversal logistics, too. As in research, they researched different aspects which include ecological factors, economic and legal aspects, and their positives from a supply chain and CLSC perspective, Özceylan, Paksoy & Bektaş (2014) suggesting a nonlinear model for constructing a closed-loop supply system in the accurate economic climate for balance separating and supply lines. The levels in the forward process include supplier, plant, retailer, and customer also the same number of levels are in the design of the reverse process which are the refurbishing center, waste center, disassembly center, and collecting center. To fulfill the need for product returns, certain writers, including Keyvanshokooh, Ryan & Kabir (2016), examined and improved the CLSC and presented models. The model used Benders’ approach with hybrid probabilistic and resilient programming to get the best solution to the issue. An increasing tendency towards networks of sustainability practices may be seen throughout the development of supply chain management. Rahemi et al. (2020), Li, Manier & Manier (2019), Balcombe et al. (2018) and Soleimani et al. (2017) included sustainability, optimization, and new distinct networks. They used multi-objective programming to construct a trustworthy and fuzzy CLSC that takes into account uncertainty. Suppliers, manufacturers, distributors, warehouses, return centers, and disposal facilities are all part of this supply chain. The maximization of overall revenue, taking sustainability principles, and enhancing customer response are examples of objective functions. To enhance the outcomes, a genetic algorithm (GA) was utilized. Devika, Jafarian & Nourbakhsh (2014) suggesting a linear model creates a closed-loop supply chain that is sustainable; to concurrently research the impacts of the environment, society, and economy in a specific context. The design contains six stages of collecting, recycled, reconstruction, refinement, waste, and supplementary consumers in the reverse direction, and four stages of supplies, manufacturing, distributing, and consumers in a forward manner. To cut carbon emissions, Talaei et al. (2016) created a closed-loop, environmentally friendly model of the supply chain. The suggested model tackles the position, in addition to optimal allocation problems, and has two objectives. The authors combined fuzzy sets and robust programming to deal with the uncertainty.

To solve the model, the epsilon constraint approach is employed. By creating a CLSC network and taking into account both optimizing profitability and social capacities by recycling items Pishvaee & Torabi (2010) and Panda, Modak & Cárdenas-Barrón (2017) attempted to maximize their suggested model by offering a bi-objective design framework about a closed-loop supply chain in an ambiguous setting. Although the secondary purpose seeks in minimizing the time needed to receive goods and services, the first goal seeks to decrease expenses. The model contains four layers of collecting, rebuilding, recycled, and third consumers in both forward and reverse directions, three phases of manufacturing, distribution, and consumers are also included. Soleimani & Kannan (2015) constructed a CLSC system on large systems. The suggested system is comprehensive, multi-level, multi-period, and multi-product. The model is solved using the GA, particle swarm optimization (PSO) technique, hybrids genetic algorithm particle swarm optimization (GAPSO), and the exact approach (CPLEX). Kaya & Urek (2016) created a CLSC network to place stocks and set prices in the best possible places. To find services and inventory, a nonlinear mixed integer programming approach is suggested in this study. The goal of the problem was to increase the suggested supply chain network’s overall profit.

While exploring the design of reversible logistics networks, taking into account multiple layers including warehouse, transit, logistic, manufacturer, distributor, and client, it was seen how each stage interacted with the supply chain (Fulconis & Philipp, 2018). Sadeghi, Mina & Bahrami (2020) used mixed integer linear programming (MILP) to create a reverse chain for industrial logistics. The levels of production, distribution, separation, and destruction were among those that were taken into account. To further clarify the problem, they developed a multi-period, multi-product model. Reverse supply networks’ environmental properties and their supply chain were analyzed and categorized by Islam et al. (2021). Govindan & Soleimani (2017) presented a review of reverse logistics and other earlier research in this field. Kazemi, Modak & Govindan (2019) developed evaluations of the work on reverse logistics and CLSC as well as a detailed analysis of content. For armed hardware including replacement parts, a closed-loop system was developed (Wang & Shi, 2019). Since it is essential to get to the military supplies, the target function of the minimal durations is used in this research’s design. To handle the suggested system, the Lion meta-heuristic method is employed. Torabi et al. (2016) created a reverse supply chain while taking the disrupting risk and operating limits into account. The suggested model is then solved using based on fuzzy logic programming, and even the goal function of minimizing cost is then optimized. Gholizadeh et al. (2022) added the creation of a renewable system to the list of priorities while utilizing reverse logistics to address both environmental and economic concerns. In this study, reproduction and recycling levels are used to reduce the manufacturing cost’s objectives functions. Zhang et al. (2021) proposes a problem with logistics management and choice in the production chain with an emphasis on the criteria for choosing the best RL providers. The best RL logistics provider for remanufacturing is suggested in this scheme, and the strategy is resolved using a multi-criteria decision-making (MCDM) approach that includes an analytical hierarchy process (AHP) also a technique for order of preference by similarity to the ideal solution (TOPSIS). Pourjavad & Mayorga (2019), Pourmehdi, Paydar & Asadi-Gangraj (2020), Gholizadeh & Fazlollahtabar (2020), Pahlevan, Hosseini & Goli (2021), Akbari-Kasgari et al. (2022), Amirian, Khalili & Mehrabian (2022) are some investigations on CLSC or RL in industrial manufacturing systems. Even if some earlier research created industrial manufacturing systems, there have been relatively few investigations on CLSC or RL for such items. Table 2 contains some studies related to sustainable industrialization.

Table 2:
Literature reivew.
Study Model Objective functions Solution approach
Akbari-Kasgari et al. (2022) Mixed integer linear programming Maximize profit, minimize water consumption and air pollutants, and maximize social desirability ε—constraint method, weighted sum method
Goodarzian et al. (2023) Mixed integer linear programming Minimize cost, minimize environmental impacts and Optimize social factors ε—constraint method
Lotfi et al. (2024) A two-stage robust stochastic multiobjective programming mode Minimize cost, CO2 emissions, energy consumption and maximising employment Lp-metric method
Khalili-Fard et al. (2024) A novel two-stage stochastic model Minimize cost, maximize the remaining fuel and maximize total job opportunities Multi-objective particle swarm optimization
Gholipour et al. (2024) Multi-objective optimization method Minimize cost, minimize risk and maximize accountability A meta-heuristic algorithm
Momeni, Jain & Bagheri (2024) Multi-objective mixed-integer programming Minimize cost, minimize environmental impact and social impact ε—constraint method
Present study Multi-objective mixed-integer programming Minimize total investment cost, minimize energy consumption, and minimize carbon emission Four-valued refined neutrosophic optimization technique
DOI: 10.7717/peerj-cs.2591/table-2

The previous studies on CLSC and RL have mostly focused on developing new systems for selling rebuilt goods, recycling product returns, and converting waste into renewable energy. Although these efforts have expanded global attention, more comprehensive approaches toward closed and reverse supply chains, beyond just economic objectives, are necessary. In addition to addressing climate change, developing a holistic CLSC system presents wide-ranging opportunities for investors to maximize the value of products across their entire life cycle. Over the past two decades, significant developments in CLSC systems have simplified policymaking for addressing complex issues such as re-manufacturing used goods, meeting customer demand, and identifying marketing strategies for re-manufactured products. However, much of the existing research has focused on the forward supply chain and its role in meeting client demands, resulting in higher revenues and lower expenses for suggested systems. While such insights are valuable, they overlook critical aspects of closed-loop and backward logistics systems that are key to achieving a circular economy.

Despite the growing research on supply chain management, there is a noticeable gap in research addressing the full spectrum of CLSC systems. Specifically, the integration of circular economy principles, the development of multi-objective optimization models that consider both economic and environmental sustainability, and the need for stakeholder collaboration remain under-explored. Additionally, while most research has concentrated on the manufacturing industry, other sectors such as electronics, textiles, and consumer goods lack detailed studies on implementing efficient backward logistics systems. Additionally, the advent of digital transformation, AI, and IoT technologies presents new challenges and opportunities for enhancing the efficiency of CLSC networks. However, the role of these technologies in building resilient CLSC systems remains insufficiently explored in the current literature. Addressing these gaps could provide a more comprehensive framework for the implementation of CLSC, promoting sustainability and long-term profitability in various industries. Therefore, the following represent the key goals and suggestions to address the study gaps throughout the research:

  • 1)

    Literature indicated that there are plenty of research studies on aforementioned topic in industrial supply chains (ISC), therefore, this study begins by examining the existing research on supply chain system architecture for the industrial sector.

  • 2)

    Here on basis of its creative system, the present investigation has created revolutionary modeling. Although the majority of earlier studies recommended forward supply chain systems.

  • 3)

    The suggested framework not only includes a tire backward route to preparing the gathered produced goods again for the secondary market, but also minimizes the total costs, overall energy, and CO2 emissions in both forward and backward reverse directions.

  • 4)

    Because of its distinctive system, the research presents a new mixed integer linear programming (MILP) approach from a mathematical standpoint.

  • 5)

    The objective will be achieved using a four-valued refined neutrosophic optimization approach.

Closed-loop manufacturing system model formulation

The imperative for returned products from the manufacturing process is one of the emerging concerns in CLMS design for the industrial system. These materials can not only be obtained as raw materials in a range of industries and enterprises but they can also be recycled and used again. In this study, a CLMS (see Fig. 3) is developed to minimize total costs, energy use, and CO2 emissions in the manufacturing system’s supply chain while taking into account returned goods on a reverse flow. Figure 3 encapsulates a schematic design of a CLMS. It moves ahead with production from the facility directed through the retailer to customers before going backward with the collection, disassembly, and refurbishment centers. During these phases, the recycling process is carried out to change the commodities into new ones. Our suggested model incorporates tactical decisions, such as production, inventory planning, capacity management, and transportation system, as well as strategic considerations, such as the best choice of facility sites made up of disassembly centers, recycling facilities, and collection centers. The locations of the factory and customers are taken for granted.

Structure of the proposed closed-loop manufacturing system.

Figure 3: Structure of the proposed closed-loop manufacturing system.

The following are the fundamental presumptions of the model:

  • 1)

    The prospective sites of each center (Plant, Retailer, Collection center, Disassembly center, waste center, and Refurbishing center) are known and each has a limited capacity.

  • 2)

    The predicted customer demand rate for each time.

  • 3)

    There cannot be a shortage.

  • 4)

    The inventory cost for the remaining products is calculated at the end of each period.

  • 5)

    All costs are taken into account after each session.

Decision variables and parameters

Notations

SU

supplier

PL

plant

RE

retailer

CU

customer

CC

collection center

DC

disassembly center

WC

waste center

RC

refurbishing center

kPLj

number of processes involved in the plant; j{1,2,,kPL}

Parameters

CCLFixed

Fixed cost ($) of closed loop system

CSU.PLRaw

Raw material cost ($)

CSURaw

Unit raw material cost ($) in SU

CPL.REProduct

Manufactured product cost ($)

CPLProduct

Unit manufactured product cost ($)

CPL.REI

Inventory cost ($) from a PL to a RE

CREI

Unit inventory cost ($)in RE

CCLT

Transportation cost ($) of closed loop system

CAT.B

Unit transportation cost ($) (per mile) of material from A to B

DA.B

Distance in miles from A to B

GA

Quantity (in kg) of the material produced from A (per month)

PLj

Machine manufacturing rate (%) involved in process j

μPLj

Machine efficiency (%) involved in process j

ρPLjaircompressor

Compressed air’s capacity ( m3/h) of compressor in process j

NAjmachine

Installed power (Kw) of machine in A involved in process j

NAjcompressor

Installed power (Kw) of compressor in A involved in process j

NAjairconditioner

Installed power (Kw) of air conditioner in A involved in process j

NAjilluminations

Installed power (Kw) of illuminations in A involved in process j

emPL

Quantity of CO2 (in kg) released from PL

emT

Quantity of CO2 (in kg) released from vehicles during the transportation

emA.BT.Raw

Quantity of CO2 (in kg) emitted during transportation of raw materials from A to B

emA.BT.Product

Quantity of CO2 (in kg) emitted during transportation of manufactured product from A to B

emPLjmachine

Quantity of CO2 (in kg) released from machines involved in process j

emPLjaircompressor

Quantity of CO2 (in kg) released from air compressor involved in process j

emPLjairconditioner

Quantity of CO2 (in kg) released from air conditioner involved in process j

emPLjilluminations

Quantity of CO2 (in kg) released from illuminations involved in process j

αPLj

CO2 emission ( kg/kWh) depending upon the energy resources that are used in a PL

αT

CO2 emission ( kg/mile) depending upon the transportation of material

βPLjmachine

The total waste ratio (in %) for machines involved in process j

γPLjairconditioner

Covering rate (%) for each air conditioner used in process j

γPLjaircompressor

Covering rate (%) for each air compressor used in process j

γPLjilluminations

Covering rate (%) for each illumination used in process j

γAairconditioner

Covering rate for each air conditioner used at A (other than PL)

γAilluminations

Covering rate (%) for each illumination used at A (other than PL)

DA

Demand (per month) of material (in Kg) at A

CaA

Capacity (in Kg) of the A

V

Capacity in kg per vehicle

R1Rate

Return rate (%) from customer to CC

R2Rate

Return rate (%) from CC to DC

R3Rate

Return rate (%) from CC to PL

R4Rate

Return rate (%) from DC to RC

R5Rate

Return rate (%) from DC to WC

Decision variables

qA.BProduct

Quantity of produced material (kg) transported from A to B

qA.BRaw

Quantity of raw material (kg) transported from A to B

qPLjRaw

Quantity of raw material (kg) used in process j in PL

qPLj+1Raw

Quantity of material (kg) shipped from process j in PL

nPLjmachine

Number of machinery (unit) in the plant involved in process j

nAairconditioner

Number of air conditioners (unit) installed at A

nAilluminations

Number of illumination (unit) installed at A

Note: In the description, instead of writing the names of each center, the change in one center to another has been written as “from A to B”.

Objective functions

Total cost invested f1

The total investment cost in a circular system is considered by combining the all fixed costs of a closed loop system which include costs of land, building, equipment, services, and wages; cost of raw material; cost of transportation of raw materials; costs of manufacturing and inventory. Here the total cost f1 to minimize is given in Eq. (1).

Minimizef1=CCLFixed+CSU.PLRaw+CPL.REProduct+CIPL.RE+CTCLwhere, the fixed cost CCLFixed of the circular supply chain model is obtained by the combination of costs of the plant CCLPL, retailer CCLRE, collection center CCLCC, disassembly center CCLDC, waste center CCLWC, and refurbishing center CCLRC. These costs are given in Eqs. (2)(6).

CCLPL=CPLLand+CPLBuilding+CPLEquipment+CPLservices+CPLWages

CCLRE=CRELand+CREBuilding+CREEquipment+CRservices+CREWages

CCLCC=CCCLand+CCCBuilding+CCCEquipment+CCCservices+CCCWages

CCLDC=CDCLand+CDCBuilding+CDCEquipment+CDCservices+CDCWages

CCLWC=CWCLand+CWCBuilding+CWCEquipment+CWCservices+CWCWages

The cost of one unit of raw material CSU.PLRaw from supplier to plant and the cost of one unit of product from plant to the retailer is defined in Eqs. (7) and (8).

CSU.PLRaw=CSURawqSU.PLRaw

CPL.REProduct=CPLProductqPL.REProduct

Inventory cost CCL.REI at the retailer is calculated by using Eq. (9).

CPL.REI=CREIqPL.REProduct

The transportation cost of material CCLT in the circular supply chain model is obtained by the combination of the transportation costs of material transported between various segments of CLSC as in Eq. (10).

CCLT=CCLT.PL+CCLT.RE+CCLT.CU+CCLT.CC+CCLT.DC+CCLT.WC+CCLT.RC

The transportation costs involved in Eq. (10) are derived from Eqs. (11)(17).

CCLT.PL=CST.PLqSU.PLRawVDSU.PL+CCCT.PLqCC.PLRawVDCC.PL+CCLT.PL=CRCT.PLqRC.PLRawVDRC.PL

CCLT.RE=CPlantT.REqPL.REProductVDPL.RE

CCLT.CU=CRET.CUqRE.CUProductVDRE.CU

CCLT.CC=CCUT.CCqCU.CCRawVDCU.CC

CCLT.DC=CCCT.DCqCC.DCRawVDCC.DC

CCLT.WC=CPlantT.WCqPL.WCProduct.WCVDPL.WC+CDCT.WCqDC.WCRawVDDC.WC

CCLT.RC=CDCT.RCqDC.RCRawVDDC.RC

Total energy consumption f2

The second objective is to minimize the total energy consumption during the production process. The objective is mathematically expressed in Eq. (18).

Minimizef2=EPL+ERE+ECC+EDC+EWC+ERC

In the above equation, the total energy consumed by the plant EPL is computed in Eq. (19).

EPL=j=1kPLEPLjmachine+j=1kPLEPLjaircopmressor+j=1kPLEPLjairconditioner+j=1kPLEPLjilluminations=j=1kPL(qPLjRawRPLjμPLjNPLjmachinenPLjmachine)+j=1kPL(qPLjRawRPLjμPLjNPLjaircompressorρPLjaircompressorνPLjaircompressornPLjmachine)+j=1kPL(NPLjairconditionernPLjairconditionerqPLj+1RawGPL)+j=1kPL(NPLjilluminationsnPLjilluminationsqPLj+1RawGPL)

Equation (20) provides the energy consumed at the retailer’s end.

ERE=NREairconditionernREairconditionerqPL.REproductGRE+NREilluminationsnREilluminationsqPL.REProductGRE

Equation (21) provides the energy consumed at the collection center.

ECC=NCCairconditionernCCairconditionerqCU.CCRawGCC+NCCilluminationsnCCilluminationsqCU.CCRawGCC

Equation (22) provides the energy consumed at the disassembly center.

EDC=NDCairconditionernDCairconditionerqCC.DCRawGDC+NDCilluminationsnDCilluminationsqCC.DCRawGDC

Equation (23) provides the energy consumed at the waste center.

EWC=NWCairconditionernWCairconditionerqDC.WCRawGWC+NWCilluminationsnWCilluminationsqDC.WCRawGWC

Equation (24) provides the energy consumed at the refurbishing center.

ERC=NRCairconditionernRCairconditionerqDC.RCRawGRC+NRCilluminationsnRCilluminationsqDC.RCRawGRC

Total carbon emission f3

The third objective is to minimize carbon emission. The total emission is comprised of two components, that are, emission from plants and emission during transportation expressed in Eq. (25):

Minimizef3=emPL+emT

The quantity of CO2 which is emitted from the plant during production can be computed by using Eq. (26).

emPL=αPLjEPLjmachineqPLjRaw+αPLjEPLjaircopmressor+αPLjEPLjairconditioner+αPLjEPLjilluminations

The quantity of CO emitted from vehicles during the transportation of material from supplier to plant and shipment of the product to the retailer, customer, collection center, disassembly center waste center, and refurbishing center can be calculated by using Eq. (27).

emT=emSU.PLT.Raw+emPL.RET.Product+emPL.WCT.Product+emRE.CUT.Product+emCU.CCT.Raw+emCC.DCT.Raw+emCC.PLT.Raw+emDC.WCT.Raw+emDC.RCT.Raw+emRC.PLT.Raw=αTqSU.PLRawVDSU.PL+αTqPL.REProductVDPL.RE+αTqPL.WCProductVDPL.WC+αTqRE.CUProductVDRE.CU+αTqCU.CCRawVDCU.CC+αTqCC.DCRawVDCC.DC+αTqCC.PLRawVDCC.PL+αTqDC.WCRawVDDC.WC+αTqDC.RCRawVDDC.RC+αTqRC.PLRawVDRC.PL

Constraints

Supply constraints

Equation (28) ensures that the quantity of material shipped to the plant, retailer, customer, collection center and refurbishing center, or disassembly center cannot be greater than their capacity.

PLqSU.PLRawCaSUREqPL.REProduct+WCqPL.WCProductCaPLCCqRE.CCRawCaREPLqCC.PLRaw+DCqPL.DCRawCaDCRCqDC.RCRaw+WCqDC.WCRawCaDCPLqRC.PLRawCaRC

Demand constraints

Equation (29) ensure that the demands of the plant, retailer, and customer are fulfilled, respectively.

SUqSU.PLRawDPLPLqPL.REProductDREREqRE.CUProductDCU

Balance equation

Equation (31) represent the balanced constrained equation.

SUqSU.PLRaw=REqPL.REProduct+WCqPL.WCProductPLqPL.REProduct=CUqRE.CUProductREqRE.CUProduct=CCqCU.CCProductCUqCU.CCRaw=PLqCC.PLRaw+DCqCC.DCRawCCqCC.DCRaw=WCqDC.WCRaw+RCqDC.RCRawDCqDC.RCRaw=PLqRC.PLRaw

Informational constraints

(1βPLjmachineqPLjRaw)qPLj+1RawγPLjairconditionernPLjairconditionernPLjmachineγPLjaircompressornPLjaircompressornPLjmachinenPLjilluminationsγPLjilluminationsnPLjmachineRPLjmachinenPLjmachineqPLj+1RawCaREnREairconditionerγREairconditionerqPL.REProductCaREnREilluminationsγREilluminationsqPL.REProductCaCCnCCairconditionerγCCairconditionerqCU.CCRawCaCCnCCilluminationsγCCilluminationsqCU.CCRawCaDCnDCairconditionerγDCairconditionerqCC.DCRawCaDCnDCilluminationsγDCilluminationsqCC.DCRawCaRCnRCairconditionerγRCairconditionerqDC.RCRawCaRCnRCilluminationsγRCilluminationsqDC.RCRaw

Transshipment constraints

CUqCU.CCRawR1RateREqRE.CUProductCUqCU.CCRawR2RateCCqCC.DCRawCUqCU.CCRawR3RateCCqCC.PLRawCCqCC.DCRawR4RateDCqDC.RCRawCCqCC.DCRawR5RateDCqDC.WCRaw

Solution methodology

The developed multi-objective linear programming model is solved by using a four-valued refined neutrosophic optimization technique (FVRNOT) proposed by Freen et al. (2020). A four-valued refined neutrosophic set is one of several types of neutrosophic sets in which indeterminacy is divided into δ= uncertainty and η= contradiction, where η=μν. The values of μ,δ,ηandν fall into the range [0,1], and 0μ+δ+η+ν4. Mathematical expression of the set in given in Eq. (33).

ZRN={(s,μZRN(s),δZRN(s),ηZRN(s),νZRN(s)):sS}.

The objective functions are transformed into neutrosophic fuzzy constraints when optimization is executed using FVRNOT. The uncertainty is effectively handled by this four-valued refined neutrosophic model that has been established and enhanced. The multi-objective four-valued refined neutrosophic linear programming model that is being described can handle uncertain data, which avoids unsustainable modeling. The computational algorithm of the stated approach includes the following steps:

1. One objective function from a set of N objectives should be used to solve the first objective function. Under the stated constraints, the values of the decision variables and objective functions will be determined.

2. Now, through using decision variables from step 1, calculate the values of the undetermined objective, i.e., (N1).

3. Steps 1 and 2 must be repeated for the rest (N1) objective functions to get the pay-off matrix 34.

[f1(s1)f2(s1)fn(s1)f1(s2)f2(s2)fn(s2)f1(sn)f2(sn)fn(sn)]

4. For every objective function fp determine the lower bound Lp and the upper bound Up such that

Lp=maxr=1n{fp(sr)}andUp=minr=1n{fp(sr)}.

5. The lower and upper bounds for truth membership function μ, falsity membership function ν, uncertainty membership function δ and contradiction membership function η of the objective function fp are computed in Eqs. (35)(38).

Lpμ=Lp;Upμ=Up

Lpν=Lpμ+t1(UpμLpμ);Upν=Upμ

Lpδ=Lpμ+t2(UpμLpμ);Upδ=Upμ

Lpη=LpμLpν;Upη=UpμUpν+t3(UpμUpνLpμLpν)where t1,t2,t3(0,1).

6. Define the truth, uncertain, false, and contradiction membership functions as in Eqs. (39)(42)

μp(fp(s))={1iffp(s)Lpμ,Upμfp(s)UpμLpμifLpμfp(s)Upμ,0iffp(s)Upμ,

δp(fp(s))={1iffp(s)Lpδ,Upδfp(s)UpδLpδifLpδfp(s)Upδ,0iffp(s)Upδ,

νp(fp(s))={0iffp(s)Lpν,fp(s)LpνUpνLpνifLpνfp(s)Upν,1iffp(s)Upν,

ηp(fp(s))={1iffp(s)Lpη,Upηfp(s)UpηLpηifLpηfp(s)Upη,0iffp(s)Upη,

7. Therefore, after finding the membership and non-membership functions, an equivalent crisp optimization problem is formed in Equation (43).

max(ϕσ+ψ+θ),such that

μp(fp(s))ϕ,νp(fp(s))σ,δp(fp(s))ψ,ηp(fp(s))θwith

0ϕ+σ+ψ+θ4andϕσ,ψ,θand subject to the original constraints.

Assessment: a case study

A real case study is employed to determine the effectiveness and validity of the previously mentioned designed fuzzy multi-objective optimization model. For this purpose data were collected from already conducted research studies and minor changes were made to make that data suitable for practice. The plant uses several distinct procedures to make its products, thus different phases require varied numbers of machines, air conditioners, air compressors, and illumination. Additionally, each center in the suggested model is made up of such equipment, and each piece of equipment has a varied mass, energy consumption, and CO2 emission quantity.

The data for a manufacturing system that produces, stores, and transports the product is shown in Table 3. The supply chain in this instance is facilitated by two different renewable energy sources, particularly oil, which is used as direct sources to generate thermal energy and electricity, respectively; illustrating how to create a sustainable manufacturing system. The solution methodology defined in the previous section is carried out in the following manner.

Table 3:
Data table.
Amount of C02 emission factor per kWh using oil (per mile)
αPLj 0.04 kg/kWh αT 0.41 kg/mile
Cost
CPLFixed 6,000,000 $ CSU.PLRaw 2 $
CPL.REProduct 3 $ CPL.REI 2 $
CCLT 2 $
Distance in miles
DSU.PL 50 DPL.RE 30
DPL.WC 05 DRE.CU 15
DCU.CC 10 DCC.PL 15
DCC.DC 20 DDC.WC 08
DDC.RC 17 DRC.PL 20
Installed power in kW
NPLjmachine 180 NPLjaircompressor 150
NPLjairconditioner 2 NPLjilluminations 0.5
NREairconditioner 0.5 NREilluminations 0.1
NCCairconditioner 2 NCCilluminations 0.3
NDCairconditioner 1.5 NDCilluminations 0.2
NRCairconditioner 1 NRCilluminations 0.1
Capacity
V (per vehicle) 20,000 kg ρPLjaircompressor 666 m3/h
CaPL 16,000 kg/month CaSU 30,000 kg/month
CaRE 14,000 kg/month CaCC 5,000 kg/month
CaDC 3,000 kg/month CaRC 2,000 kg/month
Covering rate (units) used in each center
γPLjairconditioner 1.5 γPLjaircompressor 03
γPLjilluminations 15 γREairconditioner 1.5
γREilluminations 10 γCCairconditioner 1.5
γCCilluminations 10 γDCairconditioner 1.6
γDCilluminations 9 γRCairconditioner 1.4
γRCilluminations 10
Parameters involved in manufacturing
βPLjmachine 0.001 PLj 1840
μPLj 0.6 R1Rate 0.07
R2Rate 0.02 R3Rate 0.05
R4Rate 0.10 R5Rate 0.20
Quantity of material/product produced and the demand for each center (kg/month)
GPL 751,569 GRE 632,650
GCC 53,164 GDC 514,500
GWC 50,198 GRC 495,045
DPL 15,000 DRE 12,000
DOI: 10.7717/peerj-cs.2591/table-3

Step 1–3: The three objective functions are simultaneously solved, and the values of decision variables are substituted into the objective functions that provide the pay-off matrix presented in Eq. (46).

[6200100148060006160406200100137462435062001001374624350]

Step 4–5: The lower and upper bounds for truth membership function μ, falsity membership function ν, uncertainty membership function δ and contradiction membership function η of the objective function fp are computed by using Eqs. (35)(38).

For total cost investment f1

U1ν=U1μ=6200100,L1ν=6200100+0.30=6200100,U1δ=U1μ=6200100,L1δ=6200100+0.30=6200100,L1η=62001006200100=6200100,U1η=6200100.

For total energy consumption f2

L2μ=13746,U2μ=14806000,U2ν=U2μ=14806000,L2ν=4451400,U2δ=5930600,L2δ=13746,L2η=14806000,U2η=8889000.

For total carbon emission f3

L3μ=24350,U3μ=616040,U3ν=616040,L3ν=201900,U3δ=261020,L3δ=24350,L3η=616040,U3η=379360.

Step 6: Truth, uncertain, falsity, and contradiction functions are defined by using the above bounds and Eqs. (39)(42).

Step 7: Single objective optimization problem is defined as in Eq. (43) subject to constraints in Eqs. (28)(32), (44) and (45).

Results

The specified multi-objective optimization model for a sustainable manufacturing system was solved and subsequently optimized using FVRNOT. Solutions appear to have a maximum range for the number of machines nPLjmachine=7, number of air conditioner units nPLjairconditioner=5, number of illuminations nPLjilluminations=105 and the number of air compressor nPLjairconditioner=3 that are involved in the plant for manufacturing products; where j = 1 is fixed. In the given system for the retailer center, the number of air conditioner units is found nREairconditioner=196 and the number of illumination is nREilluminations=1,305. In the collection center nCCairconditioner=13 and nCCilluminations=86. The result also depicts the ideal supply chain strategy for the raw material, that is, qSU.PLRaw=12,000kg/month and the quantity of manufactured product transported from the plant to the retailer center qPL.REProduct=12,000kg/month. This amount fulfilled the consumer demands, which is qRE.CUProduct=12,000kg/month. The amount of material that the customer has supplied to the collection center is evaluated as qCU.CCRaw=840kg/month. By using the lower and upper bound of each function the aspiration level is determined. Considering the objective function of CO2 emission whose optimal value is 30,646 kg, the aspiration level can be calculated as:

Aspirationlevel=U3μf3(s)U3μL3μ×100=6160403064661604024350×100=98.9%

Comparably, we will determine the aspiration level of the first two functions. The ultimate satisfaction level is attained in the minimization scenario when the resulting value approaches the lower bounds. When the optimized value is identical to the lower bounds, the level of satisfaction is 100%.

The result indicates that the number of machines and equipment should be carefully chosen to maximize production efficiency while minimizing energy consumption and operational costs. These quantities are indicative of a system that balances production capacity with sustainability goals, avoiding overuse of resources and potential environmental impact. The high number of air conditioners and illuminations suggests a focus on maintaining product quality and visibility. However, this also implies substantial energy consumption. The design should consider energy-efficient technologies and practices to reduce the environmental footprint while meeting operational needs. This configuration supports effective raw material handling while managing energy use. The focus on minimizing energy consumption in the collection center aligns with sustainability objectives and reduces overall operational costs. The efficient flow of materials and products ensures that consumer demand is met without excess or shortage, contributing to a balanced and responsive supply chain. This strategy minimizes waste and enhances the overall efficiency of the closed-loop system. The high aspiration level reflects the system’s effectiveness in minimizing CO emissions, which is crucial for sustainability. This result underscores the importance of integrating emission reduction strategies into the design and operation of manufacturing systems.

Comparison

The under-consideration optimization problem is solved by using the intuitionistic fuzzy optimization technique (IFOT) and neutrosophic optimization technique (NOT) method.

An intuitionistic fuzzy set is defined as ZIFS={(s,μZIFS(s),νZIFS(s)):sS}. The intuitionistic fuzzy optimization technique (IFOT) is based on pay-off matrix, Eqs. (35), (36), (39) and (41). Where max(ϕσ) is fuzzy objective and μp(fp(s))ϕ,νp(fp(s))σ,0ϕ+σ1 and ϕσ are fuzzy constraints for IFOT.

Similarly, a neutrosophic set is defined as ZNS={(s,μZNS(s),δZNS(s),νZNS(s)):sS}. The neutrosophic optimization technique (NOT) is based on pay-off matrix, Eqs. (35)(37), (39), (40) and (41). Where max(ϕσ+ψ) is neutrospohic objective function and μp(fp(s))ϕ,νp(fp(s))σ,δp(fp(s))ψ, 0ϕ+σ+ψ3 and ϕσ,ψ are neutrosophic constraints. The values of the objective functions with the three techniques are presented in Table 4. It is clear that the FVRNOT method produces the best overall results compared to the other two methods because it requires less overall investment in terms of costs, energy use by the machines, and overall CO2 emissions from all sources.

Table 4:
Comparison of techniques.
Techniques Cost Energy consumption CO2 emission
IFOT 6,200,100 $ 2,554,000 kWh 125,960 kg
NOT 6,200,100 $ 2,554,500 kWh 125,980 kg
FVRNOT 6,200,100 $ 171,140 kWh 30,646 kg
DOI: 10.7717/peerj-cs.2591/table-4

Sensitivity analysis

The parameters such as inventory, production, and capacity of the multi-objective optimization model can directly be influenced and controlled by either operations or overall supply chain management. In addition, parameters like market demand and fixed return rate are mostly governed by market factors with significant impacts on the solution approach in the model. In this research, the parameters dependent on the market are selected to perform sensitivity analysis to have an insight into their influence on the objective functions. A custom experimental design is constructed based on the chosen factors and levels to analyze, and to run the model for evaluation of all scenarios, and the respective responses are illustrated in Table 5. The sensitivity analysis is performed by keeping one parameter constant and assessing the impact of one factor at all levels. For each level from 10% to +10%, this approach is repeated systematically. Further experiments repeatedly revealed similar trends of increasing or decreasing function values, resulting in five levels being added for explanatory purposes.

Table 5:
Experimental data for sensitivity analysis.
Levels Demand Return rate Total cost Energy consumption CO2 Emission
−10% −10% −10% 6,180,090 1.55E+05 2.76E+04
−10% −5% 6,180,090 1.55E+05 2.76E+04
−10% 0% 6,180,090 1.55E+05 2.76E+04
−10% 5% 6,180,090 1.55E+05 2.76E+04
−10% 10% 6,180,090 1.55E+05 2.76E+04
−10% −10% 6,180,090 1.55E+05 2.76E+04
−5% −10% 6,190,095 1.63E+05 2.91E+04
0% −10% 6,200,100 1.71E+05 3.06E+04
5% −10% 6,210,105 1.79E+05 3.22E+04
10% −10% 6.22E+06 1.89E+05 3.37E+04
−5% −5% −10% 6,190,095 1.63E+05 2.91E+04
−5% −5% 6,190,095 1.63E+05 2.91E+04
−5% 0% 6,190,095 1.63E+05 2.91E+04
−5% 5% 6,190,095 1.63E+05 2.91E+04
−5% 10% 6,190,095 1.63E+05 2.91E+04
−10% −5% 6,180,090 1.55E+05 2.76E+04
−5% −5% 6,190,095 1.63E+05 2.91E+04
0% −5% 6,200,100 1.71E+05 3.06E+04
5% −5% 6,210,105 1.79E+05 3.22E+04
10% −5% 6,220,110 1.89E+05 3.37E+04
0% 0% −10% 6,200,100 1.71E+05 3.06E+04
0% −5% 6,200,100 1.71E+05 3.06E+04
0% 0% 6,220,110 1.71E+05 3.06E+04
0% 5% 6,200,100 1.71E+05 3.06E+04
0% 10% 6200,100 1.71E+05 3.06E+04
−10% 0% 6,180,090 1.55E+05 2.76E+04
−5% 0% 6,190,095 1.63E+05 2.91E+04
0% 0% 6,200,100 1.71E+05 3.06E+04
5% 0% 6,210,105 1.79E+05 3.22E+04
10% 0% 6,220,110 1.89E+05 3.37E+04
5% 5% −10% 6,210,105 1.79E+05 3.22E+04
5% −5% 6,210,105 1.79E+05 3.22E+04
5% 0% 6,210,105 1.79E+05 3.22E+04
5% 5% 6,210,105 1.79E+05 3.22E+04
5% 10% 6,210,105 1.79E+05 3.22E+04
−10% 5% 6,180,090 1.55E+05 2.76E+04
−5% 5% 6,190,095 1.63E+05 2.91E+04
0% 5% 6,200,100 1.71E+05 3.06E+04
5% 5% 6,210,105 1.79E+05 3.22E+04
10% 5% 6,220,110 1.89E+05 3.37E+04
10% 10% −10% 6,220,110 1.89E+05 3.37E+04
10% −5% 6,220,110 1.89E+05 3.37E+04
10% 0% 6,220,110 1.89E+05 3.37E+04
10% 5% 6,220,110 1.89E+05 3.37E+04
10% 10% 6,220,110 1.89E+05 3.37E+04
−10% 10% 6,180,090 1.55E+05 2.76E+04
−5% 10% 6,190,095 1.63E+05 2.91E+04
0% 10% 6,200,100 1.71E+05 3.06E+04
5% 10% 6,210,105 1.79E+05 3.22E+04
10% 10% 6,220,110 1.89E+05 3.37E+04
DOI: 10.7717/peerj-cs.2591/table-5

In Figs. 4 and 5, a 5% increase and 10% decrease in demand respectively, with the percentage change in demand used as a changing parameter the objective function’s values are calculated at all levels. Due to increased demand, the additional use of machinery, transportation, and refurbishing is increasing energy consumption and CO2 emission with a slight change in cost. Reducing demand leads to decreases in both energy consumption and CO2 emissions. However, these changes in demand do not substantially affect cost, which remains nearly constant, highlighting a potential challenge for cost-sensitive strategies that aim to influence demand fluctuations. Similarly, the same trend at levels such as 10%, −5%, and +10% is clearly shown in Table 5.

Effect of 5% increase in demand on objective functions values.

Figure 4: Effect of 5% increase in demand on objective functions values.

Effect of 10% decrease in demand on objective functions values.

Figure 5: Effect of 10% decrease in demand on objective functions values.

Figures 6 and 7 demonstrate that 10% an increase and 10% decrease in return rate respectively, with changes in the return rate from −10% to +10% having a minimal impact on the performance of all objective functions. These objective functions exhibit a similar, stable response in the range of changes in the return rate, which indicates that these factors (cost, energy, and CO2 emission) are relatively insensitive to changes in return. Energy and Carbon Emissions show slight deviations, with minor decreases at lower return rates and small increases at higher rates, though the changes are minimal.

Effect of 10% increase in return rate on objective functions values.

Figure 6: Effect of 10% increase in return rate on objective functions values.

Effect of 10% decrease in return rate on objective functions values.

Figure 7: Effect of 10% decrease in return rate on objective functions values.

Managerial insights

Professionals unfamiliar with the complicated mathematical framework need to concentrate on understanding the primary decisions optimized in the CLMS design. The model considers economic and environmental factors when optimizing production, distribution, inventory, recycling, and location decisions. By minimizing the system’s total cost, and energy consumption, and lowering CO2 emissions, the model seeks to advance both environmental and economic sustainability within the CLMS.

In this article, manufacturing, transportation, recycling, and location problems for CLMS are introduced. The managers of supply chains can plan the placement of their facilities to maximize consumer demand with the aid of the model. The methodology minimizes overall supply chain costs, and decisions about production, distribution, inventory holding, recycling, and location will increase profits for those involved in the supply chain in achieving financial sustainability for their companies. Supply chain managers can also use this model to minimize the quantity of CO2 emissions emitted during transportation and to guarantee that products are moved properly across the CLMS. As a result, an industrial system can use the planned sustainable CLMS network.

The fuzzy multiobjective optimization model provides valuable results for managers to increase the stability of production. It assists managers in making better decisions despite uncertainties in parameters such as demand and regulations. This research shows how closed-loop systems need to be adopted in industries because of their benefits in reducing waste and cost savings due to recycling and reusable materials implementation. Management is urged to adopt the principles of circular economy to better resource utilization and supply chain resilience. Meanwhile, this study also demonstrates that the impact of sustainability on manufacturing processes changes competitive advantages in the long perspective, enhances customer loyalty, and helps to meet environmental requirements. For solving the issues of sustainable development, resource use, and waste reduction along the principles of the circular economy, the study presents practitioners and managers with operational tools and tactics to make their activities sustainable and cost-effective even in challenging environments.

Conclusion

In designing supply chains, manufacturers increasingly prioritize analytics that balance system efficiency with sustainability. This research contributes to this field by exploring sustainable development in closed-loop production processes and product transportation across all stages of the supply chain. It introduces a refined fuzzy multi-objective model for optimizing a green production framework, carefully considering various constraints such as cost, energy consumption, and emissions. By addressing these competing objectives, the proposed model offers an inclusive tool for manufacturers to achieve sustainable production while minimizing overall investment costs, energy usage, and carbon emissions. The model was tested using data from a real-world scenario, successfully resolving the inherent uncertainties in demand and returns within a closed-loop manufacturing system. The model achieved a notable aspiration level of 98.9%, demonstrating its effectiveness in navigating the complexities of closed-loop supply chains. These findings provide supply chain designers with actionable insights and a benchmark for evaluating and implementing sustainability measures across various industries. Additionally, this research lays the foundation for optimizing production levels, determining optimal market pricing, and managing resource allocation in a circular economy framework.

However, this model has certain limitations. Firstly, its complexity, involving advanced mathematical terminologies and fuzzy logic, may limit its application to experts in the field. Industries looking to adopt this model may need to hire skilled personnel or invest in training their existing staff, which can be time-consuming and costly. Moreover, while the model’s theoretical robustness is evident, the study relied on hypothetical data based on previous research due to the unavailability of real-time data from existing industries. This reliance on secondary data highlights a gap in the direct application of the model to real-world industrial settings.

To further enhance the model’s applicability and relevance, future research should focus on determining innovative strategies that incorporate environmental and social components into the optimization process more effectively. Practical evaluations of this model in real-world settings would provide valuable insights into its strengths and potential areas for refinement. Furthermore, expanding the scope of the study by integrating more objectives, such as waste management, customer satisfaction, and supply chain resilience, could yield a more complete framework for closed-loop supply chain management. Another critical area for future research is to explore the obstacles industries may face when implementing this model. For example, identifying barriers related to data collection and processing and the need for cross-sectoral collaboration can provide a clearer road map for operationalizing such models in practice. Overcoming these challenges will validate the model’s effectiveness and ensure that it can be adapted to different industrial contexts, making a significant contribution to the advancement of sustainable supply chain management.

Supplemental Information

MATLAB code to compute bounds.

DOI: 10.7717/peerj-cs.2591/supp-2

MATLAB code for the first objective.

DOI: 10.7717/peerj-cs.2591/supp-3

MATLAB code for the objective function.

DOI: 10.7717/peerj-cs.2591/supp-4

MATLAB code of second objective function.

DOI: 10.7717/peerj-cs.2591/supp-5

All the objective functions.

DOI: 10.7717/peerj-cs.2591/supp-6

MATLAB code for third objective function.

DOI: 10.7717/peerj-cs.2591/supp-7
1 Citation   Views   Downloads