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:
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 (N
2O–He), was performed for the proposed recompression Brayton cycle [
13] to convert thermal energy into electricity. In the N
2O–He mixture’s components, the N
2O part has thermodynamic properties such as molecular weight, critical pressure, and temperature, causing almost comparable behaviors under operating conditions that may improve the N
2O system’s performance and can use the existing CO
2 configuration. Sarkar J. [
14] utilized supercritical N
2O (S-N
2O) 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-CO
2 system. Characteristics of several critical working fluids at the critical point are shown in
Table 1. Using these mixtures (CO
2-He [
15], CO
2-C
2H
6/SF
6 [
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-CO
2 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 (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%N
2O–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-(N
2O-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:
where subscript
i represents the value of the entry state point, and
e represents the value of the exit state point.
De is the exergy destruction rate in the component.
Stream exergy consists of four main parts, which are
ph (stream physical exergy),
ch (stream chemical exergy),
k (stream kinetic exergy), and
p (stream potential exergy). Regardless of the
k and
p of the working medium, total exergy of the working medium
total is modified as follows [
20]:
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]:
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].
The thermodynamic model for the HTR and the LTR are given by [
21]:
The efficiency expressions for the TB, MC, and RC can be calculated by [
14]:
The subscript “s” is the isentropic process.
The thermodynamic calculation relations of the TB, MC, and RC are as follows [
14]:
Specific heat input of the reactor (
R) is expressed as:
The net power of circulating,
Wnet, is defined as:
The system thermal efficiency can be calculated by:
The system exergy efficiency is calculated by:
where
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]:
In Equation (20), Q,j is the cost rate related to component input heat, and 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]:
In Equations (21)–(23), cj, cW, and cQ are the per unit of cost exergy, $/GJ.
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]:
The component
j’s annual levelized CI can be calculated by [
27,
29]:
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]:
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.
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 (
total) is given in
Table A4, which is the sum of these values, including the CI cost rate, the total exergy destruction cost rate (
D), and the OM cost rate.
3.3. Performance Evaluation Index of the Cycle
net, ηth, ηex, total, Atotal, and LOCE are the evaluation metric for improving the operating system performance.
Net power (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 (
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]:
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
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 N
2O and He can be calculated using the NIST REFPROP [
35] software. As the NIST REFPROP software shows, the temperature of N
2O is limited to 525 K, and to estimate the thermodynamic properties of N
2O, 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 Symbol | Select Value |
---|
T0 (K) | 298 |
P0 (MPa) | 0.101 |
P6 (MPa) | 7.4 [36] |
m (kg/s) | 1 |
Fuel cost ($/MWh) | 7.4 [31] |
ηTB | 0.88 [13] |
ηMC | 0.85 [9] |
ηRC | 0.85 [9] |
εHTR | 0.86 [9] |
εLTR | 0.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,
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
j+
D,j +
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
j+
D,j +
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,
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 (
total) first decreased and then increased slightly; the whole variation trend decreases gently. An increase in the split ratio leads to the reduction of
D,total. The minimum value of
total can be obtained at
x = 0.33. With the increase of
T6, the
total decreases, but the variation trend is not apparent. When
T6 = 315 K and
x > 0.37,
total increases gradually with
x increases. The main reason for the change in
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
total of the entire cycle, mainly due to an increase in
D,reactor. In addition, with the increase of
π, the
D,HTR decreases and
D,LTR increases, and the total exergy destruction cost rate of recuperators increases due to the reduction of recuperators’ heat transfer area.
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
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 (
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
total decrease as
T6 increases. However,
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 N
2O-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
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
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,
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 (
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-(N
2O-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 (
in) increase greater than net power increases, resulting in reduced exergy efficiency, as shown in
Figure 9b.
total is decreased, mainly caused by the decrease in
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
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
total.
D,Reactor and
D,MC are also reduced. and the two components with the more significant proportion. All decreases in
total and
D,total are a result of decreasing
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. net (net power) decreases more than 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 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),
total =
f3(
x),
Atotal =
f4(
x), and
LCOE =
f5(
x) are expressed as:
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 total
Figure 12 and
Table 7 present the three objectives optimization results of
ηth—
f1(
x),
ηex—
f2(
x), and
total—
f3(
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
total, and
ηex and
total are the optimization objectives, respectively.
Figure 12d considers the Pareto solution of
ηth,
ηex, and
total.
4.3.3. Three Objectives Optimization of ηth, ηex, and Atotal
The three objectives optimization between
ηth–
f1(
x),
ηex–
f2(
x), and
Atotal–
f4(
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
ηth–
f1(
x),
ηex–
f2(
x), and
LCOE–
f5(
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-(N
2O-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), 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, 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 j + D,j + 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 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.