Next Article in Journal
Experimental Study of the Feasibility of In-Situ Hydrogen Generation from Gas Reservoir
Previous Article in Journal
A Comparative Study of Fuzzy SMC with Adaptive Fuzzy PID for Sensorless Speed Control of Six-Phase Induction Motor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic, Exergoeconomic and Multi-Objective Analyses of Supercritical N2O-He Recompression Brayton Cycle for a Nuclear Spacecraft Application

1
School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
2
Institute of Spacecraft System Engineering, China Academy of Space Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(21), 8184; https://doi.org/10.3390/en15218184
Submission received: 8 October 2022 / Revised: 28 October 2022 / Accepted: 31 October 2022 / Published: 2 November 2022
(This article belongs to the Section J2: Thermodynamics)

Abstract

:
Detailed thermodynamic, exergoeconomic, and multi-objective analysis are performed for a supercritical recompression Brayton cycle in which the advanced working medium mixture of nitrous oxide and helium (N2O–He) is utilized for power generation. The thermodynamic and exergoeconomic models are propitious based on the standard components’ mass and energy conservation, exergy balance equation, and exergy cost calculation equation. An investigation of the sensitivity parametric is considered for judging the impact of crucial decision variable parameters on the performance of the proposed Brayton cycle. The proposed cycle’s performance is evaluated by systematic analysis of the thermal efficiency (ηth), exergy efficiency (ηex), total cost rate ( C . total), levelized cost of electricity (LCOE), and the total heat transfer area (Atotal). Furthermore, multi-objective optimization is adopted from the viewpoint of the first and second laws of exergoeconomics to find the optimum operating parameters and to improve the circular’s exergoeconomic performance. The final results illustrate that the optimization calculation is based on the fact of the exergoeconomics method; the whole system produces electrical power of 0.277 MW with C . total of USD 18.37/h, while the ηth, ηex, Atotal, and LCOE are 49.14%, 67.29%, 165.55 m2 and USD 0.0196/kWh, respectively. It is concluded that the work exergy destruction for the reactor and turbine is higher than that of other components; then, after the multi-objective optimization analysis, the ηth and ηex improved by 2.08% and 5.07%, respectively, and the C . total, Atotal, and LCOE decreased by 13.99%, 0.01%, and 5.13%, respectively.

1. Introduction

Concerns over restrictions of mass, mission time, and long-distance interstellar flight suggest an essential role for the nuclear electric propulsion (NEP) system in future deep space, even extrasolar exploration. The NEP system is simple electric propulsion in which the electricity is derived from a nuclear reactor. The thrusting fluid (gas, ion, plasma) receives its energy to accelerate from the application of electric fields, either directly heating the gas as in an arc jet, using an electric field to accelerate ions, or a magnetic field as in the magnetoplasmadynamic (MPD) approach. In order to successfully utilize nuclear power, many studies have been carried out to improve thermoelectric conversion efficiency, be lightweight, and reduce the power generation cost of spacecraft. High-efficiency thermoelectric conversion devices include static and dynamic conversion methods [1].
Static methods:
  • Thermocouple conversion system;
  • Thermionic conversion system.
Dynamic methods:
  • Brayton conversion system;
  • Rankine conversion system;
  • Stirling conversion system.
The Brayton cycle has higher thermoelectric conversion efficiency than other conversion methods. It can convert approximately 20% to 35% of the whole thermal energy transported by the reactor part into electrical energy through dynamic mechanical processes [2]. Recently, the recompression Brayton cycle power generation has become a research frontier due to its advantages of compact and higher thermal efficiency [3] with a simple structure [4]. The recompression Brayton cycle is shown in Figure 1. To select suitable supercritical cycle structures, a lot of research was carried out during the period, with over 40 structures studied [5]. In 2015, Sulaiman’s research found that a system with a re-compressor has advanced performance compared to the simple system, regeneration system, pre-compression system, partial cooling system, segmented expansion system, and other structures [6]. The cycle with a re-compressor can weaken the influence of the “pinch point” of the cycle that may occur in the recuperator and regenerator, and the supercritical cycle is found to be more efficacious [7,8].
For the supercritical Brayton cycle, Sarkar [9] studied the exergy and optimization of the cycle and investigated the influence of pressure ratio, cycle energy, and thermal efficiency. Sarkar and Bhattacharyya [10] optimized the supercritical recompression Brayton cycle, achieving the maximum efficiency and the thermal efficiency improved by 3.5% when used reheating. Wang et al. [11] adopted the genetic algorithm (GA) method for the optimization of cycle parameters to obtain the maximum efficiency (thermal efficiency and exergy efficiency). P. Serrano et al. [12] designed a new arrangement scheme in which a bypass low-temperature cladding source replaced one of the regenerators. In the new arrangement scheme, high system efficiency can be achieved at the high inlet temperature of the turbine with medium pressures.
A revolutionary supercritical working medium, nitrous oxide and helium (N2O–He), was performed for the proposed recompression Brayton cycle [13] to convert thermal energy into electricity. In the N2O–He mixture’s components, the N2O part has thermodynamic properties such as molecular weight, critical pressure, and temperature, causing almost comparable behaviors under operating conditions that may improve the N2O system’s performance and can use the existing CO2 configuration. Sarkar J. [14] utilized supercritical N2O (S-N2O) as the working medium to combine the Brayton cycle with the new generation nuclear reactors. Its thermal efficiency is higher than that of the S-CO2 system. Characteristics of several critical working fluids at the critical point are shown in Table 1. Using these mixtures (CO2-He [15], CO2-C2H6/SF6 [16], binary mixture [17]) can further improve the thermodynamic cycling performance of different Brayton cycles. The mixed working medium Brayton cycle research has apparent advantages over the pure S-CO2 Brayton cycle.
The research mentioned above is devoted primarily to thermodynamic analysis, and we know that comprehensive exergoeconomic analyses work on the recompression supercritical N2O–He mixture Brayton cycle has not been reported yet. Moreover, some literature merely concerns the exergoeconomic analyses of the S-CO2 Brayton system. A detailed exergoeconomic investigation and thermoeconomic analysis needs to also be conducted.
In this study, energy, exergy, and exergoeconomic theory of the recompression Brayton cycle are analyzed using the novelty working medium N2O-He mixture with the mole fractions of N2O and He of 85% and 15% (85%N2O–15%He). Firstly, comprehensive energy and exergy studies are analyzed, and the exergy cost theory is applied to the cycle. Moreover, a parameter sensitivity analysis is presented to reveal the influence of the proposed cycle on the thermal and exergy efficiency and total cost of the nuclear-powered spacecraft. Finally, from the perspective of exergoeconomics, the supercritical N2O–He Brayton cycle is optimized, and the appropriate parameters of various components’ thermal efficiency (ηth), exergy efficiency (ηex), total cost rate ( C . total), total heat transfer area (Atotal), and levelized cost of electricity (LCOE) are selected by the GA method.

2. System Description and Assumptions

An NEP system with a typical Brayton cycle unit can be separated into six subsystems, depicted in Figure 2. The NEP system’s six subsystems are the reactor core, shielding, power conversion unit, heat dissipation, power management and distribution (PMAD), and electric propulsion (EP). The supercritical 85%N2O–15%He recompression Brayton conversion unit consists of a heat source reactor (R), a TB (turbine), an MC (main compressor), an RC (recompression compressor), a G (generator), an HTR (high temperature recuperator), an LTR (low temperature recuperator), a cooler, and a radiator. The recompression Brayton cycle schematic diagram is presented in Figure 3. The T-s diagram for the proposed corresponding S-(N2O-He) recompression Brayton cycle is shown in Figure 4.
The supercritical N2O-He stream leaving the high temperature reactor enters the TB as a working fluid (1–2), where it expands and produces power (2–3). Then, the expansion flow successfully flows through the HTR and LTR to heat the recompression flow (3–4–5). The air flows out of the LTR and is split into two streams (state 5). One stream,5a, has a high mass flow rate, and the other stream, 5b, has a low mass flow rate. The partial flow of 5a is cooled in the cooler component (5a–6) and then compressed in the MC to higher pressure (6–7). The stream flowing out of the MC absorbs heat from the LTR (7–8a) and is mixed with the steam, which is compressed by RC to the same pressure as MC (5b–8b). The two partial split streams are mixed at state point 8 and then enter into the HTR and are heated before flowing into the reactor core part (8–1). In this way, a complete thermodynamic cycle is formed.
The space nuclear reactor system of the actual recompression Brayton cycle includes the nuclear reactor core, high temperature recuperator, low temperature recuperator, heat accumulators, panel-type radiation radiators, turbine machinery, main compressor and re-compressor machinery, heat pipe and other complex and auxiliary equipment. The equipment structure is complex in the actual operation, and the operating condition is variable, so the system needs to be simplified. After the system is simplified, the cycle process can be better modeled to explore its cycle performance.
The study of this paper is conducted on the general assumptions as below:
  • All processes reach the stable state.
  • Turbine and compressor processes are adiabatic.
  • The pressure loss of the heat exchanger is negligible.
  • Kinetic energy and potential energy and heat transfer to the environment of the system are neglected.
  • The efficiency of the turbine (TB) and compressors (MC and RC) is isentropic efficiencies.
  • The effectiveness of the HTR and the LTR is considered.
