Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application

Page created by Yvonne Wise
 
CONTINUE READING
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
Received: 2 May 2020      Revised: 14 July 2020     Accepted: 15 July 2020
DOI: 10.1002/er.5857

RESEARCH ARTICLE

Investigation of off-design characteristics of an improved
recompression supercritical carbon dioxide cycle
for concentrated solar power application

Yuegeng Ma1,2                     | Tatiana Morosuk2 |                           Ming Liu3 |                 Jiping Liu1

1
 MOE Key Laboratory of Thermal Fluid
Science and Engineering, Xi'an Jiaotong
                                                      Summary
University, Xi'an, PR China                           The off-design characteristics of an improved recompression supercritical carbon
2
 Institute for Energy Engineering,                    dioxide cycle integrated with a two-stage intercooled main compressor are inves-
Technische Universität Berlin, Berlin,
                                                      tigated with a focus on the concentrated solar power application. An off-design
Germany
3
 State Key Laboratory of Multiphase Flow
                                                      model is established for each crucial component of the cycle system of
in Power Engineering, Xi'an Jiaotong                  100-megawatt scale. Four cycle control schemes with different main compressor
University, Xi'an, PR China                           configurations or/and cycle maximum pressure modes are evaluated and com-
Correspondence
                                                      pared. A sensitivity analysis is performed on the parameters related to the cycle
Jiping Liu, MOE Key Laboratory of                     thermal input and ambient condition to predict the off-design characteristics
Thermal Fluid Science and Engineering,                due to the plant dispatch and ambient condition change in a solar power plant.
Xi'an Jiaotong University, Xianning West
Road. 28, Xi'an 710049, PR China.                     The off-design results regarding the cycle thermodynamic performance and
Email: liujp@xjtu.edu.cn                              operational issue prevention are presented. The effect of the design-point value
Yuegeng Ma, Institute for Energy
                                                      of the main compressor inlet temperature on the off-design characteristics is
Engineering, Technische Universität
Berlin, Berlin, Germany.                              evaluated with the comparison among the results at three design points. The
Email: 1977497362@qq.com                              results reveal that the compressor surge may occur to the main compressor with

[Correction added on 16 January 2021,
                                                      basic configuration as the main compressor inlet temperature decreases to a cer-
after first online publication: Yuegeng Ma            tain value beneath the corresponding design point. By contrast, the surge risk
was designated as corresponding author.]              can be prevented with the modified main compressor configuration by activat-
Funding information                                   ing the recirculation system and the cycle can thus operate normally in the
Innovation capability support program of              entire off-design range of main compressor inlet temperature. The off-design
Shaanxi, Grant/Award Number: 2018TD-
                                                      change in thermal input has overall limited effects on the cycle system control.
014; China Scholarship Council, Grant/
Award Number: 201806280078; National                  No operational compressor issues occur for the main compressor with either
Key Research and Development Program                  basic or modified configuration as the thermal input deviates from the design
of China, Grant/Award Number:
                                                      points and varies in the studied ranges. The cycle maximum pressure mode has
2016YFB0600105; Science and Technology
on Thermal Energy and Power Laboratory                slight effects on the cycle thermodynamic performance as the thermal input
Open Foundation of China, Grant/Award                 deviates from the design point. The flexible cycle maximum pressure mode has
Number: TPL2017AA001
                                                      slightly lower sensitivity to the thermal input variation in net output power due

Abbreviations: B-FP, basic configuration with fixed pressure; B-VP, basic configuration with variable pressure; CSP, concentrated solar power;
DNI, direct normal irradiance, W/m2; HTF, heat transfer media; HTR, high-temperature recuperator; ICMC, intercooled main compressor; LTR,
low-temperature recuperator; MC, main compressor; M-FP, modified configuration with fixed pressure; M-VP, modified configuration with variable
pressure; PHX, primary heat exchanger; TES, thermal energy storage.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided
the original work is properly cited.
© 2020 The Authors. International Journal of Energy Research published by John Wiley & Sons Ltd.

1818      wileyonlinelibrary.com/journal/er                                                                                        Int J Energy Res. 2021;45:1818–1835.
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
MA ET AL.                                                                                                                       1819

                                           to the counteraction of mass flowrate variation. The selection of the design-point
                                           value of the main compressor inlet temperature has significant effects on the
                                           off-design characteristics of the cycle. A low design-point value of the main com-
                                           pressor inlet temperature leads to less demanding control action for operational
                                           issue prevention whereas a high design-point main compressor inlet tempera-
                                           ture results in overall more stable thermodynamic performance under off-design
                                           conditions. Among other schemes, the fixed maximum pressure mode with the
                                           modified main compressor configuration is found to be the most satisfactory
                                           one due to the consequent superior efficiency, steady net output power, and free
                                           of hazardous operating issues despite the relatively demanding task of compres-
                                           sor surge prevention. The developed control scheme can be further improved by
                                           implementing parametric optimization during the off-design operation.

                                           KEYWORDS
                                           compressor surge prevention, concentrated solar power plant, control scheme,
                                           off-design modeling, supercritical carbon dioxide cycle

1 | INTRODUCTION                                                       cycle with an emphasis on the CSP application. Among dif-
                                                                       ferent cycle designs, the recompression with an intercooled
As a potential solution to the climate change and energy               main compressor (ICMC) was considered as an competi-
shortage, concentrated solar power (CSP) is regarded as an             tive option due to the outstanding cycle performance and
important renewable energy technologies in the future                  large temperature differential in the cycle hot end for the
energy market due to its capability of utility-scale electricity       integration of the thermal storage system.10,11 Padilla
production and dispatchable power supply with the inte-                et al12 carried out an exergetic comparison between several
gration of thermal energy storage.1,2 Heat transfer fluids,            S-CO2 cycle configurations with the integrated exergetic
such as molten salts, synthetic oil or steam are adopted as            model of the solar receiver and S-CO2 cycle system. The
the medium to transfer the solar energy into thermal                   authors indicated that the recompression cycle with an
energy to activate the power cycle system for power genera-            ICMC was the cycle configuration with the highest overall
tion or store the thermal energy to facilitate the continuous          exergetic efficiency. Ma et al13 conducted a superstructure-
operation when the desired solar power is unavailable.                 based optimization on the S-CO2 cycle system for CSP
Limited by the degradation temperatures of different heat              application with the integrated thermo-economic model of
transfer media, the power block can operate at a maximum               the entire CSP plant employing an S-CO2 cycle-based
temperature level of around 350 C-700 C.3 The commer-                power block. This work also pointed out that the adoption
cialized CSP plants employ a traditional subcritical steam             of the recompression cycle with an ICMC can minimize
cycle system for the power block, which requires a complex             the levelized cost of electricity for the CSP plant.
system configuration and large footprint and does not                      An S-CO2 cycle system is expected to frequently operate
exhibit desirable limited thermal efficiency. The upgrade              under off-design conditions in a CSP plant due to the varia-
of the power cycle system is an important approach to                  tions in ambient conditions and cycle thermal input. The
maximizing the system performance of the CSP power and                 studies regarding the cycle off-design characteristics and
achieving a competitive cost of the CSP power generation.              operational strategy development are therefore essential
Among many other candidates, supercritical carbon diox-                for the further deployment of S-CO2 cycle systems in
ide (S-CO2) cycle systems have been seen as a competitive              CSP plants. Dyreby14 developed an off-design model for a
option for the CSP plant due to the superior efficiency and            10 MW recompression cycle which allows for the design
compact configuration.4,5                                              analysis of the cycle system. The cycle performance is
    Numerous cycle designs of S-CO2 cycle were proposed                predicted with the model under various design and bound-
