Ecological Function Analysis and Optimization of a Recompression S-CO2 Cycle for Gas Turbine Waste Heat Recovery

In this paper, a recompression S-CO2 Brayton cycle model that considers the finite-temperature difference heat transfer between the heat source and the working fluid, irreversible compression, expansion, and other irreversibility is established. First, the ecological function is analyzed. Then the mass flow rate, pressure ratio, diversion coefficient, and the heat conductance distribution ratios (HCDRs) of four heat exchangers (HEXs) are chosen as variables to optimize cycle performance, and the problem of long optimization time is solved by building a neural network prediction model. The results show that when the mass flow rate is small, the pressure ratio, the HCDRs of heater, and high temperature regenerator are the main influencing factors of the ecological function; when the mass flow rate is large, the influences of the re-compressor, the HCDRs of low temperature regenerator, and cooler on the ecological function increase; reasonable adjustment of the HCDRs of four HEXs can make the cycle performance better, but mass flow rate plays a more important role; the ecological function can be increased by 12.13%, 31.52%, 52.2%, 93.26%, and 96.99% compared with the initial design point after one-, two-, three-, four- and five-time optimizations, respectively.


Introduction
A gas turbine (GT) has the advantages of a compact structure, small volume, light weight, high power density, and high thermal efficiency and is mostly used as a driving machine in surface ships and civilian ships. However, its waste heat temperature is high and has the possibility of re-utilization. Strengthening the utilization of this part of the waste heat can effectively improve the energy utilization rate [1]. The characteristics of high efficiency, better stability, economy, compactness, and simplicity make the Brayton cycle (BC) system based on supercritical carbon dioxide (S-CO 2 ) have broad application prospects [2][3][4][5][6]. Moreover, the optimal heat source temperature range of the S-CO 2 Brayton cycle (SCBC) is 450~700 • C [1], which matches the exhaust gas temperature of GT. Therefore, it is an ideal cycle for recovering the waste heat of GTs [7][8][9][10][11].
Sulaiman et al. [12] analyzed the performances of five different types of SCBCs integrated with solar thermal power plants, including regenerative, recompression, etc. The results show that the recompression S-CO 2 Brayton cycle (RCSCBC) has the highest thermal efficiency. Vasquez et al. [13] analyzed the regenerative, recompressed, partially cooled, and intercooled SCBC integrated with the solar receiver and obtained the optimal operating conditions through multi-objective thermodynamic optimization. Anton et al. [14] used thermodynamic modeling and analysis based on the combination of the RCSCBC with a sodium-cooled fast reactor to explore the effect of different structural layouts on the thermal Entropy 2022, 24, 732 2 of 19 efficiency of the cycle. Alharbi et al. [15] proposed a multi-effect desalination system driven by RCSCBC waste heat, which can be used for electricity and fresh water production.
Liu et al. [16] carried out the design and analysis of the RCSCBC for ships with an output power of 40MW. The research results showed that the designed RCSCBC reached an efficiency of 45.06% at 823.15 K, which meets the design requirements and improves by 8.28% over the simple regenerative cycle. Li et al. [17] conducted a comparative study on SCBC for waste heat recovery from GTs, where selected cycles were optimized and analyzed in terms of system efficiency, etc. Khadse et al. [18] conducted a thermodynamic optimization study of the RCSCBC for waste heat recovery applications; the results show that higher pressures can result in a higher cycle power output.
Mohammadi et al. [19] proposed a combined cycle consisting of a GT cycle, an RC-SCBC, and an organic Rankine cycle. The research results showed that the cycle could effectively improve the utilization rate of heat recovery. Saeed and Kim [20] proposed a novel SCBC, which consisted of four compression processes and one expansion process, with a better ability to integrate with heat sources by comparing with the regenerative, recompression, and intercooling and regenerative cycles, and the results showed that the thermal efficiency (TEF) was comparable to that of the intercooling and regenerative cycle. Khatoon et al. [21] investigated the thermodynamic performance of an RCSCBC integrated with a direct air-cooled heat exchanger (HEX), and the research results showed that the cooling process of S-CO 2 was sensible heat transfer, and the thermophysical properties changed nonlinearly.
Although the above work has carried out analysis and research on different types of SCBCs, these studies have not considered in detail the effects of irreversible factors such as irreversible compression and irreversible expansion. The overall optimum performances of the SCBC device have not been obtained under the condition of a constant overall size of the heat exchange devices.
Finite-time thermodynamics (FTT) [22][23][24][25][26][27] realizes the cross-integration of multiple basic discipline theories including thermodynamics, heat transfer, and fluid mechanics and is an important branch of modern thermodynamic theory. FTT seeks more practical and useful performance bounds and obtains the optimal performance under limited size constraints and the optimal way to achieve the goal [28][29][30][31], whose results have important engineering value. FTT has been used in the studies of a wide variety of processes and cycles, such as the organic Rankine cycle [32][33][34], Carnot cycle [35], thermoelectric generator [36], Kalina cycle [37,38], and Stirling heat engine [39,40]. Combining the analysis method of FTT theory with the SCBC, on the one hand, makes up for the deficiency of the traditional pure thermodynamic analysis method in the irreversibility analysis of the SCBC. On the other hand, it can further expand the scope of application research objects of FTT, which has important theoretical and engineering significances.
In Brayton cycle research based on the "ideal gas hypothesis", Cheng and Chen [41][42][43] established a Brayton cycle model by considering the irreversible losses of the cyclic processes; derived the functional expressions of the cyclic power, efficiency, and ecological functions; and optimized the cycle with these performance objectives. Ust et al. [44] analyzed and studied the irreversible closed Brayton cycle with a constant-temperature heat source and the ecological performance coefficient. Kaushik et al. [45] studied the closed variable-temperature irreversible regenerative Brayton cycle model, aiming at cycle power and efficiency. Tyagi et al. [46][47][48][49] established a closed irreversible intercooling regenerative Brayton cycle model with constant-temperature and variable-temperature heat sources, respectively, and analyzed the power, efficiency, and ecological performance coefficient of the cycle.
In the FTT study of the SCBC, Na et al. [50] studied a preheating SCBC and investigated the influences of related parameters on thermal efficiency and net power output. Jin et al. [51] established an FTT model of the regenerative S-CO 2 Brayton cycle and carried out multi-objective optimization. An RCSCBC model that considers the finite temperature difference heat transfer between the heat source and the working fluid (WF), irreversible compression, and irreversible expansion will be established in this paper. Firstly, the ecological functions are analyzed; then, under the constraint of the total HEX inventory, taking the mass flow rate, pressure ratio, diversion coefficient, and the heat conductance of each HEXs as optimization variables, performance optimization for ecological function will be performed; and the optimal design results with practical engineering application value will be given.

