A thermal heat-pump grid storage technology is described based on closed-cycle Brayton engine transfers of heat from a cryogenic storage fluid to molten solar salt. Round-trip efficiency, computed as a function of turbomachinery polytropic efficiency and total heat exchanger steel mass, is found to be competitive with that of pumped hydroelectric storage. The cost per engine watt and cost per stored joule based are estimated based on the present-day prices of power gas turbines and market prices of steel and nitrate salt. Comparison is made with electrochemical and mechanical grid storage technologies.
I. INTRODUCTION
It is widely understood that facilities for storing very large amounts of electric energy are key to any long-term plan to supplant fossil fuel with renewable energy.1–3 The desirability of switching to renewable sources is, of course, controversial, and it is not clear that eliminating fossil fuel is economically feasible at this time.4–7 But without storage, it is also impossible.8 The intermittencies of energy sources like wind and sun are fundamentally incompatible with the electricity industry's need to supply power to customers—the instant demands by customers. The overriding importance of timing is demonstrated by the frequent development of negative spot prices for electric energy in markets with large wind deployments.9–12
The purpose of this paper is to discuss a specific storage technology and make the case that it has the right engineering compromises to prevail when the need for storage eventually becomes acute. It is an implementation of pumped thermal storage, an idea already in the literature and under development in industry, and differs from others chiefly in the substitution of heat exchangers for thermoclines13–17 (Fig. 1). Instead of pumping water uphill from a low reservoir to a higher one to store energy, as occurs in pumped hydroelectricity, one pumps heat from a cold body to a hot one by means of a heat engine. In either case, the process is reversible so that energy banked can be withdrawn later to satisfy demand.
The reasoning leading to the “Brayton Battery,” as it might be called, applies the metrics of safety, low cost, and high efficiency, in this order. The average power delivered to a large metropolitan area such as Los Angeles or New York is about 1.4 × 1010 W.18 Storing this power for only one hour gives 5.04 × 1013 J or one Hiroshima-sized atomic bomb.19 It is absolutely essential that explosive release of this stored energy be physically impossible. Once this safety criterion is met, capital and maintenance costs must be brutally minimized, even at the price of a small hit in round-trip efficiency, because storage of electricity is fundamentally about value, not about conserving energy.
The round-trip storage efficiency ηstore of the configuration in Fig. 1 satisfies
where ηc and ηt are the polytropic efficiencies of the compressor and turbine, respectively, and Tdump is the (Kelvin) temperature at which waste heat is sloughed. For the prototype parameters in Table I, this limit is 0.75. The factor of 2 in Eq. (1), which comes from the need to pass the working fluid through the turbomachinery twice in each storage cycle, is the efficiency disadvantage vis-à-vis pumped hydroelectricity or thermoclines.
. | Ar . | N2 . |
---|---|---|
T0 | 180 K | 180 K |
300 K | 300 K | |
T1 | 495 K | 495 K |
823 K | 823 K | |
Tdump | 300 K | 300 K |
ξ | 1.66 | 1.66 |
ηc | 0.91 | 0.91 |
ηt | 0.93 | 0.93 |
p0 | 1.00 × 105 Pa | 1.00 × 105 Pa |
pl | 1.00 × 106 Pa | 1.00 × 106 Pa |
ph | 3.55 × 106 Pa | 6.52 × 105 Pa |
. | Ar . | N2 . |
---|---|---|
T0 | 180 K | 180 K |
300 K | 300 K | |
T1 | 495 K | 495 K |
823 K | 823 K | |
Tdump | 300 K | 300 K |
ξ | 1.66 | 1.66 |
ηc | 0.91 | 0.91 |
ηt | 0.93 | 0.93 |
p0 | 1.00 × 105 Pa | 1.00 × 105 Pa |
pl | 1.00 × 106 Pa | 1.00 × 106 Pa |
ph | 3.55 × 106 Pa | 6.52 × 105 Pa |
The structure of this paper is as follows. Section II discusses materials, with emphasis on the properties of steels and storage fluids and the constraints these impose on operating temperatures. Section III deals with regeneration, the need for which follows from the material constraints. Section IV presents a detailed round-trip efficiency computation. Section V discusses costs. Section VI gives the summary.
II. MATERIALS
A. Steel limitations
Closing the Brayton cycle as shown in Fig. 1 enables the background working fluid pressure to be raised, thus greatly reducing the cost per engine watt. The turbine rotation speed and blade angles fix the working fluid velocity to first approximation, and thus, increasing the pressure simply increases the number of moles of working fluid passing a given point per second and thus the total power. Background pressures as high as 7.7 × 106 Pa (77 atmospheres) compressed to 1.38 × 107 Pa have been experimentally benchmarked in closed-cycle Brayton engines using supercritical CO2.20 Pressures approaching 3.0 × 107 Pa (300 atm) are routine in modern supercritical steam power plants.21
However, the use of high pressure severely limits the temperatures that one can employ. As shown in Fig. 2, raising the temperature of a steel eventually causes it to exhibit creep, a slow plastic deformation that presages full mechanical failure at higher temperatures. This effect is irrelevant on short time scales but is a major design constraint on the scale of 40 years.22 Creep is what prevents conventional carbon steels from being used in high-pressure thermal applications at temperatures above about 700 K (427 °C). This limit can be raised to about 800 K (527 °C) by adding impurities in percent amounts that block grain boundary motion, but maximum resistance to creep requires full alloying with Cr and Ni to make stainless steels.22 The creep limitations of alloy steels may be seen from Fig. 2 to apply universally, including to Inconels. Thus, pressure vessels constructed from steel become problematic for temperatures much above 873 K (600 °C). This is a major consideration in present-day supercritical steam plants.21
B. Solar salt
The temperature limitations of steels make solar salt, the elementary NaNO3/KNO3 eutectic storage medium of the concentrating solar industry, to be a particularly good choice for the high-temperature storage medium.23–25 All solid and liquid substances have heat capacities of approximately 3R per mole of atoms, where R is the ideal gas constant, in the temperature range of interest, so virtually, anything will act as a thermal storage medium, including rocks. However, the cost of either gravel or salt is much less than the cost of the engine for storage times less than 1 day, and salt has the great advantage of being liquid, thus enabling heat transfer by counterflow, minimizing entropy creation. Solar salt has a well-known list of other advantages such as low vapor pressure, high compatibility with steels, and environmental friendliness. It does not disintegrate in response to thermal cycling stresses the way a solid would. It creates no explosion hazard.
The NaNO3/KNO3 phase diagram is shown in Fig. 3.26–28 At the eutectic composition (0.5 M fraction of NaNO3), the melting temperature is 495 K (220 °C), about the same as an ordinary Pb/Sn solder. There is also a decomposition boundary at approximately 823 K (550 °C), a temperature conveniently near the steel creep cliff shown in Fig. 1.29–32 The decomposition boundary is the point where the performance becomes difficult to predict, not where the salt begins to fail. The first stages of NO3 decomposition are reversible and occur via the nitrate/nitrite reaction . Approximately 3% of the nitrate has converted into nitrite at 823 K (550 °C).29 The kinetics of outgassing and re-absorption of O2 have been measured and found not to be unduly rapid.30,31 Nitrogen oxide and N2 evolve irreversibly at high rates above 923 K (700 °C).32 The true upper operating temperature of solar salt exposed to air is not presently known.32 The corrosion of stainless steels by these salts is mild and amounts to approximately 10 μm per year at 823 K (550 °C), chiefly due to oxidation.33,34
C. Hydrocarbon cryogen
Equation (1) dictates that T1 − T0 is as large as possible and ξ is as small as possible. For reasons explained in detail in Sec. III, the properties of solar salt determine and T1 = 495 K, and thus, . It remains to maintain T0 as low as possible. However, once T0 gets so low that , these two temperatures effectively switch roles. must then be greater than Tdump since there would otherwise be no way to re-initialize the engine after it had stopped and equilibrated with the environment. Thus, we have T0 = 300 K/ξ = 180 K.
For the purposes of this paper, a second storage fluid will be prototyped as n-hexane. All fluids that remain liquid between 180 K and 300 K with vapor pressure less than 1 atm are hydrocarbons or derivatives of them, so their physical properties are similar; they all have health and safety issues that must be weighed. Hexane, a component of gasoline, is pervasive, biodegradable, and cheap, and it is widely used in the food industry for extracting oils from seeds, but it is also highly flammable and neurotoxic.38–42 Anhydrous methanol, to which hexane is physically similar and which is an acceptable substitute, is vastly more toxic. All the lower alcohols—methanol, ethanol, propanol, and propylene glycol—are hygroscopic and become viscous, like honey, at cryogenic temperatures when mixed with water.43,44 Hydrocarbons can be rendered flame-resistant through halogenation, and a fluid such as CH2Cl2 would work perfectly well although it would increase the stress on stratospheric ozone.45,46
D. Working fluid
The working fluid in this application is limited to gases that are extremely stable at high temperatures and far from liquefaction or solidification phase transitions at mildly cryogenic ones. The mechanical advantages of working near a critical point are outweighed in this case by danger of fluid raining or snowing out and damaging the turbine blades. This eliminates, in particular, CO2, which has both liquefaction and freezing transitions in the range of 200–300 K. These considerations plus requirements of chemical inertness, cheapness, and environmental friendliness restrict the possibilities to Ar and N2. The advantages of a lowered compression ratio and, potentially, higher adiabatic efficiency by using a gas with no internal degrees of freedom point to Ar as the preferred choice. However, N2 has the advantage of requiring only minor modifications of rotor and stator airfoil shapes already optimized for use in air-breathing jet engines.54
Table II shows the thermophysical properties of Ar and N2 over the temperature and pressure range of interest.53 Both gases are functionally ideal (and thus polytropic) except at the lowest temperatures. The correct low-temperature equation of state must be used in fine-tuning the turbomachinery, but the ideal equation of state is sufficient for making efficiency and cost estimates.
Ar . | p = 1.0 × 106 Pa . | p = 3.55 × 106 Pa . | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
T . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . |
(K) . | (Pa s) . | (W m−1 K−1) . | (Pa s) . | (W m−1 K−1) . | ||||||
180 | 1.482 × 10−5 | 1.178 × 10−2 | 1.041 | 1.769 | 0.715 | 1.620 × 10−5 | 1.406 × 10−2 | 1.180 | 2.167 | 0.872 |
240 | 1.898 × 10−5 | 1.501 × 10−2 | 1.014 | 1.716 | 0.685 | 1.981 × 10−5 | 1.634 × 10−2 | 1.055 | 1.857 | 0.731 |
300 | 2.286 × 10−5 | 1.801 × 10−2 | 1.005 | 1.695 | 0.676 | 2.347 × 10−5 | 1.903 × 10−2 | 1.019 | 1.771 | 0.696 |
398 | 2.865 × 10−5 | 2.251 × 10−2 | 0.999 | 1.681 | 0.670 | 2.908 × 10−5 | 2.329 × 10−2 | 0.999 | 1.718 | 0.676 |
495 | 3.385 × 10−5 | 2.656 × 10−2 | 0.997 | 1.675 | 0.668 | 3.418 × 10−5 | 2.720 × 10−2 | 0.992 | 1.697 | 0.670 |
659 | 4.174 × 10−5 | 3.269 × 10−2 | 0.996 | 1.671 | 0.667 | 4.198 × 10−5 | 3.321 × 10−2 | 0.990 | 1.681 | 0.666 |
823 | 4.724 × 10−5 | 3.700 × 10−2 | 0.996 | 1.667 | 0.666 | 4.749 × 10−5 | 3.753 × 10−2 | 0.989 | 1.673 | 0.665 |
Ar . | p = 1.0 × 106 Pa . | p = 3.55 × 106 Pa . | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
T . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . |
(K) . | (Pa s) . | (W m−1 K−1) . | (Pa s) . | (W m−1 K−1) . | ||||||
180 | 1.482 × 10−5 | 1.178 × 10−2 | 1.041 | 1.769 | 0.715 | 1.620 × 10−5 | 1.406 × 10−2 | 1.180 | 2.167 | 0.872 |
240 | 1.898 × 10−5 | 1.501 × 10−2 | 1.014 | 1.716 | 0.685 | 1.981 × 10−5 | 1.634 × 10−2 | 1.055 | 1.857 | 0.731 |
300 | 2.286 × 10−5 | 1.801 × 10−2 | 1.005 | 1.695 | 0.676 | 2.347 × 10−5 | 1.903 × 10−2 | 1.019 | 1.771 | 0.696 |
398 | 2.865 × 10−5 | 2.251 × 10−2 | 0.999 | 1.681 | 0.670 | 2.908 × 10−5 | 2.329 × 10−2 | 0.999 | 1.718 | 0.676 |
495 | 3.385 × 10−5 | 2.656 × 10−2 | 0.997 | 1.675 | 0.668 | 3.418 × 10−5 | 2.720 × 10−2 | 0.992 | 1.697 | 0.670 |
659 | 4.174 × 10−5 | 3.269 × 10−2 | 0.996 | 1.671 | 0.667 | 4.198 × 10−5 | 3.321 × 10−2 | 0.990 | 1.681 | 0.666 |
823 | 4.724 × 10−5 | 3.700 × 10−2 | 0.996 | 1.667 | 0.666 | 4.749 × 10−5 | 3.753 × 10−2 | 0.989 | 1.673 | 0.665 |
N2 . | p = 1.0 × 106 Pa . | p = 6.52 × 106 Pa . | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
T . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . |
(K) . | (Pa s) . | (W m−1 K−1) . | (Pa s) . | (W m−1 K−1) . | ||||||
180 | 1.208 × 10−5 | 1.752 × 10−2 | 1.033 | 1.462 | 0.759 | 1.489 × 10−5 | 2.465 × 10−2 | 1.238 | 2.010 | 0.990 |
240 | 1.519 × 10−5 | 2.205 × 10−2 | 1.010 | 1.429 | 0.735 | 1.667 × 10−5 | 2.597 × 10−2 | 1.051 | 1.607 | 0.794 |
300 | 1.804 × 10−5 | 2.617 × 10−2 | 1.002 | 1.417 | 0.728 | 1.909 × 10−5 | 2.900 × 10−2 | 1.002 | 1.509 | 0.753 |
398 | 2.224 × 10−5 | 3.232 × 10−2 | 0.998 | 1.406 | 0.724 | 2.297 × 10−5 | 3.427 × 10−2 | 0.979 | 1.450 | 0.734 |
495 | 2.598 × 10−5 | 3.804 × 10−2 | 0.996 | 1.396 | 0.724 | 2.654 × 10−5 | 3.952 × 10−2 | 0.973 | 1.422 | 0.729 |
659 | 3.162 × 10−5 | 4.742 × 10−2 | 0.996 | 1.377 | 0.727 | 3.203 × 10−5 | 4.847 × 10−2 | 0.973 | 1.389 | 0.730 |
823 | 3.666 × 10−5 | 5.658 × 10−2 | 0.997 | 1.358 | 0.732 | 3.699 × 10−5 | 5.738 × 10−2 | 0.975 | 1.365 | 0.733 |
N2 . | p = 1.0 × 106 Pa . | p = 6.52 × 106 Pa . | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
T . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . | μ . | κ . | ρRT/(mνp) . | γ . | Pr . |
(K) . | (Pa s) . | (W m−1 K−1) . | (Pa s) . | (W m−1 K−1) . | ||||||
180 | 1.208 × 10−5 | 1.752 × 10−2 | 1.033 | 1.462 | 0.759 | 1.489 × 10−5 | 2.465 × 10−2 | 1.238 | 2.010 | 0.990 |
240 | 1.519 × 10−5 | 2.205 × 10−2 | 1.010 | 1.429 | 0.735 | 1.667 × 10−5 | 2.597 × 10−2 | 1.051 | 1.607 | 0.794 |
300 | 1.804 × 10−5 | 2.617 × 10−2 | 1.002 | 1.417 | 0.728 | 1.909 × 10−5 | 2.900 × 10−2 | 1.002 | 1.509 | 0.753 |
398 | 2.224 × 10−5 | 3.232 × 10−2 | 0.998 | 1.406 | 0.724 | 2.297 × 10−5 | 3.427 × 10−2 | 0.979 | 1.450 | 0.734 |
495 | 2.598 × 10−5 | 3.804 × 10−2 | 0.996 | 1.396 | 0.724 | 2.654 × 10−5 | 3.952 × 10−2 | 0.973 | 1.422 | 0.729 |
659 | 3.162 × 10−5 | 4.742 × 10−2 | 0.996 | 1.377 | 0.727 | 3.203 × 10−5 | 4.847 × 10−2 | 0.973 | 1.389 | 0.730 |
823 | 3.666 × 10−5 | 5.658 × 10−2 | 0.997 | 1.358 | 0.732 | 3.699 × 10−5 | 5.738 × 10−2 | 0.975 | 1.365 | 0.733 |
III. REGENERATION
The practical material limitations described in Sec. II require . As shown in Fig. 6, this condition causes heat transfers on the high-pressure and low-pressure sides of the circuit to overlap, thus eliminating the need to actually transfer heat to or from the storage fluids over this temperature range. Instead, heat may simply be transferred directly from one side of the circuit to the other through a gas-gas heat exchanger, referred to as a regenerator, or a recuperator.58 In the limit that the entropy generation by the heat exchangers is zero, regeneration has no effect at all on the round-trip efficiency but simply reduces the amount of heat exchanger steel required. It also reduces the temperature ranges over which the storage fluids are required to be liquid. A modified version of Fig. 1 with regeneration included is shown in Fig. 7.
Figures 6 and 7 clarify why cannot exceed Tdump. An engine powered down and equilibrated to Tdump could easily be re-initialized by heating the (frozen) solar salt from Tdump to T1, but it could not be easily re-initialized if had to be cooled to a value less than Tdump.
Figure 6 justifies choosing T1 to be the solar salt melting point. With and fixed at the solar salt decomposition temperature and ambient, respectively, Eq. (1) is maximized when is minimized. However, it may be seen from the plot of this function in Fig. 8 that the minimum occurs when ξ → 1. Other costs not included in the calculation run away in this limit, so the specific value ξ = 1 is not meaningful. However, the convergence of the function itself to 1 in this limit has an important implication that there is no significant round-trip efficiency penalty for 1 < ξ < 2. Virtually, any value of ξ in this range will do. Moreover, lowering T1 degrades the efficiency rather than improving it. There is thus no efficiency advantage in employing specialty salts with lower melting temperatures.
IV. ROUND-TRIP EFFICIENCY
A. Fictive temperature
The round-trip storage efficiency ηstore is most simply computed through the entropy budget. The system entropy S must be the same before and after a storage cycle because it is a property of state, so any entropy ΔS generated during the cycle must be discarded as waste heat TdumpΔS, where Tdump is the slough temperature, roughly ambient. This heat represents a loss from the stored energy Estore that cannot be re-transmitted as grid power. We thus have
The main contributions to the fictive temperature Tf come from the turbomachinery and three heat exchangers,
Other losses such as generator/motor inefficiency, unmanaged turbulence, and viscous drag in the conduits are percent-level corrections.
B. Turbomachinery
In an adiabatic (entropy-conserving) compression or expansion, the temperature change dT resulting from a pressure change dp satisfies
where γ = cp/cv. In the limit that the turbomachinery is perfectly adiabatic, the temperature ratio ξ in Fig. 1 is thus related to the ratio of the upper and lower pressures ph/pl by
In a real turbocompressor, Eq. (4) is modified to
where ηc is the compressor polytropic efficiency. This slightly elevates the discharge temperature with respect to the ideal one for a given compression ratio. Integrating between the ideal temperature and the actual one at a constant pressure gives the entropy created when ν moles of gas pass through the compressor
From the version of Eq. (6) appropriate for the turbine
one obtains similar entropy created when ν moles of gas expand through the turbine
Since the energy stored by cycling ν moles of gas is , the inverse fictive temperature of the turbomachinery is
The factor of 2 accounts for summing the entropy generated during the storage and extraction steps. Equation (1) is obtained by substituting this expression into Eq. (2).
The parameters of Table I give
The values of ηc and ηt leading to this result lie at the extreme high end of present-day gas turbine technology.59–63 They are benchmarked in the GE 90 aircraft engine. They are not characteristics of turbomachinery generally, however, so custom design is required to achieve them in practice. Turbomachinery bought off the shelf is a product fine-tuned to specific market needs that require a different set of compromises. There are theoretical indications that each ηc and ηt could be improved a further 1%.64–67
C. Heat exchangers
The fictive temperature of a heat exchanger differs from that of the turbomachinery depending on both the size and the power. The physical principles are elementary and equally applicable to all heat exchanger designs, so the simple shell-and-tube prototype of Fig. 9 suffices for estimating the amount of steel required to achieve a given fictive temperature, even if more advanced designs are actually employed.70,71
The laminar and turbulent cases are both important. The Reynolds number inside a heat exchanger tube is
where a is the tube's inner radius, N is the number of tubes, and is the number of moles of working fluid passing through the circuit per unit time. The flow is laminar if Re < 2000, turbulent if Re > 3000, and intermittent otherwise. For this particular application, the effects of turbulence are adequately accounted for by replacing μ and κ in the laminar expressions by the Darcy-Weisbach formula55,56
with the Swamee-Jain approximation for the Darcy friction factor
and the Gnielinski correlation
where Pr = μcp/(mνκ) is the Prandtl number. Figure 10 shows these modifications to μ and κ for various values of the tube surface roughness parameter ϵ.
The laminar case follows from elementary considerations. Assuming a pressure gradient ∂p/∂z along the tube and solving the Navier-Stokes equation
with a no-slip boundary condition, one obtains Hagen-Poiseuille flow
and an entropy generation due to viscous drag inside the tubes of
Assuming a temperature gradient ∂T/∂z along the tube and similarly solving the heat flow equation
one obtains the Graetz solution
from which the entropy generation due to the thermal resistance inside the tubes is computed to be
These two contributions to become equal when the heat exchanger length is L0, defined by
where ΔT is the temperature difference between the two ends of the heat exchanger and
is the working fluid scattering mean free path.57 The entropy creation outside the tube is a multiple of obtained by solving Eqs. (16) and (19) numerically for a given value of d. The result is summarized in Table III. The fictive temperature is then
where is the engine power. The factor of 2 accounts for entropy creation during both the storage and extraction steps. Entropy creation due to heat flow in the steel is too small to matter in this application.
Ar . | Low . | Regen. . | High . |
---|---|---|---|
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 3.55 × 106 Pa | 3.55 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 4.259 | 2.500 |
0.0182 | 0.6429 | 0.0119 | |
ℓ | 8.27 × 10−9 m | 4.62 × 10−9 m | 8.54 × 10−9 m |
a | 0.0015 m | 0.0015 m | 0.0015 m |
L0 | 55.2 m | 96.5 m | 54.6 m |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 2000 | 2000 | 2000 |
0.0997 W−1 | 0.0651 W−1 | 0.0463 W−1 | |
0.0551 kg W−1 | 0.0540 kg W−1 | 0.0629 kg W−1 | |
30 520 K | 30 880 K | 31 990 K | |
N2 | Low | Regen. | High |
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 6.52 × 106 Pa | 6.52 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 5.047 | 2.500 |
0.0268 | 0.7943 | 0.0209 | |
ℓ | 8.97 × 10−9 m | 2.72 × 10−9 m | 4.91 × 10−9 m |
a | 0.0015 m | 0.0015 m | 0.0015 m |
L0 | 70.7 m | 221.4 m | 137.5 m |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 2000 | 2000 | 2000 |
0.0632 W−1 | 0.0407 W−1 | 0.0290 W−1 | |
0.0350 kg W−1 | 0.0338 kg W−1 | 0.0393 kg w−1 | |
29 540 K | 28 050 K | 32 150 K |
Ar . | Low . | Regen. . | High . |
---|---|---|---|
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 3.55 × 106 Pa | 3.55 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 4.259 | 2.500 |
0.0182 | 0.6429 | 0.0119 | |
ℓ | 8.27 × 10−9 m | 4.62 × 10−9 m | 8.54 × 10−9 m |
a | 0.0015 m | 0.0015 m | 0.0015 m |
L0 | 55.2 m | 96.5 m | 54.6 m |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 2000 | 2000 | 2000 |
0.0997 W−1 | 0.0651 W−1 | 0.0463 W−1 | |
0.0551 kg W−1 | 0.0540 kg W−1 | 0.0629 kg W−1 | |
30 520 K | 30 880 K | 31 990 K | |
N2 | Low | Regen. | High |
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 6.52 × 106 Pa | 6.52 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 5.047 | 2.500 |
0.0268 | 0.7943 | 0.0209 | |
ℓ | 8.97 × 10−9 m | 2.72 × 10−9 m | 4.91 × 10−9 m |
a | 0.0015 m | 0.0015 m | 0.0015 m |
L0 | 70.7 m | 221.4 m | 137.5 m |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 2000 | 2000 | 2000 |
0.0632 W−1 | 0.0407 W−1 | 0.0290 W−1 | |
0.0350 kg W−1 | 0.0338 kg W−1 | 0.0393 kg w−1 | |
29 540 K | 28 050 K | 32 150 K |
There is no need to make the heat exchanger fictive temperature more than about 30 times . A set of parameters that give in this range for laminar flow is shown Table III. The important quantity for costing purposes is the steel mass per engine watt
Here, is the number of tubes per engine watt, given by
per Eq. (12). Both and remain unchanged if L is made smaller by keeping NL fixed (which requires lowering Re). Thus, Table III actually describes a family of designs with different aspect ratios but with similar performance characteristics.
Making L longer and keeping NL fixed push Re over the turbulence threshold. As shown in Table IV, this causes a 50% reduction in for the same value of Tf. This occurs because the turbulent enhancement of the thermal conductivity matters more than the turbulent enhancement of the viscosity when L < L0. Thus, the optimal design with respect to steel mass has Re just above the turbulence threshold. Further, cost compromise at this value of Re may result in shortening of L keeping N fixed, reducing both Tf and proportionately.
Ar . | Low . | Regen. . | High . |
---|---|---|---|
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 3.55 × 106 Pa | 3.55 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 4.258 | 2.500 |
0.0404 | 0.6429 | 0.0253 | |
1.78 × 10−8 m | 9.91 × 10−9 m | 1.82 × 10−8 m | |
a | 0.0015 m | 0.0015 m | 0.0015 m |
25.7 m | 45.0 m | 25.6 m | |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 3000 | 3000 | 3000 |
0.0665 W−1 | 0.0434 W−1 | 0.0309 W−1 | |
0.0367 kg W−1 | 0.0360 kg W−1 | 0.0419 kg W−1 | |
31 040 K | 34 440 K | 32 420 K | |
N2 | Low | Regen. | High |
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 6.52 × 106 Pa | 6.52 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 5.047 | 2.500 |
0.0610 | 0.7942 | 0.0399 | |
1.95 × 10−8 m | 5.91 × 10−9 m | 1.07 × 10−8 m | |
a | 0.0015 m | 0.0015 m | 0.0015 m |
32.5 m | 101.7 m | 63.3 m | |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 3000 | 3000 | 3000 |
0.0421 W−1 | 0.0272 W−1 | 0.0193 W−1 | |
0.0233 kg W−1 | 0.0225 kg W−1 | 0.0262 kg W−1 | |
33 920 K | 39 800 K | 44 280 K |
Ar . | Low . | Regen. . | High . |
---|---|---|---|
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 3.55 × 106 Pa | 3.55 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 4.258 | 2.500 |
0.0404 | 0.6429 | 0.0253 | |
1.78 × 10−8 m | 9.91 × 10−9 m | 1.82 × 10−8 m | |
a | 0.0015 m | 0.0015 m | 0.0015 m |
25.7 m | 45.0 m | 25.6 m | |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 3000 | 3000 | 3000 |
0.0665 W−1 | 0.0434 W−1 | 0.0309 W−1 | |
0.0367 kg W−1 | 0.0360 kg W−1 | 0.0419 kg W−1 | |
31 040 K | 34 440 K | 32 420 K | |
N2 | Low | Regen. | High |
T | 240 K | 398 K | 659 K |
ΔT | 120 K | 195 K | 327 K |
p(in) | 1.00 × 106 Pa | 6.52 × 106 Pa | 6.52 × 106 Pa |
p(out) | 1.00 × 105 Pa | 1.00 × 106 Pa | 1.00 × 105 Pa |
σmax | 1.30 × 108 Pa | 1.30 × 108 Pa | 1.00 × 108 Pa |
b/a | 1.224 | 1.224 | 1.491 |
d/b | 2.500 | 5.047 | 2.500 |
0.0610 | 0.7942 | 0.0399 | |
1.95 × 10−8 m | 5.91 × 10−9 m | 1.07 × 10−8 m | |
a | 0.0015 m | 0.0015 m | 0.0015 m |
32.5 m | 101.7 m | 63.3 m | |
L | 20.0 m | 30.0 m | 20.0 m |
Re | 3000 | 3000 | 3000 |
0.0421 W−1 | 0.0272 W−1 | 0.0193 W−1 | |
0.0233 kg W−1 | 0.0225 kg W−1 | 0.0262 kg W−1 | |
33 920 K | 39 800 K | 44 280 K |
The choice of tube inner and outer radii a and b has no effect on Tf if Re and L are held fixed (and if L < L0), but it has a large effect on . Accordingly, the total cost of steel is minimized by making a as small as possible. The minimum value of b/a required for creep resistance is determined by
where Δp is the pressure difference between the inside and outside of the tube and σmax is the maximum stress allowed in the tube steel, shown in Fig. 2. To this minimal, b/a must be added a mill tolerance margin, the default for which is a factor of 1.25 (12.5% of diameter), and a salt corrosion margin, which is 4.0 × 10−4 m for a 40-year lifespan.68,69 Thus, it is impractical to make a much smaller than 1.5 × 10−3 m in the high-temperature heat exchanger. The regenerator and low-temperature heat exchanger may involve microchannel designs of a = 5.0 × 10−4 m, thus potentially lowering in those cases.70,71
V. COST
The cost of grid storage has two distinct metrics: the cost per engine watt and the cost per stored joule . The first is characteristic of the engine and depends only on how fast energy is transferred to and from the grid, not on how much energy is stored or for how long. The second characterizes the storage medium and is completely independent of how fast energy is transferred in or out. One imagines first building the engine at a certain power rating (and cost) and then adding as much storage capacity as circumstances warrant.
A crude estimate of the cost per engine watt is
The definitions of these parameters and their values are summarized in Tables I and V.72–77 The idea behind this expression is that the turbomachinery should cost about the same as a present-day commercial gas turbine of the same size. Thus, one scales the present-day gas turbine cost by the number of moles per second required to give the given power. The marginal cost of the (very large) heat exchangers required is estimated at 2 times the cost of the steel used to make them. This is unrealistically low in the context of present-day heat exchanger markets, but it is a reasonable expectation for large-volume production. This is consistent with $2 kg−1 implicit in the heat-exchanger price figures in the study by Loh et al., assuming a tube width of 0.125 times the tube diameter, the standard ASME milling margin.78
. | Ar . | N2 . |
---|---|---|
Tf | 1089 K | 1107 K |
0.1136 kg W−1 | 0.0720 kg W−1 | |
873 K | 823 K | |
ηGT | 0.38 | 0.038 |
$0.25 W−1 | $0.25 W−1 | |
$1.00 kg−1 | $1.00 kg−1 | |
$0.27 W−1 | $0.20 W−1 | |
$0.61 kg−1 | $0.61 kg−1 | |
$0.70 kg−1 | $0.70 kg−1 | |
$3.54 × 10−6 J−1 | $3.54 × 10−6 J−1 |
. | Ar . | N2 . |
---|---|---|
Tf | 1089 K | 1107 K |
0.1136 kg W−1 | 0.0720 kg W−1 | |
873 K | 823 K | |
ηGT | 0.38 | 0.038 |
$0.25 W−1 | $0.25 W−1 | |
$1.00 kg−1 | $1.00 kg−1 | |
$0.27 W−1 | $0.20 W−1 | |
$0.61 kg−1 | $0.61 kg−1 | |
$0.70 kg−1 | $0.70 kg−1 | |
$3.54 × 10−6 J−1 | $3.54 × 10−6 J−1 |
A crude estimate for the marginal cost per stored joule is
The definitions of these parameters and their values are summarized in Tables I and V. The idea behind this expression is that the asymptotic cost per stored joule is simply the cost of the medium in which it is stored. Specifically omitted from this estimate because they are too small to matter are the costs of large storage tanks ($50 m−3) and excavation ($2.4 m−3).78–80
Equations (28) and (29) are oversimplified, and they leave out many obvious costs—flow handling, cooling structures, tanks, insulation, pumps, site preparation, and small loss accounting. However, they are sufficiently accurate to reveal the broad-brush picture: the cost of the turbomachinery is lowered by the elevated pressure in the closed Brayton loop so much that the cost per engine watt becomes dominated by the cost of the heat exchangers. The latter are conceptually trivial and scale up easily to arbitrarily large sizes. They become arbitrarily efficient when they do. Heat exchangers large enough to contribute negligibly to the total entropy budget can be built for a total cost per engine watt comparable to that of a present-day simple-cycle gas turbine. The marginal cost per stored joule, dominated by the cost of the storage fluids, is about $3.54 × 10−6 J−1 ($12.7 per kWh).
The plant cost curves associated with Table V are shown in Fig. 11. With the understanding that these are only illustrations, on account of the large error bars, one can see that the cost of the storage is negligible for storage times less than a day and thus that only the costs per engine watt and plant infrastructure matter. There is, in particular, no significant economic advantage in substituting rocks for solar salt and hexane.
VI. DISCUSSION
A. Physical constraints
Which storage technology actually prevails in the end will depend ultimately on cost, and this is something difficult to assess correctly without actually building machines and testing them. Thus, it is not possible to make a purely economic case that the technology described here, which exists only as a concept, will win out. Rather the argument rests partly on such cost analysis as one can do combined with a little common sense.
It is highly reasonable, for example, from a physics perspective that the mechanical parts of thermal storage with heat exchange should cost less than pumped hydroelectric storage, the technology with which it is most closely related. The turbines are smaller by virtue of turning faster and having greater blade and fluid velocities. They require no burners or blade cooling. Salt and hexane store energy more compactly than water does when pumped uphill: one kg of water lifted 380 m, a typical pumped storage elevation drop, stores 0.7% of the energy that one kg of nitrate salt does when heated from T1 to . One kg of water falling 380 m transmits 3.4% of the energy to the turbine blades that 1 kg of Ar working fluid does when it travels around the Brayton circuit. In the case of N2 working fluid, it is 1.7%. Thermal storage also uses less land than pumped hydroelectricity does—and, of course, requires no mountains or water supplies. This is shown explicitly in Fig. 12, where the footprint of a prototypical storage facility based on oil depot storage tanks is sketched out.79 Depending on construction details, the thermal storage land requirement can be 10% or less than the equivalent pumped hydroelectric storage requirement, reckoned from the upper reservoir area.81
It is also reasonable that the heat exchangers should cost only slightly more than the steel used to make them. This is so even though large heat exchangers with the specifications in Table IV do not presently exist as products at any price. Heat exchangers are the mechanical engineering equivalent of semiconductor memory chips or flat-panel displays: one makes them by repeating the same simple step over and over again millions of times. This means that their manufacture can be automated. Their manufacture has not yet been automated, but this is only because there is no market for such products. It would be quite straightforward to program industrial robots to perform this task at an extremely low cost. A facility with the size shown in Fig. 12 would require about 3.0 × 107 tubes of length 20 m or a total length of 6.0 × 108 m, enough to circle the earth 15 times.
B. Heat versus electrochemistry
Were batteries ever to beat the marginal cost per stored joule of pumped hydroelectric storage, the comparison with the latter would become moot. However, they have not done so yet. The continued construction of new pumped storage facilities around the world demonstrates this.82–85
There are two reasons why the battery cost problem remains so intractable. One is that all batteries require an expensive internal infrastructure that cannot be eliminated without causing the battery to short and explode. The stored energy per unit mass of a battery varies between 1.4 × 105 J kg−1 for cheap lead-acid car batteries to 9.5 × 105 J kg−1 for high-end Li-ion batteries.86–88 This is not significantly different from of nitrate salt, notwithstanding the fact that an active atom in a battery stores 10 times more energy than the atom of nitrate salt does. The extra mass is due to atoms that do not actually store any energy but guide electrons to the atoms that do. The second is that battery electrodes mechanically degrade, particularly if the battery is deeply cycled, thus creating long-term maintenance issues.89 This occurs because the electrode surface, the place where electron motion converts into ion motion, is a scene of terrible violence on the scale of the atoms. Flow batteries and liquid metal batteries reduce the collateral damage of this violence, but they do so by means of engineering compromises that increase other costs.90–93 It is not clear whether they will be cheaper than mass-produced Li-ion batteries.94
Batteries also have environmental issues associated with their metal ions that have led to mandatory recycling and the banning of household battery disposal in landfills.95–99 This issue does not exist with pumped thermal storage with heat exchange. If the facility shown in Fig. 12 had a catastrophic tank breach (and no fire), the stored energy would dissipate harmlessly as heat, and the result would be a patch of cold nitrate fertilizer that could be easily cleaned up and re-used. If, on the other hand, a vanadium flow battery of the same capacity (two such tanks are required) had such a breach, 8.0 × 106 kg of vanadium ions would be dumped on the ground along with a comparable mass of sulfuric acid.
Thus, while batteries have an advantage over all other forms of storage at small scales in having no significant entry cost per engine watt, this advantage disappears once the scale becomes large enough that the cost entry barrier is surmounted.
C. Explosion danger
With the exception of hydrogen electrolysis, which has cost and electrode issues similar to batteries, all other methods for storing grid-scale energy have nuclear-scale explosion dangers.100 This includes, in particular, flywheels, high pressure tanks, and all purely electrical media such as supercapacitors and superconducting magnetic coils.101,102 These media also have lower energy storage densities, which is a secondary concern. In the case of flywheels and pressure vessels made out of steel, the maximum energy stored per unit mass of steel is σmax/ρsteel or about 2.5% of the energy stored thermally per unit mass of nitrate salt. For supercapacitors, this factor is about 10%.
Deliberately excluded from the list of explosive technologies is compressed air storage in underground caverns.103,104 This is a special case both because it is underground (and thus not explosive) and because it is physically equivalent to thermal storage with heat exchange. It is well known that the energy expended in compressing any gas is stored in its heat. This is why adiabatically compressing N2 from 1.0 × 105 Pa to 7.0 × 106 Pa (70 atmospheres), the typical underground storage pressure, raises its temperature from 300 K to 1101 K. Placing such hot gas underground just to cool off there would make no sense. Thus, all compressed air storage technologies with high round-trip efficiencies employ above-ground heat exchangers like those in Figs. 1 and 7 to cool the gas (extract energy from it) before pumping it underground.105–109 Heat is then added back to the gas as it expands in the extraction step. The physical difference between underground storage and pure thermal storage with heat exchange is that the latter sends the working fluid through the Brayton cycle twice so as to eliminate the need to store working fluid at pressure at all—and thus to eliminate the need for the cavern.
Pumped thermal storage with heat exchange does have explosion danger. It is associated with the heat exchangers, however, not the storage media, so it scales with engine power rather than total stored energy. The parameters in Table IV give a total explosive energy of 90 s worth of generation at the design power, whatever it is, for Ar and 166 s for N2. Thus, for the configuration of Fig. 12, the explosive power for Ar is 2.25 × 1010 J or 4.8 tons of TNT. Heat exchange units of this size and pressure are common in the petrochemical industry, and it is known that they explode rarely, but that when they do the result is catastrophic.110 Thus, precautions must be taken to make sure that any tank explosion that might occur does not cascade, for example, by siting the units underground.
Another serious difficulty is the cost associated with managing the working fluid inventory in the case of breach. Both Ar and N2 are asphyxiating gases. They are quite deadly until they dissipate in the atmosphere. The total working fluid inventory in the case of Fig. 12 is 2.7 × 105 kg of Ar. For comparison, the total amount of CO2 released in the Lake Nyos disaster is estimated to be 109 kg.111
D. Additional technical issues
The power of pumped thermal storage with heat exchange is governed by adjusting the working fluid inventory up and down. The large heat capacity of the heat exchanger steel and corresponding slow thermal response call for regulating the storage fluid flow so as to keep the temperatures fixed. The power the working fluid delivers to the grid or absorbs from it is then directly proportional to the number of moles passing a given point per second. Since a motor/generator connected to the grid is phase-locked with other generators through electric forces propagating through the grid at the speed of light, the working fluid flow velocity is essentially fixed, and this means that one governs the power by reducing or increasing the background pressure of the working fluid in the circuit.
Two sets of turbomachinery may be required, one for storage and the other for extraction. For the parameters in Tables IV and V, this would increase the cost per engine watt by about $0.05 W−1. In contrast to the situation in pumped hydroelectricity, the turbomachinery in this case cannot be automatically reversed because the blades are air foils carefully crafted for maximum efficiency under specific operating conditions, notably flow direction and Mach number. It is possible to design reversible airfoils, but it is not presently known how much efficiency compromise would be required to make machinery that worked equally well in both directions. The worst-case scenario is that no set of air foils would perform this task adequately, in which case two sets of turbomachinery would be required. The cost of doubling the turbomachinery becomes less and less of an issue as the pressure is raised.
The Brayton cycle requires closing. This is most conveniently accomplished using slightly different compression ratios for storage and generation, as shown in Fig. 13. Choosing these to make match causes the heat to be dumped at the lowest possible temperature. In place of Eq. (1), we then have
The parameters in Table I give ηstore < 0.75, in agreement with Eq. (1). This number is too high, however, because cannot lie below ambient. Taking ambient to be 300 K, we obtain T0 = 195 K and thus ηstore < 0.72.
In contrast to pumped hydroelectric and electrochemical storage, pumped thermal storage does not actually lose any energy but simply degrades it into heat that can be used for other things. For example, the rejection temperature (82 °C) implicit in Fig. 13 is appropriate for preparing hot water. Cooling towers, such as those required for steam plants, are unnecessary in this case because there is no need to make a Rankine vacuum.21 Simple thermal cooling ponds will do.
E. 19th century science
Insofar as the numbers I have presented in this paper are correct, they demonstrate that energy storage is a problem of 19th century science. No future laboratory breakthroughs or discoveries are required for solving it. All that is needed is fine engineering and assiduous attention to detail. Said poetically, this is 21st century rocket science.
Moreover, it is clear from Fig. 11 that the storage capacity of months becomes feasible once the engine (including the heat exchangers) exists as a product one can purchase at a known cost, particularly if the heat is further transferred into cheaper media for longer-term storage, such as rocks underground. Thus, pumped thermal storage with heat exchange is not a niche solution to the energy storage problem but a global one. This is the reason I think it will prevail.
ACKNOWLEDGMENTS
I wish to thank P. Larochelle, N. Cizek, J. Kesseli, T. Wolf, R. Apte, and T. Geballe for helpful discussions.