A Review of Multi-Objective Optimization in Organic Rankine Cycle (ORC) System Design

: Organic Rankine cycle (ORC) is considered a promising heat-to-power technology to utilize waste heat and renewable energy, including solar, biomass and geothermal. However, since the thermodynamic, economic and environmental performance is usually conflict, the single objective design could no longer meet the requirements of the ORC system, putting forward urgent requirements for multi-objective optimization, which has attracted increasing attention with lots of papers published. However, due to these different decision variables, optimization objectives and approaches, existing research is significantly different from each other and is difficult to compare without a systematic summary. Therefore, this paper provides an overview of ORC multi-objective research from three perspectives: optimization objective, method and optimization parameters. Based on the classification of different objectives, this work summarizes the involved variables and provides a recommendation for selecting appropriate objectives in different scenarios. For the optimization method, this work compares different approaches and reveals their advantages and disadvantages. Finally, the decision variables are reviewed and classified into four levels. Then the integrated design approach considering “system-process-component-fluid” is proposed and recommended for further development.


Introduction
Energy and climate crises are the common challenges facing mankind. Reducing fossil fuel consumption, increasing the proportion of renewable energy and energy efficiency are beneficial to realizing carbon neutrality in the middle of this century, which has become the consensus of the world. Organic Rankine cycle is a widely used power system in utilizing medium-to-low temperature thermal energy, which could effectively use renewable energy, including solar energy, geothermal energy, biomass energy and ocean energy [1]. Furthermore, ORC could also recover waste heat resources, including industrial waste heat and engine waste heat [2]. Therefore, ORC could improve the proportion of renewable energy, increase energy efficiency and reduce carbon emissions, thereby showing huge application potential in a low-carbon energy system in the future [3].
In the past decade, ORC has attracted widespread interest from worldwide scholars with an amount of research conducted. These studies mainly focus on system optimization design, working fluid selection, cycle configuration improvement, component design and control strategies, which have made significant progress [4]. Particularly, as the basis of fluid selection, component design and configuration comparison, system optimization has attracted the most attention [5]. Early research only aims to optimize a single objective, usually the maximum thermal efficiency or net output power. As the research develops, the economic and environmental performance of ORC systems such as the total cost, payback period (PBP), levelized cost of energy (LCOE), net present value (NPV), total carbon emission and sustainability indicator (SI) have gradually attracted increasing attention.
Recently, there have been many researches discussing the MOO in ORC system design. For instance, Wang et al. focused on maximizing the thermal efficiency and minimizing the heat exchanger area per power output of the subcritical ORC using the simulated annealing algorithm [9]. Yang et al. optimized the evaporator/condenser pressure and superheat degree to minimize the total investment cost and maximize the W net , using the multi-objective Genetic algorithm [10]. Furthermore, the MOO method has also been used in working fluid selection [11]. Feng et al. carried out a thermo-economic comparison between pure and mixture fluids in ORC, in which the exergy efficiency and Levelized energy cost (LEC) are selected as the objective variables [12,13]. Xi et al. proposed a graphical criterion, namely the Pareto frontier, to distinguish the optimal fluid for waste heat recovery, in which the exergy efficiency and annual cash flow were considered [14]. Ghasemian et al. compared eight fluids from the perspective of thermal, exergy and cost. Results indicated that dry fluids have a better performance of costs and exergy efficiency than wet fluids [15]. Moreover, the MOO method could be used to compare different ORC architectures. Lecompte et al. explored the specific investment cost (SIC) and total power output (W net ) of subcritical and transcritical ORC using the Non-dominated Sorting Genetic Algorithm II (NSGA-II). Results indicated that subcritical ORC shows lower W net but with better economics [6,16]. Sadeghi et al. compared the simple ORC, parallel two-stage ORC and series two-stage ORC from the viewpoints of turbine size and W net. Results indicated that STORC outperforms other configurations under the same conditions [17]. Song et al. compared the subcritical, transcritical, superheated and recuperated ORC considering the exergy efficiency and payback period Results indicated that the superheating is more suitable for the fluids with lower critical parameters, while the recuperation is not attractive under most operating conditions [18].
In general, MOO has been widely applied in ORC, from system design to cycle configuration and to fluid selection. The involved objectives are also numerous, including the conventional thermodynamic indicators such as exergy efficiency, power output and also the economic indicators such as the LCOE and SIC. Moreover, there are also multiple optimization methods such as the NSGA-II and weighted sum method (WSM). Due to these different decision variables, optimization objectives and approaches, existing researches are significantly different with each other and are difficult to compare or evaluate. Therefore, this work is carried out with two main purposes. The first is to review the optimization objectives, methods and variables involved in the existing research and summarize the main findings. The second is to refine further the important issues that deserve additional attention to provide valuable references for future work.
are significantly different with each other and are difficult to compare or evaluate. Therefore, this work is carried out with two main purposes. The first is to review the optimization objectives, methods and variables involved in the existing research and summarize the main findings. The second is to refine further the important issues that deserve additional attention to provide valuable references for future work.
Authors searched on the Web of Science using the keywords of "multi-objective" and "organic Rankine cycle (Including other expressions: multi objective, multi criteria, bi-objective, tri-objective, Pareto, ORC.)". Then, based on the searched results, the literature is further judged whether its focus is on ORC, removing the combined systems such as CCHP or CHP. Finally, a total of 122 papers are obtained and analyzed in this work, which are directly related to multi-objective optimization of ORC system design. As Figure 1 shows, most papers are published after 2015, revealing the popular trend of MOO in system design. There are mainly six chapters in this work. The second section discusses the optimization objectives from thermodynamic, economic, environmental and safe indicators. The third section summarizes the popular optimization methods and further compares their advantages/disadvantages. The fourth section classifies main decision variables from the system, process, component and fluid level. This work has two main contributions: (1) Summarize the main progress made in the existing multi-objective research and provide specific suggestions for selecting optimization objectives and methods in ORC optimization. (2) Propose a four-level design approach considering "system-processcomponent-fluid" aspects, which could design ORC architecture, component structure and working fluid molecules simultaneously. This approach is expected to improve the overall performance further and is worthy of future discussion.

Optimization Objective
Optimization objective is also called the target variable in MOO, referring to the variables to be maximized or minimized. The earliest and most studied evaluation indicator is the thermodynamic criterion, including the power output, thermal efficiency and exergy efficiency. Gradually, the economics attract increasing attention, including the total costs, area power ratio (APR), electricity production cost (EPC) and the specific investment cost (SIC) [19,20]. As the climate crisis deepens, the environmental performance has also been quantified and evaluated in ORC, including carbon emission and sustainability indicators. Moreover, due to the diversification of ORC application scenarios, some special needs should be considered. For instance, ORCs deployed on mobile devices to recover exhaust heat should simultaneously consider the space limitations, equipment weight and safety. In this section, different indicators in ORC will be sorted out and analyzed.

Optimization Objective
Optimization objective is also called the target variable in MOO, referring to the variables to be maximized or minimized. The earliest and most studied evaluation indicator is the thermodynamic criterion, including the power output, thermal efficiency and exergy efficiency. Gradually, the economics attract increasing attention, including the total costs, area power ratio (APR), electricity production cost (EPC) and the specific investment cost (SIC) [19,20]. As the climate crisis deepens, the environmental performance has also been quantified and evaluated in ORC, including carbon emission and sustainability indicators. Moreover, due to the diversification of ORC application scenarios, some special needs should be considered. For instance, ORCs deployed on mobile devices to recover exhaust heat should simultaneously consider the space limitations, equipment weight and safety. In this section, different indicators in ORC will be sorted out and analyzed.
According to the optimization objective, all 122 papers are counted and classified. A list of all reviewed papers could be found in the Supplementary Material, and this section mainly summarizes the key findings. Figure 2a shows the paper number involving different indicators. Results indicate that the thermodynamic criterion is the most popular indicator, especially the exergy efficiency, accounting for almost half of all papers. The economic criterion is followed, including the LCOE, SIC and total cost. In contrast, there are very few studies on environmental performance, volume, weight and safety. Only nine papers calculate the total carbon emission with four papers discussing the sustainability index. Figure 2b shows the number of objective indicators. The majority of studies (75%) only focus on bi-objective optimization. 16% of studies focus on tri-objective optimization.
In contrast, only one study discusses the optimization with over five objectives, due partly to the immaturity of high-dimensional optimization methods, which will be discussed in Section 3. list of all reviewed papers could be found in the Supplementary Material, and this section mainly summarizes the key findings. Figure 2(a) shows the paper number involving different indicators. Results indicate that the thermodynamic criterion is the most popular indicator, especially the exergy efficiency, accounting for almost half of all papers. The economic criterion is followed, including the LCOE, SIC and total cost. In contrast, there are very few studies on environmental performance, volume, weight and safety. Only nine papers calculate the total carbon emission with four papers discussing the sustainability index. Figure 2(b) shows the number of objective indicators. The majority of studies (75%) only focus on bi-objective optimization. 16% of studies focus on tri-objective optimization. In contrast, only one study discusses the optimization with over five objectives, due partly to the immaturity of high-dimensional optimization methods, which will be discussed in Section 3.

