Parametrized Analysis and Multi-Objective Optimization of Supercritical CO 2 (S-CO 2 ) Power Cycles Coupled with Parabolic Trough Collectors

: Supercritical CO 2 (S-CO 2 ) Brayton cycles have become an e ﬀ ective way in utilizing solar energy, considering their advantages. The presented research discusses a parametrized analysis and systematic comparison of three S-CO 2 power cycles coupled with parabolic trough collectors. The e ﬀ ects of turbine inlet temperature and pressure, compressor inlet temperature, and pressure on speciﬁc work, overall e ﬃ ciency, and cost of core equipment of di ﬀ erent S-CO 2 Brayton cycles are discussed. Then, the two performance criteria, including speciﬁc work and cost of core equipment, are compared, simultaneously, between di ﬀ erent S-CO 2 cycle layouts after gaining the Pareto sets from multi-objective optimizations using genetic algorithm. The results suggest that the simple recuperation cycle layout shows more excellent performance than the intercooling cycle layout and the recompression cycle layout in terms of cost, while the advantage in speciﬁc work of the intercooling cycle layout and the recompression cycle layout is not obvious. This study can be useful in selecting cycle layout using solar energy by the parabolic trough solar collector when there are requirements for the speciﬁc work and the cost of core equipment. Moreover, high turbine inlet temperature is recommended for the S-CO 2 Brayton cycle using solar energy.


Introduction
Solar energy has played an important role in sustainability. The main routes to harvest solar energy are by concentrating solar thermal power (CSP) and solar photovoltaics (PV) [1]. There is a lot of research on PV, including economics and environmental performance [2,3]. Compared with the PV generation, CSP has certain advantages in the power generation scale. At present, multiple kinds of CSP technology, including the parabolic trough collector (PTC) [4], the linear Fresnel reflector [5], and the solar power tower (SPT) [6] have been reported. Among them, parabolic trough (PT) can be considered the most mature technology for CSP plants due to projects conducted by the United States [7] and Spain [8].
At present, the Rankine cycle is the most common power generation cycle in CSP technology, which uses parabolic trough collectors, while organic working fluids and water are commonly adopted. However, these forms of power generation face barriers, including the use of working fluids that are easy to aggravate global warming, and high investment costs [9]. The performance of the Organic Rankine Cycle (ORC) is limited by chemical degradation and auto-ignition at high temperatures [10]. Thus, ORC is a better power cycle choice only when the temperature of the heat source is lower than 1.
Three S-CO 2 power cycles coupled with PTC (namely, simple regeneration, intercooling and re-compression) are established and their parameters are analyzed for preliminary comparison.

2.
A multi-objective optimization method using genetic algorithm (GA) is presented to compare the performance, including thermodynamic performance and preliminary economic cost of cycles. The advantages and disadvantages of the three cycle kinds are thoroughly compared. 3.
The sensitivity analysis of a typical optimization result of a simple regeneration cycle is carried out; the influence of different parameters on system performance is obtained.
This work can provide a solid foundation for the fundamental design and construction of S-CO 2 Brayton cycles. We hope that this research will ignite interest in the application of S-CO 2 cycles coupled with PTC and the performance comparison of cycles considering thermodynamics and economics.