and evaluated for applications in different power generation           ary conditions. Celle15 improved Dyreby's heat exchanger
scenarios, such as nuclear power plant,6 coal-fired plant,7            model by considering the effects of variations in CO2 prop-
waste heat recovery8 and CSP plant.9 Many previous works               erties and investigate the impacts of yearly variation of the
were carried out on the parametric studies of the S-CO2                ambient condition on the thermodynamic performance of
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
1820                                                                                                                  MA ET AL.

the S-CO2 cycle. Wang et al16 performed a performance           temperature, to simulate the cycle performance under vari-
evaluation of a solar thermal power plant employing a           ous off-design conditions encountered in a real CSP plant.
direct air-cooled supercritical carbon dioxide Brayton cycle.   The cycle performance and the prevention of potential haz-
The focus has been put on the design of the air-cooled heat     ardous issues of the main compressor (MC) under off-
sink and the effects of it under off-design conditions. The     design conditions are highlighted in the off-design study.
results revealed that the directed air-cooled system can        Four operation modes featuring different MC configuration
accommodate the lower solar intensities without deteriora-      and/or cycle maximum pressure control are proposed and
tion in electricity. Son et al17 developed an one-dimension     compared regarding the cycle thermodynamic performance
mean-line off-design model for S-CO2 turbomachinery with        and operational issue control for the MC under off-design
the assistance of Deep Neural Network. With the devel-          conditions. Based on the comparison among the cycle per-
oped models, it was claimed that more accurate off-design       formance with different operation modes, the optimal off-
performance prediction could be achieved in even shorter        design control scheme is finally recommended.
time than that with the traditional correction model,
which would greatly facilitate the off-design performance
evaluation of the entire cycle system. Luu18 investigated       2 | MODEL ESTAB LISHMENT
the approach of flexible operation of S-CO2 cycle in a CSP
plant with the assistance of a fossil-fuel heat source.         Figure 1 is the schematic of the S-CO2 recompression cycle
Two operational modes regarding the turbine inlet temper-       with an ICMC and the corresponding T-S diagram. The
ature control were proposed and compared. According to          detailed on-design model of this cycle was reported in our
the comparison, flexible temperature mode was reported          previous work.11 In comparison to the basic recompression
to outperform constant temperature mode in terms of fossil      cycle, a higher cycle efficiency (ηcyc) and a larger tempera-
fuel saving. Wright et al19 conducted an experimental study     ture differential for thermal input can be achieved with the
on a small-scale Brayton cycle loop and compared the            introduction of an ICMC, which can lead to more power
results with their developed models. The experimental           output from the cycle and a lower capital cost of the ther-
results exhibited controllable and stable operation in the      mal energy storage (TES) system.13,24 The levelized cost of
critical region, and the measured results were claimed to       electricity of the entire CSP system can be reduced as a
show good coherence with their modeling data. Carsten20         result. Table 1 gives the values of cycle design parameters
developed control strategies for the dynamic operation of       to initialize the off-design calculation. As shown in Table 1,
a recompression S-CO2 cycle. Several control strategies,        a 100-megawatt S-CO2 power cycle system is considered in
namely, high and low-temperature control, turbine bypass        this work. To accommodate this large power production
control, and inventory control, were compared under differ-     scale, a split shaft configuration is applied for the turboma-
ent operating modes including part-load operation, loss-of-     chinery. A synchronous generator is tied to the turbine and
load, loss of heat sink, over-power, and start-up/shut-down.    variable-speed drive motors are used for the compressors
Yang et al21 carried out part-load performance analysis and     for the sake of system efficiency and control flexibility.25
comparisons among various S-CO2 cycle configurations,           The main compressor is a two-stage internally geared
and concluded that the modifications on the cycle configu-      compressor with an intercooler in between. A CO2-Salt
ration can lead to different effects in different scenarios.    counter-flow shell-and-tube heat exchanger is chosen for
     The cycle configurations involved in the previous          the primary heat exchanger (PHX). Counter-flow printed
off-design studies were mostly recompression S-CO2 cycle        circuit heat exchangers with flow channels that are 5 mm
and simple-recuperated cycle. The results regarding the         wide and 2.5 mm deep, which is deemed as a representa-
off-design characteristics of more promising S-CO2 cycle        tive design for using in S-CO2 cycles,14,26 is adopted for the
designs for the CSP application, for example, the recom-        high-temperature recuperator (HTR) and low-temperature
pression cycle with an ICMC, were rarely reported. Besides,     recuperator (LTR). The air-cooled heat sink is applied for
very few studies focused on the effects of the potential        both the precooler and intercooler. A buffer tank is inte-
hazardous operational issues and the associated control         grated at the inlet of the MC to impose active control on
schemes on off-design operation. Operating in the critical      the inlet pressure of the MC with inventory control.14 The
region, the variation in inlet temperature of the main com-     inlet pressure of the MC is assumed to remain the on-
pressor can lead to drastic changes in CO2 density, which       design value and the split ratio of the stream flowing
may cause surge or choke to the compressor.14,22,23             through the recompressor is also assumed to remain the
     This article develops an off-design model for the          on-design value with the active valve control.
recompression cycle with an ICMC. Sensitivity analyses are          An objected-oriented approach is employed for the
conducted on the parameters related to different boundary       off-design model development of each main component
conditions, namely, the thermal input and the ambient           in the cycle in MATLAB 2018a. The calculation of the
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
MA ET AL.                                                                                                                       1821

F I G U R E 1 Diagrams of the S-CO2 recompression cycle with an ICMC. (A) Schematic diagram, (B) T-S diagram [Colour figure can be
viewed at wileyonlinelibrary.com]

required thermal properties of the CO2 is achieved by               calculation.29,30 Applying the Stodola's ellipse method, the
calling REFPROP database from MATLAB.27 The inlet                   following relationship between the on-design (ϕd) and off-
pressure of the MC is assumed to remain the design value            design (ϕod) mass flow coefficients is obtained:
using the inventory control under off-design conditions14
                                                                                                  rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
                                                                                                                    