Thermodynamic
As the most basic index in ORC, the thermodynamic criterion indicates the system's ability to utilize heat sources, which is usually quantified by the power output, thermal efficiency and exergy efficiency. The basic concepts and equations are listed in Appendix A and will not be discussed in the main part.
Usually, the designer only chooses one of the power output and thermal efficiency as the optimization objective, according to the specific application scenario. In terms of thermal efficiency, it is more popular in solar or biomass ORC since the thermal oil cycle is usually included. This oil cycle belongs to the closed system defined by Zhai et al. [21], which is more suitable for efficiency evaluation. In contrast, the power output is more suitable for an open system, including the industrial waste heat, engine waste heat, geothermal and other scenarios without oil cycle, as shown in Figure 3. In these open systems, the outlet heat is directly discharged without being recycled.

Thermodynamic
As the most basic index in ORC, the thermodynamic criterion indicates the system's ability to utilize heat sources, which is usually quantified by the power output, thermal efficiency and exergy efficiency. The basic concepts and equations are listed in Appendix A and will not be discussed in the main part.
Usually, the designer only chooses one of the power output and thermal efficiency as the optimization objective, according to the specific application scenario. In terms of thermal efficiency, it is more popular in solar or biomass ORC since the thermal oil cycle is usually included. This oil cycle belongs to the closed system defined by Zhai et al. [21], which is more suitable for efficiency evaluation. In contrast, the power output is more suitable for an open system, including the industrial waste heat, engine waste heat, geothermal and other scenarios without oil cycle, as shown in Figure 3. In these open systems, the outlet heat is directly discharged without being recycled. The traditional exergy loss analysis for each component helps facilitate the optimal design of the component and ORC system [22]. Furthermore, the advanced exergy analysis methods have become increasingly popular in recent years,which decomposes the exergy loss into endogenous and exogenous parts, thereby determining the possibility of reducing the loss through technological improvement [23]. Therefore, if the endogenous loss of all components is summed up, the efficiency improvement limit of the ORC system could be obtained. For instance, Wang et al. used advanced exergy analysis to estimate The traditional exergy loss analysis for each component helps facilitate the optimal design of the component and ORC system [22]. Furthermore, the advanced exergy analysis methods have become increasingly popular in recent years, which decomposes the exergy loss into endogenous and exogenous parts, thereby determining the possibility of reducing the loss through technological improvement [23]. Therefore, if the endogenous loss of all components is summed up, the efficiency improvement limit of the ORC system could be obtained. For instance, Wang et al. used advanced exergy analysis to estimate the improvement potential of dual-loop ORC. Results indicated that the high-temperature turbine is the first component that needs technical modifications [23].

Economic
In addition to the thermodynamic performance, the economic criterion is also an important indicator to evaluate ORC performance, which has become a necessary option for ORC optimization. Up to date, the popular indicators include the indirect indicator, and direct indicators (total cost, SIC, PBP, etc). Particularly, exergoeconomic analysis is different from traditional economic analysis and deserves additional discussion.

Indirect Indicator
Indirect indicators do not directly measure the cost or economics of the whole ORC sysem, but use the parameters such as heat exchanger area to characterize the component cost qualitatively [24]. The commonly used index includes heat exchanger internal heat transfer requirement (UA) [25,26], heat exchanger area, turbine size, area to power ratio (APR) [27] and other indexes.
Specifically, the UA is simple to use since only the heat capacity and temperature difference are needed without complicated iterative calculation for U. Feng et al. used UA as an economic evaluation index to evaluate the cost of a supercritical-subcritical ORC in waste heat utilization [28]. Tiwari et al. optimized the exergy efficiency and the total UA of the evaporator and condenser. For a known working fluid, the total UA will exhibit the required area [29]. Further, the heat transfer coefficient U could be calculated to determine the heat exchanger area based on various heat exchange correlations [30]. Gotelip et al. designed the heat exchanger to minimize the total heat exchanger area and maximize the power of ORC in floating product storage offload [31]. The turbine diameter is also selected as the optimization objective to obtain a smaller turbine size [32]. Bahadormanesh et al. defined an objective function by dividing the turbine inlet radius by power output [33]. Rahbara et al. minimized the overall turbine size (d max ) for a highly compact ORC layout applied in diesel truck engines [34]. However, the smaller area or diameter does not necessarily represent lower costs since the impacts of material, pressure level, component structure are not considered [35].
These indirect indicators could only characterize the cost of a certain component. However, they could not represent the overall economy of the overall ORC system. For different system sizes, the results are difficult to compare, either. Therefore, more direct indicators are required to represent the costs of the whole ORC system.

Direct Indicator
Direct indicators are much more frequently used in existing research since they are strictly defined in engineering or economics, including the Total cost, Specific investment cost (SIC), Payback period (PBP), Levelized cost of electricity (LCOE), Net present value (NPV) and Internal Rate of Return (IRR), as shown in Table 1. Total cost is the most basic index to evaluate ORC economics [36], which are mainly calculated using empirical correlations, fitted from the equipment cost on the market [37]. The two most popular correlations are from Turton [38] and Smith [39]. In Sun's work, the carbon emission costs and water consumption costs were also considered [40]. Based on total cost, SIC describes the unit cost per power output [41] and has the advantage of easy use and intuitive comparison between different cases. Sun et al. maximized the exergy efficiency and minimized the specific power cost, revealing the better performance of double pressure ORC over traditional ORC [42]. However, SIC does not consider the depreciation, operation costs or discount rate [43]. PBP measures the number of years required to recover the total cost, including the static and dynamic ones. Generally, the dynamic PBP is more usual since it considers the time value over multiple years. For instance, Wang et al. used dynamic PBP to evaluate the superiority of the ORC investment plan in engine waste heat recovery [44]. LCOE denotes the cost of unit electricity and could be directly compared with the local electricity price to represent the profitability. In addition, the comparison between different energy generating technologies is also convenient, as suggested by Lecompte [6]. It's worth noting that many existing researches also use other abbreviations such as EPC or LEC to represent the cost per power generation, especially in the research of industrial/engine waste heat [45][46][47]. In this review, these abbreviations are unified as LCOE. The net present value (NPV) intuitively shows the profitability during the lifetime operation. Hu et al. used NPV to evaluate the total profits of a hybrid geothermal-solar power system during the 30-year operation [48]. Pirerobon et al. used NPV to evaluate the profitability of three waste heat recovery units for offshore platforms [49]. However, the expected discount rate should be defined in advance, which is related to the company's lowest investment rate of return and is usually difficult to determine. In contrast, IRR has an outstanding advantage in that it depends entirely on investment cash flow and is not affected by external parameters, including the discount rate [50]. However, IRR only represents a ratio rather than an absolute value. When the project scale is large, a low IRR may also result in large profits. Therefore, it is suggested that IRR and NPV should be simultaneously considered when comparing different schemes. In general, authors recommend the direct indicators in ORC design, especially the LCOE, NPV and IRR, which describe the lifetime economy of the entire system, rather than individual component or static periods.

Exergoeconomic Analysis
Exergoeconomic analysis combines advanced exergy analysis and economic analysis, aiming to explore the relative importance of each component in the design of the ORC system [51]. Different from the normal economic analysis, exergoeconomic analysis also calculates the exergy cost rate when considering the equipment cost and operation cost [52,53]. Özahi et al. calculated the capital cost rate of components using the SPECO method to determine the optimal fluid from toluene, MDM, D4 and n-decane [54]. Behzadi et al. identified the relative cost importance of each component and designed a cost-effective ORC system for waste recovery. Results indicated that R123 is the best working fluid from exergy efficiency and total production cost [55].