Physical Model
Figures 1 and 2 are the RCSCBC device diagram and T-s diagram, respectively. The RCSCBC device is mainly composed of the main compressor, re-compressor, low temperature regenerator (LTR), high temperature regenerator (HTR), heater, turbine, and cooler. In Figure 2, 1-2 S and 1-2 denote the ideal reversible and the actual irreversible adiabatic compression processes in the main compressor, respectively; 2-8 and 8-3 denote the heat absorption process in the LTR and HTR, respectively; 3-4 denotes the heat absorption process from the high-temperature heat reservoir; 4-5 and 4-5 S denote the actual irreversible and the ideal reversible adiabatic expansion processes, respectively; 5-6 and 6-7 denote the heat release processes in the LTR and HTR. At state point 7, the WF is divided by a flow divider, and a part of the WF enters the re-compressor; 8-7 is the actual irreversible adiabatic compression process in the re-compressor. After this, part of the WF leaves the re-compressor; it reaches the temperature of state point 8 after an isenthalpic mixing process with the WF after endothermic heating in the LTR, and 7-1 indicates the exothermic process. lyzed; then, under the constraint of the total HEX inventory, taking the mass flow rate, pressure ratio, diversion coefficient, and the heat conductance of each HEXs as optimization variables, performance optimization for ecological function will be performed; and the optimal design results with practical engineering application value will be given.

Physical Model
Figures 1 and 2 are the RCSCBC device diagram and T-s diagram, respectively. The RCSCBC device is mainly composed of the main compressor, re-compressor, low temperature regenerator (LTR), high temperature regenerator (HTR), heater, turbine, and cooler. In Figure 2, 1-2S and 1-2 denote the ideal reversible and the actual irreversible adiabatic compression processes in the main compressor, respectively; 2-8 and 8-3 denote the heat absorption process in the LTR and HTR, respectively; 3-4 denotes the heat absorption process from the high-temperature heat reservoir; 4-5 and 4-5S denote the actual irreversible and the ideal reversible adiabatic expansion processes, respectively; 5-6 and 6-7 denote the heat release processes in the LTR and HTR. At state point 7, the WF is divided by a flow divider, and a part of the WF enters the re-compressor; 8-7 is the actual irreversible adiabatic compression process in the re-compressor. After this, part of the WF leaves the re-compressor; it reaches the temperature of state point 8 after an isenthalpic mixing process with the WF after endothermic heating in the LTR, and 7-1 indicates the exothermic process.
We regarded wf m as the mass flow rate of the WF;   We regarded m wf as the mass flow rate of the WF; m wf,1 as the mass flow rate of WF entering the cooler through the splitter after the S-CO 2 comes out of the LTR; and m wf,2 as the mass flow rate of WF entering the re-compressor through the splitter. Defining the ratio of m wf,1 to m wf as the diversion coefficient x p , there is the following formula:  The efficiencies of main compressor, re-compressor, and turbine are represented by η c , η c,2 , and η t :    The efficiencies of main compressor, re-compressor, and turbine are represented by η c , η c,2 , and η t : where h 1 , h 2s , h 8,LTR , h 8,RC , h 8 , h 2 , h 3 , h 4 , h 5s , h 5 , h 6 , h 8s , and h 7 are the specific enthalpies of the corresponding state points. π is the pressure ratio. T H,in , T H,out , T L,in , and T L,out are the inlet and outlet temperatures of the hot and cold source, respectively. T 8,LTR is the temperature of the WF after being heated by LTR, T 8,RC is the temperature of S-CO 2 after leaving the re-compressor, and T 8 is the temperature of WF after the isenthalpic mixing process; T 0 is the ambient temperature. m H and m L are the mass flow rates of the heat source and the cold source. The heat conductances of the heater, the cooler, the LTR, and the HTR are expressed as U H , U L , U LTR , and U HTR , respectively. According to the HEX theory, the heat absorption rate Q H , heat release rate Q L , and heat recovery rates Q LTR and Q HTR are, respectively, given by [52][53][54][55][56]: Entropy 2022, 24, 732