and the split ratio of the stream flowing through the rec-                                                              2
                                                                                                1 − PPout,od
ompressor is also assumed to remain the on-design value                                  ϕod               in,od
                                                                                             = rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
                                                                                                                2             ð1Þ
with the active valve control. The model for each compo-                                 ϕd
nent is presented as follows.                                                                   1 − PPout,d in,d

2.1 | Turbine                                                       where ϕ is defined as
                                                                                                      pffiffiffiffiffiffiffi
A multi-axial flow type turbine is selected for the cycle for                                           T in
                                                                                               _ in 
                                                                                             ϕ=m                                ð2Þ
the sake of high efficiency and steady flow considering the                                            pin
100-megawatt power output.28 The off-design turbine inlet
pressure is calculated according to the Stodola's ellipse
method assuming the sliding mode and a fixed nozzle                     pin,od is then obtained as follow by substituting Equa-
area are also assumed for the turbine for the off-design            tions (1) and (2) as
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
1822                                                                                                                                         MA ET AL.

                                                                                                           T A B L E 1 On-design parameters
  Parameter                                                 Symbol            Unit         Value
                                                                                                           for the initialization of the off-design
  Cycle input heat                                          Q_ th             MW           222.2           modeling of the S-CO2 recompression
                                                                              
  Molten salt inlet temperature to the PHX                  TTES,hot              C        565             cycle with an ICMC11
                                                                              
  Turbine inlet temperature                                 TT,in                 C        550
  Main Compressor outlet pressure                           pMC,out           MPa           25
                                                                              
  Terminal temperature difference at                        ΔTC,cold-end          C         15
    the cooler cold ends
  Isentropic efficiency of turbine                          ηT                %             93
  Isentropic efficiency of compressor                       ηC                %             89
  Pressure drop in the PHX                                  ΔpPHX             kPa           50
  Pressure drop in the hot side of the HTR                  ΔpHTR,HT          kPa           60
  Pressure drop in the cold side of the HTR                 ΔpHTR,LT          kPa           30
  Pressure drop in the hot side of the LTR                  ΔpLTR,HT          kPa           20
  Pressure drop in the cold side of the LTR                 ΔpLTR,LT          kPa           40
  Pressure drop in the cooler                               ΔpC               kPa           20
  Minimal pinch point temperature difference                ΔTPinch           K                5

           qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi p2 −p2out,d    where Uin is the tip velocity at the inlet as defined in
pin,od =    m_ 2in,od  T in,od  Y d + p2out,od with Y d = in,d2
                                                                pin,d  ϕ2d   Equation (7) and Din is the inlet diameter of the compres-
                                                                              sor. hin and hout stand for the inlet and outlet enthalpies
                                                                       ð3Þ
                                                                              of CO2, respectively.
     The off-design isentropic efficiency of turbine (ηT,od)
is finally obtained as                                                                                             Din
                                                                                                          U in =       N                         ð7Þ
                                                                                                                    2
                               "                        0:1 #
                                      _ in,od  ρin,od
                                      m
            ηT,od = ηT,d  sin 0:5π                                   ð4Þ        Then, the nondimensionalization on the shaft speed
                                      m_ in,d  ρin,d
                                                                              is done for the flow coefficient, ideal head coefficient,
                                                                              and the compressor isentropic efficiency, by introducing
where ρin,d and ρin,od are the density of CO2 under on-                       the corresponding item relating to the shaft speed to the
design and off-design, respectively.                                          expressions of these parameters as displayed in Equa-
                                                                              tions (8) through (10). It should be noted that the modi-
                                                                              fied flow coefficient should be within the range of
2.2 | Compressor                                                              0.02-0.05 to avoid operational issues.14 The lower bound
                                                                              and upper bound are deemed as the threshold to trigger
The off-design performance characterization of the com-                       surge and choke. Nevertheless, a more accurate perfor-
pressors in the cycle is achieved by regressing the experi-                   mance map is necessary for the purpose of capturing the
mental data of a prototype compressor at the Sandia                           abnormal operating conditions accurately. Finally, the
National Laboratory with further nondimensionalization                        modified ideal head coefficient (ψ C ) and modified com-
on the shaft speed (N).14 The dimensionless flow coeffi-                      pressor isentropic efficiency (ηC ) are characterized as the
cient (ϕC) and ideal head coefficient (ψ C) of the compres-                   function of modified flow coefficient ϕC with polynomial
sor are first defined as shown in Equations (5) and (6).                      regression as presented in Equations (11) and (12)
                                                                                                                   
                                          _
                                          m                                                             _
                                                                                                        m          N
                           ϕC =                                        ð5Þ               ϕC   =                       ½0:02,0:05               ð8Þ
                                   ρin  U in  D2in                                             ρin  U in  D2in N d

                                                                                                                     
                                                                                                        ðhout −hin Þ N ð20ϕC Þ
                                                                                                                            3
                                 ðhout −hin Þ
                            ψC =                                       ð6Þ                     ψ C   =                                           ð9Þ
                                    U 2in                                                                  U 2in     Nd
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
MA ET AL.                                                                                                                                       1823

                                                  
                                                  N ð20ϕC Þ
                                                         5                  thermal property variations as shown in Equations (15)
                                    ηC,D,0
                   ηC   = ηC                                         ð10Þ   and (16).15
                                     ηC,D         Nd

                                                                                UAod   α − 1 + αh,d
                                                                                                −1

   ηC   =   −0:7069 + 168:6ϕC −8089ϕC 2            + 182725ϕC 3                  = −c,d                  _ 0:8  cnp  k ð1 − nÞ  μðn − 0:8Þ
                                                                                          1 + α − 1 with α = m
                                                                                UAd αc,od       h,od
             −1638000ϕC 4                                            ð11Þ                                                                     ð15Þ

     ψ C = −0:4049 + 54:6ϕC −2505ϕC 2 + 53224ϕC 3                                                          
                                                                                        Δpod     m_ d 7=4 μd 1=4 ρd − 1
                                                                                             =                                               ð16Þ
              − 498626ϕC 4                                           ð12Þ              Δpd      _ od
                                                                                                 m          μod    ρod

                                                                             where cp stands for the specific heat capacity, k stands
2.3 | Heat exchanger                                                         for the coefficient of heat conductivity. μ stands for the
                                                                             dynamic coefficient of viscosity. n equals to 0.4 and 0.3
