The technical analysis included building a high-fidelity first principles model of a Reference power plant to determine baseline performance (plant performance without an integrated ESS) over the operating load range of the plant. The model results were used to develop baseline performance curves: Input-Output (I-O), Heat Rate (HR), Incremental heat rate (IHR) and Incremental Cost (IC) curves. The model of the Reference plant was then modified by adding/integrating the ESSs selected for the analysis.
As stated earlier, the TES systems described in this paper include the low-pressure (LP) condensate storage, two-tank molten solar salt storge, and fixed bed energy storage where heat is stored in a solid medium (crushed rocks). The selected models and their performance are described in the sections below.
2.1. Low-Pressure (LP) Condensate Thermal Energy Storage
The charging of the LP condensate storage system integrated with the Reference power plant is presented schematically in
Figure 6. To decrease the plant power output, the storage tanks are filled with hot condensate taken from the outlet stream of the feedwater storage tank (FWT). To accomplish this, the steam extraction for the deaerator (D) has to be increased. The charging hot condensate displaces cold condensate from the storage tanks forcing it to merge with the main condensate flow downstream of the condensate pump. As a result, the LP condensate flow through the LP FWHs and steam extractions for the LP FWHs and deaerator increases. The increase in steam extractions decreases steam flow through the LP steam turbine and its corresponding power output. Consequently, the overall plant power output decreases. Similar system configurations involving storage of LP condensate were proposed in [
50,
51]. The red arrows in
Figure 6 indicate the hot condensate flow, while the blue arrows denote the cold and main condensate flows. During system discharge, the LP hot condensate stored in the condensate tanks is discharged into the main condensate line upstream of the deaerator (D). The discharged hot condensate is replaced by the cold condensate, reducing the main condensate flow through the LP FWHs and D, i.e., the direction of the red and blue arrows in
Figure 6 is reversed. This decrease in the condensate flow results in a decrease in steam extractions. As a result, the steam flow through the LP steam turbine increases, increasing power output of the LP turbine and the power plant.
The storage tanks were modeled as stratified storage, where a thermocline is established between the hot condensate occupying the upper portion of the tank and the cold condensate at the bottom of the tank. To reduce heat conduction between the hot and cold condensate, floating baffles were used to separate the two zones.
The effect of the condensate tank charging on the plant power output is presented in
Figure 7 where the plant power output is plotted as a function of the charging time over a range of condensate extraction flows. The results show that the power output decreases almost instantaneously when a portion of the hot condensate flow leaving the feedwater storage tank is diverted to the LP tanks and stays constant during the charging process. During discharging, the plant power output increases almost instantaneously. These results show that the LP condensate storage system provides fast response.
The effect of the LP condensate charging and discharging flows on the plant power output, presented in
Figure 8 for the minimum and full load operation, shows that the effect of the charging/discharging flow on plant power output is linear and increases with the condensate flow. The full load charging/discharging flows presented in
Figure 8 correspond to 25 to 40% of the Boiler Feed Pump (BFP) flow, and for the analyzed Reference plant result in a power output change in the ± 10 MW (± 1.5% of full load) to ± 15 MW (± 2.3%) range. In other words, integration of the LP condensate storage system would, at full load and condensate flow corresponding to 40% of BFP flow allow the Reference plant to follow load in the ± 15 MW range (2.3%) and improve its participation in frequency regulation or increase the maximum power output by 15 MW when power prices are favorable. Charging the LP condensate storage system at the minimum load would reduce plant power output by 6.2 MW (1.8%).
Since the locations where the condensate is taken from and returned to the steam turbine cycle are very close to each other, the exergy of the condensate flows at these locations is also similar, resulting in a high power-to-power and energy-to-energy RTE, as presented in
Figure 9. Also, since changes in the power output ΔP
chg and ΔP
dchg during the charging and discharging of the storage tanks remain constant, for the same charging and discharging times, according to Eqn. 3, the roundtrip energy-to-energy efficiency η
RTE is the same as the power-to-power efficiency η
PP.
For the highest analyzed charging flow rate of 40 kg/s and full power plant output, 318 MWhth of heat can be stored in three hours. The amount of stored heat at lower load is lower due to the lower flow and enthalpy of the condensate. At the minimum analyzed load of 332 MW and 40 kg/s, 144 MWhth can be stored.
Simulations were also performed for the scenario where the condensate storage tanks are charged at the minimum load and discharged at the full load to simulate the load shift described in
Figure 2. For the analyzed example, the maximum power increase at full load during the load shift is 10 MW, i.e., lower compared to the case where storage tanks are charged and discharged at full load. The change in net unit efficiency during the LP condensate tank charging / discharging at full load is approximately equal to 1%-point.
In summary, the main advantages of the LP condensate storage system include low cost (storage tanks operate at low pressure, less than 10 bar) and fast response. The change in the plant power output is proportional to the flow rate of the charging/discharging condensate flow, while the heat storage capacity is proportional to the tank volume. The maximum charging and discharging times are proportional to the tank storage volume and condensate charging / discharging flow rates.
2.2. Two-Tank Molten Solar Salt Thermal Energy Storage
The next ESS analyzed consists of two storage tanks for the cold and hot solar molten solar salt (eutectic mixture of the sodium and potassium nitrate, NaNO
3 and KNO
3) and two heat exchangers (HEX
1 and HEX
2). Charging of a two-tank molten solar salt thermal energy storage system (2-tank MSS) integrated with the Reference power plant is shown schematically in
Figure 10, where the red and blue arrows denote the flows of the cold reheat (CRHT) steam used for system charging, while the yellow lines show flows of the molten solar salt. To charge the 2-tank MSS system and decrease the plant power output, a portion of the CRHT steam exhausted from the high-pressure (HP) steam turbine, instead of expanding in the intermediate-pressure (IP) turbine, is diverted to the heat exchanger HXE
1 and used to increase the temperature of the cold molten salt flowing from the Cold Tank. The hot molten solar salt is then stored in the Hot Tank. The tanks operate at a low, near-ambient pressure due to the low vapor pressure of the molten solar salt. Application and integration of the molten salt thermal energy storage with fossil power plants is also described in [
39,
40,
52,
53,
54].
The steam leaving HXE
1 is split into two streams: one merging with the steam extraction line for the first high-pressure feedwater heater (HP FWH1), the other one merging with the steam extraction line for the deaerator (DEA). For the system configuration shown in
Figure 10, the maximum flow rate of charging steam is limited by the sum of these two steam extractions, 40 kg/s for the analyzed system. Thus, tank storage capacity is determined by the charging time. For a charging time of 3 hrs., 1,080 tons of molten solar salt at 350
oC can be stored in the Hot Tank, corresponding to a heat storage capacity of 157 MWh
th. This small storage volume results in a compact system. Heat tracing of piping and thermal insulation of the tanks is needed to maintain the minimum salt temperature above the freezing temperature of 248
oC and reduce thermal losses. Due to lower flows and enthalpies of the CRHT steam and feedwater at the minimum power plant load (336 MW), 123 MW
th of heat can be stored in three hours.
During system discharging, the hot molten salt discharged from the Hot Tank is used in HEX2 to increase the temperature of the feedwater (FW) bypass flow. The plant power output increases due to the reduced steam extractions for the HP FWH1 and/or HP FWH2. The heated FW bypass flow merges with the main FW flow between the HP FWHs 1 and 2, while the cold molten salt leaving the HXE2 is pumped back to the Cold Tank.
The effect of the 2-tank MSS system charging and discharging on the plant power output at full load is presented in
Figure 11 where the plant power output is plotted as a function of time. As shown in
Figure 11, the system is charged for 3 hours using the CRHT steam at a flow rate of 40 kg/s, placed on hold for 30 minutes, and then discharged using a range of the FW bypass flows. The results presented in
Figure 11 show that the plant power output decreases almost instantaneously when the charging steam is diverted from the HP turbine discharge and stays constant during the charging process. The power output increases almost instantaneously to the baseline level when charging stops.
When the stored heat is returned to the cycle during system discharging, the plant power output increases instantaneously. A shown in
Figure 12, the Reference plant power output decreases by 23.2 MW (3.6% of full load) during the system charging. For the range of the FW bypass flows analyzed (30 to 90 kg/s or 5.5 to 16.5% of the BFP discharge flow), the plant power output increases from 5 to 20 MW (0.8 to 3.3%). The plant power increase during discharging is a linear function of the FW bypass flow. Regarding plant flexibility, integration of a 2-tank MSS storage system with the Reference power plant, would at full load, allow the plant to follow load in the approximately ± 20 MW range and increase its participation in the frequency regulation, or to increase the maximum power output by 20 MW when power prices are favorable.
The roundtrip power-to-power and energy-to-energy efficiencies for the full load and maximum charging CRHT steam flow of 40 kg/s are shown in
Figure 13 as function of a dimensionless FW bypass flow. As the results show, the power-to-power roundtrip efficiency
PP increases linearly with the FW bypass flow since the increase in the plant power output with the FW bypass flow during discharging is linear. For the highest bypass flow of 90 kg/s analyzed in this work, η
PP is high, approximately 90%. For low FW bypass flows, η
PP is much lower since, as presented in
Figure 12, the power output increase, ΔP
dchg, during discharge is much lower compared to the power output decrease, ΔP
chg, during charging. This is due to a much higher exergy of the CRHT steam extraction flow used for charging as compared to the exergy of the FW bypass flow and steam extraction flows affected by discharging.
In contrast to η
PP, the roundtrip energy-to-energy efficiency, η
RTE, is relatively low and virtually constant. For the analyzed range of FW bypass flows and selected tank size, the value of η
RTE is in the 25 to 28% range. The roundtrip energy-to-energy efficiency is low due to the significantly shorter tank discharging time compared to the charging time. For example, as shown in
Figure 14, for the low FW bypass flow of 7% of the BFP flow, the Hot Tank is fully discharged in 135 minutes, while for the highest analyzed dimensionless FW bypass flow of 16.5%, the Hot Tank discharging time is 52 minutes, more than three times shorter compared to the charging time of 180 minutes. To satisfy the energy balance for HXE
2, a higher flow of the hot molten salt is needed for the higher FW bypass flow, which results in a shorter discharge time of the Hot Tank, larger value of parameter f
z (Eqn. 5) and a lower value of η
RTE (Eqn. 3). For the constant tank volume, the discharging time is inversely related to the hot molten salt flow from the Hot Tank and thus to the FW bypass flow. Therefore, the energy efficiency should be constant and independent of the FW bypass flow (Eqns. 3 and 4).
The results shown in
Figure 12,
Figure 13 and
Figure 14 are typical of a constant volume storage tank system where the discharging flow of the hot molten solar salt and discharging time are inversely related resulting in a linear variation of P
dchg and
PP with the discharging flow, and a constant value of
RTE. The
RTE value could be increased by increasing the volume of the storage tanks.
The charging and discharging of a 2-tank MSS system also affects cycle efficiency, which for the Reference plant decreases by 1.3%-points during charging and increases from 0.5 to 1.5%-points during discharging. Similar to the LP condensate tank storage system, the change in plant power output for the 2-tank MSS system is instantaneous during system charging and discharging.
The analysis was also performed for the scenario where a 2-tank MSS system was charged at the minimum load and discharged at the full load to simulate the load shift. For the analyzed operating conditions, the power increase of the Reference plant at the maximum load achievable by the load shift (
Figure 2) is 11.3 MW (1.8%) for a FW bypass flow of 90 kg/s, i.e., lower compared to the power increase achieved when charging and discharging is performed at full load. The
PP in this case is 96.7%. However, since the discharging flow is high compared to the charging flow, the discharging time is shorter compared to the charging time, which increases the value of
fz and reduces the value of η
RTE. The results obtained for the 2-tank MSS storage systems are similar to the results reported by Kruger et al. [
40].
2.3. Fixed Bed Thermal Energy Storage
The cost of a TES system is significantly affected by the choice of heat storage medium. A fixed bed (FB) TES system is an attractive option because of the low cost of the heat storage medium, simple operation, and easy scaleup. In the FB TES system, heat is stored in solid media, such as sand, rocks, ceramics, and other solid materials. Air, flue gas (in gas turbine cycle applications), and molten solar salt may be used as the heat transfer fluids (HTFs). As the HTF flows through the fixed bed of the heat storage medium, heat is exchanged between the solid and the HTF. During the charging phase, the hot HTF enters from the top of the bed, flows through the stationary bed of solids, transfers heat to the heat storage medium, and exits at the bottom of the bed at a lower temperature. During the discharging phase, the process is reversed; the cold HTF enters FB from the bottom and flows upward through the stationary bed. The heat stored in the heat storage medium is transferred to the HTF, increasing its temperature. The hot HTF exits FB at the top. The fixed bed energy storage system in stand-alone and integrated configuration is described in [
40,
55,
56,
57].
This mode of operation creates a hot temperature zone in the upper part of the FB and a low temperature zone at the bottom part, thus establishing a thermocline, i.e., temperature difference between the top and bottom FB sections, allowing one storage tank to operate as two storage tanks, one at a high temperature, the other one at a low temperature. According to the results published in the literature, Pacheo, [
11], a single tank system employing thermocline offers cost savings of the order of 30 to 40% compared to the two-tank system.
The integration of a FB TES system with the Reference plant analyzed in this work is presented in
Figure 15 where the FB TES system is charged by using the hot reheat (HRHT) steam extracted from the steam turbine cycle at 538
oC. Since using the HRHT steam does not affect operation of the steam reheater, boiler and HP turbine, such integration arrangement allows for high HRHT steam extraction flows and high storage capacity of a FB TES system. However, operation and power output of the intermediate-pressure (IP) and low-pressure (LP) turbines is affected due to the reduced steam flow rate. Since the heat storage medium (crushed rock, ceramics, etc. [
58]) can operate at high temperature, using high temperature charging steam results in a smaller TES system because heat is stored at a higher enthalpy.
The FB TES system may also be charged by using the main (live) steam (MST) or cold reheat (CRHT) steam. Based on the results published by Krueger et al. [
40], using the HRHT steam gives higher overall efficiency compared to the other two options. The main steam has higher exergy compared to the CRHT steam since it did not do produce work in the HP steam turbine. Using the MST as a source of heat for the FB TES charging and throttling it to a lower pressure to reduce the capital cost of the top heat exchanger (HEX
1) would result in the exergy loss and lower overall efficiency. The use of CRHT steam would, on the other hand, increase the size of the FB TES system due to the lower temperature of the charging steam.
As shown in
Figure 15, during system charging, the charging HTRH steam transfers heat to the HTF in the top heat exchanger HEX
1. The cold charging steam leaving HEX
1 is returned to the LP turbine. To reduce the capital and operating cost of HEX
1 and match the cold steam pressure to the LP turbine inlet pressure, throttling is used. The hot HTF leaving the HEX
1 enters the FB TES and exchanges heat with the heat storage medium. The cold HTF leaving FB TES is returned to HEX
1 inlet for reheating. The HTF flow loop for the charging cycle is presented in red in
Figure 15. During discharging, the HTF flow through the FB TES is reversed. The cold HTF exiting the bottom heat exchanger HEX
2 passes through the FB TES and is heated by the heat stored in the heat storage medium. The hot HTF exiting the FB TES flows through the HEX
2, exchanges heat with the FW bypass flow stream, increasing its temperature, and is returned to the FB TES for reheat. The HTF flow loop for the discharging cycle is presented in light blue in
Figure 15.
A FW bypass configuration was selected since it offers several advantages compared to heating the full feedwater flow, such as: (1) smaller heat exchanger, (2) increased operational flexibility, i.e., the FW bypass flow can be varied to match the desired power output increase during discharging, (3) the HEX2 can be installed and maintained while the power plant is in operation (the same applies to the HEX1).
To model the dynamic performance of the FB TES charging and discharging, a two equation non-equilibrium transient model of the flow and heat transfer in a fixed bed published in [
59] was used. The heat transfer between the solid and HTF was modeled using experimental correlations for the heat transfer coefficient from the literature [
60]. The model was verified against experimental data published by Hänchen at al [
61]. To reduce the individual tank storage volume to a manageable size and improve operating flexibility of the system, a modular approach was selected where six identical beds in a parallel flow arrangement, containing a total of 2,000 metric tons of quartz rock with a total thermal capacity of 175 MWh
th, was used for integration with the Reference power plant. Quartz was selected as the heat storage medium based on the results published in [
61]. This modular approach also allows thermal storage capacity of the FB system to be increased as conditions in the energy market change in the future. In this study, air was selected as the HFT since it greatly simplifies system design and operation.
The modeling results show that the HTF flow rate has a significant effect on the FB TES system performance, both during the charging and discharging. The variation of the HTF temperature at the bed exit with the charging time and HTF flow rate is presented in
Figure 16. As the results show, the charging time and the HTF flow rate have a significant effect on the HTF temperature leaving the bed. For the selected FB size and a reference HTF flow rate of 28 kg/s, the HTF temperature at the bed outlet remains constant for approximately 2.5 hours. For a longer charging time, the HTF bed outlet temperature sharply increases and reaches a maximum value in approximately 4.5 hours when the bed is fully charged. During bed discharging, the HTF temperature is initially constant but sharply decreases for longer discharging times approaching the minimum temperature when the bed is fully discharged. A lower HTF flow increases duration of the constant temperature plateau and bed charging/discharging time, while higher HTF flow has the opposite effect.
As a result of the transient performance of the FB TES, the flow rates of the charging steam and FW bypass flow are not constant during bed charging / discharging but vary with time, as presented in
Figure 17. These results were obtained by integrating the FB TES model with the Reference power plant model. As presented earlier, during the bed charging with a constant HFT flow, the HTF outlet temperature is initially constant but increases sharply as the thermocline effect propagates to the bed top. This increase in the HTF bed outlet temperature (HEX
1 inlet temperature) results in a decrease in the charging HRHT steam flow needed to maintain constant HFT temperature at the HEX
1 outlet, as shown in
Figure 17.
During the bed discharging with a constant HFT flow, the HTF bed outlet temperature is initially constant but decreases sharply as the thermocline propagates to the bed bottom resulting in less heat being available in HEX2 for the FW bypass heating. Thus, to maintain a constant temperature of the FB bypass, the FB bypass flow decreases. The charging steam and FB bypass flow reach their minimum values for the fully charged and discharged bed, respectively. For six fixed beds in service, the maximum total hot reheat steam extraction flow is 92.64 kg/s (6 x 15.44 kg/s).
The variation of the charging HRHT steam flow during bed charging and FW bypass flow during bed discharge affects the operation of the steam turbine cycle and the power plant power output, which is thus not constant as it was the case with the two previously described ESS systems but varies with the charging/discharging time. The impact of the FB TES system charging and discharging on power plant net power output is presented in
Figure 18.
As the results show, changes in plant net power output mirror changes in the charging steam and FB bypass flow. During bed charging, the plant net power output is 43 MW (6.8%) lower compared to the baseline value and remains constant for as long as the charging steam flow is constant. As the charging steam flow decreases, the net power output increases reaching the baseline value when the bed is fully charged in approximately 4.5 hours. During bed discharging, the net power output is approximately 17 MW (2.7%) higher compared to the baseline value as heat is returned to the cycle and remains constant for as long the FW bypass flow is constant. As the FW bypass flow decreases, the net power plant output decreases and returns to the baseline value when the bed is fully discharged.
Variations in plant power output during bed charging and discharging result in variations in net unit efficiency which decreases by 3-% points during charging and increases by 1.2%-points during discharging. It is worth noting that a 1%-point change in net unit efficiency is significant.
Variations in plant power output during charging and discharging also affect the roundtrip power-to-power efficiency. The value of η
PP remains close to the maximum value for as long as the ΔP
chg and ΔP
dchg values corresponding to the FB TES system charging and discharging are approximately constant. As shown in
Figure 19, as the bed is getting closer to the fully charged or fully discharged condition, and flow rates of the charging steam and FW begin to decrease, the value of
PP rapidly decreases approaching zero for the fully charged/discharged bed.
The roundtrip efficiency for the FB TES system charged with the HRHT steam is lower compared to the other energy storage systems analyzed in this work. One of the reasons for this difference is the choice of the charging and discharging sources/locations, i.e., HRHT steam and feedwater flow, with exergy of the charging steam being significantly higher compared to the exergy of the FW bypass flow. However, storing heat at high temperature reduces the size of the heat storage system and capital investment. The value of ηPP is close to the maximum value of 36.4% for as long as ΔPchg and ΔPdchg values, corresponding to charging and discharging, are approximately constant. The average ηPP value is 34.12%. As the bed is getting closer to the fully discharged state and flow rates of the charging steam and FW condensate bypass begin to decrease, the value of ηPP rapidly decreases approaching zero for the fully charged/discharged bed. The average value of the energy-to-energy roundtrip efficiency ηRTE is 31.45% at the full load operating conditions.
This analysis was also performed for the minimum plant power output and for the load shift (Figure 2). Since the changes in the plant power output during the FB TES charging and discharging are approximately the same for the full and minimum load, the roundtrip efficiencies for these operating conditions are also approximately the same. Therefore, for the load shift, the power increase at the minimum load is approximately the same as for the full load. Also, the reduction in the power output during charging at the minimum load is the same as for the full load operation.
As reported in [
62] and from the results of this analysis, low cost of the heat storage medium and simple operation are main advantages of the FB TES. The power output of a FB TES (charge / discharge rate) can be controlled by the HTF flow. However, since the FB TES is a constant volume device, the HTF flow and charging/discharging times are inversely related. Except for the initial period, power output of the FB TES varies during charging/discharging.