According to the above assumptions, the actual recompressed Brayton cycle space nuclear energy system is simplified to the recompressed closed Brayton cycle system, which becomes the most basic thermal cycle process composition.

3. Methods and Model Design

3.1. Thermodynamic Analysis

This work extends previous research [13,19], particularly concerning each state point’s thermodynamic and exergy analysis. Table 2 lists calculation equations for each component’s energy, exergy, and mass balance.
The study of the Brayton cycle is based on the theory of mass, energy, and exergy equilibrium, which is shown as follows:
m i . = m e .
Q . W . = m e . h e m i . h i
E Q . + E . i = E . W + E . e + E . D e
where subscript i represents the value of the entry state point, and e represents the value of the exit state point. E . De is the exergy destruction rate in the component.
Stream exergy consists of four main parts, which are E . ph (stream physical exergy), E . ch (stream chemical exergy), E . k (stream kinetic exergy), and E . p (stream potential exergy). Regardless of the E . k and E . p of the working medium, total exergy of the working medium E . total is modified as follows [20]:
E . t o t a l = E . p h + E . c h
In the work at this stage, there is a chemical reaction process of the working medium in the cycle, and the working fluid passes through different system components, so only consider the working medium’s physical exergy. The calculation formula can be defined as [20]:
E . p h = m . × [ ( h h 0 ) T 0 ( s s 0 ) ]
In the equation, the subscript “0” indicates the value of environmental conditions. T0 is space environment temperature. The environmental conditions for the proposed cycle are T0 = 298 K and P0 = 0.101 MPa.
The effectiveness for each recuperator is defined and can be estimated as [19].
ε H T R = T 3 T 4 T 3 T 8
ε L T R = T 4 T 5 T 4 T 7
The thermodynamic model for the HTR and the LTR are given by [21]:
h 3 h 4 = h 1 h 8
m . 4 ( h 8 a h 7 ) = m . 7 ( h 4 h 5 )
The efficiency expressions for the TB, MC, and RC can be calculated by [14]:
η T B = h 2 h 3 h 2 h 3 s
η M C = h 7 s h 6 h 7 h 6
η R C = h 8 s h 5 h 8 h 5
The subscript “s” is the isentropic process.
The thermodynamic calculation relations of the TB, MC, and RC are as follows [14]:
W . T B = m . 2 ( h 3 h 2 )
W . M C = m . 7 ( h 7 h 6 )
W . R C = m . 5 b ( h 8 h 5 )
Specific heat input of the reactor ( Q . R) is expressed as:
Q . R = m . 1 ( h 2 h 1 )
The net power of circulating, Wnet, is defined as:
W . n e t = W . T B ( W . M C + W . R C )
The system thermal efficiency can be calculated by:
η t h = W . n e t Q R .
The system exergy efficiency is calculated by:
η e x = W . n e t E . i n = W . n e t Q . R ( 1 T 0 / T R )
where E . in is the exergy entering the reactor, and TR is the reactor temperature.
The details of the relations for energy and exergy balance equations for each component are provided in Appendix A. Energy and exergy balance calculation equations are listed in Table A1.

3.2. Exergoeconomic Analysis

A cost-effective energy conversion system is vital for the cycle, so economic consideration should be considered. Exergoeconomic analysis needs to consider the theory of the exergy cost [22], the average cost analysis method [23], specific exergy cost calculation [24,25,26], and basic economic theory principles. The purpose of exergy economic analysis is to study and calculate the exergy cost of unit product flow by revealing the system’s cost formation process. For system components that receive thermal energy and generate electricity, the equation of cost balance can be expressed as follows [27]:
C . o u t , j + C . W , j = C . i n , j + C . Q , j + Z . j
In Equation (20), C . Q,j is the cost rate related to component input heat, and C . W,j is the cost rate of component output power. It can be seen from Equation (20) that the total cost rate of working medium exiting exergy is the sum of the total cost rate of working medium entering exergy and the total cost expenditure rate.
The calculation formula of the unit exergy cost rate is as follows [28]:
C . j = c j × E . j
C . W = c W × E . W
C . Q = c Q × E . Q
In Equations (21)–(23), cj, cW, and cQ are the per unit of cost exergy, $/GJ.
Z . j in calculation Formula (20) is the total cost rate related to the capital investment (CI) cost and operation and maintenance (OM) cost of the jth system component [20]:
Z . j = Z . j C I + Z . j O M
The component j’s annual levelized CI can be calculated by [27,29]:
Z . j C I = ( C R F τ ) × Z j
In Equation (25), the CRF represents the capital recovery factor, and τ is the service life of the power system.
The annual leveling of the OM cost is calculated according to [27]:
Z . j O M = ( γ j τ ) × Z j
Table 3 lists the factors CRF, τ, γk, and Zj of the system’s different components.
We update the Chemical Engineering Plant Cost Index (CEPCI) of equipment costs to 2021 values.
Z 2021 = Z o r i g i n a l × C E P C I 2021 C E P C I o r i g i n a l
where CEPCI2021 = 596.2 [30], and Zoriginal is each component’s original cost. The Zk calculation details for each component of the proposed cycle are listed in Appendix B.
The cost balance for the different components of the proposed cycle and the auxiliary equations required are all presented in Table A2. The Gauss–Seidel method is used to solve the various linear equations in Table A3. As shown in Table A4, several essential indicators are used to identify economic behavior. The total cost rate ( C . total) is given in Table A4, which is the sum of these values, including the CI cost rate, the total exergy destruction cost rate ( C . D), and the OM cost rate.

3.3. Performance Evaluation Index of the Cycle

W . net, ηth, ηex, C . total, Atotal, and LOCE are the evaluation metric for improving the operating system performance.
  • Net power ( W . net) is shown as Equation (17).
  • Thermal efficiency (ηth) is calculated by Equation (18).
  • Exergy efficiency (ηex) can be calculated by Equation (19).
  • Total cost rate ( C . total) is calculated by the equation in Table A4.
  • The price of electricity can be estimated using the levelized cost of electricity (LCOE) parameter. It is given as [31]:
L C O E = Z . t o t a l W . n e t ( $ kW h )

3.4. Parameter Selection of Sensitivity Analysis

The parameter selection of sensitivity analysis aims to analyze the influence of various working conditions and design parameters on system performance. In this research work, significant parameters are selected to study their impact on the system performance, including split ratio—x, pressure ratio—π, MC inlet temperature—T6, TB inlet temperature—T2, and MC inlet pressure—P6. The constraints are presented in Table 2. For the MW supercritical recompression Brayton cycle system, the split ratio, the pressure ratio, the MC inlet temperature, the TB inlet temperature, and the MC inlet pressure are important factors determining the cycle process.
For the compressor, the inlet temperature and pressure are two main factors directly affecting its work. The compressor pressure ratio also has a significant influence on the system cycle. When the pressure ratio is too high, the power consumption of the compressor will increase, and the cycle efficiency will decrease. When the pressure ratio is too small, the heat absorption of the working medium will be negligible, and the output work of the turbine will be further affected. Considering that the critical pressure of the mixed working medium is about 7.4 MPa, the data point near the critical point is selected. Therefore, the inlet pressure of the compressor is set as 7.4–8.4 MPa, and the outlet pressure of the compressor is required to be 20–30 MPa, so the pressure ratio is selected as 2.5–3.5. The inlet temperature of the turbine is related to the output work of the turbine and has a great influence on the system. In recent years, with material technology development, the turbine inlet temperature can be raised to 1423 K [32], so the turbine inlet temperature is selected as 1000–1200 K. According to relevant references, the split ratio is generally 0.25–0.35 [33].