of 19
From the thermal properties of the WF, Q H , Q L , Q HTR , and Q LTR are, respectively, The total HEX inventory U T is the sum of the U H , U L , U LTR , and U HTR .
Define the heat conductance distribution ratio ψ. The ψ of the heater, the cooler, LTR, and HTR are expressed as ψ H = U H /U T ,ψ L = U L /U T , ψ LTR = U LTR /U T , and ψ HTR = U HTR /U T , respectively. The ψ H , ψ L , ψ LTR , and ψ HTR have the following relationship: For the RCSCBC, the W net is the difference between the turbine power W t and the compressor power consumption, W c , and the compressor power consumption is the sum of the main compressor power consumption W c,1 and re-compressor power consumption W c,2 . They are given by The entropy production rate and ecological function E of the cycle are: Since S-CO 2 is an actual gas, the p-V-T (V is the volume) relationship does not satisfy the ideal gas equation of state. In addition, the relationship between the specific heat capacity of S-CO 2 , pressure, and temperature is very complicated, so it is impossible to obtain the explanatory formula for the temperature and specific enthalpy at each state point of the cycle. Numerical solutions for temperature and specific enthalpy can only be obtained programmatically. Under the given boundary conditions (T H,in , T L,in , U H , U L , U LTR , U HTR , η c , η c,2 , and η t ), the nonlinear equation system composed of the Equations (1)-(16) is solved by the refpropm function and the @fsolve function, and the temperature and specific enthalpy of each state point can be obtained. The calculated specific enthalpy is further imported into Equations (17)-(23) to solve the NPO and TEF under the given boundary conditions. The specific calculation flow chart is shown in Figure 3. The cycle calculation program is written by MATLAB software, and the physical properties of S-CO 2 were calculated by REFPROP [57]. temperature and specific enthalpy of each state point can be obtained. The calculated specific enthalpy is further imported into Equations (17)-(23) to solve the NPO and TEF under the given boundary conditions. The specific calculation flow chart is shown in Figure 3. The cycle calculation program is written by MATLAB software, and the physical properties of S-CO2 were calculated by REFPROP [57]. According to the Ref. [50], Table 1 gives the initial design point (IDP) parameters. Firstly, the cycle performance is analyzed, and the relationship between the E, xp, and π are studied and analyzed under the conditions of different wf m , η c , and η t . Then, with the goal of maximizing E, under the constraint that the UT is a fixed value, the optimization will be carried out with the wf m , xp, π , ψ H , ψ L , ψ LTR , and ψ HTR as the optimization variables.  According to the Ref. [50], Table 1 gives the initial design point (IDP) parameters.
Firstly, the cycle performance is analyzed, and the relationship between the E, x p , and π are studied and analyzed under the conditions of different m wf , η c , and η t . Then, with the goal of maximizing E, under the constraint that the U T is a fixed value, the optimization will be carried out with the m wf , x p , π, ψ H , ψ L , ψ LTR , and ψ HTR as the optimization variables.

Ecological Function Analysis of RCSCBC
Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and π = 3, respectively. Figure 4 shows the variation trend of E with the x p . Since the physical properties of CO 2 change drastically near the critical point, it is best to keep the T 1 above the critical point. In the case of given other parameters, the effect of the x p on the T 1 is very obvious. Figure 4 shows the minimum x p that can be selected according to the criterion that the T 1 must be higher than the critical temperature under different m wf . When the m wf is 120 kg·s −1 , the selectable minimum x p is 0.718, and when the m wf is 130 kg·s −1 , the selectable minimum x p is 0.76. When the m wf is 140 kg·s −1 , the selectable minimum x p is 0.808. Within the value range of 0.45 < x p < 1, the E increases first and then decreases with the x p . After exceeding the marked minimum x p , the E shows a decreasing trend. The x p and m wf are not just isolated variables; there is also an interaction between them. Further analysis of the regulation of the curve in Figure 4 showed that the smaller m wf , the smaller the minimum x p of the T 1 exceeding the critical temperature. This indicates that the WF split affects the temperature of each state point of the cycle, and the increase of m wf can expand the selection range of the x p .

