Power, Efficiency, Power Density and Ecological Function Optimization for an Irreversible Modified Closed Variable-Temperature Reservoir Regenerative Brayton Cycle with One Isothermal Heating Process

One or more isothermal heating process was introduced to modify single and regenerative Brayton cycles by some scholars, which effectively improved the thermal efficiency and significantly reduced the emissions. To analyze and optimize the performance of this type of Brayton cycle, a regenerative modified Brayton cycle with an isothermal heating process is established in this paper based on finite time thermodynamics. The isothermal pressure drop ratio is variable. The irreversibilities of the compressor, turbine and all heat exchangers are considered in the cycle, and the heat reservoirs are variable-temperature ones. The function expressions of four performance indexes; that is, dimensionless power output, thermal efficiency, dimensionless power density and dimensionless ecological function are obtained. With the dimensionless power density as the optimization objective, the heat conductance distributions among all heat exchangers and the thermal capacitance rate matching among the working fluid and heat reservoir are optimized. Based on the NSGA-II algorithm, the cycle’s double-, triple- and quadruple-objective optimization are conducted with the total pressure ratio and the heat conductance distributions among heat exchangers as design variables. The optimal value is chosen from the Pareto frontier by applying the LINMAP, TOPSIS and Shannon entropy methods. The results show that when the pressure ratio in the compressor is less than 12.0, it is beneficial to add the regenerator to improve the cycle performance; when the pressure ratio is greater than 12.0, adding the regenerator will reduce the cycle performance. For single-objective optimization, the four performance indexes could be maximized under the optimal pressure ratios, respectively. When the pressure ratio is greater than 9.2, the cycle is simplified to a closed irreversible simple modified Brayton cycle with one isothermal heating process and coupled to variable-temperature heat reservoirs. Therefore, when the regenerator is used, the range of pressure ratio is limited, and a suitable pressure ratio should be selected. The triple objective (dimensionless power output, dimensionless power density and dimensionless ecological function) optimization’ deviation index gained by LINMAP or TOPSIS method is the smallest. The optimization results gained in this paper could offer some new pointers for the regenerative Brayton cycles’ optimal designs.


Introduction
Gas turbine plant (Brayton cycle (BCY)) has the characteristics of high safety, high efficiency and high energy density, and is widely used in the fields of thermal power generation, aviation propulsion, pipeline transportation, distributed energy and ship power plant, etc. However, NOx in gas turbine emission pollutants causes serious environmental pollution. The concentration of NOx in pollutants is closely related to the temperature of the combustion. For this reason, Vecchiarelli et al. [1] introduced an isothermal heating process (IHP) into the BCY, which effectively improved the thermal efficiency (TEF) and significantly reduced the emissions. The regenerative BCYs with IHP were studied by using the classical thermodynamics in [2][3][4]. El-Maksoud [5] added two IHPs in a binary BCY, and the cycle performance was significantly improved.
The temperature difference between the working fluid (WF) and the heat reservoir is not considered in the classical thermodynamic analysis. Compared to the classical thermodynamics, the finite size device and the finite time process of the finite rate heat exchange between the system and the environment are considered in the finite time thermodynamics (FTT) [6] and see the books by Berry et al. [7], Dincer et al. [8], Perescu et al. [9], Kaushik et al. [10] and Sieniutycz [11], as well as the review papers by Chen et al. [12], Andresen [13], Roach et al. [14] and Feidt et al. [15]. Berry et al. [7] comprehensively summarized the previous nearly 20 years' achievements on the FTT and optimal control. Dincer et al. [8] studied all kinds of advanced multi-output thermodynamic cycles and found that compared with traditional thermodynamic cycles, the advanced cycles had higher sustainability and energy utilization efficiency. Perescu et al. [9] studied how to meet the electricity and hot water needed for the personal housing by designing suitable mirrors. By using the FTT, Kaushik et al. [10] studied Carnot, Rankine, Brayton, Stirling, Ericsson Cycles, etc. Sieniutycz [11] described the properties of complexity and complex thermo-economic systems. Chen et al. [12] discussed the performance optimization of the thermodynamic processes and devices under the constraint of finite time or finite size. Andresen [13] introduced the development of the FTT in the past ten years and focused on the methods that he considered to be the most promising. Roach et al. [14] discussed the application of the principles of the FTT and control thermodynamics to biological processes. Feidt et al. [15] suggested discussing the heat transfer entropy and the production of entropy to study the external or internal irreversibility of the converter. The FTT has been widely used in thermodynamic cycle researches, the power output [16][17][18], TEF [16][17][18], power density [19][20][21][22][23] and ecological function [24][25][26][27][28] are often targeted when applying FTT to analyze or optimize the cycle performances.
Kaushik et al. [29] applied FTT in the study of the regenerative BCY with an IHP and found that compared with the traditional BCY, the TEF of the new BCY was significantly improved. Based on this, Tyagi et al. [30][31][32][33][34][35] and Kumar et al. [36] analyzed and optimized the regenerative [30,32,36], intercooled regenerative [31,33] and complex [34,35] BCYs with IHPs. With dimensionless power output (DPO), TEF, dimensionless power density (DPD) [37] and dimensionless ecological function (DEF) [38] as the objectives, Wang et al. [37,38] analyzed an endoreversible BCY modified by isothermal heating and found that the power density and ecological function are the trade-offs between power output and TEF. Based on this, Tang et al. [39,40] further applied the NSGA-II algorithm to optimize an endoreversible cycle and carried out sensitivity analysis. Arora et al. [41,42] analyzed and optimized a regenerative BCY with an IHP. Based on [5], Qi et al. [43] analyzed a binary BCY with two IHPs, derived the expression of DPO and TEF, and found that the DPO can be maximized by choosing an optimal pressure ratio. Tang et al. [44] took the DPO as the objective and the heat conductance distributions (HCDs) as design variables to optimize a closed binary Brayton with IHPs.
FTT will be applied to establish a model of closed regenerative BCY with an IHP and variable-temperature heat reservoirs in this paper. The isothermal pressure drop ratio (IPDR) is variable. The DPO, TEF, DPD and DEF will be obtained, and the influences of the regenerator's effectiveness on the performance indexes will be analyzed. Taking the DPD as an example, the optimal HCDs among heat exchangers under maximum DPD will be gained. Based on the NSGA-II algorithm, the DPO, TEF, DPD and DEF will be maximized and the corresponding Pareto frontier will be obtained. optimal value will be chosen from the Pareto frontier by applying the LINMAP, TOPSIS and Shannon entropy methods.