Environmental
The research on environmental performance started later than thermodynamic and economic performance. The most studied environmental performance is the carbon emission. Moreover, there are also some studies on sustainability indicators (SI) [56].

Carbon Emission
Carbon emission is the most intuitive and commonly used environmental indicator. The emissions of ORC system involve all processes through the whole life cycle, including production, transportation, operation and recycling. Specifically, the emission sources include: (1) the process of manufacturing the components, (2) leakage of organic fluids in operating process, (3) replacement of coal/natural gas using ORC. Some researches also discuss the impact of recycling water on carbon emission [57]. Moreover, the methane also results in the greenhouse effect, which could be converted into carbon emissions using coefficients. Common analysis method includes the Total equivalent warming impact (TEWI) method, Life cycle climate performance (LCCP), life cycle analysis (LCA) and carbon emission coefficient method, which are described in detail: (1) Total equivalent warming impact (TEWI) TEWI method is a measure of the greenhouse gas emissions, with the advantages of simple calculation and easy application. TEWI is mainly composed of two parts [58]: (1) direct emission: emissions caused by fluid leakage and end-of-life recovery. (2) Indirect emission: energy saving due to ORC power. Yang et al. introduced an extended TEWI method to evaluate the impact of refrigerants on climate change, considering that the generated power of ORC could compensate for the warming relief [58].
(2) Life cycle climate performance (LCCP) LCCP is developed based on the TEWI method, which considers chemical manufacturing and end-of-life disposal compared with TEWI, as shown in Zhang's work [36]. The diagram of LCCP is shown in Figure 4.

(3) Carbon coefficient
Energies 2021, 14, x FOR PEER REVIEW performance of the waste heat recovery systems with simple and regenerative ORC sults indicated that toluene presented the best efficiency and lower climate effects 0.00181 kg CO2 eq·kWh −1 [62].     The carbon coefficient method is a general method for direct multiplication coefficient to calculate carbon emissions. For instance, power generation could multiply a coefficient to estimate carbon emissions [23]. Xia et al. calculated the warming value considering the conversion factor of CO 2 , CO and CH 4 based on coal-fired power plants [59]. eration, maintenance and disposal [61]. Herrera et al. assessed the exergo-environmental performance of the waste heat recovery systems with simple and regenerative ORC. Results indicated that toluene presented the best efficiency and lower climate effects with 0.00181 kg CO 2 eq·kWh −1 [62]. Table 2 summarizes related researches on carbon emission. Results show that existing studies usually calculate emission from two or three categories. However, no research has considered all factors simultaneously, which may underestimate the ORC's environmental performance.

Exergoenvironmenal Analysis
Exergoenvironmenal analysis has a similar principle to advanced exergy analysis by linking the environmental impact with the energy flow and internal inefficiency, which is different from the direct calculation of carbon emission [53]. Fergani et al. compared three fluids in cement industry ORC by taking the exergy efficiency, cost and environmental impact per exergy unit of the net produced power as optimization objectives. Results indicated that benzene is the best from an exergo-environmental point of view [65].

Sustainability Index (SI)
The definition of SI is the ratio of exergy loss rate to the exergy decline of heat source (water, exhaust gas), which is used to evaluate the impact of ORC on the environment [66].
The key idea is to improve the exergy efficiency, reduce the emitted heat and alleviate the damage to the environment [67]. A lower SI means a smaller impact [68].

Other
In addition to the thermodynamic, economic and environmental indicators, there are also other evaluation indicators, including system volume, weight and safety.

Volume
ORC could be used on mobile vehicles and marine scenarios. Due to the limited space, ORC volume is a very important factor in these scenarios. Many researchers focus on the volume of the heat exchanger since it is the largest component in the ORC system. For instance, Xu et al. [69], Baldasso et al. [70], Liu et al. [71], Barbazza et al. [72] optimized the system parameters and structure by taking the heat exchanger volume as the optimization objective. Moreover, some researchers also consider the volume of turbines, feed pumps and auxiliary equipment [73]. Based on each component's structure, related empirical correlations could also be used to calculate the volume [74].

Weight
Similar to the volume, the weight of the ORC system is also an important factor in mobile vehicles. Imran et al. [73] and Pierobon et al. [49] used correlation equations to calculate the weights of turbines, heat exchangers, feed pumps, pipes, instrumentation/control systems and other auxiliary components. Then an additional mass of 10% is supplemented to represent the mass of storage tanks. Finally, the total weight is selected as one of the optimization objectives. More detailed equations could refer to these references.

Safety
Safety is a very important index in all application scenarios but is seldom discussed. Safety mainly refers to the risk of workers when exposed to the ORC environment, which is usually measured by the individual risk [57]. These risks are mainly caused by the leakage of flammable and toxic working fluids and high-pressure pipes. There are special methods for quantitative safety assessment, such as the Quantitative Risk Analysis (QRA).
Lee et al. simplified the QRA method and focused on the riskiest part, namely the pipeline between the feed pump and evaporation, since the pressure is the highest in this process [75]. The authors focused on source modeling, dispersion modeling and fire event modeling, and used the concept of risk distance to characterize the ORC safety. In addition, some scholars directly optimize the ORC system from the perspective of fluid design to minimize the flammability, thereby improving ORC safety. Papadopoulos et al. researched the molecular design of mixed fluids by minimizing the flammability of two pure fluids using the group-contribution method [76].

Stability
Most of the aforementioned studies only discuss a single working condition. During the practical operation, the ORC system always deviates from the design condition due to the influence of solar radiation, wind speed and ambient temperature. Therefore, some scholars take the power fluctuation as the optimization objective, and minimize the fluctuation range as small as possible by parameter optimization.
Li et al. focused on solar-based ORC and selected the fluctuation of output (W·min −1 ) as the optimization objective [77]. Results indicated that a larger energy storage capacity could reduce power fluctuation, but will significantly increase the costs. Bufi et al. focused on maximizing the thermal efficiency and minimizing its variance [78]. Zhang et al. proposed a multi-objective estimation of distribution algorithm to keep superheat following a target value by controlling the pump speed [79].

Optimization Method
Multi-objective optimization method is essentially different from single-objective optimization. A single optimal solution could be obtained in single-objective optimization. However, different indicators compete with each other, and there is no unique optimal solution in multi-objective optimization (MOO), which is also more complex and timeconsuming to converge. MOO is usually divided into the Priori method and No preference method. Further, the Priori method could be divided into the Apriori method, interactive method and Aposteriori method, according to whether the preference information is determined before, during or after the optimization process, as shown in Figure 5. At present, the Apriori method and evolutionary algorithm method are widely used in ORC, including the linear weighted sum method (WSM), ε-constraint method and smart algorithms such as NSGA-II, MOPSO and etc. netic algorithm and are not distinguished in many previous researches. Therefore, this review uses NSGA-II to represent these two methods. Results show that NSGA-II is the most popular algorithm, accounting for about 66% of all existing studies. The second popular method is WSM, which accounts for 16%. Other methods such as MOPSO and εconstraint method only account for 18%. Therefore, this work will take WSM, ε-constraint and intelligent algorithm as examples to introduce the principle and application in detail, and compare the advantages and disadvantages of each method.   This work has summarized the application of these methods in the ORC MOO application, as shown in Figure 6 [7,80]. Results show that, from the perspective of optimization methods, many interesting methods have not been applied in ORC, including the interactive methods that could feedback the decision makers' preferences during the design process. Applying these methods may make the system design more in line with the needs of designers and engineering projects, thus worth future exploration.

Weighted Sum Method (WSM)
3.1.1. Principle WSM combines multiple target variables into a single one according to a certain weight ratio, thereby transforming the multi-objective optimization problem into a simpler single-objective optimization problem [81], as shown in Equation (2).
where w represents the weight factor, ranging from 0 to 1. The sum of all factors is 1. WSM has a simple principle and is easy to use. There is no theoretical upper limit to the number of optimization objectives. Thus 2, 3, 5 or even more than 10 objectives could be combined into one [82]. For instance, Arasteh et al. combined the thermal efficiency and exergy efficiency into one objective function with each factor's contribution of 0.5. Then the Genetic Algorithm is used to solve this optimization problem [83]. Zhu et al. combined the exergy efficiency and the heat exchanger area per power output into one function. Then the optimization is conducted to determine the optimal evaporation temperature, condensation temperature and working fluid [84]. In addition to the Genetic algorithm, the PSO could also be used to solve this single-objective problem [71].