3.5. Calculation Process of Optimization

To build a high-efficiency and economical energy system, we must comprehensively consider the balance between efficiency (energy efficiency and exergy efficiency) and total cost rate. The objective of the present work is to maximize the efficiency (ηth and ηex) and minimize C . total, Atotal, and LCOE. The genetic algorithm (GA) [34] method optimizes the decision parameters selected to change the system performance. Figure 5 presents the optimization flowchart of the proposed system.

4. Results and Discussion

In this section, detailed thermodynamic and exergoeconomic analyses are presented. The thermodynamic properties of N2O and He can be calculated using the NIST REFPROP [35] software. As the NIST REFPROP software shows, the temperature of N2O is limited to 525 K, and to estimate the thermodynamic properties of N2O, a thermodynamic property subroutine “PUREN2OPROP” was developed according to the previous literature report [36]. The performance of the system under different working conditions is evaluated using Excel and Python software. The basic parameters of the cycle used are listed in Table 3.
Table 3. Basic parameters selected in simulation.
Table 3. Basic parameters selected in simulation.
Parameter SymbolSelect Value
T0 (K)298
P0 (MPa)0.101
P6 (MPa)7.4 [36]
m (kg/s)1
Fuel cost ($/MWh)7.4 [31]
ηTB0.88 [13]
ηMC0.85 [9]
ηRC0.85 [9]
εHTR0.86 [9]
εLTR0.86 [9]

4.1. Energy, Exergy, and Exergoeconomic Investigation on Basic Condition

Table 4 shows the energy and exergy analysis results of the proposed cycle under typical working conditions, and the parameters are shown in Table 3. In the typical cycle, W . net = 0.277 MW, ηth = 48.14%, and ηex = 64.04%. The thermodynamic characteristics, mass flow rate, and cost values of state points under the proposed cycle based on the state point numbers in Figure 2 are provided in Appendix C.
The exergy and exergoeconomic analyses of components are listed in Table 5. Table 5 indicates that the reactor component has a higher Z . j+ C . D,j + C . L,j than other components, indicating that it should improve the reactor’s exergoeconomic performance. The exergoeconomic factor value for the reactor is 56.94%, meaning that the sum of CI and OM in the reactor has an excellent ratio. The exergoeconomic factor of the turbine is 8.04%, and the value of Z . j+ C . D,j + C . L,j is the second highest value. The cooler has a higher relative cost difference value than other cycle components. The cooler’s exergoeconomic factor is 31.54%, which reveals that the majority cost is the exergy destruction cost in the system’s total cost. Comprehensively considering overall cost values and reducing the exergy destruction of the cooler is an excellent method to change the cycle performance. The total f value of the whole cycle is 24.21%. It means that the exergy destruction reached 75.79% of the total system cost.

4.2. Sensitivity Analysis of Main Operating Parameters

4.2.1. Influence of Split Ratio—x

Figure 6 presents the influence of x (split ratio) on ηth, ηex, C . total, Atotol, and LCOE. The cooling effect of the cooler and compressor power is changed by changing the split ratio with the mass flow rate change.
Figure 6a,b depict the cycle thermal and exergy efficiencies, respectively. From Figure 6a,b, under different T6, the thermal and exergy efficiencies have the same variation trend with the increase of x, which increases slowly at first and then decreases. The ηth and ηex are between 45.0–48.3% and 59.9–64.3%. With the increase of T6, the ηth and ηex decrease gradually. When x = 0.35, the efficiency (ηth and ηex) reaches maximum with ηth = 48.34% and ηex = 64.3%.
Referring to Figure 6c, the total cost rate ( C . total) first decreased and then increased slightly; the whole variation trend decreases gently. An increase in the split ratio leads to the reduction of C . D,total. The minimum value of C . total can be obtained at x = 0.33. With the increase of T6, the C . total decreases, but the variation trend is not apparent. When T6 = 315 K and x > 0.37, C . total increases gradually with x increases. The main reason for the change in C . total is that the working temperature of 315 K is near the critical point area, and the thermodynamic characteristics of the working medium at this temperature change dramatically.
The impact of the split ratio on Atotal (the sum of the HTR heat exchange area, LTR heat exchange area, and radiator heat exchange area) is presented in Figure 6d, and the central part is the radiator area. If we change the split ratio, the total area of the heat exchange components increases with the increase of x, and the area also tends to increase with the increase of T6. Nevertheless, as the x < 0.31, the total area has little difference. When T6 is at the temperature of 315 K, the total heat exchanger area at a high split ratio x > 0.37 slightly decreases, which is associated with the physical properties of the working medium.
Figure 6e displays the variation trend of the LCOE. Under the different MC inlet temperature-T6, with the x increasing, the LCOE decreased gently and at x = 0.33 achieved the lowest LCOE; after that split ratio, the value of LCOE increased slightly.
When x is small, a “pinch point” occurs at the low temperature inlet side of the HTR; when x increases, RC mass flow rate and T5 increase, reducing the total work consumed by MC and RC. The rate of increase in net power is faster than the rate of decrease in heat absorption; then, ηth and ηex are increased. With the increase of x, the position of the “pinch point” changes, which appears at the low temperature inlet side of LTR, and the power consumption of the RC increases alongside the increase of x. As the Wnet decreases, ηth and ηex decrease. With the increase of T6 (MC inlet temperature), the compression power consumption of the whole cycle increases, resulting in the decrease of ηth and ηex. Comprehensively considering the whole cycle, when x < 0.33, the whole Wnet increases, and when x > 0.33, the whole Wnet decreases. The changing trend of Wnet results in the LCOE decreasing at first and increasing modestly.

4.2.2. Influence of Pressure Ratio—π

Figure 7 presents the cycle primary performance under the operating condition of x = 2.5–3.5, T2 = 1200 K, x = 0.3, P6 = 7.4 MPa, and T6 = 315–340 K.
The MC inlet temperature, the compressor power consumption, and the turbine’s work increase with a π increase. When the TB inlet temperature is constant at 1200 K, as the π increases, Wnet first increases and then decreases. Moreover, the LTR outlet temperature was increased due to reduced turbine outlet temperature, increased outlet temperature in the MC, and reduced heat transfer in the recuperators. As the π increases, so does Wnet. Therefore, ηth and ηex first increase and then decrease with the increase of π, as presented in Figure 7a,b.
The results in Figure 7c show that an increase in π leads to an increase in C . total of the entire cycle, mainly due to an increase in C . D,reactor. In addition, with the increase of π, the C . D,HTR decreases and C . D,LTR increases, and the total exergy destruction cost rate of recuperators increases due to the reduction of recuperators’ heat transfer area. Z . j (CI and OM costs) of TB, MC, and RC increase by π. When π increases, the total system cost rate increases, mainly because of the rising cost of CI and OM, and C . total. With the increase of π, the reactor and recuperators’ exergy destruction cost rate increases with T6.
According to Figure 7d, increasing π can reduce the heat transfer area of recuperators and radiators due to a decrease in heat exchange. The heat transfer area of all three components decreases, so the whole cycle heat transfer area is reduced by increasing π.
The change of π also affects the LCOE of the system, as plotted in Figure 7e. The output power of TB and the whole cycle net power are increasing by increasing the π. Increasing π up to around 3.0 reduces the LCOE of the system, which is mainly because the system’s net power value increases significantly. When π increases beyond 3.0, the net power value increase is less than the value of π below 3.0. Total investment and OM cost ( Z . j) of the whole cycle is increased by π increasing. However, the total cycle net power is the dominant reason for LCOE reduction when π is below 3.0, and LCOE increases in the period of π beyond 3.0, mainly due to CI and OM cost increase.