Description of the Research Systems
The S-CO 2 Brayton cycle based on solar energy has the same structure as that of a conventional heat source. Core components of the power cycle include the turbine, generator, heat regenerator, cooler, and compressor. The thermal cycles studied in this paper include three types, simple recuperated cycle, intercooling cycle, and the recompression cycle. Solar receivers are adopted to absorb solar radiation heat to warm S-CO 2 .
The simple recuperated cycle is shown in Figure 1, which has been presented by Feher [28]. The working fluid S-CO 2 is cooled before compressing to an applicable pressure for the compressor. In the recuperator, cold S-CO 2 absorbs heat from hot S-CO 2 exiting from the turbine outlet. Moreover, before expanding in the turbine, the pre-heated S-CO 2 further absorbs heat from the solar receiver.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 20 cooler, and compressor. The thermal cycles studied in this paper include three types, simple recuperated cycle, intercooling cycle, and the recompression cycle. Solar receivers are adopted to absorb solar radiation heat to warm S-CO2. The simple recuperated cycle is shown in Figure 1, which has been presented by Feher [28]. The working fluid S-CO2 is cooled before compressing to an applicable pressure for the compressor. In the recuperator, cold S-CO2 absorbs heat from hot S-CO2 exiting from the turbine outlet. Moreover, before expanding in the turbine, the pre-heated S-CO2 further absorbs heat from the solar receiver.  Figure 2 shows the layout of the intercooling cycle, which has been previously presented by Ahn et al. [29]. Compared to the simple recuperated cycle, the intercooling cycle adopts another cooler to help reduce the compression work. It can be accomplished since the temperature of the compressed S-CO2 is further reduced. This arrangement can also generate a larger pressure ratio. Moreover, other components in the intercooling cycle have high similarity with those in the above cycle.  Figure 3 presents the recompression cycle, which has also been adopted by Moisseytsev et al. [30]. This kind of cycle has high similarity with the simple recuperated cycle. Nevertheless, the  Figure 2 shows the layout of the intercooling cycle, which has been previously presented by Ahn et al. [29]. Compared to the simple recuperated cycle, the intercooling cycle adopts another cooler to help reduce the compression work. It can be accomplished since the temperature of the compressed S-CO 2 is further reduced. This arrangement can also generate a larger pressure ratio. Moreover, other components in the intercooling cycle have high similarity with those in the above cycle.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 20 cooler, and compressor. The thermal cycles studied in this paper include three types, simple recuperated cycle, intercooling cycle, and the recompression cycle. Solar receivers are adopted to absorb solar radiation heat to warm S-CO2. The simple recuperated cycle is shown in Figure 1, which has been presented by Feher [28]. The working fluid S-CO2 is cooled before compressing to an applicable pressure for the compressor. In the recuperator, cold S-CO2 absorbs heat from hot S-CO2 exiting from the turbine outlet. Moreover, before expanding in the turbine, the pre-heated S-CO2 further absorbs heat from the solar receiver.  Figure 2 shows the layout of the intercooling cycle, which has been previously presented by Ahn et al. [29]. Compared to the simple recuperated cycle, the intercooling cycle adopts another cooler to help reduce the compression work. It can be accomplished since the temperature of the compressed S-CO2 is further reduced. This arrangement can also generate a larger pressure ratio. Moreover, other components in the intercooling cycle have high similarity with those in the above cycle.  Figure 3 presents the recompression cycle, which has also been adopted by Moisseytsev et al. [30]. This kind of cycle has high similarity with the simple recuperated cycle. Nevertheless, the components are of different numbers and arrangements. Before the cooler, there are two kinds of  Figure 3 presents the recompression cycle, which has also been adopted by Moisseytsev et al. [30]. This kind of cycle has high similarity with the simple recuperated cycle. Nevertheless, the components are of different numbers and arrangements. Before the cooler, there are two kinds of recuperators, which are used in high-temperature and low-temperature, respectively. After regeneration, S-CO 2 splits into two parts. Part of the fluid flows through the cooler and then enters the main-compressor while the other part flows into the re-compressor. The fluid from the main-compressor outlet flows into the low-temperature recuperator and the high-temperature recuperator. The fluid from the outlet of the re-compressor flows only into the high-temperature recuperator. In addition, other components are similar with the above two configurations.
Appl. Sci. 2019, 9, x FOR PEER REVIEW  4 of 20 recuperators, which are used in high-temperature and low-temperature, respectively. After regeneration, S-CO2 splits into two parts. Part of the fluid flows through the cooler and then enters the main-compressor while the other part flows into the re-compressor. The fluid from the maincompressor outlet flows into the low-temperature recuperator and the high-temperature recuperator. The fluid from the outlet of the re-compressor flows only into the high-temperature recuperator. In addition, other components are similar with the above two configurations. It should be noticed that due to the inexistence of a thermal energy storage device, the systems would mostly work under off-design conditions. As shown in previous studies, the performance of systems change dramatically [31,32]. Therefore, we can add a thermal energy storage device to the systems to ensure that the systems operate in a steady state. However, the discussion of a thermal energy storage device is not the focus of this paper, so the relevant description is omitted.

Modeling of the System
The purpose is to compare the thermodynamic characteristics and preliminary performance of the three different power cycles after respective optimization. The systems are modeled in a steady state with the software MATLAB; the properties of CO2 are obtained from the software NIST REFPROP. In the meantime, the modeling can be adopted to analyze preliminary performance and realize the optimizations.
The modeling process is based on typical simplified assumptions, which are commonly utilized for a power system evaluation problem. They are listed as follows: 1. The cycle is considered steady during calculation.
2. The effects of kinetic and gravitational terms are neglected in the calculations. 3. The pressure drops and heat loss are neglected. Figure 4 shows the energy balance at the cross section of the parabolic trough solar collector used in this study. The outer blue ring represents the glass envelope while the inner yellow ring represents the metal heat-absorbing tube in the figure. The equations of energy can be written as [33]: It should be noticed that due to the inexistence of a thermal energy storage device, the systems would mostly work under off-design conditions. As shown in previous studies, the performance of systems change dramatically [31,32]. Therefore, we can add a thermal energy storage device to the systems to ensure that the systems operate in a steady state. However, the discussion of a thermal energy storage device is not the focus of this paper, so the relevant description is omitted.

Modeling of the System
The purpose is to compare the thermodynamic characteristics and preliminary performance of the three different power cycles after respective optimization. The systems are modeled in a steady state with the software MATLAB; the properties of CO 2 are obtained from the software NIST REFPROP. In the meantime, the modeling can be adopted to analyze preliminary performance and realize the optimizations.
The modeling process is based on typical simplified assumptions, which are commonly utilized for a power system evaluation problem. They are listed as follows: 1.
The cycle is considered steady during calculation.

2.
The effects of kinetic and gravitational terms are neglected in the calculations.

3.
The pressure drops and heat loss are neglected. Figure 4 shows the energy balance at the cross section of the parabolic trough solar collector used in this study. The outer blue ring represents the glass envelope while the inner yellow ring represents the metal heat-absorbing tube in the figure. The equations of energy can be written as [33]:

Solar Collector
where Q is quantity of heat, subscript fluid indicates the working fluid, subscript mo-mi the transmitting process from the outer surface of the metal absorber to the inner surface, subscript mo indicates the outer surface of absorber. Subscript conv,mo-gi indicates the convection heat transfer process from absorber to the inner surface of glass envelope, subscript rad,mo-gi indicates the radiation heat transfer process between absorber to the inner surface of glass envelope, subscript mo-loss indicates the loss from absorber to the environment through the bracket, which is at the end of absorber. Subscript gi-go indicates the conducting process from the inner surface to the outer surface of the glass cover, subscript conv,go-atm indicates the loss converted from the outer wall of glass cover, subscript rad,go-atm indicates the loss caused by radiation from the outer wall of glass cover, and subscript go indicates the outer wall of the glass cover.
gi-go conv,mo-gi rad,mo-gi conv,go-atm rad,go-atm gi-go go where Q is quantity of heat, subscript fluid indicates the working fluid, subscript mo-mi the transmitting process from the outer surface of the metal absorber to the inner surface, subscript mo indicates the outer surface of absorber. Subscript conv,mo-gi indicates the convection heat transfer process from absorber to the inner surface of glass envelope, subscript rad,mo-gi indicates the radiation heat transfer process between absorber to the inner surface of glass envelope, subscript moloss indicates the loss from absorber to the environment through the bracket, which is at the end of absorber. Subscript gi-go indicates the conducting process from the inner surface to the outer surface of the glass cover, subscript conv,go-atm indicates the loss converted from the outer wall of glass cover, subscript rad,go-atm indicates the loss caused by radiation from the outer wall of glass cover, and subscript go indicates the outer wall of the glass cover. The solar radiation absorbed by the metal absorber is not only related to the properties of the absorber, but also the properties of the trough reflector and glass envelope. It can be calculated by: Moreover, the solar radiation absorber by the glass envelope can be calculated by: where G is solar radiation intensity, Dmo is the outer diameter of absorber, l is the length of the collector, CR is the concentration ratio of the collector, ρ is the reflectivity of trough reflector, τg is the transmittance of glass envelope, αm is the absorptivity of absorber, Dgo is the outer diameter of glass envelope, αg is the absorptivity of the glass envelope. The heat transfer between the fluid and the inner surface of the absorber can be regarded as a forced convective heat transfer, which can be calculated by: The solar radiation absorbed by the metal absorber is not only related to the properties of the absorber, but also the properties of the trough reflector and glass envelope. It can be calculated by: Moreover, the solar radiation absorber by the glass envelope can be calculated by: where G is solar radiation intensity, D mo is the outer diameter of absorber, l is the length of the collector, CR is the concentration ratio of the collector, ρ is the reflectivity of trough reflector, τ g is the Appl. Sci. 2020, 10, 3123 6 of 20 transmittance of glass envelope, α m is the absorptivity of absorber, D go is the outer diameter of glass envelope, α g is the absorptivity of the glass envelope. The heat transfer between the fluid and the inner surface of the absorber can be regarded as a forced convective heat transfer, which can be calculated by: where D mi is the inner diameter of absorber, hc fluid is the convective heat transfer coefficient of fluid, T mi is the average temperature of the inner surface of the absorber, T fluid is the average temperature of the working fluid, Nu fluid is the Nusselt number of the working fluid, k fluid is the thermal conductivity of the working fluid. The Nusselt number is calculated by the Gnielinski correlation [34]: 0.11 (9) where Re fluid is the Reynolds number of working fluid, Pr fluid is the Prandtl number of fluid, and Pr mi is the Prandtl number of the inner surface of the absorber, f is the friction coefficient of fluid: The area between the glass envelope and the absorber keeps a vacuum to reduce heat loss. Therefore, the convective heat transfer between the absorber and glass cover caused by free movement of the molecules should be little. The heat transfer between the absorber and glass envelope can be given as: (12) where D mo is the outer diameter of absorber, T mo is the average temperature of outer surface of absorber, T gi is the average temperature of the inner surface of the glass envelope, σ is Boltzmann constant, which is 5.67 × 10 −8 W/(m 2 ·K 4 ), ε mo is the emissivity of outer surface of absorber, ε gi is the emissivity of inner surface of glass envelope, and hc mo-gi is the convection heat transfer coefficient of the air in the interlayer, which can be calculate by [35]: where k std is the thermal conductivity of air under standard conditions, b is the interaction coefficient, a is the thermal adaptation coefficient, which is usually taken as 1, γ is the heat capacity ratio of air which is 1.4, λ is the mean free path of molecular, T mo-gi is the average temperature of absorber and glass cover, P mo-gi is the pressure of the area approximated as 0.0001 mmHg, δ is molecular diameter of air taken as 3.53 × 10 −8 cm.
The heat transfer from the outer surface of the glass envelope to the environment can be calculated as: Q conv,go-atm = π · D go · l · hc air · (T go − T atm ) (16) where D go is the outer diameter of glass envelope, hc air is the convection heat transfer coefficient of air in the atmosphere, T go is the average temperature of outer surface of glass envelope, T atm is the temperature of atmosphere, k air is the thermal conductivity of air when temperature is 0.5(T go + T atm ), Nu air is the Nusselt number of air when temperature is 0.5(T go + T atm ), ε go is the emissivity of outer surface of glass envelope, T sky is sky temperature, which can be calculated by [36]: The Nusselt number can be calculated by the correlation [34]: Nu air = C · Re m go Pr n atm · ( Pr atm Pr go ) 0.25 (19) with: And n = 0.37, for Pr atm ≤ 10 n = 0.36, for Pr atm > 10 where subscript atm is the atmosphere and subscript go is the outer surface of glass envelope. The parameters used in this the parabolic trough solar collector are shown in Table 1. According to the parameters in the table and the above formula, the required collector area under different conditions can be calculated.