Methods to Determine the Weight
WSM is a priori method with the weight and preference being determined before In particular, MOGA and NSGA-II are both developed from the single-objective Genetic algorithm and are not distinguished in many previous researches. Therefore, this review uses NSGA-II to represent these two methods. Results show that NSGA-II is the most popular algorithm, accounting for about 66% of all existing studies. The second popular method is WSM, which accounts for 16%. Other methods such as MOPSO and ε-constraint method only account for 18%. Therefore, this work will take WSM, ε-constraint and intelligent algorithm as examples to introduce the principle and application in detail, and compare the advantages and disadvantages of each method.

Weighted Sum Method (WSM)
3.1.1. Principle WSM combines multiple target variables into a single one according to a certain weight ratio, thereby transforming the multi-objective optimization problem into a simpler single-objective optimization problem [81], as shown in Equation (2). where w represents the weight factor, ranging from 0 to 1. The sum of all factors is 1. WSM has a simple principle and is easy to use. There is no theoretical upper limit to the number of optimization objectives. Thus 2, 3, 5 or even more than 10 objectives could be combined into one [82]. For instance, Arasteh et al. combined the thermal efficiency and exergy efficiency into one objective function with each factor's contribution of 0.5. Then the Genetic Algorithm is used to solve this optimization problem [83]. Zhu et al. combined the exergy efficiency and the heat exchanger area per power output into one function. Then the optimization is conducted to determine the optimal evaporation temperature, condensation temperature and working fluid [84]. In addition to the Genetic algorithm, the PSO could also be used to solve this single-objective problem [71].

Methods to Determine the Weight
WSM is a priori method with the weight and preference being determined before optimization. Thus a consequent problem is: how to determine the weight factor of each target variable? In many previous studies, the weight factor is directly assumed. For instance, the weight is usually set as 0.5:0.5 [83,85] or 0.6:0.4 [56] when two target variables are used. When four target variables are considered, the weight is usually set as 0.1:0.2:0.3:0.4 [86]. This direct assumption usually only considers typical weight scenarios but is too subjective and does not consider the characteristics of each target variable. Therefore, in some researches, the weight factor is not subjectively determined by the designer, but using other mathematical models to calculate the coefficient according to each variable's characteristics. The typical methods include: the α-method, Analytic Hierarchy Process (AHP), and Grey relational analysis (GRA) [87].
(1). α-method α-method calculates the weight factor by considering the magnitude and the range of each optimization objective. A variable with a larger range will have a larger weight to ensure that two indicators have similar magnitudes after multiplying the weight [68,88]. This method is usually used for bi-objective optimization, in which the weight factor could be calculated by: , α-method is simple and easy to use. However, the maximum and minimum value of two optimization objectives should be determined in advance, which could be calculated through four single-objective optimizations, thereby significantly increasing the calculation complexity. For instance, Kazemi et al. firstly determines the weight coefficients of exergy efficiency and SIC using the α-method, and then explores the effects of system parameters [89].
(2). Analytic Hierarchy Process method (AHP) The AHP method is a kind of multi-criteria decision-making method, which divides complex problems into orderly levels to make them organized. Then the importance of each element is quantitatively compared and described, which could be used to calculate the weight factor of each variable, as shown in Figure 7 [90].
Zhang et al. discussed four target variables, including power, thermal efficiency, exergy efficiency and carbon emission by dividing these variables into four levels: energetic, exergetic, economic and environmental criteria. Then the author used nine absolute numbers (1-9) to indicate the importance intensity (equal, moderate, strong, extreme importance and etc.) and construct the judgment matrix [36]. It is worth noting that different indicators should be dimensionless. exergy efficiency and SIC using the α-method, and then explores the effects of system parameters [89]. (2). Analytic Hierarchy Process method (AHP) The AHP method is a kind of multi-criteria decision-making method, which divides complex problems into orderly levels to make them organized. Then the importance of each element is quantitatively compared and described, which could be used to calculate the weight factor of each variable, as shown in Figure 7 [90]. Zhang et al. discussed four target variables, including power, thermal efficiency, exergy efficiency and carbon emission by dividing these variables into four levels: energetic, exergetic, economic and environmental criteria. Then the author used nine absolute numbers (1-9) to indicate the importance intensity (equal, moderate, strong, extreme importance and etc.) and construct the judgment matrix [36]. It is worth noting that different indicators should be dimensionless. (

3). Taguchi method
Taguchi method is a statistical method to obtain the importance of different factors to the objective function, which could be used in the experiment, theory and numerical simulation [91]. Bademlioglu used this method to explore the impact of changes in different ORC parameters on system thermal efficiency and exergy efficiency. The author first used this method to obtain an orthogonal array of optimization objectives (thermal and exergy efficiency). Then this array is converted into a signal-to-noise (S/N) ratio in 27 cases. Consequently, the weight factors could be calculated as 52.61% and 47.38% [91].

ε-constraint
The core idea of the ε-constraint method is to regard a single indicator as the optimization objective and convert other objectives into constraints, thereby transforming the multi-objective optimization problem into a single-objective optimization problem. Different from the WSM method, ε-constraint could deal with the non-convex problems [63].
Multiple calculations could result in a relatively complete Pareto frontier. However, Figure 8 shows that if the ε is too small, f1 may exceed the feasible region, resulting in nonconvergence. If ε is too large, there will be more feasible solutions, resulting in lower convergence speed. Therefore, some prior knowledge is often needed to determine a reasonable range of ε. (

3). Taguchi method
Taguchi method is a statistical method to obtain the importance of different factors to the objective function, which could be used in the experiment, theory and numerical simulation [91]. Bademlioglu used this method to explore the impact of changes in different ORC parameters on system thermal efficiency and exergy efficiency. The author first used this method to obtain an orthogonal array of optimization objectives (thermal and exergy efficiency). Then this array is converted into a signal-to-noise (S/N) ratio in 27 cases. Consequently, the weight factors could be calculated as 52.61% and 47.38% [91].

ε-Constraint
The core idea of the ε-constraint method is to regard a single indicator as the optimization objective and convert other objectives into constraints, thereby transforming the multi-objective optimization problem into a single-objective optimization problem. Different from the WSM method, ε-constraint could deal with the non-convex problems [63].
Multiple calculations could result in a relatively complete Pareto frontier. However, Figure 8 shows that if the ε is too small, f 1 may exceed the feasible region, resulting in non-convergence. If ε is too large, there will be more feasible solutions, resulting in lower convergence speed. Therefore, some prior knowledge is often needed to determine a reasonable range of ε. For instance, Yi et al. firstly conducted two single-objective optimizations to determine the lower and upper bounds. Then authors implemented the ε-constrained method to solve this multi-objective optimization problem in GAMS 23.6 [63].

Intelligent Algorithm
Regardless of the WSM or the ε-constraint method, there is either the invalidity of the non-convex problem or the requirement for prior information, resulting in limitations to practical application. As the algorithm develops, some intelligent optimization algorithms with wider applicability have been gradually developed and improved, which

Intelligent Algorithm
Regardless of the WSM or the ε-constraint method, there is either the invalidity of the non-convex problem or the requirement for prior information, resulting in limitations to practical application. As the algorithm develops, some intelligent optimization algorithms with wider applicability have been gradually developed and improved, which have been widely used in different fields.
Popular intelligent algorithms include the NSGA-II [33], MOPSO [92], MOEA [93]. Non-dominated Sorting Genetic Algorithm II (NSGA-II) is an improved algorithm for NSGA based on GA's selection, crossover and mutation ideas, which was proposed by Deb in 2001 [94] It is worth mentioning that the gamultiobj function embedded in the Matlab toolbox is also a modified version of NSGA-II. Therefore, this review uses NSGA-II to simultaneously characterize the method of self-programming or calling the Matlab toolbox. The multi-objective particle swarm optimization (MOPSO) algorithm was proposed by Carlos A. Coello in 2004 for multi-objective optimization based on the PSO algorithm [95], which simplifies the crossover and mutation process and shortens the convergence time. The disadvantage of PSO is that it is easy to fall into local optimization, resulting in low convergence accuracy and poor solution diversity. Multiobjective Evolutionary Algorithm Based on Decomposition (MOEA/D) transforms the multi-objective optimization into a single-objective problem with the advantage of lower computational complexity [96]. The disadvantage is that the weight vectors need to be set artificially, which will determine the quality of the final solution [96].
In addition to the intelligent algorithms mentioned above, there are also other algorithms applied in ORC, including the multi-objective heat transfer search (MOHTS) [97], Artificial Cooperative Search (ACS) [98], multi-objective grey wolf optimizer (MOGWO) [99], multi-objective firefly algorithm (MOFA) [33], artificial bee colony algorithm (ABC) [100] and simulated annealing (SA) [101]. Even though these methods are rarely used, it will still be a very interesting topic to compare these different methods. However, for highdimensional optimization with 4 or more objectives, these intelligent algorithms are currently ineffective since the calculation time will increase significantly and the solution is not accurate, either. Therefore, WSM method is recommended for three or more optimization objectives, as shown in Table 3.

Decision Making
The multi-criteria decision-making method (MCDM) develops from scheme sorting in operations research and is also used to select the optimal solution from the Pareto in the MOO of ORC. This decision-making process selects the optimal solution for design guidance according to the decision maker's preference [103]. In ORC optimization, the typical MCDM methods include TOPSIS, LINMAP, Shannon entropy, GRA, fuzzy set theory, etc. The differences between these methods lie in the definition of the optimal solution. For instance, LINMAP only requires the solution closest to the ideal one [104], while TOPSIS requires the solution closest to the ideal solution and the farthest from the non-ideal one at the same time [105,106]. The Shannon entropy method could measure the uncertainty with information sources using the probability theory [107]. The key point is that the objective with a sharp distribution will have lower importance compared with that following the biased distribution [13]. As a part of grey system theory, grey relational analysis defines the black area, white area and grey area [108]. Theoretically, this analysis proposes a dependence to measure the correlation degree of factors, meaning that more similarity leads to more factor correlation [40].
Statistical results in Figure 9 show that nearly half of the multi-objective optimization studies use the MCDM method to determine the optimal solution. TOPSIS is the most popular approach, accounting for over 60%. LINMAP comes next, accounting for about 35%, while Shannon entropy and GRA methods are relatively less used. Since these MCDM methods have various concepts and usually lead to different final solutions, some researchers propose to apply multiple methods simultaneously and then determine the final solution using the aggregation method, which may improve the robustness of the decision-making process [13,90,103]. Detailed descriptions are shown in Table 4.

Optimization Parameter
During the MOO process, many ORC parameters could be optimized. The most popular ones are evaporation pressure, superheat, condensation pressure and other parameters, which all belong to the system level. In addition, there are also some parameters at Closest to the ideal solution.
Furthest to the non-ideal solution. [13,107] Shannon entropy A sharp distribution leads to lower importance [40,44,59,108,115] Grey relational analysis More similarity leads to more factor correlation

Optimization Parameter
During the MOO process, many ORC parameters could be optimized. The most popular ones are evaporation pressure, superheat, condensation pressure and other parameters, which all belong to the system level. In addition, there are also some parameters at the component level and working fluid level that could be optimized, which will be discussed in detail in this section.

System Level
System-level optimization parameters mainly include the evaporating pressure/ temperature, condensing pressure/temperature, subcooling and superheating. Under given conditions of heat source, the efficiency and net power output of ORC could be calculated according to the above system-level parameters. For the transcritical ORC, the superheat degree is not required, but the evaporation pressure and turbine inlet temperature should be determined at the same time [6,52]. For other new architectures such as dual-pressure evaporation ORC and two-stage ORC, the optimized parameters are more but are similar to the simple cycle [116]. System-level parameters are the most basic parameters of ORC, which are involved in almost all ORC optimization researches and will not be discussed in detail.

Process Level
Process-level design mainly refers to the design of cycle processes and system configurations, such as the conventional subcritical cycles, transcritical cycles, two-stage cycles, multi-pressure evaporation cycles. Most of the existing researches select the configuration by directly comparing the Pareto frontier of different forms through multi-objective optimization. However, this comparison could only study simple and several configurations. When there are multiple possible configurations, the computational complexity will increase sharply. Superstructure optimization could discuss various alternative configurations by analyzing the process stream, thereby parameterizing the ORC process design. Then the intelligent algorithms could be used to quickly solve the problem and obtain the best system structure and process, as shown in Figure 10. Kermani et al. [117] conducted a superstructure modeling for ORC systems driven by industrial waste heat, including regenerative, superheating, turbine-bleeding, reheating, multi-stage and transcritical cycles, etc. The multi-objective optimization is carried out with the net power output and total cost as the objective. Results indicate that the coupling of multiple forms results in better economics and thermodynamic performance.

Process Level
Process-level design mainly refers to the design of cycle processes and system configurations, such as the conventional subcritical cycles, transcritical cycles, two-stage cycles, multi-pressure evaporation cycles. Most of the existing researches select the configuration by directly comparing the Pareto frontier of different forms through multi-objective optimization. However, this comparison could only study simple and several configurations. When there are multiple possible configurations, the computational complexity will increase sharply. Superstructure optimization could discuss various alternative configurations by analyzing the process stream, thereby parameterizing the ORC process design. Then the intelligent algorithms could be used to quickly solve the problem and obtain the best system structure and process, as shown in Figure 10. Kermani et al. [117] conducted a superstructure modeling for ORC systems driven by industrial waste heat, including regenerative, superheating, turbine-bleeding, reheating, multi-stage and transcritical cycles, etc. The multi-objective optimization is carried out with the net power output and total cost as the objective. Results indicate that the coupling of multiple forms results in better economics and thermodynamic performance.

Component Level
In addition to the system-level and process-level parameters, the parameters at the component level also have a very significant impact on ORC performance. Key ORC components include the turbine, evaporator, condenser and feed pump. Up to date, most existing studies focus on the heat exchanger or the expansion device. The heat exchangers

Component Level
In addition to the system-level and process-level parameters, the parameters at the component level also have a very significant impact on ORC performance. Key ORC components include the turbine, evaporator, condenser and feed pump. Up to date, most existing studies focus on the heat exchanger or the expansion device. The heat exchangers mainly include the shell-and-tube heat exchanger, plate heat exchanger and tube-fin heat exchanger. The expansion device mainly includes the radial turbine, axial turbine and scroll expander [118]. It is worth mentioning that no multi-objective optimization on the feed pump is discovered in this literature review. Part of the reason may be that the cost of the feed pump is much lower than that of the turbine or heat exchanger. However, as one of the key components in ORC, the efficiency of the feed pump has a significant impact on system efficiency and operational stability. Therefore, the structural design of the feed pump could be further discussed in the future.
(1). Heat exchanger The heat exchanger is an important and also the largest component in ORC, including the evaporator, condenser, regenerator, preheater and superheater. The structure design of the heat exchanger is beneficial to reducing the exergy loss, improving system efficiency and increasing the safety and stability of ORC operation. Generally, the heat exchanger includes shell-and-tube, plate, tube-fin and plate-fin type [106,119,120]. The structure and optimization parameters of different heat exchanger types are also different, as shown in Table 5. In terms of the shell-and-tube heat exchanger, the optimization parameters mainly include the tube diameter, shell diameter, baffles spacing, tube number and tube length. In terms of the plate exchanger, the popular decision variables include the plate number, plate length, plate thickness, plate spacing, channel length, channel width and pass number. More detailed parameters are shown in Table 6.
(2). Expansion device The expansion device is the key component in ORC, which could be divided into the turbine, scroll expander and screw expander according to the power range. Ref. [118] introduces the categories and applications of different expansion devices in ORC. This work aims to reveal the parameters and turbine types in multi-objective optimization. According to the literature, most existing researches focus on radial turbine but pay little attention to other turbine types. The popular decision variables include the expansion ratio, specific speed, blade angle, blade height, radius ratio and flow coefficient. For an axial turbine, only Schilling et al. [121] have optimized the turbine stages but not discuss other structural parameters such as inlet angle and turbine size. For the positive displacement shown in Figure 11, both scroll expander and screw expander are more suitable for medium and small-scale ORC systems. However, there is currently no published work on incorporating expander structure parameters into the multi-objective optimization of ORC, which is worthy of further research. turbine, scroll expander and screw expander according to the power range. Ref [118] introduces the categories and applications of different expansion devices in ORC. This work aims to reveal the parameters and turbine types in multi-objective optimization. According to the literature, most existing researches focus on radial turbine but pay little attention to other turbine types. The popular decision variables include the expansion ratio, specific speed, blade angle, blade height, radius ratio and flow coefficient. For an axial turbine, only Schilling et al. [121] have optimized the turbine stages but not discuss other structural parameters such as inlet angle and turbine size. For the positive displacement shown in Figure 11, both scroll expander and screw expander are more suitable for medium and small-scale ORC systems. However, there is currently no published work on incorporating expander structure parameters into the multi-objective optimization of ORC, which is worthy of further research.

Fluid Level
Different organic working fluids significantly affect ORC performance. Therefore, determining a suitable working fluid from diverse types of fluids is particularly important for the practical application of ORC. In general, the current research on the selection of working fluids is divided into two levels. The first level is to compare the ORC performance when using different working fluids directly. The second level is to design the working fluid, that is, the working fluid is parameterized using the equation of state or the molecular group. Therefore, the second level is not limited to existing organic fluids, but can also calculate the ORC performance using potential working fluids in the future. This section will introduce the research status of these two fluid levels.

Selection from Limited Fluids
Most researches are based on comparing existing working fluids, including the pure fluids or mixtures. The typical approach is to construct the ORC model and calculate system performance, in which the physical property is obtained from software such as REFPROP or CoolProp. This method has the advantages of simplicity and high accuracy since most of the physical parameters in REFPROP have been verified and fitted by experiments.
(1) Pure fluids The optimization of pure fluid is mainly reflected in the comparison of ORC performance using different fluids but does not discuss or design the fluid itself. Therefore, the pure fluid will not be expanded in detail in this work. For more discussion, please refer to the existing literature [132].
(2) Mixture Due to the temperature slide during the phase change process, zeotropic mixtures could obtain a better temperature matching with heat/cold source, thereby reducing the heat loss in the heat exchanger and improving system performance. Different composition ratios could result in different temperature slides. Thus the design of zeotropic mixture mainly focuses on optimizing the component ratio, especially the ratio of binary mixtures [133]. However, the application potential of ternary or quaternary mixture is less discussed and could be explored in the near future.

Fluid Design
Even though the comparison of existing fluids is easy to conduct, it is usually difficult to determine the optimal fluid from diverse working fluids since the calculation process needs to be repeated many times, which significantly increases calculation time. Moreover, the direct selection of working fluids can only be based on existing working fluids. These diverse fluids usually have different critical temperatures, critical pressures, flammability and toxicity, which are difficult to quantify and compare with each other. Therefore, some researchers propose to incorporate the fluid characteristics into the design optimization of the ORC system, which is usually called computer-aided molecular design (CAMD) [122]. Currently, there are two popular molecular design methods. One is the contribution group design method, which calculates the physical properties of fluid by determining its chemical structure through the contribution group design. The other approach is to use the general equation of state (EOS) to parameterize the working fluid by characterizing the physical properties of the working fluid with critical temperature, critical pressure and eccentricity factor [134]. Consequently, the fluid characteristics could be included in the optimization process. In this section, the two methods will be introduced.
(1) Group-contribution method (GCM) The principle of GCM is that through the combination of different structural groups (-CH 3 , =CH-), any kind of organic fluid can be constructed theoretically. As shown in Figure 12, popular organic fluids in ORC such as R601a and R134a can be represented using the combination of common structural groups [135]. Using this principle, the physical properties of the working fluid are determined by the number of each group and its contribution ratio. In this GCM method, the contribution of each group is considered the same, regardless of the specific spatial structure and location. Therefore, the GCM method can predict the physical properties of both existing working fluids and also unknown fluids. GCM has a much faster calculation speed than CCM but has less accuracy. can predict the physical properties of both existing working fluids and also unknown fluids. GCM has a much faster calculation speed than CCM but has less accuracy. Papadopoulos et al. used the GCM method to generate a series of ORC candidate fluids, including existing fluids and also novel fluids. Then the optimal fluid is selected considering the requirements of safety and environmental performance [76,136,137]. Van Kleef et al. used the GCM method to estimate the critical physical properties and transport properties of the working fluid. Then the authors used a multi-objective optimization method to optimize the thermodynamic and economic performance simultaneously, and explored the optimal working fluid design in different heat source conditions [122]. Moreover, Su et al. reviewed the application of the GCM method in ORC and summarized the commonly used equations for predicting different macroscopic parameters of working fluids, such as boiling point, critical physical properties, transport properties, flammability and toxicity [135].
(2) Equation-of-State method (EOS) The equation of state method uses a few parameters, including the critical parameters and eccentricity factors, to calculate the basic and transport properties of a certain working fluid. In this way, the fluid properties and the thermodynamic and economic performance of the ORC system could be calculated. In traditional physical property software such as REFPROP, although high-precision multi-parameter equations of state can be used to ac-  Papadopoulos et al. used the GCM method to generate a series of ORC candidate fluids, including existing fluids and also novel fluids. Then the optimal fluid is selected considering the requirements of safety and environmental performance [76,136,137]. Van Kleef et al. used the GCM method to estimate the critical physical properties and transport properties of the working fluid. Then the authors used a multi-objective optimization method to optimize the thermodynamic and economic performance simultaneously, and explored the optimal working fluid design in different heat source conditions [122]. Moreover, Su et al. reviewed the application of the GCM method in ORC and summarized the commonly used equations for predicting different macroscopic parameters of working fluids, such as boiling point, critical physical properties, transport properties, flammability and toxicity [135]. The equation of state method uses a few parameters, including the critical parameters and eccentricity factors, to calculate the basic and transport properties of a certain working fluid. In this way, the fluid properties and the thermodynamic and economic performance of the ORC system could be calculated. In traditional physical property software such as REFPROP, although high-precision multi-parameter equations of state can be used to accurately obtain the thermophysical properties of the working fluid, these equations have many adjustable parameters and require a large number of experimental fittings. The forecast for novel fluids is not precise enough for ORC performance prediction. In contrast, the general-purpose state function is simpler and has fewer parameters, which can be more reliably extended to unknown working fluids, as shown in Figure 13 [134,138]. (3) Computational chemistry method (CCM) method In addition, there is a computational chemistry method (CCM) from the view of autom, which has a too long convergence time (usually from hours to days for a molecular structure). Therefore, this CCM method is not currently suitable for fluid design in the ORC system [135]. Comparison of three CAMD methods is shown in Table 7.

Discussion and Future Outlook
The single-objective optimization could no longer meet the requirements from the increasingly serious climate crisis, different technical requirements and rich application scenarios. Instead, the thermodynamic, economic, environmental and other performance should be considered simultaneously, which, however, usually conflict with each other. Thus it is difficult to obtain a unique optimal solution, proposing challenges to ORC multiobjective optimization. In recent years, there have been increasing studies involved in multi-objective optimization and discussing the ORC design. Correspondingly, multi-objective optimization has also become increasingly popular and gradually become a typical  Yang et al. proposed a corresponding states model to predict the fluid properties using five parameters: critical temperature, critical pressure, acentric factor, molar ideal gas isobaric heat capacity and temperature gradient [134]. In this way, different working fluids could be parameterized. The authors compared 13 common existing fluids and obtained the thermodynamic performance limit, which characterizes the best performance that can be obtained by improving the working fluid under ideal conditions.

(3) Computational chemistry method (CCM) method
In addition, there is a computational chemistry method (CCM) from the view of autom, which has a too long convergence time (usually from hours to days for a molecular structure). Therefore, this CCM method is not currently suitable for fluid design in the ORC system [135]. Comparison of three CAMD methods is shown in Table 7.

Discussion and Future Outlook
The single-objective optimization could no longer meet the requirements from the increasingly serious climate crisis, different technical requirements and rich application scenarios. Instead, the thermodynamic, economic, environmental and other performance should be considered simultaneously, which, however, usually conflict with each other. Thus it is difficult to obtain a unique optimal solution, proposing challenges to ORC multiobjective optimization. In recent years, there have been increasing studies involved in multi-objective optimization and discussing the ORC design. Correspondingly, multiobjective optimization has also become increasingly popular and gradually become a typical method for system design.

Discussion
In general, existing research has made significant progress, which could promote the development of many new indicators or methods and lay the basis for the practical application of the ORC system. The progress could be divided into the following areas:

Optimization Objective
The optimization objectives of the ORC system have been greatly enriched, including the thermodynamic, economic, environmental indicators, as well as the system weight, volume, safety and stability, which provide valuable guidance for the design and operation of ORC system. Different indicators have their advantages, disadvantages and application scopes, which are also summarized in this article. However, most of the existing research focuses on thermodynamic performance, followed by economic performance, and less attention to other properties. Especially when the climate issue becomes increasingly important, the research on the environmental protection performance of the system, the role of emission reduction, and the positioning of the ORC system in low-carbon energy scenarios are still insufficient.
In terms of selecting the optimization objective, unless specific preferences are required, authors are recommended to consider the thermodynamic, economic and environmental performance simultaneously. Other indicators such as safety could be represented by the working fluids. For each category, a single indicator could be selected according to the practical application scenario, as shown in Table 8. The power output and thermal efficiency are recommended for thermodynamic performance. Particularly, the power output is more suitable for geothermal, waste heat scenarios, which does not include the thermal oil cycle. In contrast, the thermal efficiency is more suitable for solar, biomass scenarios with oil cycles. In terms of the economic criterion, LCOE and IRR are recommended since these two indexes are independent of system size. The calculation of IRR is based on electricity price, thereby being more suitable for the region with a benchmark price. In contrast, LCOE describes the electricity price and is more suitable for the region with the electricity market. Particularly, LCOE is applicable for engine waste heat recovery in cars or marine applications since the electricity is usually consumed locally and not sold to the grid. Moreover, the system volume and weight should also be considered in these waste heat scenarios since the volume is very limited. carbon emission during production and transportation, and also consider the emission reduction by replacing other plants using ORC, thereby improving the accuracy over separate LCA or carbon coefficient method.

Optimization Method
As the algorithm develops, ORC optimization approaches are also gradually improved. The number of optimization objectives increases from a single one to two, three and multiple. The approach varies from the linear weighting method to the ε-constraint method, and the intelligent algorithms, which could obtain the Pareto frontier more conveniently. Furthermore, based on the Pareto frontier, many studies introduce the multi-criteria decision theory to select the desired solution with corresponding weights. Thus the unique solution could be determined to guide engineering design according to the decision maker's preference. However, most of the studies are based on optimization with three or fewer objectives. When the objective number exceeds three, the current intelligent algorithms could not converge very well. Up to date, only the traditional weighted method could be used to deal with high-dimensional optimization (>3) by combining multiple indicators into one. However, this method only considers a special scenario, and the weight is difficult to determine.
Based on this review, the authors recommend the intelligent algorithms in the optimization with 2-3 objectives, such as the NSGA-II, which is relatively simple to use and could effectively solve non-convex problems. These non-convex problems are very common in the design of configuration and system planning in complex thermal systems [140], which cannot be effectively solved by the WSM method. However, when the optimization objectives are more than or equal to 4, these intelligent algorithms such as NSGA-II could not work effectively with an explosive computation time and poor convergence at this stage. Therefore, the authors recommend the WSM method for this high-dimensional optimization. The specific weights can be calculated according to the AHP or α-method discussed in Section 3.1.
In addition, when using intelligent algorithms for optimization, it is usually necessary to adopt a decision-making method to obtain a unique solution for engineering applications. The authors recommend using multiple decision-making methods simultaneously by comparing the decision-making results to determine the final solution, which could effectively improve the robustness of decision-making.

Optimization Parameter
According to Chapter 4, the optimization parameters involved in existing research cover four levels: working fluid, components, processes and systems. Particularly, the system level is the most studied, focusing on the optimization design of evaporation/condensing pressure, superheat degree and etc. The research on process mainly aims to design the system structure and architecture by establishing a superstructure model. The research on components focuses on heat exchangers (tube and shell type, plate type and etc.) and turbines (mainly radial), while little attention is paid to feed pumps. For the working fluid, the general approach of using the group contribution method or the equation of state method to determine the optimal working fluid has been widely applied. However, the current research only focuses on one of these levels, and the research on multi-objective optimization of ORC systems considering multiple (four) levels is still scarce.
Authors recommend the collaborative design approach by considering multiple levels, especially the system, fluid and component levels. At the fluid level, the GCM or the EOS method could be used to parameterize the physical properties of the fluid and incorporate them into decision variables. At the component level, component size and structural parameters could also be included in the optimization parameters. Finally, a comprehensive design method considering four levels of system-process-component-fluid will be developed.

Future Outlook
In view of the above research status and problems, the following directions need to be developed.

Discussion on Carbon Emissions
The climate issue has become a major challenge threatening human survival. To achieve carbon neutrality, it is necessary to make concessions or even sacrifice certain economics. Therefore, future attention to environmental performance needs to be further strengthened.
Even though ORC could use turn heat resources into treasure (power), it is not completely carbon-free. Existing studies have used various methods such as LCA, LCCP, TEWI to measure its carbon emissions during the manufacturing process. However, there are some important issues that need further exploration, including: (1) Under the existing energy structure, how much room for emission reduction could be achieved by large-scale promotion of ORC? (2) In the future energy structure with a high proportion of renewable energy, how does ORC's emission reduction compare to its lifetime carbon emissions? Should ORC be promoted on a large scale? (3) With the continuous improvement of the carbon trading market, could ORC independently or in conjunction with other units benefit from the carbon trading market?

High-Dimension Optimization
When the number of optimization objectives exceeds 3, the results and convergence of existing intelligent algorithms, including NSGA-II and MOPSO, cannot meet the requirements. At present, only a few studies have focused on high-dimensional multi-objective optimization [36,64,86,141], and the methods used are mainly linear weighted method or multi-criteria decision-making method. However, as mentioned above, the weights should be determined in advance, which makes their application limited. Therefore, the following issues need to be further discussed: (1) How to simultaneously consider multiple indicators such as thermodynamic, economic, environmental and safety index for high-dimension Pareto optimization? (2) How to reasonably convert the preferences of decision-makers and the constraints of practical projects such as limited funds and limited equipment space into weights and reflect them in the objective function?

Consider the Practical Operation
Most of the existing multi-objective optimizations are based on design conditions. That is, the heat source and environmental conditions are assumed the same. However, this assumption is inconsistent with practical operating conditions, resulting in significant differences between expected operation and practical operation [111], or even leading to completely opposite conclusions [142,143]. Therefore, some research proposed to integrate long-term off-design operation into the system design process. Specifically, the typical heat source conditions and ambient temperature are selected as boundary conditions to guide the multi-objective design, resulting in a more robust design scheme [116,144]. Even though this method is closer to practical engineering, it also has limitations such as time-consuming calculations and uncertain weather and heat source conditions. Even though Hu et al. proposed an approach to reduce the calculation time using the ANN method, a duration of 1-2 days is still required [48]. Yang [145,146], but the mapping between environmental fluctuation and ORC output needs to be further established. To improve the practicability of this off-design approach, some novel methods need to be improved to reduce the calculation time and improve precision. The following aspects could be explored in the future: (1) How to integrate the dynamic response of ORC into system design? Especially the unit response speed (for heat exchanger or heat storage) and the loss during the start-stop process.
(2) How much impact does the fluctuation of user-side load have on ORC design? (3) How much influence will participation in the carbon trading market have on system design?  (4) How to improve the design method to speed up the optimization process based on a large number of off-design conditions and dynamic responses? Possible methods may include deep learning methods such as CNN, and the specific approach still needs to be developed.

Multi-Level Optimization
At present, most researches focus on the optimization of the system level, or the coupling of the two levels such as "system-component" and "system-fluid". There is still a lack of collaborative design of "system-process-component-fluid", which may further improve the results over conventional system optimization or component design, as shown in Figure 14. In this integrated design process, there are numerous decision variables, including the physical parameters of the working fluid, the structural parameters of components, the process parameters and the system parameters. Therefore, the optimization method and calculation capabilities need certain requirements. In general, the following aspects need further exploration: (1) Existing studies paid little attention to fluid pumps. How does the structural design of the fluid pump affect the thermodynamic and economic performance of the ORC system? (2) How much could the proposed "system-process-component-fluid" method improve system performance compared with conventional methods? What is the computational cost? (3) What's the difference between optimization results of four levels, three levels and two levels? Which level of optimization is most appropriate considering the improvements and computation cost? does the structural design of the fluid pump affect the thermodynamic and economic performance of the ORC system? (2) How much could the proposed "system-process-component-fluid" method improve system performance compared with conventional methods? What is the computational cost? (3) What's the difference between optimization results of four levels, three levels and two levels? Which level of optimization is most appropriate considering the improvements and computation cost? Figure 14. The multi-level optimization approach considering "system-process-component-fluid".

Conclusions
A single evaluation indicator could no longer meet the development needs of power systems such as ORC. Collaborative optimization considering multiple indicators such as thermodynamic, economic and environmental indexes has become an inevitable trend in system design. This article summarizes the research on multi-objective optimization of ORC system design. From the perspective of optimization objective, optimization methods and optimization parameters, this work classifies and statistically summarizes the existing research, and prospects the future development. Main conclusions are as follows 1. For the optimization objectives, multiple aspects should be considered, including the thermodynamic, economic and environmental indicators, which should be selected according to specific application scenarios and the local electricity market. For the thermodynamic indicator, the output power is recommended for geothermal and waste heat. Thermal efficiency is recommended for solar energy and biomass. For economics, LCOE is recommended for benchmark electricity price, and IRR is recommended for the electricity market. Furthermore, for the waste heat recovery in mobile transportation, additional attention should be paid to the ORC volume and weight. 2. In terms of the optimization method, intelligent algorithms including NSGA-II and MOPSO are recommended for low-dimension optimization (objective is less than 4).  Figure 14. The multi-level optimization approach considering "system-process-component-fluid".

Conclusions
A single evaluation indicator could no longer meet the development needs of power systems such as ORC. Collaborative optimization considering multiple indicators such as thermodynamic, economic and environmental indexes has become an inevitable trend in system design. This article summarizes the research on multi-objective optimization of ORC system design. From the perspective of optimization objective, optimization methods and optimization parameters, this work classifies and statistically summarizes the existing research, and prospects the future development. Main conclusions are as follows 1.
For the optimization objectives, multiple aspects should be considered, including the thermodynamic, economic and environmental indicators, which should be selected according to specific application scenarios and the local electricity market. For the thermodynamic indicator, the output power is recommended for geothermal and waste heat. Thermal efficiency is recommended for solar energy and biomass. For economics, LCOE is recommended for benchmark electricity price, and IRR is recommended for the electricity market. Furthermore, for the waste heat recovery in mobile transportation, additional attention should be paid to the ORC volume and weight.

2.
In terms of the optimization method, intelligent algorithms including NSGA-II and MOPSO are recommended for low-dimension optimization (objective is less than 4). In contrast, the WSM method is recommended for high-dimension optimization (the objective number is more than 3), which should determine the weight in advance and then convert the complex problem into single-objective optimization.

3.
For optimization parameter, the researches are mainly carried out from four levels: system, process, component and fluid. However, the existing researches are relatively independent. A comprehensive design method that couples multiple levels simultaneously is expected to improve the system performance further.

4.
For future development, the following aspects could be further explored: (1) Exploring the emission reduction potential of ORC in the future high-proportion renewable energy system and its economic value in the carbon trading market. (2) Transform the preferences and engineering constraints into boundary constraints in multi-objective optimization.
(3) Propose and expand the comprehensive optimization methods considering four levels "system-process-component-fluid". The product of the overall heat transfer coefficient and the total area WSM Weighted sum method the heat source is solar heat transfer oil, the power consumption of the oil pump should be considered, as shown in Equation (A1).
Some studies further use the total power generation to represent electricity generation over a period (usually a single year or 20 years). The difference lies in considering the generator and motor efficiency, although they are usually 95% [20].
where η is the efficiency of turbine or pump. E is the total electricity generation.
Appendix A.1.2. Thermal Efficiency Different from the power output, thermal efficiency measures the system's external output and input energy simultaneously, characterizing the utilization level of the input heat by ORC. Thermal efficiency could be calculated by Equation (A3).
where Q in denotes the absorbed heat by ORC, m H denotes the mass flow rate of the heat source. h H,in and h H,out represents the inlet and outlet enthalpy of the heat source.

Exergy Efficiency
Exergy efficiency further considers the energy grade and describes the effective utilization of exergy by ORC. In particular, the exergy loss analysis for each component helps facilitate the optimal design of the component and ORC system. The calculation of exergy efficiency is as follows: where the Ex in , Ex out , Ex loss represent the exergy at the heat source inlet, outlet and exergy loss in ORC, respectively.

Appendix A.2. Economic Index
Appendix A.2.1. UA UA could evaluate the heat exchanger cost according to the log mean temperature difference (LMTD) method [28,29]. A lower UA indicates lower costs and better economic performance [30], which could be calculated by: UA has the advantage of simple calculation. However, UA does not consider the impact of different working fluids and heat transfer capabilities, resulting in a relatively large cost deviation.

Appendix A.2.2. Total Cost
Total cost is the most basic index to evaluate ORC economics. Almost all direct economic indicators are calculated based on the total cost. The component costs are mainly calculated using empirical correlations, fitted from the cost of different types and sizes of equipment on the market. The two most popular correlations are from Turton [38] and Smith [39]: where C i denotes the total investment cost of each component, including the turbine, heat exchanger, pump. CEPCI denotes the correction to inflation or deflation. m denotes the benchmark year when fitting the correlations [37].
Appendix A.2.3. Specific Investment Cost (SIC) SIC is a very common index to evaluate the thermo-economic performance of ORC, which describes the unit cost per power output and could be calculated by [41]: SIC has the advantage of easy use and intuitive comparison between different cases. The disadvantage is that SIC is too simplified and does not consider the depreciation, operation costs or discount rate [43]. PBP measures the number of years required to recover the total cost, mainly including the static and dynamic PBP [44]. The calculation processes are shown in Eqs. A8 and A9, respectively. The dynamic PBP is more frequently used since it considers the time value and has higher accuracy than static PBP.
Appendix A.2.5. Levelized Cost of Electricity (LCOE) LCOE denotes the cost of unit electricity considering the project construction, operation and maintenance, depreciation and residual value [6]. This indicator could be directly compared with the local electricity price to represent the profitability. If LCOE is lower than the electricity price, then this project will be economically feasible. The calculation process is: where C resudual represents the residual value of the asset at the end. The denominator represents the discount of the annual power generation.
Appendix A.2.6. Net Present Value (NPV) The net present value (NPV) represents the difference between the discounted cash flow of an investment in the future and the total cost [48]. The expected discount rate is determined according to the company's lowest investment rate of return, which is the lowest acceptable limit. A positive NPV represents a feasible project, and a larger NPV represents a better return on investment. NPV could be calculated by: where the pc represents the electricity price. E yr represents annual electricity generation.
LT is the operating years of ORC. C OM is the operation and maintenance cost. DEP is the annual depreciation with linear or accelerated depreciation. tax is the tax rate of electricity sales, which varies in different countries or regions. r denotes the discount rate, which could be taken as the company's expected return on investment. IRR is defined as the discount rate when the total present values of capital inflow and the total outflow are equal to the net present value of 0. IRR could also represent the ability to resist inflation, the calculation of which usually requires relatively complex iteration. A higher IRR indicates better economic performance. Generally, the project is considered feasible if IRR is higher than the benchmark value of the return rate [50]. where C k denotes the exergy cost of component k. Z k represents the total investment and operation cost of component k.

Appendix A.3. Environmental Index
Appendix A.3.1. Total Equivalent Warming Impact (TEWI) TEWI = GWP · M · L · n + GWP · M · (1 − a recovery ) + E annual · β · n (A14) where L is the annual leakage rate, kg [58]. n is operating life, year. M is the refrigerant charge, kg. a recovery is the recycling factor. E annual is the annual energy consumption, kWh. β is the indirect emission factor, kg·kWh −1 . where n is the lifetime [36]. MM is the CO 2 production of material, kg·kg −1 . m is the mass of unit material, kg. RFM is the refrigerant manufacturing emissions, kg·kg −1 . The environmental balance in LCA is shown below, which helps in exploring the origin of pollution in each process [62]: where M i denotes the mass input. E i is the energy input. P i denotes the outlet stream. W i represents the residues. v is the synergy matrix of mass and energy.

Appendix A.3.4. Exergoenvironmenal Analysis
Main process is divided into three steps [53]. The first step is to perform exergy analysis on the overall system to identify the exergy loss in each component. The second step is to use the LCA method to analyze each component and corresponding energy input. The third step is to allocate the environmental impact to exergy flow in ORC using LCA. . Y k denotes the environmental impact during the production, transportation and installation. By establishing the balance equation of each component, the evaluation index F env of exergoenvironmental performance could be calculated by: where B D,k represents the environmental impact caused by exergy loss in component k.
Appendix A.3.5. Sustainability Index (SI) SI could be calculated by: where Ex loss represents the exergy loss in ORC. m H denotes the heat source flow rate. h in , h out represents the input and outlet enthalpy of the heat source. ∆s denotes the entropy change of heat source.