Cycle Model
The schematic and T-s diagrams of a closed irreversible regenerative modified BCY with an IHP and coupled to variable-temperature heat reservoirs (VTHRs) are shown in Figures 1 and 2, respectively. A compressor (C), a regenerator, a regular combustion chamber (RCC), a converging combustion chamber (CCC), a turbine (T) and a precooler are the main components of the cycle. The cycle consists of an irreversible adiabatic compression process (process 1-2), two isobaric heat absorption processes (processes 2-3 and 3-4), an IHP (process [4][5], an irreversible adiabatic expansion process (process 5-6) and two isobaric exothermic processes (processes 6-7 and 7-1). The processes of 1-2 s and 5-6 s are the corresponding isentropic processes.
It is assumed that the WF is a compressible ideal gas. The temperature of the WF gradually increases as it flows through a heating pipeline at subsonic speeds; the temperature decreases as it flows through a pipeline with a decreasing cross-sectional area and without heating at subsonic speeds. Combined with the two points, the convergent pipeline can be precisely and appropriately adjusted to achieve the isothermal heating of the WF [1]. Thus, the WF is isothermally heated in the CCC, and the thermal energy absorbed by the WF is converted into the kinetic energy.     It is assumed that the WF is a compressible ideal gas. The temperature of the WF gradually increases as it flows through a heating pipeline at subsonic speeds; the temperature decreases as it flows through a pipeline with a decreasing cross-sectional area and without heating at subsonic speeds. Combined with the two points, the convergent pipeline can be precisely and appropriately adjusted to achieve the isothermal heating of the WF [1]. Thus, the WF is isothermally heated in the CCC, and the thermal energy absorbed by the WF is converted into the kinetic energy.
The thermal capacitance rates of the WF and outer fluids at the RCC, CCC and precooler are C w f , C H , C H1 and C L , respectively. The dimensionless maximum specific capacity is v 5 /v 1 . The heat conductances of the regenerator, RCC, CCC and precooler are U R , U H , U H1 and U L , the numbers of heat transfer units are N R , N H , N H1 and N L , and the effectivenesses of the heat exchangers' are E R , E H , E H1 and E L , respectively. The corresponding expressions are: where (3) and (5) are simplified as: Processes 1-2s and 5-6s are the isentropic processes, so there are: where m = (k − 1)/k and k is the ratio of specific heats. For the IHP of process 4-5, one has: .
where . m is the mass flow rate, h is the specific enthalpy, π t (= p 5 /p 4 ≤ 1) is the IPDR, and . Q 4−5 is heat transfer rate. The degree of isothermal heating can be measured by the value of π t . Simultaneously, π t satisfies the following constraints: The velocity of the WF must be subsonic, i.e., M 5 ≤ 1. Considering that the WF has an initial velocity, M 4 = 0.2 and M 5 = 1, thus (M 2 5 − M 2 4 ) max = 0.96 and (π t ) min = 0.5107 can be obtained. The irreversibilities in the compressor and turbine are considered. The efficiencies of them are: The heat transfer rate ( . Q 2−3 ) at the regenerator, heat-absorbing rates ( . Q 3−4 and . Q 4−5 ) at the RCC and CCC and heat-releasing rate ( . Q 5−1 ) at the precooler are, respectively: . . . .
The power output (W), TEF (η), power density (P), entropy production rate (s g ) and ecological function (E) are, respectively, calculated as: where T 0 is the ambient temperature. In combination with Equations (8)- (10) and Equations (13)- (18), the analytical expressions of the state points' temperatures are: Energies 2020, 13, 5133 6 of 23 By substituting Equations (24)- (34) into Equations (19)-(23), the DPO (W), DPD (P) and DEF (E) can be obtained, as follows: The major differences between this paper and the [40] are the following aspects: Firstly, the irreversibilities in the compressor and turbine are considered in this paper; this is also the difference between the endoreversible model and the irreversible one. Secondly, the regenerator is added in this paper, that is, this paper studies the performance of the regenerated cycle, while [40] studied the performance of a simple cycle. Thirdly, the ecological function is studied in this paper. Fourthly, the double, triple or quadruple objective optimization is considered in this paper, but only triple objective optimization is considered in the [40]. Therefore, this paper is significantly different from the [40]. Equations (35)-(38) have a certain universality. When C H , C H1 , C L , E H1 , η c and η t are different values, the cycle could be simplified to different special cycles. When η c = η t = 1, Equations (35)- (38) are simplified into W, η, P and E of a closed endoreversible regenerative modified BCY with an IHP and coupled to VTHRs: where When C H = C L → ∞ , Equations (35)- (38) are simplified into W, η, P and E of a closed irreversible regenerative modified BCY with an IHP and coupled to constant-temperature heat reservoirs (CTHRs): Energies 2020, 13, 5133 When C H = C L → ∞ and η c = η t = 1, Equations (35)- (38) are simplified into W, η, P and E of a closed endoreversible regenerative modified BCY with an IHP and coupled to CTHRs: where Energies 2020, 13, 5133 9 of 23 When E H1 = 0, Equations (35)- (38) are simplified into W, η, P and E of a closed irreversible regenerative BCY coupled to VTHRs. When E H1 = 0 and η c = η t = 1, Equations (35)- (38) are simplified into W, η, P and E of a closed endoreversible regenerative BCY coupled to VTHRs. When E H1 = 0 and C H = C L → ∞ , Equations (35)- (38) are simplified into W, η, P and E of a closed irreversible regenerative BCY coupled to CTHRs. When E H1 = 0, C H = C L → ∞ and η c = η t = 1, Equations (35)- (38) are simplified into W, η, P and E of a closed endoreversible regenerative BCY coupled to CTHRs.

Single Objective Analysis and Optimization
In numerical calculations, C H1 = 0.6 kW/K, The numerical solution of y is obtained by combining Equations (11), (17) and (29). The numerical solution of y is substituted into Equations (35)- (38) to obtain the numerical solutions of W, η, P and E.

Single Objective Analysis
The influences of the effectiveness (E R ) of regenerator on W, η, P and E, π t and v 5 /v 1 are analyzed under the given π. The relationships of W and η, P and E versus π with different E R are shown in Figures 3 and 4, respectively. From the two figures, W, η, P and E all rise firstly and then fall as π rises, and there are respective optimal values of π to make W, η, P and E reach their maximum values, respectively. When π < 12.0, W, η, P and E all increase as E R increases; when π > 12.0, W, η, P and E all decrease as E R increases. The maximums of W, η, P and E gradually move to the left as E R increases, and the influence on η among them is the greatest. When π = 12.0 (i.e., T 2 = T 6 ), increasing π will reduce T 3 , which allows the exhaust to carry away more heat.   / v v versus π with different R E . From Figure 5, R E has a great influence on t π . Given different R E , the variation rule between t π and π are different. R E has a great influence on 3 T : when π increases as R E increases, which reduces the degree of the isothermal heating; when 12.0 π > , t π decreases as R E increases, which aggrandizes the degree of the isothermal heating. 5 1 / v v decreases as π increases. When 12.0  / v v versus π with different R E . From Figure 5, R E has a great influence on t π . Given different R E , the variation rule between t π and π are different. R E has a great influence on 3 T : when π increases as R E increases, which reduces the degree of the isothermal heating; when 12.0 π > , t π decreases as R E increases, which aggrandizes the degree of the isothermal heating. 5 1 / v v decreases as π increases. When 12.0   Figure 5 shows the relationships of π t and v 5 /v 1 versus π with different E R . From Figure 5, E R has a great influence on π t . Given different E R , the variation rule between π t and π are different. E R has a great influence on T 3 : when π < 12.0, π t increases as E R increases, which reduces the degree of the isothermal heating; when π > 12.0, π t decreases as E R increases, which aggrandizes the degree of the isothermal heating. v 5 /v 1 decreases as π increases. When π < 12.0, v 5 /v 1 increases as E R increases; when π > 12.0, v 5 /v 1 decreases as E R increases.  / v v versus π with different R E . Figure 6 shows the relationships of W and η versus π with different c η and t η . Figure 7 shows the relationships of P and E versus π with different c η and t η . It can be seen from the figure that W , η , P and E increase as c η or t η increases. As π increase, c η or t η has a greater influence on W , η , P and E . When 1 c t η η = = , the cycle is endoreversible, that is, the irreversibilities in the compressor and turbine is not considered. Figure 8 shows the relationships of   Figure 6 shows the relationships of W and η versus π with different η c and η t . Figure 7 shows the relationships of P and E versus π with different η c and η t . It can be seen from the figure that W, η, P and E increase as η c or η t increases. As π increase, η c or η t has a greater influence on W, η, P and E. When η c = η t = 1, the cycle is endoreversible, that is, the irreversibilities in the compressor and turbine is not considered. Figure 8 shows the relationships of π t and v 5 /v 1 versus π with different η c and η t . It can be seen from the figure that π t and v 5 /v 1 decrease as η c or η t increases. It indicates that the degree of the cycle's isothermal heating correction is improved and the volume of the cycle device is reduced by reducing the irreversibility of the cycle. greater influence on W , η , P and E . When 1 c t η η = = , the cycle is endoreversible, that is, the irreversibilities in the compressor and turbine is not considered. Figure 8 shows the relationships of t π and 5 1 / v v versus π with different c η and t η . It can be seen from the figure that t π and 5 1 / v v decrease as c η or t η increases. It indicates that the degree of the cycle's isothermal heating correction is improved and the volume of the cycle device is reduced by reducing the irreversibility of the cycle.

Optimal Distributions of Heat Exchanger Inventory
In the actual design process, the heat exchanger inventory is limited or given, which indicates that the total heat conductance (U T ) is limited or given. According to Equations (35)- (38), W, η, P and E are related to U H , U H1 , U L and U R . When U T is given, the reasonable distribution of the heat conductances is conducive to improving the cycle performance. With fixed total heat conductance U T (= U H + U H1 + U L + U R ), the HCDs among the RCC, CCC, precooler and regenerator are: The HCDs satisfy the following relationships: With a given U T , the HCDs are optimized to maximize W, η, P and E. Figure 9 presents the optimization flowchart for single objective. The optimization calculations of different objectives are similar. Therefore, only the optimization results of P is taken as an example for discussions. The optimization results of W, η and E are not presented herein.
With a given T U , the HCDs are optimized to maximize W , η , P and E . Figure 9 presents the optimization flowchart for single objective. The optimization calculations of different objectives are similar. Therefore, only the optimization results of P is taken as an example for discussions. The optimization results of W , η and E are not presented herein.   10-12 present the relationships of the maximum power density (P opt ) and the corresponding power (W P opt ), TEF (η P opt ), ecological function (E P opt ), and other optimal parameters (π t ) P opt , (u H ) P opt , (u H1 ) P opt , (u L ) P opt and (u R ) P opt versus π. From Figure 10, W P opt and P opt increase first and then decrease as π increases. When π ≤ 9.2, η P opt first augments and then slightly decreases as π augments; When π > 9.2, η P opt first augments and then decreases as π augments.  first decreases and then augments as π augments, and the degree of the isothermal heating first augments and then decreases. In combination with the previous performance analysis, the regenerator is meaningless when π exceeds a certain value. In this case, R E should be 0, and the cycle is simplified to a closed irreversible modified simple BCY with an IHP and coupled to VTHRs. In short, the effect of the regenerator is better when π is relatively small; it is not recommended to add a regenerator to the cycle when π is relatively large.   first decreases and then augments as π augments, and the degree of the isothermal heating first augments and then decreases. In combination with the previous performance analysis, the regenerator is meaningless when π exceeds a certain value. In this case, R E should be 0, and the cycle is simplified to a closed irreversible modified simple BCY with an IHP and coupled to VTHRs. In short, the effect of the regenerator is better when π is relatively small; it is not recommended to add a regenerator to the cycle when π is relatively large.  The influences of the thermal capacitance rate matching among the WF and heat reservoir on the cycle performance are discussed. Figure 13 illustrates the relationships of opt P versus / H wf C C Figure 12. Relationships of (u H1 ) P opt , (u L ) P opt and (u R ) P opt versus π.
Energies 2020, 13, 5133 14 of 23 From Figure 11, when π ≤ 9.2, E P opt augments first and then slightly augments as π augments; when π > 9.2, E P opt augments first and then slightly decreases as π augments. (π t ) P opt first decreases and then augments as π augments, and the degree of the isothermal heating first augments and then decreases. (u H ) P opt decreases as π augments. From Figure 12, when π ≤ 9.2 and π increases, (u H1 ) P opt diminishes first and then increases slightly, (u L ) P opt increases rapidly, and (u R ) P opt drops rapidly; when π > 9.2 and π increases, (u H1 ) P opt decreases first and then increases, (u L ) P opt slowly increases, and (u R ) P opt is always 0.
In combination with the previous performance analysis, the regenerator is meaningless when π exceeds a certain value. In this case, E R should be 0, and the cycle is simplified to a closed irreversible modified simple BCY with an IHP and coupled to VTHRs. In short, the effect of the regenerator is better when π is relatively small; it is not recommended to add a regenerator to the cycle when π is relatively large.

Optimal Thermal Capacitance Rate Matching among the WF and Heat Reservoir
The influences of the thermal capacitance rate matching among the WF and heat reservoir on the cycle performance are discussed. Figure 13 illustrates the relationships of P opt versus C H /C w f and C H1 /C w f when π = 8. The figure shows that P opt first increases and then tends to remain constant as C H /C w f increases. Because when C H < C w f , C Hmax = C w f and C Hmin = C H ; when C H > C w f , C Hmax = C H and C Hmin = C w f ; P opt first increases and then tends to remain constant as C H1 /C w f increases. Because . Q 3−4 first increases and then tends to remain constant as C H1 /C w f increases.
To further increase the DPD, it is necessary to increase . Q 2−3 , that is, to increase u H or U H . But as . Q 2−3 increases, . Q 3−4 will also increase, which forms a balance between the two. Finally, the DPD tends to remain constant. By numerical calculations, W P opt , η P opt and E P opt also first increase and then tend to remain constants as C H /C w f or C H1 /C w f increases. Thus, there is an optimal C H /C w f or C H1 /C w f to make W P opt , η P opt , P opt and E P opt to reach the optimal values.    Figure 15 illustrates the relationship of opt P versus    Figure 14 illustrates the relationship of P opt versus C H /C w f with different C L /C H . Figure 15 illustrates the relationship of P opt versus C H1 /C w f with different C L /C H1 . As can be seen from the figures, when C H /C w f or C H1 /C w f is smaller, P opt increases as C L /C H or C L /C H1 increases; when C H /C w f or C H1 /C w f is larger, P opt is unaffected by C L /C H or C L /C H1 . By adjusting the value of C L /C H or C L /C H1 , the value of C H /C w f or C H1 /C w f could be changed when P opt reaches the optimal values. Under the condition that the DPD is kept constant, the requirement of the heat exchangers corresponding to the high-temperature reservoirs can be reduced by increasing C L /C H or C L /C H1 and reducing C H /C w f or C H1 /C w f .   Figure 16 shows the optimization flowchart for calculating the Pareto frontier using the NSGA-II algorithm [39][40][41][45][46][47][48][49][50][51][52][53][54]. With π and HCDs as the design variables and W , η , P and E as the optimization objectives, the cycle's double, triple or quadruple objective optimization is conducted.

Multi-Objective Optimization
The Pareto frontier contains many feasible solutions, which need to be selected by some methods. The LINMAP, TOPSIS and Shannon entropy methods are three common decision-making methods [40]. These three methods are used to get the optimal scheme from the Pareto frontier, and the results are compared. Finally, the merits and demerits of each decision method are judged by the deviation index (D). Positive ideal point (PIP) is the optimal value combination of four objective functions in the Pareto frontier while negative ideal point (NIP) is the worst value combination. The LINMAP method, TOPSIS method and D are all based on the PIP and NIP.   Figure 16 shows the optimization flowchart for calculating the Pareto frontier using the NSGA-II algorithm [39][40][41][45][46][47][48][49][50][51][52][53][54]. With π and HCDs as the design variables and W , η , P and E as the optimization objectives, the cycle's double, triple or quadruple objective optimization is conducted. The Pareto frontier contains many feasible solutions, which need to be selected by some methods. The LINMAP, TOPSIS and Shannon entropy methods are three common decision-making methods [40]. These three methods are used to get the optimal scheme from the Pareto frontier, and the results are compared. Finally, the merits and demerits of each decision method are judged by the deviation index (D). Positive ideal point (PIP) is the optimal value combination of four objective functions in the Pareto frontier while negative ideal point (NIP) is the worst value combination. The LINMAP method, TOPSIS method and D are all based on the PIP and NIP.  Figure 16 shows the optimization flowchart for calculating the Pareto frontier using the NSGA-II algorithm [39][40][41][45][46][47][48][49][50][51][52][53][54]. With π and HCDs as the design variables and W, η, P and E as the optimization objectives, the cycle's double, triple or quadruple objective optimization is conducted.

Multi-Objective Optimization
The Pareto frontier contains many feasible solutions, which need to be selected by some methods. The LINMAP, TOPSIS and Shannon entropy methods are three common decision-making methods [40]. These three methods are used to get the optimal scheme from the Pareto frontier, and the results are For the Shannon entropy method, it is assumed that the Pareto frontier contains n feasible solutions, and each feasible solution contains m objectives. F ij (1 ≤ i ≤ m, 1 ≤ j ≤ n) represents the value of the j-th optimization objective for the i-th feasible solution, and i opt is the best possible solution in Equation (63).
i opt ∈ min P ij · w j (63) where P ij represents the dimensionless value of F ij , and SE j and w j are the Shannon entropy and the weight of j-th optimization objective, respectively. For LINMAP method, i opt is the best possible solution in Equation (66).
where ED i+ represents the Euclidean distance between the i-th feasible solution and the PIP, and F positive represents the value of the PIP. For the TOPSIS method, i opt is the best possible solution in Equation (68).
where ED i− represents the distance between the i-th feasible solution and the NIP, and F negative represents the value of the NIP. The deviation index is calculated according to Equation (70): where G j represents the value of the j-th optimization objective corresponding to the decision method. Figure 17 shows the Pareto frontier with W, η, P and E as the optimization objectives. The line's color of the Pareto frontier represents the value of the DEF. Table 1 lists the comparison of double, triple and quadruple objective optimization results. For quadruple objective optimization, the results gained by LINMAP and TOPSIS methods are the same, and D is smaller than that gained by the Shannon entropy method. It indicates that the results gained by the LINMAP method are better than those gained by the Shannon entropy method. The u R gained by LINMAP or TOPSIS method is extremely small and tends to 0. The IPDRs of the three decision methods have little difference, thus the isothermal heating degrees of the three methods have little difference. W and P gained by LINMAP or TOPSIS method are larger than those gained by the Shannon entropy method, but η and E gained by LINMAP or TOPSIS method are smaller than those gained by Shannon entropy method. For triple objective optimization, the results gained by LINMAP and TOPSIS methods are the same. The triple objective (W, P and E) optimization' deviation index gained by LINMAP or TOPSIS method is the smallest. For double objective optimization, all of u R corresponding to the double objective optimization are extremely small and tend to 0. where j G represents the value of the j-th optimization objective corresponding to the decision method. Figure 17 shows the Pareto frontier with W , η , P and E as the optimization objectives. The line's color of the Pareto frontier represents the value of the DEF. Table 1

Conclusions
FTT is applied to establish a model of a closed irreversible regenerative BCY with an IHP and variable IPDR in this paper. The DPO, TEF, DPD and DEF are obtained, and the influences of E R on the performance indexes are analyzed. With P as the optimization objective and HCDs as the design variables, the single objective analysis and optimization of the cycle are conducted. Based on the NSGA-II algorithm, the DPO, TEF, DPD and DEF are maximized with π and HCDs as the design variables and the corresponding Pareto frontier is obtained. Three decision methods are applied to get the best scheme. The main conclusions manifest that: (1) For the single objective analysis, the DPO, TEF, DPD and DEF all augment first and then reduce as π augments. When π < 12.0, it is beneficial to add the regenerator to improve the cycle performance; when π > 12.0, adding the regenerator will reduce the cycle performance. For the single-objective optimization, when π > 9.2, the effectiveness at the regenerator should be 0, and the cycle will be simplified into a closed irreversible modified BCY with an IHP and coupled to VTHRs. Hence, the regenerator is added, and a suitable pressure ratio should be selected in the limited range. In the actual process, the pressure ratio should be strictly controlled in this new-type gas turbine, or the new-type gas turbine is suitable for the situations where the pressure ratio is relatively low. The DPO, TEF, DPD and DEF increase as η c or η t increases. π t and v 5 /v 1 decrease as η c or η t increases. It indicates that the cycle performance and the degree of the cycle's isothermal heating correction are improved and the volume of the cycle device is reduced by reducing the irreversibility of the cycle. (2) For the quadruple objective optimization, the results gained by LINMAP and TOPSIS methods are the same, and D is smaller than that gained by the Shannon entropy method. Therefore, the results gained by LINMAP and TOPSIS methods are better than those gained by the Shannon entropy method. For the double, triple or quadruple objective optimization, the triple objective (W, P and E) optimization' deviation index gained by LINMAP or TOPSIS method is the smallest. According to decision-makers' actual needs, they could choose the best scheme from the Pareto frontier, design the new-type gas turbine reasonably, and then determine the construction or the operation of the gas turbine.
FTT is a powerful theoretical tool and the NSGA-II algorithm is an effective mathematical tool for the performance analysis and optimization of closed irreversible regenerative BCY with an IHP and variable IPDR. The results obtained in this paper have some guiding significance for the design of a new-type gas turbine.