The off-design model is developed for the heat exchangers                    for hot and cooling fluid, respectively.
in the S-CO2 cycle system including the HTR, LTR and
PHX. For the intercooler and precooler, the off-design
models are simplified by assuming a constant cold-end                        2.4 | Model validation
temperature difference of 15 K with the active control of
the cooling air flowrate. The off-design model is an one-                    The developed model is validated with the S-CO2 off-
dimension counter-flow heat exchanger model divided into                     design code developed by Dyreby.14 The results obtained
several axial nodes as commonly applied for the off-design                   from the developed model are compared to Dyreby's under
performance prediction of the heat exchangers in the                         both on-design and the off-design conditions as shown
system-level modeling in previous works.14,26 The devel-                     in Table 2. It is found the results of the developed model
oped off-design model considers both the heat transfer and                   show good coherence with the results obtained with
fluid dynamic characteristics of the heat exchanger. For                     Dyreby's code under both on-design and off-design condi-
the purpose of capturing the effect of fluid property varia-                 tions with a maximum difference of 0.14%. The accuracy
tion, the models are discretized into 20 sub-nodes. The                      of this model is therefore validated.
pressure loss (Δp) and the conductance of heat transfer
(UA) are calculated for each sub-node following the
method developed in Patnode's thesis31 with some custom-                     3 | CONS IDERATIONS FOR THE
ized modifications. The effects of fluid property variation                  OFF-DESIGN CHARACTERISTIC
on the off-design performance are considered in different                    STUDY
way for each type of heat exchanger. For the HTR and
PHX which have insignificant thermal property variations                     The S-CO2 cycle system employed in a CSP plant fre-
of the CO2 fluid, the effects of thermal property variations                 quently operates under off-design conditions due to the
of each sub-node in the heat transfer process are neglected.                 implementation of CSP dispatch and variation of ambient
The conductance of heat transfer (UAod) and pressure loss                    conditions. The CSP dispatch leads to the variation of
(Δpod) under off-design conditions are treated as the func-                  thermal input to the cycle in the hot end, which may
tions of the mass flowrate as follows.                                       change the mass flowrate (m      _ salt) and inlet temperature
                                                                             (Tsalt,in) of the molten salt at the inlet of the primary heat
                         UAod m_ − 0:8 + m − 0:8
                                         _ h,d                               exchanger. This may then affect the performance of tur-
                             = c,d
                                 −         − 0:8                      ð13Þ   bine and other cycle components. The variation of ambi-
                         UAd m _ c,od + m
                                   0:8
                                         _ h,od
                                                                             ent temperature causes the change of CO2 temperature at
                                                                           the outlets of precooler and intercooler, which in turn
                           Δpod     m_ d 7=4
                                =                                     ð14Þ   lead to the change in inlet temperature of the MC
                           Δpd      _ od
                                    m                                        (TMC,in). The variations in the MC inlet temperature tend
                                                                             to cause drastic changes in the density of CO2 at the inlet
    The CO2 fluid flowing through the LTR displays                           of the main compressor, especially when the pressure of
significant and nonlinear variation in fluid properties                      CO2 is close to the critical pressure. Surge/choke may
relating to the heat transfer and hydraulic characteristic                   occur to the MC as the volume flowrate drastically fluctu-
due to the vicinity to the critical point. Therefore,                        ates down/up with density if no measure is taken. There-
the UAod and Δpod are calculated considering the                             fore, regarding the control of the potential operational
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
1824                                                                                                                           MA ET AL.

TABLE 2         Validation of the developed model

                                          On-design condition                                  Off-design condition

                                                                              Error                                            Error
  Parameter                Unit           Present           Dyreby            (%)              Present             Dyreby      (%)
  Input heat               [MW]           10                10                Input            5.22                5.22        −0.04
  pmin                     [kPa]          9000              9000              Input            9000                9000        Input
  pmax                     [kPa]          25 000            25 000            Input            18 933.28           18 993.20   0.32
  Split ratio              [—]            0.3               0.3               Input            0.3                 0.3         Input
  effHTR                   [—]            0.95              0.95              Input            0.94                0.94        −0.05
  effLTR                   [—]            0.95              0.95              Input            0.90                0.90        −0.03
  Cycle efficiency         [%]            44.67             44.67             0.00             40.46               40.46       0.02
  Mass flowrate            [kg/s]         51.76             51.76             0.00             35.90               35.90       0.00
  t1 (t2,t3)               [ C]          41                41                Input            50                  50          Input
  t4                       [ C]          93.10             93.10             0.00             106.82              106.82      0.00
       0                    
  t5                       [ C]           210.47            210.88            0.20             239.85              240.19      0.14
       00                   
  t5                       [ C]           197.40            197.40            0.00             196.08              196.09      0.00
  t5                       [ C]          206.51            206.80            0.14             226.41              226.65      0.10
                            
  t6                       [ C]           395.19            395.19            0.00             431.77              431.72      −0.01
                            
  t7                       [ C]           550               550               Input            550                 550         Input
                            
  t8                       [ C]           431.12            431.12            0.00             461.20              461.20      0.00
  t9                       [ C]          217.67            218.03            0.17             240.91              241.24      0.14
                            
  t10                      [ C]           98.35             98.36             0.00             118.82              118.87      0.04

FIGURE 2          Two configurations of the main compressor. (A) basic configuration, (B) modified configuration

issue to the MC under varied conditions of inlet tempera-                additional parallel compressor for each compressor stage
ture, two configurations of the MC, namely, the basic con-               besides using the shaft speed control to deal with the
figuration and the modified configuration, are evaluated                 potential operational issues. The prevention of the surge
and compared in the following discussion. The simplified                 condition is achieved by the adjustments of both shaft
diagrams of the two MC configurations are presented in                   speed and the recirculation flows in two stages (m_ rec,1 and
Figure 2. The shaft speed is the only control variable for               m_ rec,2). The choke control is achieved by the adjustment
the MC with the basic configuration for operational issue                of shaft speed and the introduction of the second parallel
prevention, whereas the MC with the modified configura-                  compressor configured for both stages if necessary. These
tion is further configured with an anti-surge valve and an               modifications are expected to expand the applicable
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
MA ET AL.                                                                                                                       1825

F I G U R E 3 Four operation schemes under off-design conditions. (A) Basic main compressor configuration with fixed cycle maximum
pressure (B-FP), (B) modified main compressor configuration with fixed cycle maximum pressure (M-FP), (C) basic main compressor
configuration with variable cycle maximum pressure (B-VP); (D) modified main compressor configuration with variable cycle maximum
pressure (M-VP)
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
1826                                                                                                                                        MA ET AL.