Ecological Function Analysis of RCSCBC
Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and π = 3, respectively. Figure 4 shows the variation trend of E with the xp. Since the physical properties of CO2 change drastically near the critical point, it is best to keep the T1 above the critical point. In the case of given other parameters, the effect of the xp on the T1 is very obvious. Figure 4 shows the minimum xp that can be selected according to the criterion that the T1 must be higher than the critical temperature under different mwf. When the mwf is 120kg·s −1 , the selectable minimum xp is 0.718, and when the mwf is 130kg·s −1 , the selectable minimum xp is 0.76. When the mwf is 140kg·s −1 , the selectable minimum xp is 0.808. Within the value range of 0.45 < xp < 1, the E increases first and then decreases with the xp. After exceeding the marked minimum xp, the E shows a decreasing trend. The xp and mwf are not just isolated variables; there is also an interaction between them. Further analysis of the regulation of the curve in Figure 4 showed that the smaller mwf, the smaller the minimum xp of the T1 exceeding the critical temperature. This indicates that the WF split affects the temperature of each state point of the cycle, and the increase of mwf can expand the selection range of the xp. Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and xp = 0.8, respectively. Figure 5 shows the variation trend of the E with the π under different wf m s. Under the same wf m , the E increases first and then decreases with the increase of the π . Additionally, there are the following characteristics, that is, with the increase of the wf m , the optimal π corresponding to the maximum E gradually decreases. This is due to that the π affects the size of the net power output in both indirect and direct processes, and this influence will change with the change of wf m and act on the specific value of E. In practical engineering, an appropriate π should be selected according to the corresponding environment to achieve better ecological performance. Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and x p = 0.8, respectively. Figure 5 shows the variation trend of the E with the π under different m wf s. Under the same m wf , the E increases first and then decreases with the increase of the π. Additionally, there are the following characteristics, that is, with the increase of the m wf , the optimal π corresponding to the maximum E gradually decreases. This is due to that the π affects the size of the net power output in both indirect and direct processes, and this influence will change with the change of m wf and act on the specific value of E. In practical engineering, an appropriate π should be selected according to the corresponding environment to achieve better ecological performance.
Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and x p = 0.8, respectively. Figure 6 shows the variation trend of the E with the π under different η t s and η c s. With the increase of η t and η c , the maximum E of the cycle increases, and the corresponding optimal π also increases gradually. This is due to that the values of η t and η c reflect the irreversibility of the expansion and compression processes. The E, as a thermodynamic index that compromises the s g and W net , can measure the irreversibility of the cycle to a certain extent; that is, there is a correlation between η t , η c , and the E. Therefore, the higher the η t and η c , the smaller the energy loss caused by the irreversibility of the cycle, and the larger the E value.  Figure 6 shows the variation trend of the E with the π under different η t s and η c s. With the increase of η t and η c , the maximum E of the cycle increases, and the corresponding optimal π also increases gradually. This is due to that the values of η t and η c reflect the irreversibility of the expansion and compression processes. The E, as a thermodynamic index that compromises the sg and Wnet, can measure the irreversibility of the cycle to a certain extent; that is, there is a correlation between η t , η c , and the E. Therefore, the higher the η t and η c , the smaller the energy loss caused by the irreversibility of the cycle, and the larger the E value.    Figure 6 shows the variation trend of the E with the π under different η t s and η c s. With the increase of η t and η c , the maximum E of the cycle increases, and the corresponding optimal π also increases gradually. This is due to that the values of η t and η c reflect the irreversibility of the expansion and compression processes. The E, as a thermodynamic index that compromises the sg and Wnet, can measure the irreversibility of the cycle to a certain extent; that is, there is a correlation between η t , η c , and the E. Therefore, the higher the η t and η c , the smaller the energy loss caused by the irreversibility of the cycle, and the larger the E value.   Figure 7 shows the three-dimensional relationship between the E, m wf , and π. From Figure 7, E decreases gradually with the increase of the m wf and increases first and then decreases with the increase of the π. The spherical point in Figure 7 is the maximum E when the T 1 exceeds the critical temperature. Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and xp = 0.8, respectively. Figure 7 shows the three-dimensional relationship between the E, wf m , and π . From Figure 7, E decreases gradually with the increase of the wf m and increases first and then decreases with the increase of the π . The spherical point in Figure 7 is the maximum E when the T1 exceeds the critical temperature. Let mwf = 140 kg·s −1 , ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, and ψ L = 0.3, respectively. Figure 8 shows the three-dimensional relationship between the E, xp, and π . The spherical point in Figure 8 is the maximum E when the T1 exceeds the critical temperature. It can be seen that the extreme points in Figures 7 and 8 are far away from the actual maximum value point. This is due to that the π , mwf, and xp have great influences on the temperature at each state point of the cycle, and the T1 will directly affect the stability of the cycle operation. Under the combined action of these factors, the value range of the E is not the entire surface, so the optimal value that can be selected is not the actual maximum value.  Let m wf = 140 kg·s −1 , ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, and ψ L = 0.3, respectively. Figure 8 shows the three-dimensional relationship between the E, x p , and π. The spherical point in Figure 8 is the maximum E when the T 1 exceeds the critical temperature. It can be seen that the extreme points in Figures 7 and 8 are far away from the actual maximum value point. This is due to that the π, m wf , and x p have great influences on the temperature at each state point of the cycle, and the T 1 will directly affect the stability of the cycle operation. Under the combined action of these factors, the value range of the E is not the entire surface, so the optimal value that can be selected is not the actual maximum value.
Let ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, ψ L = 0.3, and xp = 0.8, respectively. Figure 7 shows the three-dimensional relationship between the E, wf m , and π . From Figure 7, E decreases gradually with the increase of the wf m and increases first and then decreases with the increase of the π . The spherical point in Figure 7 is the maximum E when the T1 exceeds the critical temperature. Let mwf = 140 kg·s −1 , ψ LTR = 0.1, ψ HTR = 0.2, ψ H = 0.4, and ψ L = 0.3, respectively. Figure 8 shows the three-dimensional relationship between the E, xp, and π . The spherical point in Figure 8 is the maximum E when the T1 exceeds the critical temperature. It can be seen that the extreme points in Figures 7 and 8 are far away from the actual maximum value point. This is due to that the π , mwf, and xp have great influences on the temperature at each state point of the cycle, and the T1 will directly affect the stability of the cycle operation. Under the combined action of these factors, the value range of the E is not the entire surface, so the optimal value that can be selected is not the actual maximum value.  Let m wf = 140 kg·s −1 , x p = 0.8, π = 3, and ψ LTR = 0.1, respectively. Figure 9 shows the three-dimensional relationship between the E, ψ H , and ψ HTR . The spherical point in the Figure 9 is the maximum point of the E, and the ψ H at this point is significantly higher than that of the ψ HTR . It shows that under the condition of this parameter, the heater has a great influence on the E. The best performance of the cycle can be obtained by adjusting the values of the m wf , π, and x p and further optimizing the ψ H , ψ HTR , ψ L , and ψ LTR .
Let mwf = 140 kg·s , xp = 0.8, = 3, and ψ LTR = 0.1, respectively. Figure 9 shows the three-dimensional relationship between the E, ψ H , and ψ HTR . The spherical point in the Figure 9 is the maximum point of the E, and the ψ H at this point is significantly higher than that of the ψ HTR . It shows that under the condition of this parameter, the heater has a great influence on the E. The best performance of the cycle can be obtained by adjusting the values of the wf m , π , and xp and further optimizing the ψ H , ψ HTR , ψ L , and ψ LTR .

