Multi-Objective Optimization and Fluid Selection of Different Cogeneration of Heat and Power Systems Based on Organic Rankine Cycle

: Cogeneration of heat and power systems based on the organic Rankine cycle (ORC-CHP) has been proven to be an effective way to utilize waste heat at medium and low temperatures. In this work, three ORC-CHP (combined heat and power based on organic Rankine cycle) systems are simulated and compared, including the SS (serial system), the CS (the condensation system), and the SS/CS. The multi-objective genetic algorithm (MOGA) is used to optimize the three systems respectively to achieve higher exergy efﬁciency and proﬁt ratio of investment (PRI). The optimal thermal-economic performance is obtained. Twelve organic ﬂuids are adopted to evaluate their performance as working ﬂuids. The calculation results show that SS has the highest exergy efﬁciency, while SS/CS has the best economic performance. Compared with the highest exergy efﬁciency of SS and the best economic performance of SS/CS, CS will be the optimal solution considering these two objective functions. Under the optimal working conditions, SS has the highest thermal efﬁciency because it has the highest net power output. The components with the largest proportion of exergy destruction are the heat exchangers, which also has the highest cost.


Introduction
The energy problem is an important factor affecting the development of human society and economy. How to effectively utilize the limited energy to promote social development in a sustainable way has always been a research hotspot. Nowadays, about 90% of the world's energy use involves generating or consuming heat over a wide temperature range [1]. However, due to the limited technical level, these heat sources have not been fully used, and the utilization rate is not high, especially for some low-temperature heat [2]. For example, in the process of industrial production, at least 50% of the energy is wasted, mainly in the form of low-grade waste heat [3]. Recycling these low-temperature waste heat sources is an important way to improve energy utilization efficiency.
The organic Rankine cycle (ORC) is considered to be a promising technology for low-grade waste heat recovery [4]. The ORC system uses low boiling point organic matter as the working fluid, so it can output power at lower temperature conditions, and it has been widely studied and applied because of its simple design and operation [5]. System configuration, fluid screening, and performance optimization of ORC systems have been carried out by many researchers. Li et al. [6] compared the ORC system with the two-stage series organic Rankine cycle (TSORC) system by recovering heat from the hot flue gas. Various thermodynamic and economic indicators were set, respectively. The results showed that TSORC is superior to ORC in terms of net power output when the flue gas temperature is between 200 and 300 • C, and ORC has more advantages in terms of economic performance because the increased investment of TSORC is lower than the income. Bin et al. [7] used ORC to recover the residual heat of the passenger hybrid vehicle. The maximum efficiency of the ORC system they obtained was 5.4%. The extra power output was 2.02 kW, which improved the fuel economy by 1.0% and 1.2% in two driving cycle tests, respectively. Thermodynamic characteristics of the working fluid largely determine ORC system performance, and the working fluid should also meet environmental and safety standards [8]. This, therefore, indicates that selecting the best working fluid is an important task when considering ORC for waste heat recovery processes [9]. Multi-objective optimization combined with improved grey relational analysis (GRA) was adopted by Xia et al. [10] to screen fluids for a dual-loop organic Rankine cycle (DORC) system. The boiling temperature was set as a working fluid selection criterion for the DORC system to make the system have both high exergy efficiency and excellent economic performance. They found that the most suitable boiling point temperature of working fluid in the high-temperature (HT) loop and low-temperature (LT) loop is 330-363 and 255-305 K, respectively. Cyclohexane/butane is the optimal working fluid when flue gas temperature is 573.15 K. In our previous work, research on working fluid selection [11,12] and experimental tests [13][14][15] has also been carried out.
A number of research efforts have also focused on ORC-based combined heat and power generation (CHP) systems to improve the energy efficiency and economic benefit of ORC systems. The waste heat in the ORC-CHP system can be continued to be used to provide heat to thermal users, so the ability to provide different forms of energy (power and heat) is an important feature of the ORC-CHP system [16]. Lots of works have been published to evaluate the thermo-economic aspect of different ORC-CHP systems. Oyewunmi et al. [16] investigated the adoption of mixed working fluids in an ORC-CHP system. Three typical flue gas waste-heat streams were considered in this work, where the flue-gas streams have temperatures of 150, 250, and 330 • C, with mass flow rates of 30, 120, and 560 kg/s, respectively. The results indicated that pure fluids (especially pentane) show the best performance at low hot-water supply temperatures, and fluid mixtures deliver higher power outputs and exergy efficiencies at higher temperatures. However, the selection of the best working fluid was based on the exergy efficiency of the system, and the economic performance was not considered. Eyerer et al. [17] performed an advanced ORC architecture for geothermal combined heat and power generation, and a novel ORC-CHP test rig with regenerative preheating was constructed. The results showed that the new ORC-CHP system has a wide operating range with a minimum load of 15.3% of the nominal power output. When the system has both power output and heat output, the thermodynamic efficiency and economic performance of the system are always not optimal at the same time. Therefore, a series of optimizations are needed to make the system have both good thermal and economic performance. Zhu et al. [18] analyzed a biomass-fired ORC-CHP system integrated with MEA-based CO 2 capture. With a thermo-economic evaluation, the optimal boiler hot water temperature with a specific evaporation temperature as well as condensation temperature was obtained. At the same time, a thermo-economic comparison of eleven candidate working fluids concerning CO 2 capture was conducted. They concluded that for R245fa, the optimal boiler hot water temperature is obtained as 433 K with the optimized evaporation temperature of 390 K, and cyclopentane shows excellent performance in thermodynamics, followed by R141b, R113, R123, and Pentane. Wang et al. [19] developed a novel optimization algorithm for thermal economic optimization of an ORC-based micro-cogeneration system. R1234 and R141b were selected as candidate working fluids. The optimization results showed that the system performance with R141b as the working medium was better than that with R123 as the working medium. Under the optimal condition, the thermal efficiency and the capital cost rate were found to be 77.70% and 0.363 $/h, respectively. A trigeneration hybrid system consisting of a supercritical water oxidation technology (SCWO) process and an ORC system for power, freshwater, and heat was proposed in our previous work [20]. The energy, exergy, and economic analyses were conducted to assess the performance of the system. The results showed that by treating 5020.0 kg/h of sewage, 286.5 kW of power and 1081 kW of heat could be generated and the electricity efficiency, thermal efficiency, energy efficiency, and exergy efficiency could reach up to 9.56%, 44.93%, 53.47%, and 22.45%, respectively.
From the above literature review, it is perspicuous that for ORC-CHP systems, several system configurations have been proposed for different heat sources and purposes. However, most researchers focused their studies on parameter optimization and thermal economic analysis of a specific ORC-CHP system, and relatively few studies have been performed on the comparison of different ORC-CHP systems. Van et al. [21] studied four ORC-CHP systems, and the economic and thermodynamic optimization of the four systems were carried out, respectively. They found that tandem systems were suitable for low-temperature heating systems and that HB4 configurations were more suitable for high-temperature heating systems with low heat requirements. Finally, it was concluded that net present value (NPV) was more suitable to be optimized than power based on the purpose of investing in geothermal projects. However, it should be noted that the net power output and net present value are optimized and compared respectively, which means that the economic and thermodynamic constraints of each system are not clearly expressed. In this work, three ORC-CHP systems were simulated and compared, including the SS (serial system), the CS (the condensation system), and the SS/CS. The three systems provided domestic hot water as well as electricity, which means there will be a trade-off between thermal and economic performance. Based on the above analysis, the thermodynamic model and economic model of the three systems were established, and then the multi-objective genetic algorithm was adopted to optimize the three systems respectively to achieve higher exergy efficiency and profit ratio of investment (PRI) under different heat sources. The optimal thermal economic performance was obtained. Furthermore, twelve organic fluids were adopted to evaluate their performance as working fluids in the three ORC-CHP systems. Compared with the work already completed in the above studies, this study adopted the multi-objective genetic algorithm as the optimization method, and focused on the presentation of the system's economic performance and thermodynamic efficiency and the variation relationship. The selected objective functions' mutual restriction relationship was presented more visually through the use of graphics. In addition, compared with the optimization analysis for a specific system or working condition, this study would provide a basis for exploring the optimal system or organic working medium under different working conditions.