conditions during off-design operation. Besides, two cycle                       4 | RESULTS A ND DISCUSSION
maximum pressure modes, namely, fixed pressure mode
and variable pressure mode, are evaluated and compared                           The cycle off-design performance is evaluated with a sensi-
with each other. The main compressor outlet pressure                             tivity analysis in this section. Three design points of the
(pMC,out) is always controlled at 25 MPa in the fixed pres-                      main compressor inlet temperature (TMC,in,d), i.e., 32 C,
sure mode, whereas in the variable pressure mode, pMC,out                        41 C and 50 C, are selected to investigate the effects of the
is not controlled unless it would cross the maximum bound                        TMC,in,d on the off-design performance. The design vari-
of 30 MPa. Four operation schemes are thereby presented                          ables are optimized for three design points prior to the off-
for the analysis, namely. (1) Basic main compressor config-                      design analysis. The values of the relevant parameters at
uration with Fixed cycle maximum Pressure mode (B-FP);                           the three corresponding on-design points are presented in
(2) modified main compressor configuration with Fixed                            Table 3. The thermodynamic states of the cycle at these
cycle maximum Pressure (M-FP); (3) basic main compres-                           three design points are presented in Table 4. The paramet-
sor configuration with Variable cycle maximum pressure                           ric analyses are performed for Tsalt,in and m_ salt to consider
(B-VP); (4) modified main compressor configuration with                          the off-design changes of the thermal input and for
Variable cycle maximum Pressure (M-VP). Figure 3A-D                              TMC,in to consider the effect of the ambient temperature.
presents the realization of these four schemes.                                  The information regarding the off-design sensitivity

                                                                                                          T A B L E 3 The values of the
                                                                   Design value of TMC,in [ C]
                                                                                                          relevant parameters at three optimized
  Parameter                                              Unit      32            41            50         on-design points
  Inlet pressure of the MC                               MPa       6.63          8.06          9.82
  Split ratio                                            —         0.62          0.66          0.70
  Inlet rotor diameter of the first MC stage             m         0.50          0.49          0.49
  Inlet rotor diameter of the second MC stage            m         0.41          0.42          0.43
  Shaft speed of the first MC stage                      rpm       2361          2562          2131
  Shaft speed of the second MC stage                     rpm       12 000        10 329        8786
  Inlet rotor diameter of the first RC stage             m         0.88          0.84          0.88
  Inlet rotor diameter of the second RC stage            m         0.35          0.41          0.48
  Shaft speed of the RC                                  rpm       11 977        11 208        10 215
  Cycle mass flow rate                                   kg/s      859.18        986.49        1164.39
  Design-point cycle efficiency                          %         48.00         46.44         44.81

                                                                                                          T A B L E 4 The thermodynamic
                    tMC,in,d = 32 C           tMC,in,d = 41 C             tMC,in,d = 50 C
                                                                                                          states of the cycle at three design points
  State point       t [ C]       p [kPa]      t [ C]          p [kPa]     t [ C]        p [kPa]
  1                 32            6627.80      41               8060.06     50             9817.88
  2                 45.25         7909.19      54.90            9804.57     59.95          11 475.54
  3                 32            7889.19      41               9784.57     50             11 455.54
  4                 64.97         25 000       75.42            25 000      84.43041       25 000
  5                 201.89        24 960       196.47           24 960      191.011        24 960
  50                202.05        24 960       198.96           24 960      195.672        24 960
      00
  5                 201.62        24 960       191.72           24 960      180.6101       24 960
  6                 355.47        24 930       374.73           24 930      396.002        24 930
  7                 550           24 880       550              24 880      550            24 880
  8                 392.26        6727.80      413.57           8160.06     435.78         9917.88
  9                 211.52        6667.80      207.35           8100.06     203.08         9857.88
  10                72.97         6647.80      82.62            8080.06     91.73          9837.88
Investigation of off-design characteristics of an improved recompression supercritical carbon dioxide cycle for concentrated solar power application
MA ET AL.                                                                                                                                    1827

analysis is shown in Table 5. It should be noted that the                  issues occur to the MC in the entire studied range of Tsalt,in.
upper limit of Tsalt,in selected here is selected only for the             The MC configuration is found to have no effects on the
sake of parametric analysis. The actual operating Tsalt,in                 system control. The cycle maximum pressure is slightly dif-
for the CSP using solar salt as storage media is rec-                      ferent in two cycle maximum pressure mode as can be
ommended to remain lower than 580 C considering the                       found by comparing the results of FP and VP cases. As
chemical stability of the solar salt.32                                    Tsalt,in rises, the cycle maximum pressure increases linearly
                                                                           under the VP mode and the shaft speed decreases linearly
                                                                           under the FP mode. The maximum variation of the cycle
4.1 | Analysis on thermal input                                            maximum pressure is within ±0.15 MPa and the variation
                                                                           of the shaft speed is within ±0.3% as Tsalt,in deviates the
Figure 4 presents the results of the variables relating to the             design point by ±30 C. The value of TMC,in,d does not show
cycle control as the function of Tsalt,in. No operational                  apparent effects on the cycle control.
                                                                               Figures 5 and 6 display the results of variables relating
T A B L E 5 The settings for the sensitivity analysis parameters           to the cycle thermodynamic performance. As shown in
under off-design operation                                                 Figures 5 and 6, the variation of Tsalt,in has relatively mild
                                    Design            Range for            effects on the cycle thermodynamic performance. The max-
                                    point             sensitivity          imum variations of ηen,cycle and W   _ net are both within
  Parameter           Unit          value             analysis             ±5.5% as Tsalt,in deviates the design point by ±30 C. The
  Tsalt,in            
                          C         565               535-595              choice of control schemes has limited effects on the ther-
                                                                           modynamic characteristics under off-design conditions
  _ salt/m
  m      _ salt,d     \             1                 0.6-1.2
                      
                                                                           of Tsalt,in as evidenced by the close variation trends of
  TMC,in                  C         32/41/50          32-50
                                                                           the four control scheme cases. The MC configuration is

F I G U R E 4 The variations of relative shaft speed (RN = N/Nd) and main compressor outlet pressure (pMC,out) with molten salt
temperature (Tsalt,in) under three design points of main compressor inlet temperature (TMC,in,d). (A) the case of TMC,in,d = 32 C, (B) the case
of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]
1828                                                                                                                                       MA ET AL.

FIGURE 5         The variations of net output power (W  _ net) and cycle energetic efficiency (ηen,cycle) with molten salt temperature (Tsalt,in)
under three design points of main compressor inlet temperature (TMC,in,d). (A) The case of TMC,in,d = 32 C, (B) the case of TMC,in,d = 41 C,
(C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]

F I G U R E 6 The variations of outlet temperature of the molten salt (Tsalt,out), turbine inlet temperature (TT,in) and mass flowrate of cycle
working fluid CO2 (m_ CO2 ) with molten salt temperature (Tsalt,in) under three design points of main compressor inlet temperature (TMC,in,d).
(A) The case of TMC,in,d = 32 C, (B) the case of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at
wileyonlinelibrary.com]
MA ET AL.                                                                                                                                     1829