Performance Optimization
This section will aim to maximize E under the constraint that the UT is a fixed value. The optimization is carried out with the wf m , xp, π , ψ H , ψ HTR , and ψ LTR as the optimization variables, and the specific optimization process is shown in Figure 10.

Performance Optimization
This section will aim to maximize E under the constraint that the U T is a fixed value. The optimization is carried out with the m wf , x p , π, ψ H , ψ HTR , and ψ LTR as the optimization variables, and the specific optimization process is shown in Figure 10. LTR three-dimensional relationship between the E, ψ H , and ψ HTR . The spherical point in the Figure 9 is the maximum point of the E, and the ψ H at this point is significantly higher than that of the ψ HTR . It shows that under the condition of this parameter, the heater has a great influence on the E. The best performance of the cycle can be obtained by adjusting the values of the wf m , π , and xp and further optimizing the ψ H , ψ HTR , ψ L , and ψ LTR .

Performance Optimization
This section will aim to maximize E under the constraint that the UT is a fixed value. The optimization is carried out with the wf m , xp, π , ψ H , ψ HTR , and ψ LTR as the optimization variables, and the specific optimization process is shown in Figure 10. The parameter constraints are shown as follows: At the same time, it is best to keep the T 1 above the critical temperature, so an additional constraint is given by The performance optimization of the RCSCBC includes six optimization variables, including m wf , π, x p and ψ LTR , ψ H , ψ HTR . Through the traditional solution method, there will be a problem of a large amount of calculation, and the @fsolve function is used in the calculation process to rely heavily on the initial value. The interruption of the calculation program caused by calling the REFPROP physical property library to report an error has added more difficulties to the research work. Using neural network prediction can increase the calculation speed, eliminate the need to calculate the ecological function by solving the nonlinear equation system, and get rid of the program's dependence on the calculation of the initial value. In the case of many optimization variables, neural network prediction is an effective method to simplify the calculation steps.
Therefore, this section introduces neural network prediction in the optimization of the cycle's single-objective performance.
The method of using the neural network is relatively simple. As long as the corresponding optimization variables are input into the neural network, the corresponding performance target can be obtained. A neural network can be thought of as a very convenient function. In the optimization process, it is only necessary to specify the value range of the corresponding optimization variables and then use the global search algorithm @globalsearch to call the neural network to obtain the optimal ecological function under different conditions. Among them, constructing a neural network requires correct sample data to construct a mapping function, so its training is carried out on the basis of computational data. By constructing a neural network with a large sample amount of data, the target value of the cycle can be predicted only by entering the corresponding parameter values. Additionally, the more sample data, the closer the final trained model is to the real function. The parameter settings for training the neural network are shown in Table 2. Neural networks are not just the starting point for optimization algorithms. When there are few optimization variables, the calculation process is relatively simple, but every time one more variable is released, the calculation process becomes complicated. Therefore, a global search algorithm is used for each optimization to call the neural network for optimization.
The construction of the neural network model requires sample points, and each optimization variable is used as the input value; ecological function is the output value. The value range of each variable is shown in Equation (24). There were 6227 random sample points within the calculation range obtained by calculation, and these data are used as samples to train the neural network prediction model.
To verify the reliability of neural network prediction model, it is necessary to re-collect 36 sets of data and compare the calculated values with the predicted values. The re-collected 36 sets of test data cannot be any of the sample data. The comparison results are shown in Figure 11. According to Figure 11, neural network prediction is reliable.
Minimum training target error 1 × 10 −6 1 × 10 −7 Performance function mse mse Neural networks are not just the starting point for optimization algorithms. When there are few optimization variables, the calculation process is relatively simple, but every time one more variable is released, the calculation process becomes complicated. Therefore, a global search algorithm is used for each optimization to call the neural network for optimization.
The construction of the neural network model requires sample points, and each optimization variable is used as the input value; ecological function is the output value. The value range of each variable is shown in Equation (24). There were 6227 random sample points within the calculation range obtained by calculation, and these data are used as samples to train the neural network prediction model.
To verify the reliability of neural network prediction model, it is necessary to re-collect 36 sets of data and compare the calculated values with the predicted values. The recollected 36 sets of test data cannot be any of the sample data. The comparison results are shown in Figure 11. According to Figure 11, neural network prediction is reliable. Let π = 3, xp = 0.8, mwf = 140 kg·s −1 , and ψ LTR = 0.1. Figure 12 shows the variation trend of E and its corresponding ψ H and ψ L with ψ HTR . From Figure 12, with the increase of ψ HTR , E first increases and then decreases, ψ H gradually decreases, and ψ L increases first and then decreases. When ψ HTR = 0.38, ψ H = 0.25, and ψ L = 0.27, E reaches the maximum value. The reason for the difference from the optimal value given in Figure 9 comes from two aspects: one is that there is a certain error in the neural network itself, and the other is that the prediction of the neural network has volatility, but this part of the gap is within an acceptable range. Let π = 3, xp = 0.8, and mwf = 140. Figure 13 shows the maximum E and the corresponding optimal ψ HTR , ψ H , and ψ L under different ψ LTR s. From Figure 13, with the increase of ψ LTR , E first increases and then decreases, and its corresponding optimal ψ HTR has a very obvious downward trend. ψ L first increases and then decreases, but the change is small. Although ψ H also shows a downward trend, the change is relatively small. It shows that under the premise of given other parameters, there is a relationship Let π = 3, x p = 0.8, m wf = 140 kg·s −1 , and ψ LTR = 0.1. Figure 12 shows the variation trend of E and its corresponding ψ H and ψ L with ψ HTR . From Figure 12, with the increase of ψ HTR , E first increases and then decreases, ψ H gradually decreases, and ψ L increases first and then decreases. When ψ HTR = 0.38, ψ H = 0.25, and ψ L = 0.27, E reaches the maximum value. The reason for the difference from the optimal value given in Figure 9 comes from two aspects: one is that there is a certain error in the neural network itself, and the other is that the prediction of the neural network has volatility, but this part of the gap is within an acceptable range.
Let π = 3, x p = 0.8, and m wf = 140. Figure 13 shows the maximum E and the corresponding optimal ψ HTR , ψ H , and ψ L under different ψ LTR s. From Figure 13, with the increase of ψ LTR , E first increases and then decreases, and its corresponding optimal ψ HTR has a very obvious downward trend. ψ L first increases and then decreases, but the change is small. Although ψ H also shows a downward trend, the change is relatively small. It shows that under the premise of given other parameters, there is a relationship of mutual influence and restriction between the LTR and HTR. The reasonable distribution of ψ HTR and ψ LTR can make the ecological performance of the RCSCBC reach a better level. of mutual influence and restriction between the LTR and HTR. The reasonable distribution of ψ HTR and ψ LTR can make the ecological performance of the RCSCBC reach a better level.  Let π = 3 and xp = 0.8. Figure 14 shows the E of the RCSCBC and its corresponding ψ HTR , ψ H , ψ LTR , and ψ L under different wf m . It can be seen from Figure 14 that, given the xp and π , the optimal value of the cyclic E first increases and then decreases with the increase of the wf m . Additionally, it reaches the maximum value when the wf m is in the range of 100~110 kg·s −1 . The changing trends of the four heat conductance distribution ratios with the mwf show that ψ H gradually decreases, ψ L gradually increases, and ψ LTR and ψ HTR fluctuate up and down within a certain range. When the mwf is 70, 80, of mutual influence and restriction between the LTR and HTR. The reasonable distribution of ψ HTR and ψ LTR can make the ecological performance of the RCSCBC reach a better level.  Let π = 3 and xp = 0.8. Figure 14 shows the E of the RCSCBC and its corresponding ψ HTR , ψ H , ψ LTR , and ψ L under different wf m . It can be seen from Figure 14 that, given the xp and π , the optimal value of the cyclic E first increases and then decreases with the increase of the wf m . Additionally, it reaches the maximum value when the wf m is in the range of 100~110 kg·s −1 . The changing trends of the four heat conductance distribution ratios with the mwf show that ψ H gradually decreases, ψ L gradually increases, and ψ LTR and ψ HTR fluctuate up and down within a certain range. When the mwf is 70, 80, Let π = 3 and x p = 0.8. Figure 14 shows the E of the RCSCBC and its corresponding ψ HTR , ψ H , ψ LTR , and ψ L under different m wf . It can be seen from Figure 14 that, given the x p and π, the optimal value of the cyclic E first increases and then decreases with the increase of the m wf . Additionally, it reaches the maximum value when the m wf is in the range of 100~110 kg·s −1 . The changing trends of the four heat conductance distribution ratios with the m wf show that ψ H gradually decreases, ψ L gradually increases, and ψ LTR and ψ HTR fluctuate up and down within a certain range. When the m wf is 70, 80, and 90 kg·s −1 , there are inflection points in the ψ H and ψ LTR . The main reason for these inflection points is that the T 1 needs to exceed the critical temperature of carbon dioxide.  Taking into account the possibility of too many variables causing the curves to be unclear in the picture, the images after four-and five-time optimizations are no longer given, and the results after four-and five-time optimizations are directly given in Table 3. Table 3 shows the results obtained by performing one-, two-, three-, four-and fivetime optimizations in turn with the goal of maximizing E. It can be seen from Table 3 that after the one-time optimization, the E can be improved by 12.13% compared with the IDP. After the two-time optimization, the E can be improved by 31.52%. After the three-time optimizations, the E can be improved by 52.2%. After four-time optimization, the E can be improved by 93.26%. After five-time optimization, the E can be improved by 96.99%. From the changes of cycle parameters in the one-, two-and three-time optimizations, it can be found that when the wf m is large, the ψ LTR is the main factor affecting the E. After the three-time optimization, the wf m corresponding to the optimal value of the E decreases, the ψ LTR decreases, and the ψ HTR increases. This shows that when the wf m is small, the HTR is the main factor affecting the E. In addition, in the three-time optimization, the improvement of the E produces a sudden change. It shows that the heat conductance distribution ratio can make the distribution of cycle performance better, but the wf m plays a much more important role.
Comparing three-, four-and five-time optimizations, it can be found that after releasing the π and wf m , the parameters corresponding to the maximum E have changed significantly, which indicates that the π and wf m are the main factors affecting the E.
In the five-time optimization, the xp corresponding to the best E reaches 0.9, which is the calculation boundary of the selected parameters, which is influenced by other parameters.
Comparing the results of the five-time optimization with those of the four-time optimization, in addition to the increase of the xp, the optimal π corresponding to the optimal E of the cycle also gradually increases, the ψ H increases, the ψ HTR increases, the ψ LTR gradually decreases, and the ψ L increases. The values of all parameters have been Taking into account the possibility of too many variables causing the curves to be unclear in the picture, the images after four-and five-time optimizations are no longer given, and the results after four-and five-time optimizations are directly given in Table 3.  Table 3 shows the results obtained by performing one-, two-, three-, four-and fivetime optimizations in turn with the goal of maximizing E. It can be seen from Table 3 that after the one-time optimization, the E can be improved by 12.13% compared with the IDP. After the two-time optimization, the E can be improved by 31.52%. After the three-time optimizations, the E can be improved by 52.2%. After four-time optimization, the E can be improved by 93.26%. After five-time optimization, the E can be improved by 96.99%. From the changes of cycle parameters in the one-, two-and three-time optimizations, it can be found that when the m wf is large, the ψ LTR is the main factor affecting the E. After the three-time optimization, the m wf corresponding to the optimal value of the E decreases, the ψ LTR decreases, and the ψ HTR increases. This shows that when the m wf is small, the HTR is the main factor affecting the E. In addition, in the three-time optimization, the improvement of the E produces a sudden change. It shows that the heat conductance distribution ratio can make the distribution of cycle performance better, but the m wf plays a much more important role. Comparing three-, four-and five-time optimizations, it can be found that after releasing the π and m wf , the parameters corresponding to the maximum E have changed significantly, which indicates that the π and m wf are the main factors affecting the E.
In the five-time optimization, the x p corresponding to the best E reaches 0.9, which is the calculation boundary of the selected parameters, which is influenced by other parameters.
Comparing the results of the five-time optimization with those of the four-time optimization, in addition to the increase of the x p , the optimal π corresponding to the optimal E of the cycle also gradually increases, the ψ H increases, the ψ HTR increases, the ψ LTR gradually decreases, and the ψ L increases. The values of all parameters have been re-arranged and changed, which further illustrates that the relationship between the various cycle parameters is also mutual influential. Table 4 shows the maximum E of the RCSCBC and its corresponding optimal parameters under different m wf s. According to Table 4, when the m wf is small, the π and x p corresponding to the optimal value of the E are relatively large. When the m wf is larger, the π and x p are smaller. According to the definition of the x p , the larger the x p , the smaller m wf entering the re-compressor. When the x p is one, the RCSCBC is transformed into a RSCBC with two regenerators connected in series. This shows that when the m wf is small, the E of the RSCBC is larger and more advantageous, and when the m wf is large, the E of the RCSCBC is larger, which is more advantageous than the RSCBC. Besides, the optimal π gradually decreases with the increase of m wf and directly affects the ecological performance of the cycle. In actual engineering, appropriate parameters can be selected for different situations according to the data in Tables 3 and 4 to achieve the best ecological performance of the cycle.
Comparing the data in Table 4, it can be found that when the m wf is 65, 75 and 100 kg·s −1 , there are some inflection points in the values of the x p and π.This is due to the influence of the computational boundary, the additional constraint from the T 1 , and the influence of the volatility of the neural network prediction. Considering the functionality of the RCSCBC, the x p cannot take a value of 1, and according to Equation (24), the x p cannot take a value greater than 0.9. Under this premise, when the m wf is small, the calculation of the optimal value of the E is limited.
To sum up, when the m wf is small, the π, heater, and HTR are the main factors influencing the ecological function. At a higher m wf , the effects of the re-compressor, LTR, and cooler on the E increases. According to Table 4, after optimization, the E can be improved by at least 53.74% compared with the IDP. Figure 15 shows the T 1 , T H,out , and the T 4 at each value point in Figure 14. From Figure 15, the T 1 is maintained above 304.42K but has not changed significantly. With the increase of the m wf , both the T H,out and T 4 decrease. The temperature of each state point can provide a reference for the stable operation of the cycle and reflect the accuracy of the optimization results to a certain extent. of the RCSCBC, the xp cannot take a value of 1, and according to Equation (24), the xp cannot take a value greater than 0.9. Under this premise, when the wf m is small, the calculation of the optimal value of the E is limited. To sum up, when the wf m is small, the π , heater, and HTR are the main factors influencing the ecological function. At a higher wf m , the effects of the re-compressor, LTR, and cooler on the E increases. According to Table 4, after optimization, the E can be improved by at least 53.74% compared with the IDP. Figure 15 shows the T1, TH,out, and the T4 at each value point in Figure 14. From Figure  15, the T1 is maintained above 304.42K but has not changed significantly. With the increase of the wf m , both the TH,out and T4 decrease. The temperature of each state point can provide a reference for the stable operation of the cycle and reflect the accuracy of the optimization results to a certain extent.