Basic Organic Rankine Cycle
The cycle configuration and the corresponding T-s diagram of the simple organic Rankine cycle shown in Figure 1 have been described in detail in the literature (e.g., in [22][23][24]). The organic Rankine cycle uses organic working medium as the working fluid, which means that the organic Rankine cycle can recover heat at lower temperatures than the steam Rankine cycle, and has a simpler structure and more stable operation [25]. The organic working medium is heated in the high-temperature heat exchanger (namely the evaporator and the preheater) to the high-temperature and high-pressure overheating or saturation state (process 5a to 1), and then the gaseous organic working medium enters the expander to do work, and further drives the generator to generate electric energy (process 1 to 2a). Next, the working medium flows through the condenser and pump, and re-enters the high-temperature heat exchanger for the next round of circulation (process 2a to 5a). As can be seen in Figure 1, two parts of the heat can be further utilized in the simple organic Rankine cycle: (a) the residual heat of the hot gas after the working fluid is heated, and (b) the heat released by the organic working fluid in the condenser. Therefore, three ORC-CHP systems were simulated to make further use of these two parts of heat, namely the serial system (SS) [26], the condensation system (CS) [18], and the serial/condensation system, respectively. In addition to generating electricity, the remaining usable heat is provided to the thermal consumer during the operation phase of these three systems.

Serial System (SS)
The cycle configuration of the serial ORC-CHP system (SS) and the T-s diagram are shown in Figure 2. It can be seen that after heating the organic working medium, the heat source still has part of the waste heat, which enters the hot water heater (Hwh) to heat domestic hot water (process 9 to 15). At the same time, a cooling tower is applied to reduce the temperature of the cold source. In this study, the heat source comes from the flue gas discharged from industrial production, and water is adopted as a cold source [27]. organic Rankine cycle: (a) the residual heat of the hot gas after the working fluid is heated, and (b) the heat released by the organic working fluid in the condenser. Therefore, three ORC-CHP systems were simulated to make further use of these two parts of heat, namely the serial system (SS) [26], the condensation system (CS) [18], and the serial/condensation system, respectively. In addition to generating electricity, the remaining usable heat is provided to the thermal consumer during the operation phase of these three systems.