found to have no effects here according to the compari-                    increases under the variable pressure mode and the
son between the results of different MC configurations.                    shaft speed decreases linearly under the FP mode, both
The cycle maximum pressure control has slight effects on                   at a decreasing rate. The cycle maximum pressure is
the results, specifically, the variation rate of m   _ CO2 and W _         decreased by around 0.35-0.55 MPa for the cycle with the
net with the change of Tsalt,in. The cycle with the VP mode                VP mode and the shaft speed is increased by 1% -1.5% as
has slightly lower variation rate of m      _ CO2 than that with           the m  _ salt decreases to 60% of the design value, while these
the FP mode due to the alleviation from the change of                      two variables only increases and decreases by less than
maximum pressure with Tsalt,in. Correspondingly, the W           _         0.04 MPa and 0.1%, respectively, as the m          _ salt rises to 120%
net of the cycle with the VP mode has slightly higher vari-                of the design value. The reduction in TMC,in,d slightly
ation rate than that with the FP mode. The cycle with                      increases the sensitivity of the control variable to the
lower TMC,in,d exhibits higher ηen,cycle and lower m       _ CO2 as        m_ salt variation.
expected. The effect of TMC,in,d appears to be independent                      Figures 8 and 9 display the sensitivity analysis results
of the variation of Tsalt,in since the studied variables have              on the variables relating to the cycle thermodynamic per-
similar variation tendencies and rates as the Tsalt,in devi-               formance as the functions of m                 _ salt,d. As shown in
                                                                                                                   _ salt/m
ates from the design point.                                                Figures 8 and 9, the thermodynamic performance is more
     Figure 7 presents the sensitivity analysis results on                 sensitive to the decrease in m      _ salt than to its increase. The
the variables relating to the cycle control as the function                ηen,cycle and W   _ net decrease by 15%-22% as m        _ salt decreases
of m       _ salt,d. No operational compressor issues occur as
    _ salt/m                                                               to 60% of the design point. The control scheme has over-
m_ salt deviates from the design point. The main compres-                  all limited effects on the thermodynamic characteristics
sor configuration does not affect the system control under                 as m _ salt deviates. The MC configuration is found to have
the off-design conditions of m     _ salt. The cycle maximum               no effects here according to the comparison between the
pressure is slightly different in two cycle maximum pres-                  results of the cases with different MC configurations.
sure modes according to the comparison between the                         The cycle maximum pressure mode has slight effects
results of FP and VP cases. As m    _ salt changes from the 60%            on the variation rate of m      _ CO2 and W   _ net with the change
to 120% of the design value, the cycle maximum pressure                    of m _ salt, especially as m_ salt decreases. The VP mode leads

F I G U R E 7 The variations of relative shaft speed (RN = N/Nd) and main compressor outlet pressure (pMC,out) with the molten salt mass
                           _ salt,d) under three design points of main compressor inlet temperature (TMC,in,d). (A) The case of TMC,in,d = 32 C,
                    _ salt/m
flow rate fraction (m
(B) the case of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]
1830                                                                                                                                         MA ET AL.

FIGURE 8             The variations of net output power (W _ net) and cycle energetic efficiency (ηen,cycle) with the molten salt mass flow rate
fraction (m       _ salt,d) under three design points of main compressor inlet temperature (TMC,in,d). (A) the case of TMC,in,d = 32 C, (B) the case
           _ salt/m
of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]

F I G U R E 9 The variations of outlet temperature of the molten salt (Tsalt,out), turbine inlet temperature (TT,in) and mass flowrate of cycle
working fluid CO2 (m_ CO2 ) with the molten salt mass flow rate fraction (m _ salt/m
                                                                                   _ salt,d) under three design points of main compressor inlet
temperature (TMC,in,d). (A) the case of TMC,in,d = 32 C, (B) the case of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be
viewed at wileyonlinelibrary.com]
MA ET AL.                                                                                                                                      1831

to slightly higher W _ net than the FP mode does under the                 to significantly different results in the cycle maximum
low m _ salt conditions. The reason for this is also due to                pressure and shaft speed. When TMC,in increases to 50 C,
the variation of the m  _ CO2 as explained above in the Tsalt,in           the cycle maximum pressure is reduced to around 12 MPa
cases. The cycle with a lower TMC,in,d exhibits higher                     in the VP mode and the shaft speed is increased by around
ηen,cycle and lower m _ CO2 as expected. It is also found that             50% in the FP mode to keep the cycle maximum pressure
the decrease of TMC,in,d leads to more drastic changes in                  at 25 MPa. In the cases with a TMC,in,d of 41 or 50 C, com-
the cycle variables with m   _ salt, especially as the m
                                                       _ salt has a        pressor surge occurs to the MC with basic configuration
significant reduction relative to the design value.                        when the TMC,in reduces to be lower than a certain value
                                                                           as shown in Figure 10B,C. For the cycle with the modified
                                                                           main compressor configuration, the recirculation system is
4.2 | Analysis on ambient temperature                                      activated when the compressor surge is approaching. The
                                                                           decrease in TMC,in entails more recirculating flow to pre-
Figure 10 presents the results of the variables relating to                vent the potential surge condition. For the cycle with the
the cycle control as the function of TMC,in. In comparison                 FP mode, the shaft speed decreases as TMC,in decreases.
to the thermal input variation, the variation in ambient                   For the cycle with the VP mode, the cycle maximum pres-
temperature has more significant effects on the cycle con-                 sure increases as TMC,in decreases until it reaches the upper
trol. The cycle at different design points of the TMC,in,d                 limit of the maximum pressure, the further decrease in
exhibit different off-design characteristics. The results of               TMC,in leads to the decrease in shaft speed. Besides, the
different TMC,in,d are therefore discussed separately. No                  cycle with the VP mode is less demanding on the preven-
operational compressor issues occur as TMC,in deviates                     tion of compressor surge as the cycle with the VP mode
from the design point for the cycle with a TMC,in,d of 32 C,              can operate in a larger range of TMC,in without the use of
and the MC configuration does not affect the system con-                   recirculation system. The design under a high TMC,in,d
trol in this case. The two maximum pressure modes lead                     entails more demanding control actions for compressor

F I G U R E 1 0 The variations of relative shaft speed (RN = N/Nd) and main compressor outlet pressure (pMC,out) with the main
compressor inlet temperature (TMC,in) under three design points of main compressor inlet temperature (TMC,in,d). (A) The case of
TMC,in,d = 32 C, (B) the case of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]
1832                                                                                                                                       MA ET AL.

FIGURE 11           The variations of net output power (W  _ net) and cycle energetic efficiency (ηen,cycle) with the main compressor inlet
temperature (TMC,in) under three design points of main compressor inlet temperature (TMC,in,d). (A) the case of TMC,in,d = 32 C, (B) the case
of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be viewed at wileyonlinelibrary.com]

F I G U R E 1 2 The variations of outlet temperature of the molten salt (Tsalt,out), turbine inlet temperature (TT,in) and mass flowrate of
cycle working fluid CO2 (m_ CO2 ) with the main compressor inlet temperature (TMC,in) under three design points of main compressor inlet
temperature (TMC,in,d). (A) The case of TMC,in,d = 32 C, (B) the case of TMC,in,d = 41 C, (C) the case of TMC,in,d = 50 C [Colour figure can be
viewed at wileyonlinelibrary.com]
MA ET AL.                                                                                                                            1833

surge prevention as can be observed by comparing the                 compressor issues occur when the thermal input devi-
results for the cycle with different TMC,in,d.                       ates from the on-design conditions in the studied ranges
    Figures 11 and 12 display the results of the variables           regardless of the main compressor configuration. The
relating to the cycle thermodynamic performance as the               cycle maximum pressure mode has slight effects on the
functions of TMC,in. As shown in Figures 11 and 12, the              cycle thermodynamic performance under off-design con-
effects of TMC,in on the cycle thermodynamic perfor-                 dition. The cycle in flexible maximum pressure mode
mance are significantly different for the cycle with differ-         is less sensitive to the thermal input variation in W           _ net
ent TMC,in,d. The W  _ net reduces as TMC,in increases for all       due to the counteract of m       _ CO2 variation. The maximum
three TMC,in,d cases, with more significant reduction                variations of ηen,cycle and W     _ net are both within ±5.5% as
observed with lower TMC,in,d. The cycle maximum pres-                Tsalt,in deviates the design point by ±30 C. The ηen,cycle
sure mode has significant effects on the thermodynamic               and W   _ net decrease by 15%-22% as m      _ salt decreases to 60%
characteristics as TMC,in deviates from the design value.            of the design point and increase by 1.0%-1.2% and 1.4%-
The FP mode can lead to much more stable ηen,cycle                   1.9% as m   _ salt increases to 120% of the design point.
and W _ net. A main reason for this is that the m
                                                _ CO2 remains      • The selection of TMC,in,d significantly affect the off-
constant for the cycle with a fixed maximum pressure. It             design characteristics of the cycle under varied condi-
is found that the cycle performance exhibits lower sensi-            tions of TMC,in. A low TMC,in,d leads to less demanding
tivity to the TMC,in with a higher TMC,in,d according to the         control task regarding the operational issue prevention
comparison among the results of the cycle with different             as no surge/choke is reported at the design point of
TMC,in,d.                                                            TMC,in,d = 32 C. The cycle with a high TMC,in,d is likely
                                                                     to have stable thermodynamic performance under off-
                                                                     design conditions. The maximum reductions in ηen,cycle
5 | C ON C L U S I ON                                                and W    _ net are 12.1% and 17.2% at the design point of
                                                                     TMC,in,d = 50 C and 25.6% and 84.8% at the design
This article develops an off-design model for a                      point of TMC,in,d = 32 C.
recompression S-CO2 cycle with an ICMC with an empha-              • Despite the relatively demanding control task for com-
sis on CSP application. The off-design characteristics with          pressor surge prevention, M-FP appears to be the most
respect to the thermodynamic performance and opera-                  satisfactory control scheme for the consequent steady
tional issue prevention are highlighted. Four control                ηen,cycle and W     _ net as well as the effective prevention of
schemes with different cycle maximum pressure mode                   operational issue of the main compressor. However,
or/and main compressor configuration are evaluated and               the cycle performance may be further improved when
compared. The effects of the off-design changes in cycle             the real-time parametric optimization is applied for the
thermal input and ambient temperature are investigated               cycle under off-design conditions.
                                              _ salt and TMC,in.
through the sensitivity analyses on Tsalt,in, m
Three different TMC,in,d are selected to investigate the           ACKNOWLEDGEMENTS
effects of the choice of TMC,in,d on the cycle off-design          This work is supported by the National Key Research and
performance. The following conclusions are drawn:                  Development Program of China (No. 2016YFB0600105),
                                                                   the Science and Technology on Thermal Energy and
• The compressor surge may occur to the main compres-              Power Laboratory Open Foundation of China (No.
  sor with basic configuration when the TMC,in attains a           TPL2017AA001), and Innovation capability support pro-
  certain value lower than the corresponding TMC,in,d. By          gram of Shaanxi (No. 2018TD-014). The first author
  contrast, the surge risk can be prevented with the mod-          Yuegeng Ma would also like to thank China Scholarship
  ified MC configuration by activating the recirculation           Council (No. 201806280078) for the financial support. Open
  system and can thus operate normally in the entire off-          access funding enabled and organized by Projekt DEAL.
  design range of TMC,in. The choke condition does not
  occur to the cycle under all the off-design conditions,          NOMENCLATURE
  the parallel compressor is not activated for either stage
  of the main compressor. This indicate that the risk of           Symbol
  choke appears controllable with shaft speed control. A           A          area of heat exchanger, m2
  customized performance map is still required to pre-             Cp         specific heat at constant pressure, kJ kg−1 K−1
  dict the abnormal conditions and exert the prevention            CSP        concentrated solar power
  control for the main compressor with better accuracy.            D          diameter, m
• The off-design change in thermal input has relatively            h          specific enthalpy, kJ kg−1
  mild effects on the cycle system control. No operational         H          height, m
1834                                                                                                                                  MA ET AL.

m_        mass flow rate, kg s−1                                             Proc of the ASME 2012 6th International Conference on Energy
N         shaft speed, —                                                     Sustainability; 2012;375-383.
p         pressure, kPa                                                 5.   Crespi F, Gavagnin G, Sánchez D, Martínez GS. Supercritical

Q_
                                                                             carbon dioxide cycles for power generation: a review. Appl
          heat transfer rate, kW
                                                                             Energy. 2017;195:152-183.
S-CO2     supercritical carbon dioxide                                  6.   Yoon HJ, Ahn Y, Lee JI, Addad Y. Potential advantages of cou-
SR        split ratio                                                        pling supercritical CO2 Brayton cycle to water cooled small and
T         temperature, K/ C                                                 medium size reactor. Nucl Eng Des. 2012;245:223-232.
U         tip velocity                                                  7.   Liu M, Zhang X, Ma Y, Yan J. Thermo-economic analyses on a
UA        heat transfer conductance, MW/K                                    new conceptual system of waste heat recovery integrated with
V         volume, m3                                                         an S-CO2 cycle for coal-fired power plants. Energy Convers
W _       work, kW                                                           Manage. 2018;161:243-253.
                                                                        8.   Wang X, Dai Y. An exergoeconomic assessment of waste heat
                                                                             recovery from a gas turbine-modular helium reactor using two
                                                                             transcritical CO2 cycles. Energy Convers Manage. 2016;126:561-572.
Greek symbols                                                           9.   Ma Y, Zhang X, Liu M, Yan J, Liu J. Proposal and assessment
k    coefficient of heat conductivity, W m−1 K−1                             of a novel supercritical CO2 Brayton cycle integrated with LiBr
Δp pressure drop, kPa                                                        absorption chiller for concentrated solar power applications.
ΔT temperature difference, K/ C                                             Energy. 2018;148:839-854.
η    isentropic efficiency                                             10.   Wang K, Li M-J, Guo J-Q, Li P, Liu Z-B. A systematic compari-
                                                                             son of different S-CO2 Brayton cycle layouts based on multi-
μ    dynamic coefficient of viscosity, kg m−1 s−1
                                                                             objective optimization for applications in solar power tower
ρ    density, kg/m3                                                          plants. Appl Energy. 2018;212:109-121.
ϕ    mass flow coefficient, —                                          11.   Ma Y, Morozyuk T, Liu M, Yan J, Liu J. Optimal integration of
ϕ*   mass flow coefficient, —                                                recompression supercritical CO2 Brayton cycle with main com-
ψ    ideal head coefficient, —                                               pression intercooling in solar power tower system based on
ψ* ideal head coefficient, —                                                 exergoeconomic approach. Appl Energy. 2019;242:1134-1154.
                                                                       12.   Vasquez R, Chean Y, Too S, Benito R, Stein W. Exergetic analy-
                                                                             sis of supercritical CO2 Brayton cycles integrated with solar
                                                                             central receivers. Appl Energy. 2015;148:348-365.
Subscript
                                                                       13.   Ma Y, Morosuk T, Luo J, Liu M, Liu J. Superstructure design
1,2… state point                                                             and optimization on supercritical carbon dioxide cycle for
cyc   cycle                                                                  application in concentrated solar power plant. Energy Convers
d     design point                                                           Manage. 2020;206:112290.
in    inlet                                                            14.   Dyreby JJ. Modeling the Supercritical Carbon Dioxide Brayton
od    off-design                                                             Cycle with Recompression. Madison, Wisconsin: Univ Wisconsin-
out   outlet                                                                 Madison; 2014.
                                                                       15.   de la Calle A, Bayon A, Soo Too YC. Impact of ambient tem-
rec   recirculation
                                                                             perature on supercritical CO2 recompression Brayton cycle in
S1    stage 1
                                                                             arid locations: finding the optimal design conditions. Energy.
S2    stage 2                                                                2018;153:1016-1027.
salt  molten salt                                                      16.   Wang X, Li X, Li Q, Liu L, Liu C. Performance of a solar ther-
                                                                             mal power plant with direct air-cooled supercritical carbon
                                                                             dioxide Brayton cycle under off-design conditions. Appl Energy.
ORCID                                                                        2020;261:114359.
Yuegeng Ma       https://orcid.org/0000-0001-7452-3080                 17.   Son S, Jeong Y, Cho SK, Lee JI. Development of supercritical
                                                                             CO2 turbomachinery off-design model using 1D mean-line
                                                                             method and deep neural network. Appl Energy. 2020;263:114645.
R EF E RE N C E S                                                      18.   Luu M, Milani D, Mcnaughton R, Abbas A. Analysis for flexi-
 1. Mehos M, Turchi C, Vidal J, et al. Concentrating solar Power             ble operation of supercritical CO2 Brayton cycle integrated with
    Gen3 demonstration roadmap. NREL Tech Rep. 2017;1-140.                   solar thermal systems. Energy. 2017;124:752-771.
 2. Islam MT, Huda N, Abdullah AB, Saidur R. A comprehensive           19.   Wright SA, Radel RF, Vernon ME, Rochau GE, Pickard PS.
    review of state-of-the-art concentrating solar power (CSP) tech-         Operation and analysis of a supercritical CO2 Brayton cycle;
    nologies: current status and research trends. Renew Sustain              September 2010.
    Energy Rev. 2018;91:987-1018.                                      20.   Carstens NA. Control strategies for supercritical carbon dioxide
 3. Power CS. Technology Roadmap: CSP Int Energy Agency;                     power conversion systems. Massachusetts Institute of Technol-
    2010;5:1-52.                                                             ogy; 2007.
 4. Turchi CS, Ma Z, Neises T, Wagner M. Thermodynamic study of        21.   Yang J, Yang Z, Duan Y. Part-load performance analysis and
    advanced supercritical carbon dioxide power cycles for high per-         comparison of supercritical CO2 Brayton cycles. Energy Conver
    formance concentrating solar power systems. Paper presented at:          Manage. 2020;214:112832.
MA ET AL.                                                                                                                              1835

22. Kurz R, Marechale RK, Fowler EJ, Cave MJ. Operation of cen-          30. Cook DH. Modeling of off-design multisage turbine pressures
    trifugal compressors in choke conditions. Paper presented at:            by Stodola's ellipse. Paper presented at: Energy Incorporated
    Proc 40th Turbomach Symp; 2011.                                          PEPSE User's Group Meeting; 1983.
23. White RC, Kurz R, White RC, Engineer P. Surge avoidance for          31. Patnode AM. Simulation and Performance Evaluation of Para-
    compressor systems. Paper presented at: Proc 35th Turboma-               bolic Trough Solar Power Plants. Madison, Wisconsin: Univ
    chinery Symp; 2006.                                                      Wisconsin-Madison; 2006.
24. Ma Y, Liu M, Yan J, Liu J. Thermodynamic study of                    32. Romero M, González-Aguilar J. Next generation of liquid metal
    main compression intercooling effects on supercritical CO2               and other high-performance receiver designs for concentrating
    recompression Brayton cycle. Energy. 2017;140:746-756.                   solar thermal (CST) central tower systems. Advances in Concen-
25. Weiland N, Thimsen D. A practical look at assumptions and                trating Solar Thermal Research and Technology. Cambridge,
    constraints for steady state modeling of sCO2 Brayton Power              England: Woodhead Publishing; 2017:129-154.
    cycles. Paper presented at: 5th Int Symp—Supercrit CO2 power
    cycles; 2016.
26. Dostal V. A supercritical carbon dioxide cycle for next generation
    nuclear reactors. Massachusetts Institute of Technology; 2004.
27. Lemmon EW, Huber ML, Mclinden MO. NIST Reference Fluid
    Thermodynamic and Transport Properties—REFPROP; 2013.
28. Sienicki JJ, Moisseytsev A, Fuller RL, Wright SA, Pickard PS.           How to cite this article: Ma Y, Morosuk T,
    Scale dependencies of supercritical carbon dioxide Brayton              Liu M, Liu J. Investigation of off-design
    cycle technologies and the optimal size for a next-step super-
                                                                            characteristics of an improved recompression
    critical CO2 cycle demonstration. Paper presented at: 3rd Int
    Symp—Supercrit CO2 Power Cycles; 2011.
                                                                            supercritical carbon dioxide cycle for concentrated
29. Li H, Yang Y, Cheng Z, Sang Y, Dai Y. Study on off-design per-          solar power application. Int J Energy Res. 2021;45:
    formance of transcritical CO2 power cycle for the utilization of        1818–1835. https://doi.org/10.1002/er.5857
    geothermal energy. Geothermics. 2018;71:369-379.
You can also read