Conclusions
In this paper, an FTT model of the RCSCBC was established. First, the performance analysis of the E was carried out. Then, on the premise of a certain total heat exchanger inventory, the E as optimized with the wf m , π , xp, ψ H , ψ HTR , and ψ LTR as optimization variables, and the following conclusions were obtained.
(1) The values of η t and η c reflect the irreversibility of the expansion process and the compression process, and the E is used as a thermodynamic index for compromising entropy yield and net power. The irreversibility of the cycle can be measured to a certain extent, that is, there is a correlation between η t , η c , and E. The higher the η c and η t , the smaller the energy loss caused by the irreversibility of the cycle, and the larger the E value.
(2) When the wf m is small, the π , heater and HTR are the main influencing factors on the ecological function. When the wf m is large, the influences of the re-compressor, LTR, and cooler on the E increases. A reasonable adjustment of the distribution ratios

Conclusions
In this paper, an FTT model of the RCSCBC was established. First, the performance analysis of the E was carried out. Then, on the premise of a certain total heat exchanger inventory, the E as optimized with the m wf , π, x p , ψ H , ψ HTR , and ψ LTR as optimization variables, and the following conclusions were obtained.
(1) The values of η t and η c reflect the irreversibility of the expansion process and the compression process, and the E is used as a thermodynamic index for compromising entropy yield and net power. The irreversibility of the cycle can be measured to a certain extent, that is, there is a correlation between η t , η c , and E. The higher the η c and η t , the smaller the energy loss caused by the irreversibility of the cycle, and the larger the E value. (2) When the m wf is small, the π, heater and HTR are the main influencing factors on the ecological function. When the m wf is large, the influences of the re-compressor, LTR, and cooler on the E increases. A reasonable adjustment of the distribution ratios of ψ H , ψ HTR , ψ LTR , and ψ L can make the cycle performance better, but m wf plays a much more important role. (3) Each cycle parameter not only affects the performance of the cycle, but also has a mutual influence relationship. The ecological function can be increased by 12.13%, 31.52%, 52.2%, 93.26%, and 96.99% compared with the IDP after one-, two-, three-, four-, and five-time optimizations.