Serial System (SS)
The cycle configuration of the serial ORC-CHP system (SS) and the T-s diagram are shown in Figure 2. It can be seen that after heating the organic working medium, the heat source still has part of the waste heat, which enters the hot water heater (Hwh) to heat domestic hot water (process 9 to 15). At the same time, a cooling tower is applied to reduce the temperature of the cold source. In this study, the heat source comes from the flue gas discharged from industrial production, and water is adopted as a cold source [27]. organic Rankine cycle: (a) the residual heat of the hot gas after the working fluid is heated, and (b) the heat released by the organic working fluid in the condenser. Therefore, three ORC-CHP systems were simulated to make further use of these two parts of heat, namely the serial system (SS) [26], the condensation system (CS) [18], and the serial/condensation system, respectively. In addition to generating electricity, the remaining usable heat is provided to the thermal consumer during the operation phase of these three systems.

Serial System (SS)
The cycle configuration of the serial ORC-CHP system (SS) and the T-s diagram are shown in Figure 2. It can be seen that after heating the organic working medium, the heat source still has part of the waste heat, which enters the hot water heater (Hwh) to heat domestic hot water (process 9 to 15). At the same time, a cooling tower is applied to reduce the temperature of the cold source. In this study, the heat source comes from the flue gas discharged from industrial production, and water is adopted as a cold source [27].

Condensation System (CS)
The cycle configuration of the condensation system (CS) and the corresponding T-s diagram are shown in Figure 3. The condensation system is configured similarly to the basic ORC system, except that the domestic hot water in the condensation system is heated by the heat released by the gaseous organic working substance in the condenser 2.1.4. Serial/Condensation System (SS/CS) The cycle configuration and the corresponding T-s diagram of the serial/condensation system (SS/CS) are shown in Figure 4. Different from the two systems proposed above, the heat supplied to the heat consumer can be divided into two parts: one part is the heat released by the gaseous organic working fluid in the condenser, and the other part comes from the residual heat of the heat source after heating the organic working fluid. The cycle configuration of the condensation system (CS) and the corresponding T-s diagram are shown in Figure 3. The condensation system is configured similarly to the basic ORC system, except that the domestic hot water in the condensation system is heated by the heat released by the gaseous organic working substance in the condenser 2.1.4. Serial/Condensation System (SS/CS) The cycle configuration and the corresponding T-s diagram of the serial/condensation system (SS/CS) are shown in Figure 4. Different from the two systems proposed above, the heat supplied to the heat consumer can be divided into two parts: one part is the heat released by the gaseous organic working fluid in the condenser, and the other part comes from the residual heat of the heat source after heating the organic working fluid.

Condensation System (CS)
The cycle configuration of the condensation system (CS) and the corresponding T-s diagram are shown in Figure 3. The condensation system is configured similarly to the basic ORC system, except that the domestic hot water in the condensation system is heated by the heat released by the gaseous organic working substance in the condenser

Serial/Condensation System (SS/CS)
The cycle configuration and the corresponding T-s diagram of the serial/condensation system (SS/CS) are shown in Figure 4. Different from the two systems proposed above, the heat supplied to the heat consumer can be divided into two parts: one part is the heat released by the gaseous organic working fluid in the condenser, and the other part comes from the residual heat of the heat source after heating the organic working fluid.

Mathematical Methods
The proper mathematical model and calculation should be adopted to further accurately study the performance of the proposed systems. At the same time, it is necessary to make reasonable assumptions about some conditions to reduce the complexity of the model. Before the analysis of the system, the following assumptions were adopted in this study: (1) all the systems and components operate under stable conditions, (2) the pressure drops and heat loss of pipes and heat exchangers (evaporator, preheater, water heater, condenser) are ignored, (3) efficiencies of the pump and the turbine remain stable, and (4) the values of the parameters used in the calculation of the model are provided in Table 1.

Selection of Working Fluid
Choosing the right working fluid is one of the necessary conditions for the efficient operation of ORC systems. Different criteria need to be considered when selecting the working fluid, such as thermal efficiency, exergy efficiency, and system economy [30]. In addition, the thermodynamic properties of the organic working fluid need to be reasonably matched with the temperature of the heat source [31]. The critical temperature of the working fluids should not be much lower than the temperature of the heat source, and the reason for this is that if the system operates near the critical point, small perturbations may cause drastic changes in the state of the fluids, which may make the system performance unstable [32]. Furthermore, working fluids that have a negative impact on the environment should not be adopted according to the Montreal Agreement [33] (e.g., R11, R12, R123, R141b, and others).
Based on the above analysis, 12 organic working fluids were pre-selected, and their thermodynamic properties are shown in Table 2. All the data were obtained from NIST REFPROP 9.1.

Exergy Analysis
The subscripts of the formulas below represent the state points in Figures 1-3. Exergy analysis was performed in this study to evaluate the thermodynamic performance of the three ORC-CHP systems. Exergy, which is the part of the energy that can theoretically be converted into any other form of energy when the system changes reversibly from any state to a state in equilibrium with the given environment, represents the quality of energy contained in the working substance at each state point [34]. The exergy of the working fluid at each state point, i, can be written as follows: where h is the enthalpy of the fluid and s is the entropy. h 0 and s 0 represent enthalpy and entropy values at ambient temperature (298.15 K) and pressure (101.3 kPa), respectively.
The key exergy equation for each component can be described as: where ∑ · E in and ∑ · E out represent exergy rates of entry and exit from each component, respectively.
· I is the exergy destruction of each component. The total exergy destruction of the system is the sum of the exergy loss of all parts. The exergy destruction calculation formulas of the three systems are shown in Table 3. Table 3. Total exergy destruction of each system.