Other Thermodynamic Components
The compression and expansion process are calculated from the inlet. The calculation procedure is shown as below: In the above two equations, h is the enthalpy and η is the isentropic efficiency. For the subscripts, c stands for compressor and t stands for turbine. Moreover, 'in' refers to inlet conditions; 'out' refers to outlet conditions, and 'id' represents ideal state with no entropy increase.
In this work, printed circuit heat exchanger (PCHE) is used for the heat exchanger and the minimum end difference during heat exchange is set to 5 • C in the calculation. Then, a few assumptions are made for simplification. The mass flow through each channel and the geometry of each channel are the same. Based on the assumptions, heat exchanger is divided into several axial units, shown in Figure 5 (the number of units is i in the figure). The j-th node is taken out to establish the heat transfer equation.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 20 In the above two equations, h is the enthalpy and η is the isentropic efficiency. For the subscripts, c stands for compressor and t stands for turbine. Moreover, 'in' refers to inlet conditions; 'out' refers to outlet conditions, and 'id' represents ideal state with no entropy increase. In this work, printed circuit heat exchanger (PCHE) is used for the heat exchanger and the minimum end difference during heat exchange is set to 5°C in the calculation. Then, a few assumptions are made for simplification. The mass flow through each channel and the geometry of each channel are the same. Based on the assumptions, heat exchanger is divided into several axial units, shown in Figure 5 (the number of units is i in the figure). The j-th node is taken out to establish the heat transfer equation. Then, the calculation formula of overall heat transfer coefficient for unit j can be obtained according to [37]: Where the subscripts h and c here represent hot side and cold side respectively. They are similar during the calculation. Besides, k refers to the thermal conductivity of material and t stands for the plate thickness. Where Uj_h can be calculated as [37]: Where: In the heat exchangers, the equations above are also adopted to calculate Uj_c.

Equipment Cost Accountings
The equipment cost is used to compare the economics of the proposed three system configurations, so we make the following assumptions: 1. The cost of the core equipment is considered the total cost of the system. In actuality, the cost of the core equipment consists of only part of the total cost, with pipes, construction fee, and engineering design as additions. Nevertheless, in this design process, we assume these costs are proportional to the total amount. Hence, we can obtain sufficient analysis to conduct the comparison Then, the calculation formula of overall heat transfer coefficient for unit j can be obtained according to [37]: where the subscripts h and c here represent hot side and cold side respectively. They are similar during the calculation. Besides, k refers to the thermal conductivity of material and t stands for the plate thickness. Where U j_h can be calculated as [37]: where: In the heat exchangers, the equations above are also adopted to calculate U j_c .

Equipment Cost Accountings
The equipment cost is used to compare the economics of the proposed three system configurations, so we make the following assumptions: 1.
The cost of the core equipment is considered the total cost of the system. In actuality, the cost of the core equipment consists of only part of the total cost, with pipes, construction fee, and engineering design as additions. Nevertheless, in this design process, we assume these costs are proportional to the total amount. Hence, we can obtain sufficient analysis to conduct the comparison of the proposed three system configurations, considering thermodynamics and thermoeconomic performance.

2.
Characteristic parameters are used to calculate the cost of different core equipment. However, geometric design parameters and flow forms of turbomachinery will affect the purchasing price of the device, which will undoubtedly lead to some deviations compared to actual cost of the S-CO 2 system. However, this deviation can be acceptable for these three kinds of systems because we can consider it as consistent.
Appl. Sci. 2020, 10, 3123 9 of 20 The common cost calculation of the S-CO 2 cycle is mostly high temperature, while the temperature involved in this paper is relatively low. Hence, we use the trans-critical CO 2 cycle as the cost reference, which corresponds to more similar temperature levels. The functions built to calculate the equipment costs are presented in Table 2. They are obtained from Henchoz et al [38] and Price et al [10]. Moreover, in the calculation process of the cooler, the constant heat transfer coefficient between water and supercritical CO 2 is set to 1600 W·m −2 · • C −1 . This value is obtained by averaging the heat exchange process between water and S-CO 2 in a larger pressure range [39].

Performance Criteria
In this study, specific work, overall efficiency, and cost of core equipment are chosen as the performance criteria.
During operation, the systems can produce electric power (W net ) by the absorbed solar energy (Q solar ). Therefore, the overall efficiency of cycle can be calculated by W net /Q solar . The specific work is calculated as W net /m c (m c means the mass flow). These two parameters are used to characterize the thermodynamic properties of the system.
The economic performance of the system is mainly calculated by the cost of core equipment mentioned in Section 3.3. In fact, the thermodynamic performance, such as specific work, can be used to calculate the return of the systems. However, due to the lack of other data, the calculation of the investment recovery of the systems is not involved in this study.
In addition, for the calculation of the intercooling cycle and the convenience of comparison, we consider that the cooler between the no.1 compressor and no.2 compressor can reduce the inlet temperature of the no.2 compressor, which makes it equal with the inlet temperature of the no.1 compressor. The pressure ratios of the two compressors are also the same.
For the calculation of the recompression cycle, we assume that the system has achieved the best shunt coefficient to ensure that the minimum end difference is exactly the set temperature difference.

Model Validation
In this paper, the calculation model is verified by comparing it with Turchi [21] and Guo [40]. The cycle layout in those studies show a recompression re-heating cycle, and the parameters shown in Table 3 are used. The comparison of the results is shown in Figure 6, which can illustrate that the model in the present study is very consistent with that of Turchi and Guo.   Figure 6. Results of the present study compared to data from Turchi's and Guo's studies.

Parameter Study
It is very important to study the influence of key parameters on system performance. Thus, the parameter study, including effects of turbine inlet temperature and pressure, compressor inlet temperature, and pressure on specific work, overall efficiency, and cost of core equipment is performed. In order to make a fair comparison of the three types of cycles, and considering that the solar collector used is the parabolic trough solar collector, the basic condition is shown in Table 4.  9 show the results of the parameter study on three cycle layouts. It should be noted that the heat exchange efficiency and working temperature of the parabolic trough solar collector are low, so overall efficiency of the three cycles are low, but the actual thermal cycle efficiency of them is about 25%-30%, while the turbine inlet temperature is 350 ˚C. In contrast, due to the limitation of the inlet temperature of the turbine, the actual thermal efficiency of ORC Rankine cycle is basically 20%-24% [41], which is lower than the actual thermal efficiency of the S-CO2 cycles in this study. However, Figure 6. Results of the present study compared to data from Turchi's and Guo's studies.

Parameter Study
It is very important to study the influence of key parameters on system performance. Thus, the parameter study, including effects of turbine inlet temperature and pressure, compressor inlet temperature, and pressure on specific work, overall efficiency, and cost of core equipment is performed. In order to make a fair comparison of the three types of cycles, and considering that the solar collector used is the parabolic trough solar collector, the basic condition is shown in Table 4.  Figures 7-9 show the results of the parameter study on three cycle layouts. It should be noted that the heat exchange efficiency and working temperature of the parabolic trough solar collector are low, so overall efficiency of the three cycles are low, but the actual thermal cycle efficiency of them is about 25%-30%, while the turbine inlet temperature is 350 • C. In contrast, due to the limitation of the inlet temperature of the turbine, the actual thermal efficiency of ORC Rankine cycle is basically 20%-24% [41], which is lower than the actual thermal efficiency of the S-CO 2 cycles in this study. However, for the steam Rankine cycle, the overall efficiency is 24.5% when the turbine inlet temperature is 390 • C, which is higher than the S-CO 2 cycles in this study [42]. However, S-CO 2 cycles have fewer phase changes, so it needs no condenser arrangement and can reduce the number of valves, resulting in fewer total area and cost. Figure 7 illustrates the effects of the above parameters on a simple recuperated cycle, respectively. Specific work, overall efficiency, and cost of core equipment increase with the rising turbine inlet temperature and pressure, as shown in Figure 7a,b. Figure 7c shows that the specific work decreases and cost of core equipment increases with the rise of the compressor inlet temperature due to the increase of compression work. However, the overall efficiency changes tinily, and then decreases as the compressor inlet temperature rises. Figure 7d is the influence of compressor inlet pressure on specific work, overall efficiency, and cost of core equipment. In the range of compressor inlet pressure, there will be extreme values of specific work and overall efficiency due to pseudo critical pressure. Figure 8 illustrates the effects of the same parameters on the intercooling cycle respectively. Similar to the simple recuperated cycle, the relevant parameters generally increase as the turbine inlet temperature and pressure rises, which is shown in Figure 8a,b. It can be seen from Figure 8c that the cost of core equipment increases while specific work and overall efficiency decreases as the compressor inlet temperature rises. Figure 8d shows that overall efficiency and the specific work firstly increases as the compressor inlet pressure rises, and then decreases. Figure 9 illustrates the effects of the same parameters on the recompression cycle, respectively. The change trend shown in Figure 9a-c are, basically, similar to the above two cycles, the only difference is that with the turbine inlet pressure rising, the increase speed of overall efficiency slows down. For the recompression cycle layout, overall efficiency and the specific work have a similar trend-to intercool the cycle with the rise of the compressor inlet pressure. for the steam Rankine cycle, the overall efficiency is 24.5% when the turbine inlet temperature is 390 ˚C, which is higher than the S-CO2 cycles in this study [42]. However, S-CO2 cycles have fewer phase changes, so it needs no condenser arrangement and can reduce the number of valves, resulting in fewer total area and cost. Figure 7 illustrates the effects of the above parameters on a simple recuperated cycle, respectively. Specific work, overall efficiency, and cost of core equipment increase with the rising turbine inlet temperature and pressure, as shown in Figure 7(a) and Figure 7(b). Figure 7(c) shows that the specific work decreases and cost of core equipment increases with the rise of the compressor inlet temperature due to the increase of compression work. However, the overall efficiency changes tinily, and then decreases as the compressor inlet temperature rises. Figure 7(d) is the influence of compressor inlet pressure on specific work, overall efficiency, and cost of core equipment. In the range of compressor inlet pressure, there will be extreme values of specific work and overall efficiency due to pseudo critical pressure. Figure 8 illustrates the effects of the same parameters on the intercooling cycle respectively. Similar to the simple recuperated cycle, the relevant parameters generally increase as the turbine inlet temperature and pressure rises, which is shown in Figure 8(a) and Figure 8(b). It can be seen from Figure 8(c) that the cost of core equipment increases while specific work and overall efficiency decreases as the compressor inlet temperature rises. Figure 8(d) shows that overall efficiency and the specific work firstly increases as the compressor inlet pressure rises, and then decreases. Figure 9 illustrates the effects of the same parameters on the recompression cycle, respectively. The change trend shown in Figure 9(a), Figure 9(b), and Figure 9(c) are, basically, similar to the above two cycles, the only difference is that with the turbine inlet pressure rising, the increase speed of overall efficiency slows down. For the recompression cycle layout, overall efficiency and the specific work have a similar trend-to intercool the cycle with the rise of the compressor inlet pressure.   It can be seen from the results of the parameter study that the performance of the system is obviously affected by several parameters, and the influence of the compressor inlet temperature and pressure on performance are not simply monotone. Furthermore, using the same parameters set  It can be seen from the results of the parameter study that the performance of the system is obviously affected by several parameters, and the influence of the compressor inlet temperature and pressure on performance are not simply monotone. Furthermore, using the same parameters set It can be seen from the results of the parameter study that the performance of the system is obviously affected by several parameters, and the influence of the compressor inlet temperature and pressure on performance are not simply monotone. Furthermore, using the same parameters set cannot maximize the specific work, overall efficiency, and minimize the cost of core equipment at the same time.
Moreover, the parameter study is difficult to correct and complete comparison between the cycle layouts because the analysis is calculated on the base condition. Therefore, it is necessary to perform multi-objective optimization and then compare the various cycle layouts.

Multi-Objective Optimization
In the previous section, we analyzed three performance parameters as specific work, overall efficiency, and cost of core equipment. For the S-CO 2 power cycle using solar energy, we were concerned about the net power-that the system could provide in the operation condition and the purchase cost of the system during the construction period. Thus, we chose the specific work and cost of core equipment during the optimization.
The objective of this study is to maximize the system net power and, in the meantime, minimize the investment cost. Table 5 gives the selected decision variables and its corresponding range. The multi-objective optimization process is conducted with GA and the solver is MATLAB R2019b. The basic steps of genetic algorithm include coding, population initialization, fitness calculation, evolution (composed of selection, crossover, mutation), and iteration. In this work, the optimal setting of the GA algorithm is set to 50 for each generation of population at population initialization. Moreover, when the genetic algebra reaches 400, the calculation terminates, which can limit the number of iterations. The reason for this set is-we find that if the total number of calculations is less than this set, the results have much more randomness. When we chose these two parameters, the results of the multiple calculations had little change. Based on the three research cycles, the optimal Pareto fronts are presented in Figures 10-12.
The optimal Pareto fronts show the trade-off between the cost of core equipment and the specific work under the optimum conditions. In addition, the scatter points in the figures are for each working condition and the optimization results of the calculation process. It is worth noting that there is a significant discontinuity in the Pareto optimization results in Figures 10-12, especially when the specific work is low. The results from the algorithm and its capability of evolution can first analyze the first-generation cluster. Then, the algorithm keeps the representative individuals and proceeds. These representative individuals are selected from the clusters of best performance. When a series of results barely improve through optimization, a few elite individuals are further selected. Hence, this results in a smaller portion of the global Pareto front with fewer solutions condensing. However, this discontinuity has no effect on the analysis of the results and the comparison of the systems.
From Figure 10, the values of the optimization objectives for a simple recuperated cycle increased from 30 to 70 kW/kg for specific work, and from $155,000 to $210,000 USD for the cost of core equipment. Outside these ranges, we can conclude that-for the poorer performance solutions-the cost cannot reduce, or, due to thermodynamic limits, the system net power cannot increase.  Figure 11 presents the result of intercooling configuration. The values of the optimization objectives extend between 27 and 65 kW/kg for specific work and from $255,000 to $365,000 USD for the cost of core equipment. Compared with a simple recuperated cycle, the economic cost of the intercooling cycle is too expensive and the specific work has not been significantly improved. Moreover, it can be concluded that the optimization result of the range of cost of the core equipment in the intercooling cycle is larger than the optimization result range of the simple recuperated cycle. Hence, for this cycle, reasonable selection of system operation parameters is more important for the cost reason.    Figure 11 presents the result of intercooling configuration. The values of the optimization objectives extend between 27 and 65 kW/kg for specific work and from $255,000 to $365,000 USD for the cost of core equipment. Compared with a simple recuperated cycle, the economic cost of the intercooling cycle is too expensive and the specific work has not been significantly improved. Moreover, it can be concluded that the optimization result of the range of cost of the core equipment in the intercooling cycle is larger than the optimization result range of the simple recuperated cycle. Hence, for this cycle, reasonable selection of system operation parameters is more important for the cost reason.  Figure 11 presents the result of intercooling configuration. The values of the optimization objectives extend between 27 and 65 kW/kg for specific work and from $255,000 to $365,000 USD for the cost of core equipment. Compared with a simple recuperated cycle, the economic cost of the intercooling cycle is too expensive and the specific work has not been significantly improved. Moreover, it can be concluded that the optimization result of the range of cost of the core equipment in the intercooling cycle is larger than the optimization result range of the simple recuperated cycle. Hence, for this cycle, reasonable selection of system operation parameters is more important for the cost reason.    To keep things straight, the optimized Pareto set of three cycles are compared and shown in Figure 13. These cycle layouts are compared in the coordinate system of "cost of core equipment and specific work". In the "cost of core equipment and specific work" coordinate system, the location closer to the lower right corner means more excellent performance. Although the device cost of this study is not fully calculated and the actual system payback period cannot be calculated correctly. However, from Figure 13, we can still see that the simple recuperated cycle has obvious advantages over the recompression cycle and intercooling cycle after optimization. In detail, under the same level of specific work, the cost of the core equipment of the simple recuperated cycle is lower than that of the recompression cycle and intercooling cycle, which can be clearly seen in the figure.
However, a previous study has shown that when the central receiver tower with heliostat field is used, the intercooling cycle and recompression cycle will achieve greater efficiency [23]. The reason for this phenomenon may be that the temperature of the S-CO2 entering turbine is too low when using the parabolic trough solar collector. Because of this bottleneck, the complex system structure cannot produce an obvious advantage of specific work while the increase of its cost is significant. Moreover, complex cycle structures are often used in superposition, such as recompression, intercooling, and reheating. Thus, they are more suitable for higher temperature scenarios. To keep things straight, the optimized Pareto set of three cycles are compared and shown in Figure 13. These cycle layouts are compared in the coordinate system of "cost of core equipment and specific work". In the "cost of core equipment and specific work" coordinate system, the location closer to the lower right corner means more excellent performance. Although the device cost of this study is not fully calculated and the actual system payback period cannot be calculated correctly. However, from Figure 13, we can still see that the simple recuperated cycle has obvious advantages over the recompression cycle and intercooling cycle after optimization. In detail, under the same level of specific work, the cost of the core equipment of the simple recuperated cycle is lower than that of the recompression cycle and intercooling cycle, which can be clearly seen in the figure.
However, a previous study has shown that when the central receiver tower with heliostat field is used, the intercooling cycle and recompression cycle will achieve greater efficiency [23]. The reason for this phenomenon may be that the temperature of the S-CO 2 entering turbine is too low when using the parabolic trough solar collector. Because of this bottleneck, the complex system structure cannot produce an obvious advantage of specific work while the increase of its cost is significant. Moreover, complex cycle structures are often used in superposition, such as recompression, intercooling, and reheating. Thus, they are more suitable for higher temperature scenarios. To keep things straight, the optimized Pareto set of three cycles are compared and shown in Figure 13. These cycle layouts are compared in the coordinate system of "cost of core equipment and specific work". In the "cost of core equipment and specific work" coordinate system, the location closer to the lower right corner means more excellent performance. Although the device cost of this study is not fully calculated and the actual system payback period cannot be calculated correctly. However, from Figure 13, we can still see that the simple recuperated cycle has obvious advantages over the recompression cycle and intercooling cycle after optimization. In detail, under the same level of specific work, the cost of the core equipment of the simple recuperated cycle is lower than that of the recompression cycle and intercooling cycle, which can be clearly seen in the figure.
However, a previous study has shown that when the central receiver tower with heliostat field is used, the intercooling cycle and recompression cycle will achieve greater efficiency [23]. The reason for this phenomenon may be that the temperature of the S-CO2 entering turbine is too low when using the parabolic trough solar collector. Because of this bottleneck, the complex system structure cannot produce an obvious advantage of specific work while the increase of its cost is significant. Moreover, complex cycle structures are often used in superposition, such as recompression, intercooling, and reheating. Thus, they are more suitable for higher temperature scenarios.

Sensitivity Analysis of a Typical Result
Through the description in Section 4.2, it can be found that the simple recuperated cycle has incomparable advantages in S-CO 2 power generation cycles using solar energy by a parabolic trough solar collector. Even though it is difficult to carry out a complete investment recovery analysis through the specific work and the cost of core equipment, we can still conclude that the simple recuperated cycle has the best economy among the three cycle types. Therefore, a typical optimization result in the simple recuperated cycle is selected for sensitivity analysis in this section.
The parameters in typical optimization results selected are shown in Table 6, and the rest of the parameters of the system are shown in Table 4. First, the cost of each core component is compared with histogram and cumulative proportion, as shown in Figure 14. It can be seen that the cost of turbine machinery, including turbine and compressor, is the main component of equipment cost. This is because the technology of other components in S-CO 2 power is mature, while the manufacturing of turbine machinery using S-CO 2 is relatively difficult and the cost is high. Therefore, it can be concluded that reducing the manufacturing cost of the turbine machinery can achieve higher economic benefits.

Sensitivity Analysis of a Typical Result
Through the description in Section 4.2, it can be found that the simple recuperated cycle has incomparable advantages in S-CO2 power generation cycles using solar energy by a parabolic trough solar collector. Even though it is difficult to carry out a complete investment recovery analysis through the specific work and the cost of core equipment, we can still conclude that the simple recuperated cycle has the best economy among the three cycle types. Therefore, a typical optimization result in the simple recuperated cycle is selected for sensitivity analysis in this section.
The parameters in typical optimization results selected are shown in Table 6, and the rest of the parameters of the system are shown in Table 4.  Figure 14. It can be seen that the cost of turbine machinery, including turbine and compressor, is the main component of equipment cost. This is because the technology of other components in S-CO2 power is mature, while the manufacturing of turbine machinery using S-CO2 is relatively difficult and the cost is high. Therefore, it can be concluded that reducing the manufacturing cost of the turbine machinery can achieve higher economic benefits. Secondly, sensitivity analysis is considered on the cost of the core equipment and specific work. Since the cost of turbine machinery accounts for the main part, and the compressor inlet temperature after optimization is rather low (which resulted in less change space), the analysis variables are mainly turbine inlet temperature, turbine inlet pressure, compressor inlet pressure, isentropic efficiency of turbine, and isentropic efficiency of the compressor. The five parameters vary by -5% and +5%. The results of the sensitivity analysis are shown in Table 7.  Secondly, sensitivity analysis is considered on the cost of the core equipment and specific work. Since the cost of turbine machinery accounts for the main part, and the compressor inlet temperature after optimization is rather low (which resulted in less change space), the analysis variables are mainly turbine inlet temperature, turbine inlet pressure, compressor inlet pressure, isentropic efficiency of turbine, and isentropic efficiency of the compressor. The five parameters vary by -5% and +5%. The results of the sensitivity analysis are shown in Table 7. Through the sensitivity analysis, it can be found that most of the changes of these variables only affect the cost of a single component, so the overall cost is not sensitive to them. The specific work is still sensitive to these variables. For the cost of core equipment, the change of pressure will cause great fluctuation; this is because the change of pressure affects the performance of two types of turbomachinery (turbine and compressor), resulting in a great change of cost. As for the specific work, in addition to the compressor efficiency, the influence of other parameters is more obvious, which is also reflected in Section 4.1. The influence of the system design parameters on the cost of core equipment is less than that on the specific work of the system. Therefore, in the actual process, more consideration should be given to the thermal performance of the system regarding the selection of system parameters, and the structural optimization and layout optimization of components should be adopted to reduce the system cost.

Conclusions
In this paper, a preliminary analysis and systematic comparison of three S-CO 2 power generation cycles using solar energy by the parabolic trough solar collector are studied. Firstly, the effects of the turbine inlet temperature and pressure, compressor inlet temperature and pressure on specific work, overall efficiency, and cost of core equipment of S-CO 2 Brayton cycles are discussed. After that, multi-objective optimizations are studied, taking the specific work and the cost of core equipment as the objectives. The optimized Pareto sets obtained through multi-objective optimizations are used to compare different cycle layouts. The main conclusion follows: (1) It results in an optimal specific work ranging between 30 and 70 kW/kg, and optimal cost of core equipment ranging between $155,000 and $210,000 USD for the simple recuperated cycle. The values of the optimization objectives extend between 27 and 65 kW/kg for specific work and between $255,000 and $365,000 USD for cost of core equipment in the intercooling cycle. Moreover, the values of the optimization objectives extend between 25 and 50 kW/kg for specific work and from $270,000 and $370,000 USD for cost of core equipment in the recompression cycle. (2) The simple recuperation cycle layout shows more excellent performance than the intercooling cycle layout and the recompression cycle layout in terms of cost, while the advantage in specific work of the intercooling cycle layout and the recompression cycle layout is not obvious. (3) The sensitivity analysis of a typical optimization result of a simple regeneration cycle shows that the change of parameters have little influence on the cost of core equipment and much influence on specific work. (4) Through this study, it can be found that the turbine inlet temperature is low using a parabolic trough solar collector, and the complex SCO 2 Brayton cycles are not dominant. (5) The preliminary analysis and systematic comparison in this study can be useful in selecting cycle layout using solar energy by a parabolic trough solar collector when there are requirements for the specific work and the cost of core equipment. Moreover, it can also provide reference for the comparison and selection of other cycle types.
In future studies, a comprehensive analysis of investments profitability based on more accurate cost accounting methods will obtain more practical system comparison results. For components, improving the manufacturing level of turbine machinery and reducing its cost will make the S-CO 2 power generation cycles more competitive.