4.2.3. Influence of MC Inlet Temperature—T6

The influence of T6 on the cycle key performance indicators is shown in Figure 8. Figure 8a–c illustrate ηth, ηex, and C . total decrease as T6 increases. However, C . total is not sensitive to the change of T6. The main reason for cyclic ηth and ηex reduction is the difference between the maximum and minimum cyclic temperatures (Carnot principle). A higher value of T6 leads to a higher average temperature of N2O-He at the cooler, creating a more significant heat loss in the cooler, a higher compression work value, and a lower cycle net power. By increasing T6, the C . D,j of the reactor and the cooler is reduced. At the same time, the increase in the HTR and LTR heat transfer area will increase the Z . j value of the recuperators. Therefore, the increase of HTR and LTR cost rate is less than the decrease in the reactor and the cooler exergy destruction cost. To sum up, C . total of the circulating system is reduced.
A variation trend of the total heat transfer area has little influence on the increase of the MC inlet temperature, as presented in Figure 8d. Recuperators’ heat transfer area is increased with the increase of T6; nevertheless, the radiator heat transfer area decreases. The heat transfer area of the radiator decrease is greater than the increase of that in the recuperators. Therefore, Atotal is slightly reduced.
The change of T6 also influences the LCOE of the cycle, as depicted in Figure 8e. The total investment, OM costs ( Z . j), and whole system net power are decreased by T6 increasing. However, the whole LOCE increases with a T6 increase. The main reason for such a change is that the decrease in Wnet is more significant than that of total investment and OM cost.

4.2.4. Influence of TB Inlet Temperature—T2

The influence of T2 on the performance of the proposed S-(N2O-He) cycle with x = 0.3, P6 = 7.4 MPa, π = 3.0, and T6 = 315–340 K is shown in Figure 9.
With the increase of T2, the TB power output produced by the input rate of specific heat increases significantly due to the increase of TB input enthalpy. Insensitive to the effects of MC and RC work leads to an increase in Wnet of the system. Therefore, the temperature difference between the maximum and minimum operating temperatures increases as T2 increases, and ηth increases as shown in Figure 9a. The enthalpy of inflow and outflow of the reactor increases with the increases of T2. Cycle exergy decreases with the increase of maximum operating temperature due to thermal exergy ( E . in) increase greater than net power increases, resulting in reduced exergy efficiency, as shown in Figure 9b. C . total is decreased, mainly caused by the decrease in C . D,total, and the CI and OM costs of the reactor are decreased by the increase in maximum operating temperature, as displayed in Figure 9c.
Accordingly, the heat transfer of the radiator, HTR, and LTR increase with the TB inlet temperature, leading to a higher heat transfer area of these components in the system, as shown in Figure 9d. Referring to Figure 9e, LCOE reduces because the net power increases as T2 increases.

4.2.5. Influence of MC Inlet Pressure—P6

Figure 10 indicates the cycle performance variation at different P6, T2 = 1200 K, x = 0.3, and π = 3.0.
ηth and ηex decrease with increasing P6, as presented in Figure 10a and Figure 10b, respectively. ηth and ηex decrease slowly and then decrease rapidly. As the working medium is worked close to the critical point area, the compression factor is small, and the consumption power of the compressor is increased, resulting in a decrease in Wnet. The compression factor increases gradually as P6 increases, resulting in more pressure loss, higher power consumption, and a lower net power of the cycle, thus reducing the ηth and ηex, as shown in Figure 10a,b.
Figure 10c presents the influence of P6 on the C . total. As the CI and OM cost of the reactor is reduced, the total investment and OM cost is reduced, which is the highest value part of the Z . total. C . D,Reactor and C . D,MC are also reduced. and the two components with the more significant proportion. All decreases in Z . total and C . D,total are a result of decreasing C . total.
Atotal is increased by the increase in P6, as presented in Figure 10d. As P6 increases, heat transfer in LTR and HTR increases. As a result, the heat transfer area of recuperators is also increased. The heat transfer area variation of HTR and LTR is the dominant factor causing the increase of Atotal.
As mentioned above, Wnet and OM costs decrease with the increase in P6. W . net (net power) decreases more than Z . total, thereby allowing a high value of LCOE due to an increase in MC inlet pressure.

4.3. Exergoeconomic Optimization

Based on the sensitivity analysis discussed above, ηth and ηex improvements are always inconsistent with LCOE. However, in the actual thermodynamic system, efficiency (ηth and ηex) and cycle cost are considered simultaneously. Multi-objective optimization is used to maximize thermal performance and minimize cost. For the recompression Brayton cycle, its optimization objectives are thermal efficiency ηth = f1(x), exergy efficiency ηex = f2(x), total cost rate C . total = f3(x), total area of heat transfer Atotal = f4(x), and levelized cost of electricity LCOE = f5(x).
The response surface methodology (RSM) [21] method is utilized to make a link between the objective functions and decision parameters for optimization. The regression equation is established by experimental design and variance analysis. Finally, based on the Pareto optimality, multi-objective optimization is carried out to seek the balance between efficiency and cost.
Appropriate sample points are selected using a central composite design (CCD) [21] to approximate the response function with higher accuracy. In the CCD method, a central design point needs to be added, and the distance from the center is α. The selected decision parameters are the five test factors shown in Table 6. The Design-Expert software is used to create second-order model regression for the corresponding values of the loop.
Therefore, the multi-objective optimization problem of the whole loop can be defined as:
Max[f1(x), f2(x)]; Min[f3(x), f4(x), f5(x)].
The regression model functions of ηth = f1(x), ηex = f2(x), C . total = f3(x), Atotal = f4(x), and LCOE = f5(x) are expressed as:
f 1 ( x ) = η t h = 342.95946 1.00182 × A + 0.207314 × B + 102.6095 × C + 18.70625 × D + 233.2025 × E 0.000034 × A B + 0.0132 × A C 0.0512 × A D + 0.96 × A E 0.0115 × B C + 0.0561 × B D 0.2555 × B E 5.68 × C D 30.1 × C E 23.5 × D E + 0.001108 × A 2 0.000085 × B 2 4.3875 × C 2 3.9575 × D 2 + 83.58333 × E 2
f 2 ( x ) = η e x = 186.98794 4.05654 × A 0.025545 × B + 85.97158 × C + 67.343 × D + 1105.8 × E + 0.000136 × A B + 0.1256 × A C 0.0684 × A D 1.42 × A E 0.01515 × B C + 0.0595 × B D 0.252 × B E 8.63 × C D 63.5 × C E 31 × D E + 0.00522 × A 2 2.85417 × 10 6 × B 2 4.27417 × C 2 8.06417 × D 2 + 406.25 × E 2
f 3 ( x ) = C . t o t a l = 458.17302 1.60636 × A 0.004646 × B 40.88086 × C + 18.39035 × D 161.8025 × E 0.000275 × A B + 0.14 × A C + 0.0664 × A D + 0.792 × A E 0.001725 × B C 0.01184 × B D + 0.165 × B E 1.974 × C D 34.5 × C E 5 × D E + 0.00052 × A 2 + 0.00004 × B 2 + 0.770917 × C 2 1.62475 × D 2 15.94167 × E 2
f 4 ( x ) = A t o t a l = 120.48618 + 0.278564 × A + 0.009831 × B 2.801 × C + 6.75212 × D 46.054 × E + 0.000048 × A B + 0.0136 × A C 0.01036 × A D 0.028 × A E + 0.001272 × B C + 0.00016 × B D + 0.001 × B E 0.1989 × C D + 5.3 × C E + 38 × D E 0.000598 × A 2 0.000016 × B 2 0.24075 × C 2 2.49298 × D 2 151.70167 × E 2
f 5 ( x ) = L C O E = 0.1 * ( 33.32761 + 0.126413 × A + 0.006422 × B + 1.82567 × C + 0.731167 × D + 7.04317 × E + 2.0 × 10 9 × A B + 0.000043 × A C + 0.000034 × A D 0.00056 × A E + 4.6 × 10 6 × B C 0.000053 × B D + 0.000201 × B E + 0.000148 × C D + 0.005 × C E + 0.031 × D E 0.000193 × A 2 2.89721 × 10 6 × B 2 0.117088 × C 2 0.116345 × D 2 12.02983 × E 2 )
where A—MC inlet temperature T6, B—TB inlet temperature T2, C—MC inlet pressure P6, D—compressor pressure ratio π, and E—split ratio x.

