Exergy-Based Optimization of a CO 2 Polygeneration System: A Multi-Case Study

: A polygeneration system for power, heat, and refrigeration has been evaluated and optimized using exergy-based methods. CO 2 is the working fluid. The study considered two environmental conditions for the potential implementation of the polygeneration system: cold (Case cold ) and hot (Case hot ). Aspen HYSYS ® was used to perform steady-state simulations, Python was used for the automation of the process, and the connection of Aspen HYSYS ® with Python was successfully applied for single-objective and multi-objective optimizations. A wide range of decision variables was implemented. The minimization of the average cost of a product per unit of exergy was the goal of single-objective optimization and was included in the multi-objective optimization in addition to the maximization of the overall exergy efficiency. Single-objective and multi-objective optimization were applied. Both optimization algorithms result in the necessity to increase the pinch temperature in the heat exchanger ( ∆ T pinch,HE ), maintain the pinch temperature in the gas cooler ( ∆ T pinch,GC ), and augment this value for the evaporator ( ∆ T pinch,EVAP ). Notably, higher isentropic efficiency for turbomachinery correlates with improved optimization outcomes. These findings contribute to the applicability and performance of the polygeneration system, offering potential advancements in sustainable energy solutions.


Introduction
The World Energy Outlook 2023 [1] (International Energy Agency, IEA) reports an average annual growth rate of 0.7% in total energy demand (Stated Policies Scenario) and flattening total energy demand (Pledges Scenario) because of improvements in the energy efficiency of existing technologies and advantages in technologies powered by electricity (mainly, electricity from renewable energy sources [2]) over fossil-fuel-based alternatives.In the net-zero emissions scenario, the primary energy demand has a negative growth rate of 1.2% per year until 2030.
Regarding CO 2 emissions, the technologies still unavailable on the market, i.e., those at the prototype or demonstration phase, delivered nearly 50% of the emissions reductions (net zero emissions scenario in 2050).The actual status is around 35%. Particular attention should be given to the off-grid solutions [1].
The configuration and the application of polygeneration systems depends on the scale [2,3]: • Large-scale industrial systems find a wide application in complex manufacturing processes such as petrochemicals, textiles, and food processing, where three energetic effects (electricity, heat (in the form of steam), and cooling) are required simultaneously.The production of chemical substances is also possible.Polygeneration is a sustainable solution for urban areas as district energy systems.These systems can be designed as a Energies 2024, 17, 291 2 of 17 synergy between the energy and non-energy sectors by centralizing energy production and distribution.

•
Middle-scale systems are well-known for their application in hospitals, university campuses, and small manufacturing plants.The number of residential, business, and commercial buildings that use air conditioning intensively has increased dramatically in recent years.A large amount of cold is required for servers (equipment for IT-related technologies, etc.).This is a global problem that needs attention.• Small-scale systems often belong to off-grid solutions [4].Such systems have high design flexibility and the highest potential for integration of renewable sources (solar and wind) up to 100%.This approach addresses the challenges of renewable energy variability by providing a consistent energy supply through complementary sources.
Small-scale off-grid systems contribute positively to the net zero emissions scenario.A considerable boost in system efficiency has been observed in polygeneration systems.
The modern term "polygeneration" can be treated as a synonym for "cogeneration", "tri-generation", and "multigeneration".Not only can energetic effects be included in the generated products list, but freshwater due to water desalination, hydrogen from water electrolysis, syngas from biomass gasification, etc., can also be included in the generated products (for example, [5,6]).
There are a large number of research papers dedicated to polygeneration systems.The Scopus database (November 2023) was used to identify these publications and cluster the keywords.Such an approach can help readers obtain an overview of the already-published results and highlight the directions for research.Thousands of publications were identified using the keyword "polygeneration"; later, they were filtered using the field of research "energy" with attention to "carbon dioxide" (as emissions and as the working fluid).Only publications in English were considered.Finally, a list of around 500 research publications was generated.The dynamics of the publications are shown in Figure 1a.Only papers published during the last 10 years were further used to cluster the keywords (Figure 1b).
From Figure 1b, we can see that a large cluster corresponds to the publications addressing the combination of heat and power generation.The topic of waste heat utilization is also well-identified, with several large clusters.However, tri-generation systems (bottom of Figure 1a) and "cooling" ("refrigeration") present relatively small clusters with weak connections to others.The research methods include theoretical research (and optimization), laboratory experiments, and pilot plants.The papers related to the optimization and application of exergy-based methods (thermoeconomics = exergoeconomics) form clusters in green and yellow colors; they were published during the last 3 to 4 years.Solar-powered, as well as solar-biomass-powered or solar-geothermal-powered, polygeneration systems have been the subject of extensive study, analysis, and optimization because of their potential to contribute positively to the goals of sustainability, namely, energy efficiency, economic feasibility, and reducing emissions.For example, a comprehensive review of perspectives on solar-driven polygeneration systems of different scales is reported in [7].The so-called 4E analysis has been applied to evaluate and optimize a medium-scale system based on the combination of solar and geothermal resources [8].The multicriteria analysis and optimization of the medium-scale polygeneration hybrid solar-biomass system are reported in [9].A novel stochastic planning model for renewable distributed generation (DG) in distribution networks was proposed, considering investment in conventional assets and smart grid assets like demand-side response and coordinated voltage control.The authors also considered active power generation curtailment.The model used a node-variable formulation to alleviate computational burden.The study demonstrated smart technologies' strategic value due to flexibility in system evolution [10].
can cover all scales, from very large (for example, [12,13], where the natural-gaspolygeneration systems are evaluated) to medium-and small-sale.For the generat the refrigeration and/or heat capacity, ejector systems [14] or their combinations [ be used.Vapor-compression refrigeration systems dominate in stand-alone applic driven by electricity (including small-scale systems driven by locally generated elec from renewables).However, compression refrigeration machines can also be driv rectly by internal combustion engines or small-or medium-scale power systems bas Waste heat utilization represents a large cluster in the research on polygeneration.A review of existing technologies and their perspectives is reported in [11].Such systems can cover all scales, from very large (for example, [12,13], where the natural-gas-based polygeneration systems are evaluated) to medium-and small-sale.For the generation of the refrigeration and/or heat capacity, ejector systems [14] or their combinations [6] can be used.
Vapor-compression refrigeration systems dominate in stand-alone applications driven by electricity (including small-scale systems driven by locally generated electricity from renewables).However, compression refrigeration machines can also be driven directly by internal combustion engines or small-or medium-scale power systems based on the organic Rankine cycle.A combination has a classification name: "thermally-driven compression refrigeration systems".
This paper addresses a small-and medium-scale polygeneration system with CO 2 as the working fluid.A conceptual design is shown in Figure 2. The authors already discussed and evaluated the application of such a system in [15,16].
nergies 2024, 17, x FOR PEER REVIEW the organic Rankine cycle.A combination has a classification name: "the compression refrigeration systems".
This paper addresses a small-and medium-scale polygeneration system the working fluid.A conceptual design is shown in Figure 2. The author cussed and evaluated the application of such a system in [15,16].During the last decade, CO2 as a working fluid for refrigeration, heat p or power generation has been extensively studied.For example, the author to this research in [17][18][19].
A CO2 cogeneration system, as the combination of the vapor-compress tion system and the recompression Brayton cycle, was evaluated in [20] usin ergy, and economic analyses.The parametric sensitivity analysis was cond cover a set of decision variables.One-criteria optimization (maximum ener maximum exergy efficiency, or minimum total product cost) has been cond generation, the optimal (minimal) cost of the products is 3.5% lower than the for optimal (maximum) energy and efficiency.Under the assumption tha capacity is the sole effect, the optimal (minimum) total product cost can be 19% compared to a base case, 4% lower than the cost obtained for optima energy and efficiency.
In [21], the system for utilizing the heat from exhaust gases from the s turbine system is evaluated.A cogeneration system consists of the supercri scritical CO2 cycles.The coefficient of performance (COP) of 2.75 was calc base case and can be improved by more than 10% because of improvements (supercritical) cycle.
A supercritical CO2 cycle with preheating was proposed for the waste from an engine [22].The maximum net power output can be increased by ar results demonstrate that improving the preheating process within the S-C achieve deeper utilization of the waste heat, leading to an increase in the During the last decade, CO 2 as a working fluid for refrigeration, heat pump systems, or power generation has been extensively studied.For example, the authors contributed to this research in [17][18][19].
A CO 2 cogeneration system, as the combination of the vapor-compression refrigeration system and the recompression Brayton cycle, was evaluated in [20] using energy, exergy, and economic analyses.The parametric sensitivity analysis was conducted to discover a set of decision variables.One-criteria optimization (maximum energy efficiency, maximum exergy efficiency, or minimum total product cost) has been conducted.For cogeneration, the optimal (minimal) cost of the products is 3.5% lower than the cost obtained for optimal (maximum) energy and efficiency.Under the assumption that refrigeration capacity is the sole effect, the optimal (minimum) total product cost can be decreased by 19% compared to a base case, 4% lower than the cost obtained for optimal (maximum) energy and efficiency.
In [21], the system for utilizing the heat from exhaust gases from the shipboard gas-turbine system is evaluated.A cogeneration system consists of the supercritical and transcritical CO 2 cycles.The coefficient of performance (COP) of 2.75 was calculated for the base case and can be improved by more than 10% because of improvements in the Brayton (supercritical) cycle.
A supercritical CO 2 cycle with preheating was proposed for the waste heat recovery from an engine [22].The maximum net power output can be increased by around 7%.The results demonstrate that improving the preheating process within the S-CO 2 cycle can achieve deeper utilization of the waste heat, leading to an increase in the entire system performance by another 7%.
A novel two-stage transcritical CO 2 refrigeration cycle was proposed in [23].Two ejectors are implemented in the vapor-compression cycle.The reported results show that Energies 2024, 17, 291 5 of 17 the improvement potential is between 20% and 80% (value of COP) and heavily depends on the operation conditions.The experimental results of the transcritical CO 2 systems are reported in [24,25].The first paper reports on an ejector refrigeration system, with particular attention to the design of the ejector.The second paper reviews experimental investigations of the heat-to-upgraded-heat and heat-to-power systems (terminology used in [25]).Particular attention is given to the heat transfer characteristics and the design of the compressors and expanders.
A multi-objective optimization for a large-scale integrated energy system based on energy hubs was conducted in [26].The method considers economic analysis, energy consumption, and environmental benefits and employs the ε-constraint-fruit fly optimization algorithm.Results show that the proposed method reduces annual total cost, primary energy consumption, and CO 2 emissions by approximately 2%, 5%, and 7%, respectively.
The bilayer model which utilizes historical data and conservativeness to mitigate uncertainties for multi-energy building microgrids was presented in [27].It determines the energy storage dispatch, demand-response, and on-off hourly operation of a power generation system.Numerical case studies show the effectiveness of the evaluation approach in achieving economically effective operations with high computational performance and immunity against uncertainties.
The potential of heat-driven vapor-compression refrigeration systems is still not fully understood.This study makes a substantial contribution to the continuing investigation of heat-driven vapor-compression refrigeration systems, with a particular focus on the off-grid operation.A gap in knowledge and application is addressed by describing the optimization processes and associated algorithms.It explores the complexities of simulation and automation processes.In addition, the study stands out because it explores various optimization goals by conducting single-objective and multi-objective optimization.
This paper aims to do the following: • To describe the simulation and automation processes as the prerequisites for optimization; • To evaluate comprehensively the polygeneration system with CO 2 as the working fluid to identify the decision variables; • To describe the optimization procedures and associated algorithms; • To conduct single and multi-objective optimization; • To conduct a comparative study between evaluated and reported solar-driven and waste heat-driven polygeneration systems.

System Description
Figure 3a shows the basic design of the polygeneration system being evaluated, and Figure 3b illustrates its thermodynamic cycle.The polygeneration system consists of nine components within two sub-cycles: • A power sub-cycle is composed of a compressor for the power cycle (CM-P), a heat exchanger (HE), and an expander (EX).The "driving energy" (for example, solar energy, heat from biomass combustion, waste heat, etc.) for the entire system is supplied to the HE.• A refrigeration sub-cycle includes a throttling valve (TV), an evaporator (EVAP), and a compressor for the refrigeration cycle (CM-R).The refrigeration capacity is generated within the EVAP.
A gas cooler (GC), a mixer (MIX), and a splitter (SPLIT) link these sub-cycles.Carbon dioxide is the working fluid for the entire system.The heat capacity is generated within the gas cooler.The system can generate net power as the difference between the power from EX and the power required by both compressors (CM-P and CM-R).

System Simulation and Automation
Aspen HYSYS ® (ver.10, AspenTech, Bedford, MA, USA) has been used as a simulation tool.To apply mass and energy balances, the Span-Wagner equations of state were selected.These equations are more accurate in predicting CO2 thermodynamic properties in a wide range of temperatures and pressures, including the vicinity of the critical point.
The simulation results obtained from Aspen HYSYS ® were exported for computing the analysis of system and component levels, evaluation, and optimization.With the aid of programming software, the parameters and values being exported from the simulation software and imported into the simulation software can be automated.Python was connected to Aspen HYSYS ® to automate the entire calculation process.The communication was managed through a binary interface component object model.Automation is advantageous, particularly for optimization.
The initially developed design needs to be evaluated.For this purpose, energetic and exergy-based analyses were conducted using a large set of simulations.The economic analysis has been conducted separately.The authors report this process in detail in [15].Thermodynamic performance, cost of the system product(s), environmental aspects, safety, and reliability were considered.The obtained results are used for this research.

System Simulation and Automation
Aspen HYSYS ® (ver.10, AspenTech, Bedford, MA, USA) has been used as a simulation tool.To apply mass and energy balances, the Span-Wagner equations of state were selected.These equations are more accurate in predicting CO 2 thermodynamic properties in a wide range of temperatures and pressures, including the vicinity of the critical point.
The simulation results obtained from Aspen HYSYS ® were exported for computing the analysis of system and component levels, evaluation, and optimization.With the aid of programming software, the parameters and values being exported from the simulation software and imported into the simulation software can be automated.Python was connected to Aspen HYSYS ® to automate the entire calculation process.The communication was managed through a binary interface component object model.Automation is advantageous, particularly for optimization.
The initially developed design needs to be evaluated.For this purpose, energetic and exergy-based analyses were conducted using a large set of simulations.The economic analysis has been conducted separately.The authors report this process in detail in [15].Thermodynamic performance, cost of the system product(s), environmental aspects, safety, and reliability were considered.The obtained results are used for this research.

System Evaluation
In this study, evaluation and optimization are conducted at the system level.A component-level evaluation was already reported by the authors in [15,16].Only physical exergy is considered for the exergy analysis.Chemical, kinetic, and potential exergy are neglected.The approach "exergy of fuel/exergy of the product" is applied to conduct the exergy and exergoeconomic analysis.
The exergy balance for the overall system is written as .E F,tot = .
and the overall exergy efficiency is where for the system being evaluated (Figure 3), the exergy of fuel is . .
Here, T av HS is the average temperature of the heat source supplied to the system in the HE.
The cost balance for the overall system is .
For calculating the capital investment rate .
Z k the total revenue requirement method (the application can be found in many papers, for example, in [17]) of the economic analysis was applied with the following assumptions: (a) the plant's economic lifetime is 20 years; (b) the effective interest rate is 10%; and (c) the average general inflation rate is 2.5%.For estimating the purchased equipment cost (PEC) for each component, the following was considered (detailed explanations are reported by authors in detail [15]):

•
Heater (HE) and gas cooler (GC) are compact printed circuit heat exchangers made of stainless steel.Therefore, the cost of such a heat exchanger is estimated by its weight [28]: • Evaporator (EVAP) [30]: where C B = 32,800 USD; X B = 80 m 2 ; M = 0.68, and f P = 1.3.
Energies 2024, 17, 291 8 of 17 • The cost of the throttling valve (TV) equals PEC TV = 100 EUR [17], and the costs of the mixer and the splitter are neglected.

System Optimization
Parameter optimization and structure optimization are two methods for improving the performance of the overall system.Product (single or cumulative for a multigeneration system) cost minimization is one of the most commonly applied objective functions for commercial systems, while energy (exergy) efficiency maximization is associated with thermodynamic optimizations, i.e., minimization of the primary energy for the energy-conversion system (and associated emissions).A set of decision variables and their constraints were discussed by the authors in [15].In this study, two types of parameter optimization have been applied: single-objective and multi-objective.
To solve the design optimization problems, stochastic algorithms-individual-based and population-based-are selected.The population-based algorithms widely applied in engineering are evolutionary, physical, and swarm-based algorithms.Genetic algorithms are the earliest and most well-known evolutionary algorithms.
Two relatively new algorithms (Figure 4) were selected and applied to solve the single-optimization problem: • The differential evolution (DE) algorithm was introduced in 1997 [31].The DE al- gorithm is similar to the genetic algorithm but modified for more straightforward implementation, less computation time, reliability, and robustness.• The particle swarm optimization (PSO) algorithm was first suggested in 1995 [32].
The potential candidates in this algorithm are called particles.A group of particles work together to continuously improve their individual and collective performance on a given optimization task [33].In addition, less computational effort is required for solving moderate-dimensional problems.It is also robust and straightforward.PSO is an excellent option for solving high-dimensional optimization problems.
Energies 2024, 17, x FOR PEER REVIEW 8 of

System Optimization
Parameter optimization and structure optimization are two methods for improvi the performance of the overall system.Product (single or cumulative for a multigenerati system) cost minimization is one of the most commonly applied objective functions commercial systems, while energy (exergy) efficiency maximization is associated w thermodynamic optimizations, i.e., minimization of the primary energy for the energ conversion system (and associated emissions).A set of decision variables and their co straints were discussed by the authors in [15].In this study, two types of parameter op mization have been applied: single-objective and multi-objective.
To solve the design optimization problems, stochastic algorithms-individual-bas and population-based-are selected.The population-based algorithms widely applied engineering are evolutionary, physical, and swarm-based algorithms.Genetic algorith are the earliest and most well-known evolutionary algorithms.
Two relatively new algorithms (Figure 4) were selected and applied to solve the s gle-optimization problem:

•
The differential evolution (DE) algorithm was introduced in 1997 [31].The DE alg rithm is similar to the genetic algorithm but modified for more straightforward i plementation, less computation time, reliability, and robustness.

•
The particle swarm optimization (PSO) algorithm was first suggested in 1995 [3 The potential candidates in this algorithm are called particles.A group of partic work together to continuously improve their individual and collective performan on a given optimization task [33].In addition, less computational effort is requir for solving moderate-dimensional problems.It is also robust and straightforwa PSO is an excellent option for solving high-dimensional optimization problems.
(a) (b) For multi-objective optimization problems, the non-dominated sorting genetic alg rithm-II (NSGA-II) is chosen [31] as one of the most popular and successful optimizati tools.The NSGA-II can produce a Pareto frontier, or non-dominated set of solutio For multi-objective optimization problems, the non-dominated sorting genetic algorithm-II (NSGA-II) is chosen [31] as one of the most popular and successful optimization tools.The NSGA-II can produce a Pareto frontier, or non-dominated set of solutions, which shows the trade-off between the two objective functions without any solution being superior to the others.
Optimization problems consist of three essential components:

• Objective functions
Single-objective parametric optimization: This involves the global optimization of decision parameters.Multi-objective parametric optimization: min c av p, tot andmax ε tot .
• Set of decision variables: , where ∆T represents the temperature differential (pinch point) in a heat exchanger; and η represents the isentropic efficiency of turbomachinery.

•
Constraints of decision variables: Two operation conditions were considered to represent the performance of the proposed polygeneration system in hot and cold climates, as follows:

•
Case hot for the hot climate operation with an average environmental temperature of 35 A medium-temperature heat source suitable for the evaluated polygeneration systems is in the range of 320-590 • C [15].

Single-Objective Parametric Optimization
The results of the single-objective optimization are shown in Table 1 (Case hot ) and Table 2 (Case cold ).The comparison of the values obtained using both DE and PSO algorithms is given as well.After the optimization for Case hot , the average cost of a product decreased by 25% from its original value, and the optimum solutions between DE and PSO are within the limit of 1.0%.The differences between DE and PSO optimum solutions are in the range of 2.5% for isentropic efficiency of turbomachinery and p Merging .The optimal minimum temperature in all heat exchangers, ∆T pinch , has relatively large (in the limit of 25%) differences between the initial and optimal values obtained by the DE and PSO algorithms.However, relative values should not be the only ones to be considered.The change in the absolute values of ∆T pinch by 1 K shows that the obtained results are within the limits of technical possibility for improving a heat exchanger while keeping the same design.
The optimization results for Case cold show a decrease in the average cost of a product by 15% from its original value, and optimum solutions between DE and PSO are at the limit of 1.0%.
When comparing the optimization results using DE and PSO algorithms, the largest relative difference of 60% is observed in the pinch point temperature difference for the EVAP.Starting with the initial value of ∆T pinch,EVAP = 5 K, the DE algorithm suggests increasing to 10 K, while the PSO algorithm suggests decreasing down to 4 K. Similar dynamics have also been suggested for the value of ∆T pinch,GC : a smaller value (6 K) suggested by the PSO algorithm compared to 7 K by the DE algorithm.All values obtained by optimization are within the limits of technical possibility for heat exchangers and turbomachinery.
The execution time for Case hot varied between 1820 s (DE) and 1633 s (PSO); while for Case cold , it varied between 1944 s (DE) and 1699 s (PSO).The PSO algorithm completes the optimization process roughly 200 s faster.For example, the best values of decision variables obtained from both optimization algorithms for Case cold are shown in Figure 5.The PSO algorithm already had a better starting position in the first iteration and continued to lower the objective function until the fifth iteration.A local minimum was identified within the sixth to ninth iterations, which had further minor corrections until the tenth iteration.No better value is obtained after the 10th iteration.Another dynamic is observed for the application of the DE algorithm.Significant progress was achieved in the second iteration; a local minimum was identified within the eighth to eighth iterations, with minor corrections within the fourteenth to fifteenth iterations.

Multi-Objective Parametric Optimization
Figure 6 shows the Pareto frontier for the entire polygeneration system at warm and cold operation conditions.As expected, an absolute optimal solution does not exist, indicating that the higher the system exergy efficiency, the lower the product cost.For the Casehot, an extensive flat profile is observed, which indicates a "slow" growth in the average product cost while improving the system efficiency.Such a flat area in the Casecols is

Multi-Objective Parametric Optimization
Figure 6 shows the Pareto frontier for the entire polygeneration system at warm and cold operation conditions.As expected, an absolute optimal solution does not exist, indicating that the higher the system exergy efficiency, the lower the product cost.For the Case hot , an extensive flat profile is observed, which indicates a "slow" growth in the average product cost while improving the system efficiency.Such a flat area in the Case cols is minimal.

Comparative Analysis with Solar-Driven Polygeneration Systems
A set of six papers [25,[33][34][35][36][37] was selected to compare the results obtained by the authors.In these papers, the exergy study, economic evaluation, and exergoeconomic variables were used for the optimization.Table 3 presents the comparative insights, which In the limit of decision variables, for Case hot operation conditions, maximum and minimum overall exergy efficiencies are 0.480 and 0.414, respectively, with the trade-off of the average cost of the product per unit of exergy 0.183 USD/kW hex and 0.141 USD/kW hex .A relative difference between the limiting values is an increase in the exergy efficiency of 16%, which corresponds to a decrease in the average product cost of 23%.The optimum values are opt ε tot = 0.443 and opt c av p, tot = 0.143 USD/kW hex .For Case cold operation conditions, optε tot = 0.549 and opt c av p,tot = 0.096 USD/kW hex .Maximum and minimum overall exergy efficiencies are 0.607 and 0.528, i.e., a relative difference of 15%.In contrast, the minimum and maximum values of the average cost of the product per unit of exergy are 0.095 USD/kW hex and 0.127 USD/kW hex , respectively, which is 25%.

Comparative Analysis with Solar-Driven Polygeneration Systems
A set of six papers [25,[33][34][35][36][37] was selected to compare the results obtained by the authors.In these papers, the exergy study, economic evaluation, and exergoeconomic variables were used for the optimization.Table 3 presents the comparative insights, which include net power production ( .W net , kW), overall exergy efficiency (ε tot ), and the mean total cost of the product (c av p, tot , USD/GJ ex ).The compared systems are classified as middlescale off-grid polygeneration systems, with capacities between 50 and 370 kW, which is what this comparative analysis is all about.The differences in the reported values are caused by evaluated system configurations, selected working fluids, particular operating conditions, and the evaluation and optimization routines in each study.Nevertheless, the authors' results fit the entire map of the results and contribute to the state-of-the-art.

Conclusions
The emphasis of this research was placed on the evaluation and optimization of the off-grid polygeneration system with CO 2 as the working fluid.Exergy-based methods were applied on the system level.The simulation of the system was conducted using Aspen HYSYS ® , and to ensure a comprehensive assessment and optimization process, an automated procedure was established to connect Aspen HYSYS ® with Python.
The goal of the single-objective optimization was to minimize the average cost per unit of exergy of a product.This objective transitioned into multi-objective optimization, where the emphasis extended beyond cost considerations to maximization of overall exergy efficiency.
This comprehensive approach contributes to the scientific understanding of polygeneration systems and lays the groundwork for informed decision-making in the pursuit of sustainable and efficient energy solutions.
The main obtained results can be summarized as follows: • The polygeneration system demonstrates sustainability benefits by simultaneously producing multiple energy effects, i.e., power, heat, and refrigeration capacities.• The polygeneration system exhibits adaptability for a wide range of thermally driven applications, from solar-thermal to waste heat utilization.

•
Parametric optimization methods, encompassing single-and multi-objective optimization, were applied to optimize the polygeneration system.Two operation conditions, Case hot and Case cold , were selected for thorough evaluation and optimization.

•
Single-objective optimization was completed using alternative optimization algorithms: DE (differential evolution) and PSO (particle swarm optimization).The implementation and impact of these algorithms on results are reported and discussed.

•
Both optimization algorithms performed effectively, consistently indicating the need to increase ∆T pinch,HE , maintain ∆T pinch,GC , and augment ∆T pinch,EVAP .Higher values for the isentropic efficiency of turbomachines correlated with improved optimization outcomes, with PR EX around 3. The optimal value of p Merging slightly increased, remaining significantly distant from the critical point of CO 2 .• The PSO algorithm demonstrated a shorter optimization process duration compared to the DE algorithm.
The following study directions can be recommended for future evaluation of the polygeneration systems:

•
Structure optimization of the polygeneration system including a polytechnological approach with many smart technologies within optimization [38]; • The environmental assessment of the polygeneration system; • Implementing the polygeneration systems to the local electrical grids; • The potential of mixtures as working fluids (including the CO 2 -based mixtures) to enhance system performance by adjusting their type and application scenarios.

Conflicts of Interest:
The authors declare no conflicts of interest.

Figure 2 .
Figure 2. A conceptual design of an off-grid polygeneration system.

Figure 2 .
Figure 2. A conceptual design of an off-grid polygeneration system.

Figure 5 .
Figure 5.The optimization process using the DE and PSO algorithms (Case cold ).

AuthorFunding:
Contributions: J.L.: Software; Investigation.B.T.: Methodology; Investigation; Writing the original draft.T.M.: Reviewing and Editing; Supervision.All authors have read and agreed to the published version of the manuscript.This work was conducted during a summer visit supported by Jordan University of Science and Technology (JUST) at the Technical University in Berlin (TUB) with grant number: 20230183.Data Availability Statement: Data are contained within the article.

Table 1 .
A single-objective parametric optimization utilizing DE and PSO algorithms for the polygeneration system at Case hot operation conditions.

Table 2 .
A single-objective parametric optimization utilizing DE and PSO algorithms for the polygeneration system at Case cold operation conditions.

Table 3 .
Key system variables of the reported solar-driven polygeneration systems and current work.