System
Total Exergy Loss W is the power output of the expander. For the pumps, the value of · W is the negative form of the value of the pump power consumption. Taking the serial system, for example, the power generated by the expander and pump power consumption can be calculated as follows: where η exp and η pump are the isentropic efficiency of the expander and pump, respectively.
The net output power of the ORC system, · W net , can be computed as: According to the net output power of the system and the exergy destruction of all components, the overall exergy efficiency of the system can be obtained by: Furthermore, to calculate the net power output of the system, it is necessary to obtain the mass flow rate of the organic working fluid, which can be calculated as: where · m h and c ph are mass flow rate and the specific heat of the hot gas at constant pressure, respectively. T 8 is the temperature of the source at the end of the evaporation process, which can be expressed as: where T Pinch is the pinch point, which is the minimum temperature difference of heat exchangers, and T 6 can be determined by the saturation point of the organic working fluid.

Economic Analysis
Since the income of the three proposed systems includes power income and heating income, the profit ratio of investment (PRI) was selected as the economic evaluation index of the system, which can be defined as: where C tot signifies the initial investment cost required for the system. The cost of expanders, heat exchangers, pumps, and other components, such as cooling towers, were all taken into account. The total initial investment cost of the serial system can be calculated as follows: where CEPCI is the chemical engineering plant cost index considering the effect of time on purchased equipment cost, which is widely adopted for updating the capital costs of process engineering projects [35]. CEPCI 1996 and CEPCI 2019 represent the chemical engineering plant cost index in 1996 and 2019 respectively, and their values are 382 and 652.9, respectively. The cost of the cooling tower, C ct , accounts for 23.48% of the cost of the condenser, and C wf means the cost of the working fluid, accounting for 2% of the total investment [18]. C ORC refers to the total initial investment required for the ORC subsystem, which is described as: The calculation formulas of total initial investment costs of CS and SS/CS are shown in Equations (12) and (13), respectively: NAI means the annual net profit of the ORC-CHP system, which is expressed as: where EN E and EN H represent the annual income of the system from generating electrical power and supplying heat, respectively. The two kinds of income can be calculated as follows: where · m dhw is the mass flow rate of domestic hot water to be heated, and can be calculated according to the law of conservation of energy. For example, in the serial system, the mass flow rate of domestic hot water is defined as: where c pc is the specific heat of the hot water at constant pressure, t op is the operation time of the system, which is assumed as 7500 h, and P e and P dhw are the sales prices of electric power and domestic hot water respectively, which are supposed as 0.13 $/kWh and 3.6 $/t in this study [18]. C man is the cost of ensuring the stable operation and maintenance of the system, and it can be calculated as: where γ M is the ratio of maintenance costs for the total system investment, which is assumed as 1.5%. The module costing technique (MCT) is usually used to calculate the cost of specific components, so as to analyze the economic performance of the systems [36]. The cost of purchasing a component can usually be assessed as follows: where X represents the different components shown in Table 4, and Y means the area of the heat exchanger or the capacity of the expander and pump, respectively. C P is the basic cost of each component, which is obtained when the component runs under rated pressure. However, these components often do not work under rated pressure, which means that pressure correction of the above formulas should be carried out. The corrected cost of pumps and heat exchangers can be expressed as: where C BM is the corrected cost, and F M is the material correction factor. F P is the pressure correction factor, which is defined as: where P is the working pressure of each component. For the expander, the corrected cost is calculated as: The correction coefficients mentioned in Equations (19)- (22) can be found in the reference [11].

Optimization Method
All codes in this study were written in MATLAB R2018a and the thermal properties of the organic working fluids used were obtained by NIST REFPROP 9.1 [37].

Multi-Objectives Optimization
Many optimization problems in engineering practice are multi-objective optimization design problems. In general, multiple objectives are in a state of conflict, and there is no optimal design to make all objectives achieve the optimal states at the same time. The improvement of the performance of one target often comes at the cost of the reduction of the performance of one or more other targets [38]. The multi-objective optimization problem can be described as follows: where x is the variable space to be optimized, k is the number of target functions, and X represents the set of all the solutions that satisfy the constraint.
The Pareto optimal solution of the multi-objective optimization problem is a set of acceptable "not bad" solution sets. The solution is Pareto optimal if no other solution exists to make each objective function better. Compared with single-objective optimization, the final result converges to a specific point, and the result of multi-objective optimization exists in the form of a set of solution sets.