4.3.1. Dual Objective Optimization of ηth and ηex

Firstly, the system’s thermodynamic performance is considered by using the dual-objective optimization method. Taking the cycle’s ηth—f1(x) and ηex— f2(x) as the response targets, the functions are defined in Equations (29) and (30).
In the process of dual-objective optimization, the initial population size is 1000, the crossover factor is 0.7, and the variation factor is 0.3. The whole Pareto front of ηth and ηex has 300 data in Figure 11 and partial data are shown in Table 6. After considering the ηth and ηex, we select the final optimization result.

4.3.2. Three Objectives Optimization of ηth and ηex, and C . total

Figure 12 and Table 7 present the three objectives optimization results of ηthf1(x), ηexf2(x), and C . totalf3(x). The optimized functions are Equations (29), (30), and (31).
In the three objectives optimization, the initial population size is 1000, the crossover factor is 0.7, and the variation factor is 0.3. Figure 12a–c show the Pareto optimal boundary, which ηth and ηex, ηth and C . total, and ηex and C . total are the optimization objectives, respectively. Figure 12d considers the Pareto solution of ηth, ηex, and C . total.

4.3.3. Three Objectives Optimization of ηth, ηex, and Atotal

The three objectives optimization between ηthf1(x), ηexf2(x), and Atotalf4(x) are shown in Figure 13 and Table 8. The optimized functions are Equations (29), (30), and (32).
In the three objectives optimization, the initial population size is 1000, the crossover factor is 0.7, and the variation factor is 0.3. Figure 13a–c show the Pareto optimal boundary with ηth and ηex, ηth and Atotal, and ηex and Atotal as the optimization objectives. Figure 13d considers the Pareto solution of the three objectives of ηth, ηex, and Atotal.

4.3.4. Three Objectives Optimization of ηth, ηex, and LCOE

The result of three objectives, optimization of ηthf1(x), ηexf2(x), and LCOEf5(x), are shown in Table 9. Figure 14 depicts the Pareto frontier for the three objectives optimization. The optimized functions are Equations (29), (30), and (33).
In the three objectives optimization, the initial population size is 1000, the crossover factor is 0.7, and the variation factor is 0.3. Figure 14a–c shows the Pareto optimal boundary, with ηth and ηex, ηth and LCOE, and ηex and LCOE as the optimization objectives, respectively. Figure 14d considers the Pareto solution of three objectives: ηth, ηex, and LCOE. In practical application, suitable operating parameters can be selected according to specific conditions and requirements.
According to this part, multi-optimization and comprehensiveness consider the optimization results, and the final vital parameters of the proposed S-(N2O-He) cycle are selected and detailed in Table 10.
Table 11 shows the optimization operation parameters and performance metrics after the optimization process, while other conditions are under the same primary operating conditions.

5. Conclusions

A supercritical recompression Brayton cycle with 85% N2O–15% He as the working medium is used for detailed thermodynamic analysis, exergoeconomic analysis, and parametric optimization. The sensitivity analysis and study of cycle parameters are carried out to study the influence of critical thermodynamic parameters on efficiency (ηth and ηex), C . total, Atotal, and LCOE. Then, dual-objective and multi-objective function optimization is performed out of the loop. The multi-objective is based on the Pareto frontier, selecting ηth, ηex, C . total, Atotal, and LCOE as the objective functions. At the same time, the essential parameters studied are selected as decision variables. For this work, to evaluate the thermodynamic and exergoeconimic performance, the results can be listed as follows:
  • The exergy destruction rates in the reactor, TB, and RC are higher than the value in other components.
  • The reactor and turbine have larger values of Z . j + C . D,j + C . L,j than other components in the cycle, which can be effectively reduced with a split ratio increase, TB inlet temperature, and MC inlet pressure.
  • Compared with the initial basic design, ηth and ηex improved by 2.08% and 5.07% when optimized based on the first and second laws of thermodynamics.
  • When optimized based on the exergy economics, the sum of the system C . total is reduced by 13.99%.
  • When the cycle operates on the optimized conditions, the net power increases by 60.35%. Based on the critical parametric analysis, the total heat transfer area and LCOE decrease by 0.01% and 5.13%.
Considering the particularity of the space system, it is necessary to focus on the mass characteristics of the nuclear power system and use the radiator model, exergoeconomic model, and levelized cost of electricity model to make a quantitative evaluation of its volume characteristics, economic cost characteristics, and power generation cost characteristics. The parametric sensitivity analysis concluded that there is always a contradiction between the thermal and economic performance of the system. Therefore, a multi-objective optimization method based on the Pareto frontier is used to determine the optimal operating condition and design parameters. Minimizing the total cost and maximizing the overall thermal and exergy efficiency reduces the total product unit cost and LCOE price. The proposed 85% N2O–15% He recompression Brayton cycle could produce net power of 0.364 MW with a thermal efficiency of 49.14% at a LOCE price of USD 0.0196/kWh. The thermal efficiency is increased by 15.66%, and the LCOE price is reduced by 52.3% compared with the S-CO2 recompression Brayton cycle under the same operating conditions. Accordingly, the performance of the proposed cycle is promising and can offer a potential working fluid selection for the space nuclear power system in the near future.

Author Contributions

Conceptualization, X.M. and H.Z.; formal analysis, X.M. and Y.X.; investigation, X.M.; methodology, X.M., Q.W. and W.S.; software, X.M.; writing—original draft, X.M.; funding acquisition, H.Z. and Y.X.; supervision, H.Z. and Y.X.; data curation, Q.W. and W.S.; resources, Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China, Grant number 2020YFB1901800.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Energy and Exergy Balance Equations of the System Components

For the energy and exergy analysis, the calculation equations for each system component are listed in Table A1.
Table A1. Energy balance and exergy balance calculation equations of circulating components.
Table A1. Energy balance and exergy balance calculation equations of circulating components.
ComponentsEnergy Conservation Equation Exergy   Destruction   ( E . F E . P )
Reactor Q . R = m . 1 ( h 2 h 1 ) E . D , R = ( E . 1 + E . Q R ) E . 2
Turbine W . T = m . 2 ( h 3 h 2 )
η T = h 2 h 3 h 2 h 3 s
E . D , T = ( E . 2 E . 3 ) W T .
Main Compressor W . M C = m . 7 ( h 7 h 6 )
η M C = h 7 s h 6 h 7 h 6
E . D , M C = W . M C ( E . 7 E . 6 )
Recompression Compressor W . R C = m . 5 b ( h 8 h 5 )
W . R C = m . 5 b ( h 8 h 5 )
E . D , R C = W . R C ( E . 8 b E . 5 b )
HTR ε H T R = T 3 T 4 T 3 T 8
h 3 h 4 = h 1 h 8
E . D , H T R = ( E . 3 E . 4 ) ( E . 1 E . 8 )
LTR ε L T R = T 4 T 5 T 4 T 7
m . 4 ( h 8 a h 7 ) = m . 7 ( h 4 h 5 )
E . D , L T R = ( E . 4 E . 5 ) ( E . 8 a E . 7 )
Cooler m . 5 a ( h 5 h 6 ) = m . 9 ( h 10 h 9 ) E . D , C o o l e r = ( E . 5 a E . 6 ) ( E . 10 E . 9 )

Appendix B. Investment Cost of the System Components

