An Optimization Study to Evaluate the Impact of the Supercritical CO2 Brayton Cycle’s Components on Its Overall Performance

The rising environmental problems due to fossil fuels’ consumption have pushed researchers and technologists to develop sustainable power systems. Due to properties such as abundance and nontoxicity of the working fluid, the supercritical carbon (sCO2) dioxide Brayton cycle is considered one of the most promising technologies among the various sustainable power systems. In the current study, a mathematical model has been developed and coded in Matlab for the recompression of the supercritical carbon dioxide Brayton cycle sCO2-BC. The real gas properties of supercritical carbon dioxide (sCO2) were incorporated into the program by pairing the NIST’s Refporp with Matlab© through a subroutine. The impacts of the various designs of the cycle’s individual components have been investigated on the performance of sCO2 − BC. The impact of various sedative cycle parameters, i.e., compressor’s inlet temperature (T1), and pressure (P1), cycle pressure ratio (Pr), and split mass fraction (x), on the cycle’s performance (ηcyc) were studied and highlighted. Moreover, an optimization study using the genetic algorithm was carried out to find the abovementioned cycle’s optimized values that maximize the cycle’s per-formance under provided design constraints and boundaries.


Introduction
The concept of conversion of heat to mechanical energy is well-known using both internal and external combustion gas cycles. The external combustion gas cycles such as the Rankin cycle have an advantage in their compliance with a wide range of heat sources and many fuels. However, the component design and cycle layouts of the Rankine cycle are quite complicated due to its small molar mass and high critical temperatures that limit the operation of the cycle from low to moderate temperature sources. Closed Brayton cycles that utilize ideal gases as working fluids [1] and can be linked with high-temperature sources, up to 800 • C, have recently been proposed as an alternate to Rankin cycles.
The novel supercritical carbon dioxide Brayton cycle sCO 2 -BC has acquired much interest lately due to its higher conversion efficiency, the compact size of its components, and its more straightforward layouts. It is anticipated that, due to the considerably higher conversion efficiencies linked with the cycle and its smaller layout, it will generate new markets that are expected to be more cost-effective and greener [2,3]. Formally renowned cycles are the steam Rankine cycle and air Brayton cycle. The former is cable of operating at a higher turbine inlet temperature while the latter is known to operate with a requirement of substantially smaller pumping powers. Despite the fact that the former can operate at much higher turbine inlet temperatures, its work ratio is exceedingly small due to huge Several combinations of sCO 2 -BC were analyzed in the literature [3,4,[15][16][17][18] to furth improve its efficiency [3,4,[15][16][17][18]. Among various cycle configurations, the most naïve sCO 2 -BC layout is known as a recuperated/regenerative cycle. It comprises three heat exchangers (a regenrative, an intermediate heat exchanger, and a precooler), an expander, and a compressor. The cycle was initially suggested by Feher [16]. Other variants of the regenerative sCO 2 -BC involve a reheated compression cycle (with an additional heat exchanger between stages of the turbine) and intercooled compression (an intercooler between two stages of the compressor) [15,19].
The single recuperator in the cycle was replaced with two smaller recuperators in series to overcome the pinch-point issues linked with the regenerative cycle. The subsequent cycle is known as the recompression cycle and involves two compressors, i.e., maincompressor and re-compressor. Recompressioned sCO 2 -BC has been referred to as the most efficient layout amount other available layouts of the supercritical cycle [8,20,21]. Numerous studies have been performed to optimize the boundary conditions of sCO 2 -BC. Sarkar and Bhattacharyya [22] performed an optimization analysis for sCO 2 -BC by adding a reheater between the two stages of the turbine. They assessed the effects of various cycles, including the pressure ratio (Pr), compressor's inlet temperature, effectiveness of the hightemperature recuperator (HTR) ( HTR ) and effectiveness of low-temperature recuperator (LTR) ( HTR ). Both reported that the cycle's efficiency can be enhanced considerably by utilizing the reheating technology. Similarly, Reyes et al. [20] and Sharma et al. [23] attempted to optimize the recompression cycle utilizing a central receiver solar plant (CSP) as a heat source and for marine gas-turbine systems. Sarkar et al. [24] investigated the impact of heat exchange and pressure losses on the performance of sCO 2 -BC. Saeed et al. [3] evaluated the recompression cycle using mean line design models of the individual components and also evaluated impacts of different channel geometries of PCHEs on the performance of the cycle [25,26]. Saad et al. [27] conducted an exergy destruction assessment to explore the effects of different ambient and water temperatures on the performance of sCO 2 -BC Saeed and Kim developed an optimization study for the sCO 2 -BC's expander system.
From the studies mentioned above, it can be seen that even though numerous investigations have been conducted to analyze sCO 2 -BC, there is no optimization study available in the literature based on a genetic algorithm. Until now, only search algorithms have been employed. In this context, the current study provides a detailed analysis of sCO 2 -BC. A mathematical model for sCO 2 -BC has been developed and coded in Matlab©. Properties of sCO 2 have been integrated into the code by linking the NIST's Refprop to Matlab using an in-house code. The effects of the various components on the overall performance of sCO 2 -BC have been investigated, which has never been reported in previous studies to the authors' knowledge. Moreover, the genetic algorithm was employed to optimize the cycle parameters to maximize their efficiency.

Power Cycle Modeling
Copious arrangements of sCO 2 -BC have been suggested and analyzed [15][16][17]28]. Feher proposed the most straightforward version of the cycle known as the recuperated or regenerative cycle. Dostal [28] proposed and investigated simple recuperated and recompression cycles, while Angelino [17] proposed a partial cooling cycle. Other layouts include intercooled sCO 2 -BC, and partial recompression. For the current work, recompression sCO 2 -BC was chosen on the basis of its excellent performance [8,20,21] and more straightforward layout [28].
As explained earlier and shown in Figure 1, the cycles consist of an additional recuperator and compressor compared with the regenerative cycle. The recompression sCO 2 -BC consists of four heat exchanges (two recuperators, i.e., HTR and LTR (for heat recuperation), an external heat exchanger (to exchange heat from the source to the working fluid), and a precooler (for additional heat rejection to complete the cycle) [8,29,30].

Turbomachinery Models
As presented in [18], turbomachinery modeling is based on a single equation model of the turbine and compressor utilizing isentropic efficiencies and pressure ratio information, whereby isentropic efficiencies are approximated using the Balje chart [31].
2.1.2. Recuperator Model The cycle boundary conditions were chosen carefully through an extensive review of the literature. The selected cycle's boundary conditions are provided in Table 2, as a function of on which all simulations were evaluated through an in-house code written in the Matlab environment. As presented in [18], turbomachinery modeling is based on a single equation model of the turbine and compressor utilizing isentropic efficiencies and pressure ratio information, whereby isentropic efficiencies are approximated using the Balje chart [31].

Recuperator Model
There have been different approaches adopted in the literature to model the heat exchangers and recuperators. Most common methodologies are based on either U.A. values or effectiveness values . Recuperator [27,32] modeling techniques based on the effectiveness of the heat exchanger are computationally less costly. Thus, the NTU method that utilizes the heat exchanger's effectiveness was adopted for the current study, as given by the following equations.
The definitions of C min,LTR and C min,HTR are given by Equations (7) and (8).
C min,HTR = min(Cp h,HTR , Cp c,HTR ) C min,LTR = min(Cp h,HTR , xCp c,HTR ) It was assumed that both recuperators, i.e., the high-temperature recuperator (HTR) and low-temperature recuperator (LTR), are entirely insulated. Therefore, the resulting energy balance equations can be presented in the form provided below.
where x is the fraction of mass flowing toward the recompressor, termed the split mass fraction. It can be calculated using the following relationship: Considering ideal mixing, the mixing value was modeled using Equation (10).

of 17
To estimate the pressure losses, Equation (12) was considered. The definition of the cycle efficiency used in the current work was provided by Equation (13).
Definitions of w t , w c , and w rc are given by the following equations: The cycle's selected conditions listed in Table 2 are from the literature [8]. Cycle calculations under selected boundaries ( Table 2) were computed using a Matlab code for various values of x. The values in Table 2 were approximated from the literature [4,7,33].  As shown in Figure 2, the code initializes with known boundary conditions ( Table 2). The turbine and main compressor were solved first using Equation (1) and Equation (3), respectively, and later HTR was solved based on an assumed value of T 6 based on which T 4 was computed. Finally, T 4 was used to solve recompressor and LTR. The above-mentioned process keeps iterating in a loop until T 6 is converged.

Cycle Optimization Using a Genetic Algorithm
For the current study, four parameters were considered for the optimization process, i.e., T 1 , P 1 , Pr, and x. The most frequently adopted optimization algorithm in the literature are genetic algorithms [34][35][36][37][38][39][40] and particle swarm optimizations [41][42][43]. Both genetic algorithms (G.A.s) and particle swarm optimization (PSO) are iteration-based and start with a pool of initial values and heuristic algorithms; however, the former uses discrete data and the latter uses continuous data. It is reported in the literature that both algorithms are extremely accurate, but PSO offers a relatively higher efficiency. On the other hand, G.A.s are more useful for constraint satisfaction problems. Therefore, the method opted for optimizing the present problem is a multiobjective genetic algorithm (MOGA). The design variables mentioned above with their high and low limits are displayed in Table 3. The wellknown genetic algorithm has frequently been adopted in the literature [34][35][36][37][38][39][40] to work out global maximization and minimization problems. MOGA works on the principle known as survival of the fittest. Usually, MOGA is provided by several combinations of the design variables ( Table 3) that are generated randomly, known as populations. The population size can be decided on the basis of the best practices found in the literature [41]. The MOGA tested each individual from the populations through a fitness function table. Usually, the top 1-5% individuals in the population are moved to the next generation unchanged, while the remaining population is generated through mutations and crossovers as displayed in the Figure 3. The process continues unless stopping criteria are met. The Matlab environment and a population size of 50 were used for the current optimization problem, as recommended in the literature [18,29] for five or fewer design variables. The rank and stochastic uniform of built functions were used for the scaling and selection. Elite size, crossover, and mutation fractions were used as 5%, 80%, and 15%, respectively. A forward migration factor of 0.2 was used and augmented. The Lagrangian function was used as a nonlinear constraint algorithm. All simulations were computed on a system with a Xeon processor and 128GB Ram. Computation time consumed by an optimization iteration was found to be around 7 hours.

Results and Discussions
The current study was conducted to analyze and optimize the recompression of the supercritical carbon dioxide cycle (sCO2-BC). A mathematical model was based on abrupt variations in the thermophysical properties of the sCO2 The cycle was developed by coupling NIST's Refprop with the main cycle simulation code. Initially, the effects of the various cycle components such as the recuperator, compressor, and turbine performance were investigated on the overall performance of the sCO2-BC. Later, the effects of critical operating parameters on the cycle's performance were analyzed, i.e., , , Pr, and . Lastly, these parameters were optimized using a genetic algorithm to maximize the overall performance of the sCO2-BC.
For the validation of the code, current work was compared with the available data in the literature, as listed in the Table 4. This showed the computed results using the employed model are in close agreement with the published data.

Results and Discussions
The current study was conducted to analyze and optimize the recompression of the supercritical carbon dioxide cycle (sCO 2 -BC). A mathematical model was based on abrupt variations in the thermophysical properties of the sCO 2 The cycle was developed by coupling NIST's Refprop with the main cycle simulation code. Initially, the effects of the various cycle components such as the recuperator, compressor, and turbine performance were investigated on the overall performance of the sCO 2 -BC. Later, the effects of critical operating parameters on the cycle's performance were analyzed, i.e., T 1 , P 1 , Pr, and x. Lastly, these parameters were optimized using a genetic algorithm to maximize the overall performance of the sCO 2 -BC.
For the validation of the code, current work was compared with the available data in the literature, as listed in the Table 4. This showed the computed results using the employed model are in close agreement with the published data. It can be seen from Figure 1 that the main flow stream after the high-temperature recuperator (HTR) divided into two creeks, i.e., 10a and 10b. The flow toward 10a with a split mass fraction "x" was directed to the central compressor after passing through the cooler, while fraction "1 − x" of the flow passed through the recompression compressor. This was done to avoid early pinching within the low-temperature recuperator. Split streams, which were 3 (after gaining heat from the LTR) and 5 (from the re-compressor) combined again, forming stream 4. There is an optimized value of "x" for any set of operating conditions for maximum cycle efficacy. Figure 4 shows how split mass fractions affect the cycle's performance and its components. It can be visualized that cycle efficiency increased at first with the rise in the value of split mass fraction "x", reaching a maximum value and then decreasing. Thus, only a single value of "x" exists for which the cycle's efficiency is maximum. The split mass fraction's optimized value corresponds to the minimum temperature difference between streams 3 and 5, as shown in Figure 5. Moreover, Figure 5 shows a variation in the values of T 3 , T 5 , and cycle efficiency with split mass fraction. Cycle efficiency was maximum for the split mass fraction value of 0.82, corresponding to a minimum difference, as shown by the dotted line. The effects of the split mass fraction on the component's power can also be seen in Figure 4. It can be seen here that the primary compressor's power increased with the surge in the value of split mass fraction, while that of the recompression compressor and turbine decreased. split mass fraction "x" was directed to the central compressor after passing through the cooler, while fraction "1 − x" of the flow passed through the recompression compressor. This was done to avoid early pinching within the low-temperature recuperator. Split streams, which were 3 (after gaining heat from the LTR) and 5 (from the re-compressor) combined again, forming stream 4. There is an optimized value of "x" for any set of operating conditions for maximum cycle efficacy. Error! Reference source not found. shows how split mass fractions affect the cycle's performance and its components. It can be visualized that cycle efficiency increased at first with the rise in the value of split mass fraction "x", reaching a maximum value and then decreasing. Thus, only a single value of "x" exists for which the cycle's efficiency is maximum. The split mass fraction's optimized value corresponds to the minimum temperature difference between streams 3 and 5, as shown in Figure 5. Moreover, Figure 5 shows a variation in the values of , , and cycle efficiency with split mass fraction. Cycle efficiency was maximum for the split mass fraction value of 0.82, corresponding to a minimum difference, as shown by the dotted line. The effects of the split mass fraction on the component's power can also be seen in Error! Reference source not found.. It can be seen here that the primary compressor's power increased with the surge in the value of split mass fraction, while that of the recompression compressor and turbine decreased.

Effect of Cycle's Components' Performance on the Cycle's Overall Performance
Figures 6-9 reveal the impact of the components' performances on the performance of sCO2-BC. Figure 6 shows the variations in the cycle's efficiency ( ) with a split mass fraction for numerous combinations of effectiveness of the high-temperature recuperator and low-temperature recuperator. A significant increase in the cycle's efficiency was observed, with an increase in value effectiveness for both recuperators. Specifically, an increase in the effectiveness from 0.9 to 0.98 increases the cycle's efficiency by 22.5%. Moreover, it can be seen that an increase in the recuperator's performance affected the split mass fraction's optimal value. The value decreased marginally with a rise in the effectiveness values of the recuperator. Figure 7 shows a modification of the cycle's efficiency with a split mass fraction for various turbine efficiency values. It was observed that turbine performance significantly affects the cycle's performance; a 2% increase in the turbine's performance increased the cycle's performance by 1%. A slight decrease in the values of optimized split mass fraction can also be seen in Figure 7, with an increase in the turbine's efficiency. This reflects that load on the main compressor will decrease with an increase in the value of the turbine's efficiency. It should also be noted that the sensitivity of the cycle's performance to the effectiveness of the recuperator increases with an increase in the value of the split mass fraction as the curve diverges along with the increasing value of . Additionally, the sensitivity of to decreases with the growth in x. Error! Reference source not found. shows the variation in with the disparity in the evaluation of x for various values of and . Results suggest that curves exist with various = values. Similar trends are displayed by curves in Error! Reference source not found. for different values of . Cycle efficiency significantly increased with an upsurge in the compressor's efficiency. The value of the optimal split mass fraction shifted toward lower values of the split mass fraction with an upsurge in the value of . Moreover, the sensitivity of with deceased with an increase in the value of the split mass fraction. Error! Reference source not found. shows a comparison of how performance parameters associated with the cycle components affected the overall cycle's performance. It can be seen that the cycle's performance was most vulnerable to the recuperator's performance, followed by the performance of the turbomachin-   Figure 6 shows the variations in the cycle's efficiency (η cyc ) with a split mass fraction for numerous combinations of effectiveness of the high-temperature recuperator and low-temperature recuperator. A significant increase in the cycle's efficiency was observed, with an increase in value effectiveness for both recuperators. Specifically, an increase in the effectiveness from 0.9 to 0.98 increases the cycle's efficiency by 22.5%. Moreover, it can be seen that an increase in the recuperator's performance affected the split mass fraction's optimal value. The value x decreased marginally with a rise in the effectiveness values of the recuperator. Figure 7 shows a modification of the cycle's efficiency with a split mass fraction for various turbine efficiency values. It was observed that turbine performance significantly affects the cycle's performance; a 2% increase in the turbine's performance increased the cycle's performance by 1%. A slight decrease in the values of optimized split mass fraction can also be seen in Figure 7, with an increase in the turbine's efficiency. This reflects that load on the main compressor will decrease with an increase in the value of the turbine's efficiency. It should also be noted that the sensitivity of the cycle's performance to the effectiveness of the recuperator increases with an increase in the value of the split mass fraction as the η cyc curve diverges along with the increasing value of x. Additionally, the sensitivity of η cyc to η t decreases with the growth in x. Figure 8 shows the variation in η cyc with the disparity in the evaluation of x for various values of η c and η rc . Results suggest that η cyc curves exist with various η c = η rc values. Similar trends are displayed by η cyc curves in Figure 7 for different values of η t . Cycle efficiency significantly increased with an upsurge in the compressor's efficiency. The value of the optimal split mass fraction shifted toward lower values of the split mass fraction with an upsurge in the value of η t . Moreover, the sensitivity of η cyc with η c deceased with an increase in the value of the split mass fraction. Figure 9 shows a comparison of how performance parameters associated with the cycle components affected the overall cycle's performance. It can be seen that the cycle's performance was most vulnerable to the recuperator's performance, followed by the performance of the turbomachinery. Appl. Sci. 2021, 11, x FOR PEER REVIEW 11 of 19 Figure 6. Impact of recuperator's effectiveness on the cycle's efficiency.

Effect of the Cycle's Parameters on Its Performance
The cycle's performance is highly dependent on its selected value parameters. Cycle components are designed on the basis of these parameters. Nevertheless, the selection of cycle parameters under required conditions is tricky and requires extensive cycle simulations. The most sensitive parameters for the recompression sCO2-BC are compressor inlet temperature ( ) compressor inlet pressure ( ), split mass fraction ( ), and cycle pressure ratio. The influence of the abovementioned cycle parameters on the cycle's performance was analyzed in the current study, followed by its optimization using the genetic algorithm. Figures 10-13 show the dependence of the cycle's split mass fraction and the cycle's efficiency on the parameters of the cycle-i.e., inlet temperature ( ) compressor inlet pressure ( ) and cycle's pressure ratio. Figure 10 displays a response of the surface bounded by the two-cycle operating parameters, compressor inlet pressure, and cycle's pressure ratio. Response surface shows that the cycle's efficiency is more sensitive to the pressure ratio than the compressor inlet pressure. The cycle efficiency, which is dependent on the pressure ratio, is relatively prominent at higher values of pressure ratio when the compressor's inlet pressure is close to the critical value of the pressure for sCO2-i.e., 7.34MPa. On the other hand, the cycle's efficiency is equally sensitive to the cycle's pressure ratio for all values of compressor inlet pressure. The cycle's efficiency increases initially when the rise in the value of the reaches a maximum value and then drops with further increase in the value of the cycle's pressure ratio. Thus, there exists an optimal value of the pressure ratio corresponding to each value of the compressor's inlet pressure ( ). This value shifts towards a lower value of the cycle's pressure ratio with an increase in the value of the compressor's inlet pressure ( 1 ). Figure 11 shows variations of the optimized split mass fraction with the compressor's inlet pressure ( 1 ) and cycle's overall pressure ratio. The graph shows that the optimized split mass fraction is overly sensitive to the compressor's inlet pressure for all the pressure ratio values. The split mass fraction decreases initially with an increase in the value of ( 1 ) then it starts increasing with further increase in the value of ( 1 ). However, was also significantly affected by . Initially, decreased rapidly with an increase in the value of although this variation leveled off at higher values of .

Effect of the Cycle's Parameters on Its Performance
The cycle's performance is highly dependent on its selected value parameters. Cycle components are designed on the basis of these parameters. Nevertheless, the selection of cycle parameters under required conditions is tricky and requires extensive cycle simulations. The most sensitive parameters for the recompression sCO 2 -BC are compressor inlet temperature (T 1 ) compressor inlet pressure (P 1 ), split mass fraction (x), and cycle pressure ratio. The influence of the abovementioned cycle parameters on the cycle's performance was analyzed in the current study, followed by its optimization using the genetic algorithm. Figures 10-13 show the dependence of the cycle's split mass fraction and the cycle's efficiency on the parameters of the cycle-i.e., inlet temperature (T 1 ) compressor inlet pressure (P 1 ) and cycle's pressure ratio. Figure 10 displays a response of the surface bounded by the two-cycle operating parameters, compressor inlet pressure, and cycle's pressure ratio. Response surface shows that the cycle's efficiency is more sensitive to the pressure ratio than the compressor inlet pressure. The cycle efficiency, which is dependent on the pressure ratio, is relatively prominent at higher values of pressure ratio when the compressor's inlet pressure is close to the critical value of the pressure for sCO 2 -i.e., 7.34 MPa. On the other hand, the cycle's efficiency is equally sensitive to the cycle's pressure ratio for all values of compressor inlet pressure. The cycle's efficiency increases initially when the rise in the value of the Pr reaches a maximum value and then drops with further increase in the value of the cycle's pressure ratio. Thus, there exists an optimal value of the pressure ratio corresponding to each value of the compressor's inlet pressure (P 1 ). This value shifts towards a lower value of the cycle's pressure ratio with an increase in the value of the compressor's inlet pressure (P 1 ). Figure 11 shows variations of the optimized split mass fraction with the compressor's inlet pressure (P 1 ) and cycle's overall pressure ratio. The graph shows that the optimized split mass fraction is overly sensitive to the compressor's inlet pressure for all the pressure ratio values. The split mass fraction decreases initially with an increase in the value of (P 1 ) then it starts increasing with further increase in the value of (P 1 ). However, x opt was also significantly affected by Pr. Initially, x opt decreased rapidly with an increase in the value of Pr although this variation leveled off at higher values of Pr.  Error! Reference source not found. shows the variation in with compressor inlet temperature ( ) and cycle overall pressure ratio ( ). It dictates that the cycle's effi-

Optimization of sCO2-BC Parameter
The above-discussed results reveal that the cycle's performance was susceptible to the operating parameters associated with sCO2-BC, i.e., compressor inlet temperature ( ), compressor inlet pressure ( 1 ), cycle pressure ratio , and split mass fraction. An optimization procedure was required to compute the optimal values of the operating conditions under which cycle components would be designed in later stages. For the current study, the design variables selected for the optimization study are listed in Error! Reference source not found.. A genetic algorithm was adopted for the optimization study, the details of which are available in Section 3. The cycle simulation code explained in Section 2 of the manuscript was used as a fitness function for the optimization procedure. The

Optimization of sCO2-BC Parameter
The above-discussed results reveal that the cycle's performance was susceptible to the operating parameters associated with sCO2-BC, i.e., compressor inlet temperature ( ), compressor inlet pressure ( 1 ), cycle pressure ratio , and split mass fraction. An optimization procedure was required to compute the optimal values of the operating conditions under which cycle components would be designed in later stages. For the current study, the design variables selected for the optimization study are listed in Error! Reference source not found.. A genetic algorithm was adopted for the optimization study, the details of which are available in Section 3. The cycle simulation code explained in Section 2 of the manuscript was used as a fitness function for the optimization procedure. The Figure 13. Deviation of x with T 1 and cycle overall pressure ratio (Pr). Figure 12 shows the variation in η cyc with compressor inlet temperature (T 1 ) and cycle overall pressure ratio (Pr). It dictates that the cycle's efficiency was quite sensitive to the pressure ratio for all values of the compressor's inlet temperature (T 1 ). However, its sensitivity to the pressure ratio was more prominent when the compressor's inlet temperature was close to the critical temperature value of CO 2 . For all values of T 1 , η cyc increased initially and reached a maximum value before dropping with a further increase in the value of Pr. Thus, there exists an optimized value of the Pr corresponding to each value of T 1 that shifted toward the lower value of Pr at higher values of T 1 . On the other hand, η cyc was significantly affected by the value of T 1 . At any fixed value of r, η cyc declined with the growth in the value of T 1 . It can be seen that η cyc was sensitive to compressor inlet temperature escalations with a rise in Pr. Moreover, the cycle's efficiency was significantly higher when the compressor's inlet temperature was close to the critical temperature T c that decreased rapidly as T 1 moved away from T c Figure 13 represents a response surface bounded by T 1 and Pr. It shows that x opt was influenced by both the compressor inlet temperature and the cycle pressure ratio. The value of x opt decreased initially and then, after a certain pressure ratio for which its value was minimum, it started increasing. On the other hand, the value of x opt increased rapidly as the compressor inlet temperature moved away from T c ; however, at higher values of T 1 , the rate of increase in the value of x opt became smaller. Thus, it can be depicted from here that selection of T 1 is crucial for the main and re-compressor sizes.

Optimization of sCO 2 -BC Parameter
The above-discussed results reveal that the cycle's performance was susceptible to the operating parameters associated with sCO 2 -BC, i.e., compressor inlet temperature (T 1 ), compressor inlet pressure (P 1 ), cycle pressure ratio Pr, and split mass fraction. An optimization procedure was required to compute the optimal values of the operating conditions under which cycle components would be designed in later stages. For the current study, the design variables selected for the optimization study are listed in Table 3. A genetic algorithm was adopted for the optimization study, the details of which are available in Section 3. The cycle simulation code explained in Section 2 of the manuscript was used as a fitness function for the optimization procedure. The convergence of the genetic algorithm is shown in Figure 14. The genetic algorithm was converged in 90 iterations, as shown in Figure 14. For the current optimization procedure, two cases have been studied. Case I has been defined and optimized with no constraints, while for case II one constraint has been defined-i.e., T 1 > 310. This constraint has been defined due to abrupt variations in the properties of CO 2 near the critical point that can lead the system to instability during off-design conditions. Results of both cases have been listed in Table 5. convergence of the genetic algorithm is shown in Figure 14. The genetic algorithm was converged in 90 iterations, as shown in Figure 14. For the current optimization procedure, two cases have been studied. Case I has been defined and optimized with no constraints, while for case II one constraint has been defined-i.e., > 310. This constraint has been defined due to abrupt variations in the properties of CO2 near the critical point that can lead the system to instability during off-design conditions. Results of both cases have been listed in Table 5.

Conclusions
In the current study, a mathematical model for sCO2-BC was developed and coded in Matlab©. An in-house code was used to couple NIST's Refprop with Matlab© to evaluate the properties of sCO2. The effects of various cycle components were evaluated on the cycle's overall performance. Moreover, the sensitivity of its operating conditions was studied with respect to its performance. Lastly, an optimization study using a genetic algorithm was carried out to optimize the cycle's operating conditions and maximize the cycle efficiency. The following deductions were made from the study:

•
The performance of the cycle modules significantly affected the overall performance of the cycle. A 2% variation in the efficiency of turbines or compressors changed the cycle's performance by 1%. Simultaneously, a 2% variation in the effectiveness of both recuperators affected the cycle's performance by 3%. Thus, the cycle's efficiency was most sensitive to the performance of the recuperator.

•
The sensitivity of the cycle's performance to the recuperator's effectiveness increased with a rise in . Simultaneously, it decreased with an escalation in the values of the isentropic efficiencies of the turbomachinery components.

Conclusions
In the current study, a mathematical model for sCO 2 -BC was developed and coded in Matlab©. An in-house code was used to couple NIST's Refprop with Matlab© to evaluate the properties of sCO 2 . The effects of various cycle components were evaluated on the cycle's overall performance. Moreover, the sensitivity of its operating conditions was studied with respect to its performance. Lastly, an optimization study using a genetic algorithm was carried out to optimize the cycle's operating conditions and maximize the cycle efficiency. The following deductions were made from the study:

•
The performance of the cycle modules significantly affected the overall performance of the cycle. A 2% variation in the efficiency of turbines or compressors changed the cycle's performance by 1%. Simultaneously, a 2% variation in the effectiveness of both recuperators affected the cycle's performance by 3%. Thus, the cycle's efficiency was most sensitive to the performance of the recuperator.

•
The sensitivity of the cycle's performance to the recuperator's effectiveness increased with a rise in x. Simultaneously, it decreased with an escalation in the values of the isentropic efficiencies of the turbomachinery components.

•
The value of η cyc was susceptible to the cycle's operating parameters if the compressors' inlet conditions were close to carbon dioxide's critical point. Consequently, the optimization results suggested placing the compressor away from the critical point to ensure better off-design performance. Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the privacy policy of the organization.

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