Multi-Objective Genetic Algorithm
The multi-objective genetic algorithm has been widely used in engineering optimization technology [39]. The genetic algorithm (GA) is inspired by natural biology and speeds up convergence through selection, crossover, mutation, and other operations. It is a kind of search algorithm, through continuous search and trial and error, and when the result meets the accuracy requirement or iterates to the maximum algebra, the program will be terminated.
In this work, the temperature and pressure at the expander inlet (i.e., T 1 and P 1 ) are selected as two parameters to be optimized. A two-dimensional array [T 1 , P 1 ] is generated as a chromosome for each individual. Constraint conditions required to generate T 1 and P 1 are as follows: T 1 max = min (T crit − 10), (T gas − 8) ; where (T crit − 10) is set to prevent instability caused by the organic working fluid operating near the critical point [40]. At the same time, the maximum temperature, T 1 , cannot be greater than (T gas − 8) due to the limitation of the temperature difference between pinch points. P sat,T 1 is determined by T 1 to ensure that the organic working medium is saturated or overheated when entering the expander. The lower limits of T 1 and P 1 are also set so that the outlet temperature of the heat source (state point 15 in Figures 2 and 4, and state point 9 in Figure 3) is 15 K higher than the acid dew temperature. The exergy efficiency, η ex , and profit ratio of investment, PRI, are selected as the target functions in this work, which can be expressed as: As an important step in the genetic algorithm, the "selection" process embodies the rule of survival of the fittest. The individuals with higher fitness are selected from all the individuals to produce the next-generation population. As in natural processes, "crossover" involves the swapping of genes on the chromosome of the selected parent. The process of "mutation" causes some genes on a single chromosome to change to avoid obtaining the local optimal solution. Last but not least, the elite retention strategy is adopted to protect the genes of high-performing individuals from being damaged. The complete implementation process of the genetic algorithm is shown in Figure 5. Table 4 shows the values of each parameter of the GA in this work. The exergy efficiency, ex  , and profit ratio of investment, PRI, are selected as the target functions in this work, which can be expressed as: As an important step in the genetic algorithm, the "selection" process embodies the rule of survival of the fittest. The individuals with higher fitness are selected from all the individuals to produce the next-generation population. As in natural processes, "crossover" involves the swapping of genes on the chromosome of the selected parent. The process of "mutation" causes some genes on a single chromosome to change to avoid obtaining the local optimal solution. Last but not least, the elite retention strategy is adopted to protect the genes of high-performing individuals from being damaged. The complete implementation process of the genetic algorithm is shown in Figure 5. Table 4 shows the values of each parameter of the GA in this work.

Optimization Results and Analysis
The optimization results of different organic working fluids adopted in the three ORC-CHP systems under different heat source temperatures are shown in Figures 6-8,

Optimization Results and Analysis
The optimization results of different organic working fluids adopted in the three ORC-CHP systems under different heat source temperatures are shown in Figures 6-8, respectively. Different lines are used to distinguish between SS (dotted line), CS (short line), and SS/CS (solid line). At the same time, the different colored lines in the legend represent different organic working fluids, and this is true for the three figures. T g refers to the initial temperature of the heat source gas, which is 373. 15        As can be seen from the optimization results, the Pareto optimal solution exists in the form of a group of numbers, and the improvement of exergy efficiency is bound to lead to the decline of the economic performance of the system, and vice versa. Meanwhile, the economic benefits and thermal performance of different systems and organic working substances are distributed in different ranges, which are determined by the form of the systems and the thermodynamic properties of the organic working fluids. For example, when the temperature of the heat source is 373.15 K, it can be seen from Figure 6 that the overall exergy efficiency of SS is higher than that of CS, while the PRI value of the two systems is almost within the same range. On the other hand, the economic performance of SS/CS is unimpeachable at this time, but its exergy efficiency is generally low. Therefore, we can judge that when the heat source temperature is relatively low, SS is the best choice, SS/CS is the second, and CS should not be recommended under this condition. Then, when the temperature reaches 398.15 K, it can be seen from Figure 7 that SS and CS are neck and neck in the evaluation of the two performance indexes, and even the economic performance of CS is slightly higher than that of SS. Therefore, through the analysis of the above case, it is essential to choose the appropriate ORC-CHP system and organic working fluid according to the actual engineering requirements. The selection of system form and organic working fluid should be carried out according to specific conditions and requirements. The following examples illustrate how to select the corresponding system and working medium according to the given conditions and requirements: Assuming the temperature of the heat source is 373.15 K and the exergy efficiency of the system is required to be above 35%, it can be seen from Figure 6 that only SS can meet the above conditions, and R245ca can be considered as the best organic working fluid. The maximum profit ratio of investment the system can achieve at this point is 0.219. In the same way, the best system configurations and organic working fluids can be selected to meet the specific requirements, within the range of effects that can be achieved.
Furthermore, it can be seen from Figures 6-8 that there are some curves with very low slope, such as the optimization results of SS/CS using R1234ze when the heat source temperature is 398.15 K. In this case, when the system has higher exergy efficiency, the As can be seen from the optimization results, the Pareto optimal solution exists in the form of a group of numbers, and the improvement of exergy efficiency is bound to lead to the decline of the economic performance of the system, and vice versa. Meanwhile, the economic benefits and thermal performance of different systems and organic working substances are distributed in different ranges, which are determined by the form of the systems and the thermodynamic properties of the organic working fluids. For example, when the temperature of the heat source is 373.15 K, it can be seen from Figure 6 that the overall exergy efficiency of SS is higher than that of CS, while the PRI value of the two systems is almost within the same range. On the other hand, the economic performance of SS/CS is unimpeachable at this time, but its exergy efficiency is generally low. Therefore, we can judge that when the heat source temperature is relatively low, SS is the best choice, SS/CS is the second, and CS should not be recommended under this condition. Then, when the temperature reaches 398.15 K, it can be seen from Figure 7 that SS and CS are neck and neck in the evaluation of the two performance indexes, and even the economic performance of CS is slightly higher than that of SS. Therefore, through the analysis of the above case, it is essential to choose the appropriate ORC-CHP system and organic working fluid according to the actual engineering requirements. The selection of system form and organic working fluid should be carried out according to specific conditions and requirements. The following examples illustrate how to select the corresponding system and working medium according to the given conditions and requirements: Assuming the temperature of the heat source is 373.15 K and the exergy efficiency of the system is required to be above 35%, it can be seen from Figure 6 that only SS can meet the above conditions, and R245ca can be considered as the best organic working fluid. The maximum profit ratio of investment the system can achieve at this point is 0.219. In the same way, the best system configurations and organic working fluids can be selected to meet the specific requirements, within the range of effects that can be achieved.
Furthermore, it can be seen from Figures 6-8 that there are some curves with very low slope, such as the optimization results of SS/CS using R1234ze when the heat source temperature is 398.15 K. In this case, when the system has higher exergy efficiency, the cost of reduced economic performance is within the acceptable range. Therefore, the system should be operated under the condition of high exergy efficiency as far as possible. On the contrary, for those curves with high slope, the sharp decline in economic performance caused by higher exergy efficiency is worth considering.
Taking the heat source temperature of 398.15 K as an example [41,42], the performance of the three systems under optimal working conditions should be analyzed in detail. To facilitate the further analysis of the optimal working conditions of each system, it is necessary to find an optimal solution of each system when the heat source temperature is 398.15 K. The optimization results of SS are more clearly shown in Figure 9. When the economic performance of the CHP system is the primary goal, point A is the optimal solution of the system, but at this time, the exergy efficiency is also the lowest. On the contrary, when the exergy efficiency of the system is regarded as the greatest goal, there is no doubt that point B will be the best solution for the system, but the economic performance of point B is not very ideal, which cannot be ignored. Therefore, in this study, point O is used to find the best solution of the system [19]. It can be seen that point O is the point with the best economic performance and the highest exergy efficiency. However, since the performance of point O is beyond the reach of the actual system operation, point C, that is, the point closest to point O in the curve, is taken as the optimal solution of the system. At the same time, since point C falls on the curve representing R245ca, R245ca can be considered as the most competitive working medium applied to SS when the heat source temperature is 398.15 K. Similarly, the optimal solutions of CS (F) and SS/CS (I) are shown in Figures 10 and 11, respectively. It can be seen that when the heat source temperature is 398.15 K, the optimal working medium of SS is R245ca and R1233zd(E) (38.67%, 0.338), and that of CS is R245ca (32.53%, 0.496). For SS/CS, it is obvious that R1233zd(E) (26.85%, 0.718) will be the optimal choice. cost of reduced economic performance is within the acceptable range. Therefore, the system should be operated under the condition of high exergy efficiency as far as possible.
On the contrary, for those curves with high slope, the sharp decline in economic performance caused by higher exergy efficiency is worth considering. Taking the heat source temperature of 398.15 K as an example [41,42], the performance of the three systems under optimal working conditions should be analyzed in detail. To facilitate the further analysis of the optimal working conditions of each system, it is necessary to find an optimal solution of each system when the heat source temperature is 398.15 K. The optimization results of SS are more clearly shown in Figure 9. When the economic performance of the CHP system is the primary goal, point A is the optimal solution of the system, but at this time, the exergy efficiency is also the lowest. On the contrary, when the exergy efficiency of the system is regarded as the greatest goal, there is no doubt that point B will be the best solution for the system, but the economic performance of point B is not very ideal, which cannot be ignored. Therefore, in this study, point O is used to find the best solution of the system [19]. It can be seen that point O is the point with the best economic performance and the highest exergy efficiency. However, since the performance of point O is beyond the reach of the actual system operation, point C, that is, the point closest to point O in the curve, is taken as the optimal solution of the system. At the same time, since point C falls on the curve representing R245ca, R245ca can be considered as the most competitive working medium applied to SS when the heat source temperature is 398.15 K. Similarly, the optimal solutions of CS (F) and SS/CS (I) are shown in Figures 10 and 11, respectively. It can be seen that when the heat source temperature is 398.15 K, the optimal working medium of SS is R245ca and R1233zd(E) (38.67%, 0.338), and that of CS is R245ca (32.53%, 0.496). For SS/CS, it is obvious that R1233zd(E) (26.85%, 0.718) will be the optimal choice.      The above method can be used to analyze the working conditions with the heat source temperatures of 373.15 and 423.15 K, and the results show that when the heat source temperature is 373.15 K, the optimal working mediums corresponding to SS, CS, and SS/CS are, respectively: R245ca (38.99%, 0.178), R245fa (20.38%, 0.242), and R365mfc (15.80%, 0.468). When the heat source temperature reaches 423.15 K, the optimal working mediums corresponding to the three systems are R1233zd(E) (42.68%, 0.378), R245ca (41.22%, 0.671), and R245fa (35.81%, 0.881), respectively. Therefore, it is feasible to select the system under different heat sources according to the above optimum conditions and combined with the actual engineering requirements. Under the condition of a heat source temperature of 373.15 K, SS has the highest exergy efficiency, while SS/CS has the best economic performance, which is similar to the results when the heat source temperature is 398.15 K. When the temperature of the heat source reaches 423.15 K, compared with the highest exergy efficiency of SS and the best economic performance of SS/CS, CS will be the optimal solution considering these two objective functions.
According to the above calculation of optimal working conditions, it can be learned that R245ca and R1233zd(E) are the working mediums that can satisfy the most system and working conditions, followed by R245fa. R245fa has been considered the optimal working medium in some studies that have been carried out [43][44][45]. At the same time, it has also been mentioned in some studies that R245ca is the optimal working medium for ORC systems [46,47]. Therefore, for ORC-CHP systems, the optimal working medium is not the same, which is mainly determined by the system form and the research or optimization objectives. It is worth noting that the high GWP (Global Warming Potential) values of R245fa and R245ca will increase the potential risks for climate change mitigation. However, as a transitional working medium, R245ca's excellent performance in this study is still of important reference significance. Meanwhile, it is necessary to study some environmental protection working fluids such as R1233zd(E) from the perspective of the environment [48].