For the exergoeconomic analysis, several calculation equations are used to estimate the investment costs of each system’s main components. The investment cost is expressed as a function of the appropriate thermodynamic variables. The investment cost equations for each component of the proposed system are given in Table A2.
Table A2. CI cost equations and economic value of each component [37,38].
Table A2. CI cost equations and economic value of each component [37,38].
ComponentsCI Cost Equation (Zoriginal)
Reactor Z R = c R * Q . R ,   c R = 283   $ / k W
Turbine Z T B = 479.34 * m . 2 * [ 1 / ( 0.93 η T ) ] * l n ( π ) * [ 1 + e x p ( 0.036 * T 2 54.4 ) ]
Main Compressor Z M C = 71.1 * m . 6 * [ 1 / ( 0.92 η M C ) ] * π * l n ( π )
Recompression Compressor Z R C = 71.1 * m . 5 b * [ 1 / ( 0.92 η M C ) ] * π * l n ( π )
HTR Z H T R = 2681 * A H T R 0.59
LTR Z L T R = 2681 * A L T R 0.59
Cooler Z C o o l e r = 2143 * A C o o l e r 0.514
FactorValue
Number of system operation years (n)20
Annual system working hours (τ)8000
Interest rate of system (ir)12%
Maintenance factor (ik)0.06
Capital recovery factor (CRF) C R F = i r * ( 1 + i r ) n ( 1 + i r ) n 1
The cost analysis for the circulatory system’s different components of the cycle and auxiliary equations required for the exergy analysis are all presented in Table A3. As shown in Table A4, several vital indicators are used to identify economic behavior.
Table A3. Exergy balance for system different components and required auxiliary equations.
Table A3. Exergy balance for system different components and required auxiliary equations.
ComponentsExergy Balance EquationAuxiliary Equation
Reactor C . 2 = C . F u e l + C . 1 + Z . R non
Turbine C . 3 + C . T B = C . 2 + Z . T B C . 3 E . 3 = C . 2 E . 2 , or c 2 = c 3
Main Compressor C . 7 = C . M C + C . 6 + Z . M C C . W M C W . M C = C . W T B W . T B , or c M C = c T B
Recompression Compressor C . 8 b = C . R C + C . 5 b + Z . R C C . W R C W . R C = C . W T B W . T B , or c R C = c T B C . 5 b = x C . 5
HTR C . 1 + C . 4 = C . 8 + C . 3 + Z . H T R C . 3 E . 3 = C . 4 E . 4 , or c 3 = c 4 C . 8 = C . 8 a + C . 8 b
LTR C . 5 + C . 8 a = C . 7 + C . 4 + Z . L T R C . 4 E . 4 = C . 5 E . 5 , or c 4 = c 5
Cooler C . 6 + C . 10 = C . 9 + C . 5 a + Z . C o o l e r C . 6 E . 6 = C . 5 a E . 5 a , or c 6 = c 5 a C . 9 = 0 , C . 5 a = ( 1 x ) C . 5
Table A4. Exergoeconomic evaluation parameters [26,27].
Table A4. Exergoeconomic evaluation parameters [26,27].
Average exergy cost per unit of fuel c F , j = C . F , j / E . F , j
Average exergy cost per unit of product c P , j = C . P , j / E . P , j
Relative cost difference r j = ( c P , j c F , j ) / c F , j
Exergy destruction cost rate C . D , j = c F , j × E . D , j
Exergy economic factor f j = Z . j / ( Z . j + C . D , j + C . L , j )
Total cost rate C . t o t a l = j Z . j + j C . D , j

Appendix C. Thermodynamic Characteristics of Cycle State Point

Table A5 lists the thermodynamic characteristics, mass flow rate, and cost values of state points under the proposed cycle based on the state point numbers in Figure 2.
Table A5. Thermodynamic characteristics at different cycle state point.
Table A5. Thermodynamic characteristics at different cycle state point.
State PointT (K)P (MPa) m .   ( kg / s ) h (kJ/kg)s (kJ/kgK) E .   ( MW ) Costs
C .   ( $ / h ) c ($/GJ)
1955.722.211790.96.131.19837.8488.772
2120022.212366.26.611.63147.9888.173
3943.27.411765.46.110.95928.2148.173
4860.97.411544.15.891.02330.1078.173
5528.17.41985.25.080.70620.7858.173
63157.40.7568.24.070.4138.4995.720
7473.922.20.7889.14.890.46816.4199.745
8847.522.211569.65.921.04139.4410.518
92980.11.5629.46.26000
10437.10.11.58246.850.0252.57228.509
W . TB------20.5219.488
W . MC------7.6739.488
W . RC------5.9899.488