Thermodynamic Cycle Efficiencies of the Optimal Solutions
The thermal efficiency of the systems can be defined as: where . Q presents the total heat that the working medium obtains from the heat source. For CS, it is given by: For SS and SS/CS, the total heat absorbed by the two systems from the heat source can be expressed as: Therefore, according to the above formulas and the optimal working conditions obtained, the thermal efficiency of each system can be obtained under the optimal operating conditions. Under different heat source temperatures, with the optimal combination of exergy efficiency and PRI value, the thermal efficiency of each system is shown in Figure 12. It can be seen that SS always has higher thermal efficiency than the other two systems when all three systems operate under the optimal conditions. The main reason for this is that the SS has more net power output under the best operating conditions. At the same time, it can also be seen that the thermal efficiency of the three systems increases with the increase of the temperature of the heat source. In other words, higher reservoir temperatures lead to higher energy efficiency. However, it should be noted that such results are only found in this study; that is, within the temperature range of the heat source selected in this study, the thermal efficiency of the three systems proposed in the optimal operating conditions will increase with the temperature of the heat source. Energies 2021, 14, x FOR PEER REVIEW 17 of 22 Figure 12. Thermal efficiency of the systems under optimal operating conditions.

Exergy Analysis of the Optimal Solutions
The optimal conditions of the three ORC-CHP systems have been obtained when the heat source temperature is 398.15 K. Exergy analysis on the three optimal working conditions is necessary to carry out. The contribution of each component of the three ORC-CHP systems to the rate of overall exergy destruction at the optimal operating conditions is shown in Figure 13. It can be seen that in the SS (Figure 13a), the exergy destruction in the water heater (Hwh) (38.7%) accounted for the largest proportion of the total exergy destruction of the system, while for CS (Figure 13b), the exergy destruction was mostly caused by the evaporator (40.5%) and condenser (40.1%) under the optimal working condition. In addition, when SS/CS (Figure 13c) was adopted, the components with the largest proportion of exergy destruction were the three heat exchangers, namely: water heater (33.2%), condenser (23.1%), and evaporator (21.5%). From the above exergy analysis, it can be concluded that in the three ORC-CHP systems, the heat exchangers are always the components with the highest contribution rate to the overall exergy destruction. Furthermore, exergy loss caused by the pumps accounted for the smallest proportion of exergy loss in all three systems.

Exergy Analysis of the Optimal Solutions
The optimal conditions of the three ORC-CHP systems have been obtained when the heat source temperature is 398.15 K. Exergy analysis on the three optimal working conditions is necessary to carry out. The contribution of each component of the three ORC-CHP systems to the rate of overall exergy destruction at the optimal operating conditions is shown in Figure 13. It can be seen that in the SS (Figure 13a), the exergy destruction in the water heater (Hwh) (38.7%) accounted for the largest proportion of the total exergy destruction of the system, while for CS (Figure 13b), the exergy destruction was mostly caused by the evaporator (40.5%) and condenser (40.1%) under the optimal working condition. In addition, when SS/CS (Figure 13c) was adopted, the components with the largest proportion of exergy destruction were the three heat exchangers, namely: water heater (33.2%), condenser (23.1%), and evaporator (21.5%). From the above exergy analysis, it can be concluded that in the three ORC-CHP systems, the heat exchangers are always the components with the highest contribution rate to the overall exergy destruction. Furthermore, exergy loss caused by the pumps accounted for the smallest proportion of exergy loss in all three systems.

Economic Analysis of the Optimal Solutions
According to Equations (10)- (22), the economic performance of the three systems under optimal working conditions was evaluated. The component costs of the three systems under the optimal operating conditions described above are shown in Figures 14-16, respectively. It can be seen that the components with the highest cost are heat exchangers, which is similar to the result of the exergy analysis above. For example, the three components with the highest equipment cost in SS are the evaporator (34.2%), water heater (24.2%), and condenser (15.1%), followed by the expander (9.7%). Therefore, if the cost of the system needs to be further optimized in the future, the heat exchanger is undoubtedly the key point of optimization.

Economic Analysis of the Optimal Solutions
According to Equations (10)- (22), the economic performance of the three systems under optimal working conditions was evaluated. The component costs of the three systems under the optimal operating conditions described above are shown in Figures 14-16, respectively. It can be seen that the components with the highest cost are heat exchangers, which is similar to the result of the exergy analysis above. For example, the three components with the highest equipment cost in SS are the evaporator (34.2%), water heater (24.2%), and condenser (15.1%), followed by the expander (9.7%). Therefore, if the cost of the system needs to be further optimized in the future, the heat exchanger is undoubtedly the key point of optimization.