References

  1. Zhang, H.; Yin, D.; Chai, X.; Kong, B.; Liu, X. Thermoelectric conversion performance of combined thermoions system for space nuclear power supply. Front. Energy Res. 2020, 7, 167. [Google Scholar] [CrossRef]
  2. Moran, M.J.; Shapiro, H.N.; Boettner, D.D.; Bailey, M.B. Fundamentals of Engineering Thermodynamics, 9th ed.; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2018. [Google Scholar]
  3. Hu, H.; Liang, S.; Jiang, Y.; Guo, C.; Guo, Y.; Zhu, Y.; Cai, H. Thermodynamic and exergy analysis of 2 MW S-CO2 Brayton cycle under full/partial load operating conditions. Energy Convers. Manag. 2020, 211, 112786. [Google Scholar] [CrossRef]
  4. Iverson, B.D.; Conboy, T.M.; Pasch, J.J.; Kruizenga, A.M. Supercritical CO2 Brayton cycles for solar-thermal energy. Appl. Energy 2013, 111, 957–970. [Google Scholar] [CrossRef] [Green Version]
  5. Crespi, F.; Gavagnin, G.; S’anchez, D.; Martínez, G.S. Supercritical carbon dioxide cycles for power generation: A review. Appl. Energy 2017, 195, 152–183. [Google Scholar] [CrossRef]
  6. Al-Sulaiman, F.A.; Atif, M. Performance Comparison of Different Supercritical Carbon Dioxide Brayton Cycles Integrated with a Solar Power Tower. Energy 2015, 82, 61–71. [Google Scholar] [CrossRef]
  7. Dostal, V. A Supercritical Carbon Dioxide Cycle for Next Generation Nuclear Reactors. Ph.D. Thesis, Czech Technical University in Prague, Prague, Czech Republic, 2004. [Google Scholar]
  8. Marchionni, M.; Bianchi, G.; Tsamos, K.M.; Tassou, S.A. Techno-economic comparison of different cycle architectures for high temperature waste heat to power conversion systems using CO2 in supercritical phase. Energy Procedia 2017, 123, 305–312. [Google Scholar] [CrossRef]
  9. Sarkar, J. Second law analysis of supercritical CO2 recompression Brayton cycle. Energy 2009, 34, 1172–1178. [Google Scholar] [CrossRef]
  10. Sarkar, J.; Bhattacharyya, S. Optimization of recompression S-CO2 power cycle with reheating. Energy Convers. Manag. 2009, 50, 1939–1945. [Google Scholar] [CrossRef]
  11. Wang, J.; Sun, Z.; Dai, Y. Parametric optimization design for supercritical CO2 pow-er cycle using genetic algorithm and artificial neural network. Appl. Energy 2010, 87, 1317–1324. [Google Scholar] [CrossRef]
  12. Serrano, I.P.; Linares, J.I.; Cantizano, A.; Moratilla, B.Y. Enhanced arrangement for rec-uperators in supercritical CO2 Brayton power cycle for energy conversion in fusion reactors. Fusion Eng. Des. 2014, 89, 1909–1912. [Google Scholar] [CrossRef]
  13. Miao, X.Y.; Zhang, H.C.; Sun, W.B.; Wang, Q.; Zhang, C.X. Optimization of a recompression supercritical nitrous oxide and helium Brayton cycle for space nuclear system. Energy 2022, 242, 123023. [Google Scholar] [CrossRef]
  14. Sarkar, J. Thermodynamic analyses and optimization of a recompression N2O Brayton power cycle. Energy 2010, 35, 3422–3428. [Google Scholar] [CrossRef]
  15. Zhang, X.; Zhang, H.C.; Liu, X.T.; You, E.S.; Yin, D.Z. Performance Analysis and Optimization of a Closed Regenerative Brayton Cycle for Space Nuclear Power System. In Proceedings of the 25th National Academic Conference on Engineering Thermophysics of Colleges and Universities, Roorkee, India, 28–31 December 2019. [Google Scholar]
  16. Wright, S.A.; Pickard, P.S.; Vernon, M.E.; Radel, R.F. Enhancing Power Cycle Efficiency for a Supercritical Brayton Cycle Power System Using Tunable Supercritical Gas Mixtures. U.S. Patent 20130033044 A1, 7 February 2013. [Google Scholar]
  17. Hu, L.; Chen, D.; Huang, Y.; Li, L.; Cao, Y.; Yuan, D.; Wang, J.; Pan, L. Investigation on the performance of the supercritical Brayton cycle with CO2-based binary mixture as working fluid for an energy transportation system of a nuclear reactor. Energy 2015, 89, 874–886. [Google Scholar] [CrossRef]
  18. Lemmon, E.W.; Bell, I.H.; Huber, M.L.; Mc Linden, M.O. REFPROP database. In Reference Fluid Thermodynamic and Transport Properties. NIST Standard Reference Database 23, Version 10.0; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2013. [Google Scholar] [CrossRef]
  19. Miao, X.Y.; Zhang, H.C.; Zhang, D.; Zhang, C.X.; Huang, Z.L. Properties of Nitrous Oxide and Helium mixtures for space nuclear recompression Brayton cycle. Energy Rep. 2022, 8, 2480–2489. [Google Scholar] [CrossRef]
  20. Wang, X.R.; Yang, Y.; Zheng, Y.; Dai, Y.P. Exergy and exergoeconomic analyses of a supercritical CO2 cycle for a cogeneration application. Energy 2017, 119, 971–982. [Google Scholar] [CrossRef]
  21. Morgon, R. Comparison of High Power Nuclear Systems for Space. In Proceedings of the 18th Intersociety Energy Conversion Conference, Orlando, FL, USA, 21–26 August 1983. [Google Scholar]
  22. Lozano, M.A.; Valero, A. Theory of the exergetic cost. Energy 1993, 18, 939–960. [Google Scholar] [CrossRef]
  23. Lazzaretto, A.; Tsatsaronis, G. A General Process-Based Methodology for Exergy Costing; American Society of Mechanical Engineers: New York, NY, USA, 1996; Volume 36, pp. 413–423. [Google Scholar]
  24. Shokati, N.; Mohammadkhani, F.; Yari, M.; Mahmoudi, S.M.S.; Rosen, M.A. A comparative exergoeconomic analysis of waste heat recovery from a gas turbine-modular helium reactor via organic Rankine cycles. Sustainability 2013, 6, 2474–2489. [Google Scholar] [CrossRef] [Green Version]
  25. Mohammadkhani, F.; Shokati, N.; Mahmoudi, S.M.S.; Yari, M.; Rosen, M.A. Exergoeconomic assessment and parametric study of a gas turbine-modular helium reactor combined with two organic rankine cycles. Energy 2014, 65, 533–543. [Google Scholar] [CrossRef]
  26. Khaljani, M.; Saray, R.K.; Bahlouli, K. Comprehensive analysis of energy, exergy and exergo-economic of cogeneration of heat and power in a combined gas turbine and organic rankine cycle. Energy Convers. Manag. 2015, 97, 154–165. [Google Scholar] [CrossRef]
  27. Bejan, A.; Tsatsaronis, G.; Moran, M. Thermal Design and Optimization; John Wiley and Sons Inc.: New York, NY, USA, 1996. [Google Scholar]
  28. Ghaebi, H.; Amidpour, M.; Karimkashi, S.; Rezayan, O. Energy, exergy and thermoeconomic analysis of a combined cooling, heating and power (CCHP) system with gas turbine prime mover. Int. J. Energy Res. 2011, 35, 697–709. [Google Scholar] [CrossRef]
  29. Misra, R.D.; Sahoo, P.K.; Gupta, A. Thermoeconomic evaluation and optimization of an aqua-ammonia vapour-absorption refrigeration system. Int. J. Refrig. 2006, 29, 47–59. [Google Scholar] [CrossRef]
  30. Okoro, O.V.; Sun, Z. Desulphurisation of Biogas: A Systematic Qualitative and Economic−Based Quantitative Review of Alternative Strategies. ChemEngineering 2019, 3, 76. [Google Scholar] [CrossRef] [Green Version]
  31. Ramy, H.M.; Naef, A.A.; Syed, M.Z. Enhancing the thermal and economic performance of supercritical CO2 plant by waste heat recovery using an ejector refrigeration cycle. Energy Convers. Manag. 2020, 224, 113340. [Google Scholar] [CrossRef]
  32. Duan, C.J.; Wang, J.; Yang, X.Y. Features of supercritical carbon dioxide Brayton cycle coupled with reactor. At. Energy Sci. Technol. 2010, 44, 1341–1348. [Google Scholar]
  33. Zhu, Y. Dynamic Analysis of a Closed Brayton Cycle Using Supercritical Carbon Dioxide. J. Phys. Conf. Ser. 2021, 2030, 012010. [Google Scholar] [CrossRef]
  34. Holland, J.H. An introductory analysis with applications to biology, control and artificial intelligence. In Adaptation in Nature and Artificial Systems; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar] [CrossRef]
  35. Yari, M.; Sirousazar, M. A novel recompression S-CO2 Brayton cycle with pre-cooler exergy utilization. J. Power Energy 2010, 224, 931. [Google Scholar] [CrossRef]
  36. Lemmon, E.W.; Huber, M.L.; McLinden, M.O. NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 9.1; NIST NSRDS; NIST: Gaithersburg, MD, USA, 2010.
  37. Akbari Ata, D.; Mahmoudi Seyed, M.S. Thermoeconomic analysis & optimization of the combined supercritical CO2 (carbon dioxide) recompression Brayton/organic Rankine cycle. Energy 2014, 78, 501–512. [Google Scholar] [CrossRef]
  38. Zare, V.; Mahmoudi Seyed, M.S.; Yari, M. An exergoeconomic investigation of waste heat recovery from the gas turbine-modular helium reactor (GT-MHR) employing an ammoniaewater power/cooling cycle. Energy 2013, 61, 397–409. [Google Scholar] [CrossRef]