Economic Analysis of the Optimal Solutions
According to Equations (10)- (22), the economic performance of the three systems under optimal working conditions was evaluated. The component costs of the three systems under the optimal operating conditions described above are shown in Figures  14-16, respectively. It can be seen that the components with the highest cost are heat exchangers, which is similar to the result of the exergy analysis above. For example, the three components with the highest equipment cost in SS are the evaporator (34.2%), water heater (24.2%), and condenser (15.1%), followed by the expander (9.7%). Therefore, if the cost of the system needs to be further optimized in the future, the heat exchanger is undoubtedly the key point of optimization. Cost ratio/% Figure 14. Component costs of SS. Figure 14. Component costs of SS. Cost ratio/% Figure 16. Component costs of SS/CS.
The electricity and hot water revenues of the three systems under optimal conditions are shown in Figure 17. It can be found that when the three systems are operating in a state of excellent thermal and economic performance, the supply of hot water is the major approach for the systems to obtain revenues. It is important to note that this does not apply in every case, as the results are mainly determined by the local price of hot water and electricity.   Cost ratio/% Figure 16. Component costs of SS/CS.
The electricity and hot water revenues of the three systems under optimal conditions are shown in Figure 17. It can be found that when the three systems are operating in a state of excellent thermal and economic performance, the supply of hot water is the major approach for the systems to obtain revenues. It is important to note that this does not apply in every case, as the results are mainly determined by the local price of hot water and electricity.  The electricity and hot water revenues of the three systems under optimal conditions are shown in Figure 17. It can be found that when the three systems are operating in a state of excellent thermal and economic performance, the supply of hot water is the major approach for the systems to obtain revenues. It is important to note that this does not apply in every case, as the results are mainly determined by the local price of hot water and electricity. Cost ratio/% Figure 16. Component costs of SS/CS.
The electricity and hot water revenues of the three systems under optimal conditions are shown in Figure 17. It can be found that when the three systems are operating in a state of excellent thermal and economic performance, the supply of hot water is the major approach for the systems to obtain revenues. It is important to note that this does not apply in every case, as the results are mainly determined by the local price of hot water and electricity. Figure 17. The electricity and hot water revenues of the three systems. Figure 17. The electricity and hot water revenues of the three systems.

Conclusions and Future Work
In the present study, firstly, the thermodynamic and economic models of the three systems (SS, CS, and SS/CS) were established. Then, under the conditions that the heat source temperatures are 373.15, 398.15, and 423.15 K, the multi-objective genetic algorithm was used to optimize the three systems respectively, to obtain higher exergy efficiency and PRI. Twelve organic fluids were adopted to evaluate their performance as working fluids in ORC-CHP systems. Finally, the optimal working conditions and corresponding thermal efficiency of the three systems were obtained, and exergy analysis and economic analysis were carried out for the three systems under these optimal working conditions. The main conclusions of the study are summarized as follows: (1) When the heat source temperature was 373.15 K, the optimal working mediums corresponding to SS, CS, and SS/CS were, respectively: R245ca (38.99%, 0.178), R245fa (20.38%, 0.242), and R365mfc (15.80%, 0.468). Under the heat source temperature of 398.15 K, the optimal working mediums of SS were R245ca and R1233zd(E) (38.67%, 0.338), and that of CS was R245ca (32.53%, 0.496). For SS/CS, R1233zd(E) (26.85%, 0.718) will be the optimal choice. When the heat source temperature reached 423.15 K, the optimal working mediums corresponding to the three systems were R245ca (42.68%, 0.378), R245ca (41.22%, 0.671), and R245fa (35.81%, 0.881), respectively.
(2) Under the condition of a heat source temperature of 373.15 K, SS had the highest exergy efficiency, while SS/CS had the best economic performance, which are similar to the results when the heat source temperature was 398.15 K. When the temperature of the heat source reaches 423.15 K, compared with the highest exergy efficiency of SS and the best economic performance of SS/CS, CS will be the optimal solution considering these two objective functions.
(3) SS always had higher thermal efficiency than the other two systems when all three systems operate under the optimal conditions. The main reason for this is that the SS has more net power output under the best operating conditions. At the same time, the thermal efficiency of the three systems increases with the increase of the temperature of the heat source.
(4) The exergy analysis result under the three optimal working conditions showed that in the SS, the exergy destruction in the water heater (Hwh) (38.7%) accounted for the largest proportion of the total exergy destruction of the system, while for CS, the exergy destruction was mostly caused by the evaporator (40.5%) and condenser (40.1%) under the optimal working condition. In addition, when SS/CS was adopted, the components with the largest proportion of exergy destruction were the three heat exchangers, namely: water heater (33.2%), condenser (23.1%), and evaporator (21.5%).
(5) The components with the highest cost are heat exchangers, which is similar to the result of exergy analysis.
(6) When the three systems are operating in a state of excellent thermal and economic performance, the supply of hot water is the major approach for the systems to obtain revenues.
The results of this study can provide references for ORC-CHP system selection and organic working fluid selection under specific conditions. To cover a wider range of working conditions, it is necessary to conduct further studies on higher temperature conditions and more system forms.