Figure 1. Diagram for the recompression Brayton cycle.
Figure 1. Diagram for the recompression Brayton cycle.
Energies 15 08184 g001
Figure 2. NEP subsystems and conceptual design.
Figure 2. NEP subsystems and conceptual design.
Energies 15 08184 g002
Figure 3. Structure diagram for the proposed S-(N2O-He) cycle.
Figure 3. Structure diagram for the proposed S-(N2O-He) cycle.
Energies 15 08184 g003
Figure 4. T-s diagram representation for the proposed S-(N2O-He) cycle [13].
Figure 4. T-s diagram representation for the proposed S-(N2O-He) cycle [13].
Energies 15 08184 g004
Figure 5. Optimization flowchart of the proposed S-(N2O-He) system.
Figure 5. Optimization flowchart of the proposed S-(N2O-He) system.
Energies 15 08184 g005
Figure 6. Influence of x on the proposed S-(N2O-He) cycle.
Figure 6. Influence of x on the proposed S-(N2O-He) cycle.
Energies 15 08184 g006aEnergies 15 08184 g006b
Figure 7. Influence of π on the proposed S-(N2O-He) cycle.
Figure 7. Influence of π on the proposed S-(N2O-He) cycle.
Energies 15 08184 g007aEnergies 15 08184 g007b
Figure 8. Influence of T6 on the proposed S-(N2O-He) cycle.
Figure 8. Influence of T6 on the proposed S-(N2O-He) cycle.
Energies 15 08184 g008aEnergies 15 08184 g008b
Figure 9. Influence of T2 on the proposed S-(N2O-He) cycle.
Figure 9. Influence of T2 on the proposed S-(N2O-He) cycle.
Energies 15 08184 g009aEnergies 15 08184 g009bEnergies 15 08184 g009c
Figure 10. Influence of P6 on the proposed S-(N2O-He) cycle.
Figure 10. Influence of P6 on the proposed S-(N2O-He) cycle.
Energies 15 08184 g010aEnergies 15 08184 g010b
Figure 11. Pareto optimal solution for ηth and ηex.
Figure 11. Pareto optimal solution for ηth and ηex.
Energies 15 08184 g011
Figure 12. Three objectives Pareto optimal solution for ηth, ηex, and C . total.
Figure 12. Three objectives Pareto optimal solution for ηth, ηex, and C . total.
Energies 15 08184 g012aEnergies 15 08184 g012b
Figure 13. Three objectives Pareto optimal solution for ηth, ηex, and Atotal.
Figure 13. Three objectives Pareto optimal solution for ηth, ηex, and Atotal.
Energies 15 08184 g013aEnergies 15 08184 g013b
Figure 14. Three objectives Pareto optimal solution for ηth, ηex, and LCOE.
Figure 14. Three objectives Pareto optimal solution for ηth, ηex, and LCOE.
Energies 15 08184 g014aEnergies 15 08184 g014b
Table 1. Physical properties of the working medium at critical point [18].
Table 1. Physical properties of the working medium at critical point [18].
NameMolecular FormulaCritical
Temperature
(K)
Critical
Pressure
(MPa)
Critical
Density
(kg/m3)
Carbon DioxideCO2304.057.38467.6
EthaneC2H6305.354.87206.2
Nitrous OxideN2O309.557.24452.0
Sulfur HexafluorideSF6318.753.75724.3.
WaterH2O647.0522.13322.0
Table 2. Range of the decision operating parameters.
Table 2. Range of the decision operating parameters.
ParameterRange
Split ratio—x0.25–0.4
Pressure ratio—π2.5–3.5
MC inlet temperature—T6315–340 K
TB inlet temperature—T21000–1200 K
MC inlet pressure—P67.4–8.4 MPa
Table 4. Energy and exergy results for a basic parameters design.
Table 4. Energy and exergy results for a basic parameters design.
ItemsValue
ηth (%)48.14
ηex (%)64.04
W . net (MW)0.277
Table 5. Exergy and exergoeconomic results for a basic parameters design.
Table 5. Exergy and exergoeconomic results for a basic parameters design.
Components E . F , j ( MW ) E . P , j ( MW ) E . D , j ( MW ) cF,j
($/GJ)
cP,j
($/GJ)
C . D , j ( $ / h ) Z . j ( $ / h ) Z . j +   C . D , j +   C . L , j ( $ / h ) rj
(%)
Fj
(%)
Reactor1.9771.9010.0767.7828.1722.1302.8174.9475.0156.94
Turbine0.7210.6010.1208.1729.4883.5260.3083.83416.108.04
Main Compressor0.2250.1860.0389.48811.9081.3050.0691.37425.514.99
Recompression Compressor0.1750.1010.0749.48810.8082.5310.0292.56113.911.15
HTR0.2130.1840.0298.1728.7420.8500.0840.9336.988.98
LTR0.2610.2280.0348.1728.7420.9920.0511.0436.984.92
Cooler0.0820.0190.0638.17228.5091.8520.8532.705248.8631.54
System3.6533.2200.434--13.1864.21217.397-24.21
Table 6. The partial Pareto optimal solution for ηth and ηex.
Table 6. The partial Pareto optimal solution for ηth and ηex.
FactorT6 (K)T2 (K)P6 (MPa)πxηth (%)ηex (%)
Optimal solution 1337.751020.98.213.130.33133.0445.95
2338.991001.08.353.420.25132.9945.96
3339.991010.08.33.450.28132.0447.01
4337.991012.08.13.490.33232.0546.95
5338.651000.08.03.490.32232.8146.00
6336.891020.08.23.490.33333.0845.95
7339.991028.08.33.490.25531.6547.43
9339.991031.08.43.490.32532.8846.98
10338.941010.28.23.490.25832.6946.27
11337.901020.58.13.340.32232.7746.02
Table 7. Pareto optimal solution for ηth, ηex, and C .  total.
Table 7. Pareto optimal solution for ηth, ηex, and C .  total.
FactorT6 (K)T2 (K)P6 (MPa)πxηth (%)ηex (%)
C . total   ( $ / h )
Optimal solution 1315.01117.718.33.40.34937.7451.1116.01
2317.51086.148.43.50.34736.9750.2716.36
3315.01112.478.43.50.34537.6551.0916.01
4315.11095.638.33.40.34937.4751.2316.06
5315.11065.918.33.50.34936.7050.9516.17
6315.61091.518.23.40.34937.3150.9816.13
7319.11065.938.33.30.34636.4949.8916.69
8318.31077.118.23.30.32536.8550.2116.51
9322.51066.458.33.50.34836.2448.9617.10
10321.81016.968.33.40.34834.7948.8417.45
Table 8. Pareto optimal solution for ηth, ηex, and Atotal.
Table 8. Pareto optimal solution for ηth, ηex, and Atotal.
FactorT6 (K)T2 (K)P6 (MPa)πxηth (%)ηex (%)Atotal (m2)
Optimal solution 1337.01004.28.33.350.33333.0845.95164.83
2335.21020.38.43.450.30232.3946.34164.01
3337.51010.58.23.460.32332.8446.00164.60
4338.91078.68.23.420.33233.0445.95164.80
5339.11025.38.43.480.33133.0245.95164.78
Table 9. Pareto optimal solution for ηth, ηex, and LCOE.
Table 9. Pareto optimal solution for ηth, ηex, and LCOE.
FactorT6 (K)T2 (K)P6 (MPa)πxηth (%)ηex (%)LCOE ($/kWh)
Optimal solution 1315.041000.137.42.580.25848.7767.160.01732
2315.011000.167.42.550.25249.0867.410.01697
3315.031000.327.42.560.25549.2467.120.01792
4317.021002.897.42.520.26349.3067.560.01680
5318.031010.167.42.560.26548.9267.280.01719
6315.121000.598.42.510.25153.5172.920.01640
7315.011000.897.42.510.26849.4567.680.01654
83161000.467.42.520.25549.3667.620.01662
Table 10. Final results of the key optimization parameters.
Table 10. Final results of the key optimization parameters.
FactorT6 (K)T2 (K)P6 (MPa)πxηth (%)ηex (%) C . total   ( $ / h ) Atotal (m2)LCOE ($/kWh)
Optimal
solution
316.510107.42.580.2649.1467.2919.37165.550.0196
Table 11. Key operating parameters after cycle optimization design.
Table 11. Key operating parameters after cycle optimization design.
ParameterValue in Basic CaseValue in Optimization CaseChange from the Base Case (%)
T6 (K)315316.50.48
T2 (K)12001010−15.83
P6 (MPa)7.47.40
π3.02.58−14.00
x0.30.26−13.33
W . (MW)0.2270.364+60.35
ηth (%)48.1449.14+2.08
ηex (%)64.0467.29+5.07
C . total ($/h)22.5219.37−13.99
Atotal (m2)165.57165.55−0.01
LCOE ($/kWh)0.020660.0196−5.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miao, X.; Zhang, H.; Wang, Q.; Sun, W.; Xia, Y. Thermodynamic, Exergoeconomic and Multi-Objective Analyses of Supercritical N2O-He Recompression Brayton Cycle for a Nuclear Spacecraft Application. Energies 2022, 15, 8184. https://doi.org/10.3390/en15218184

AMA Style

Miao X, Zhang H, Wang Q, Sun W, Xia Y. Thermodynamic, Exergoeconomic and Multi-Objective Analyses of Supercritical N2O-He Recompression Brayton Cycle for a Nuclear Spacecraft Application. Energies. 2022; 15(21):8184. https://doi.org/10.3390/en15218184

Chicago/Turabian Style

Miao, Xinyu, Haochun Zhang, Qi Wang, Wenbo Sun, and Yan Xia. 2022. "Thermodynamic, Exergoeconomic and Multi-Objective Analyses of Supercritical N2O-He Recompression Brayton Cycle for a Nuclear Spacecraft Application" Energies 15, no. 21: 8184. https://doi.org/10.3390/en15218184

APA Style

Miao, X., Zhang, H., Wang, Q., Sun, W., & Xia, Y. (2022). Thermodynamic, Exergoeconomic and Multi-Objective Analyses of Supercritical N2O-He Recompression Brayton Cycle for a Nuclear Spacecraft Application. Energies, 15(21), 8184. https://doi.org/10.3390/en15218184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop