Maximizing Efficiency in Compressed Air Energy Storage: Insights from Thermal Energy Integration and Optimization

: Motivated by the suboptimal performances observed in existing compressed air energy storage (CAES) systems, this work focuses on the efficiency optimization of CAES through thermal energy storage (TES) integration. The research explores the dependence of CAES performance on power plant layout, charging time, discharging time, available power, and cavern volume. Hence, a range of solutions are examined, encompassing both solid and liquid TES options, alongside the potential utilization of external air heaters. Inefficiencies in solid TES due to significant retention of thermal power within the medium after complete discharge are identified and mitigated through optimization strategies. In addition, solutions to prevent ice formation at the low-pressure expander phase are suggested to avoid icing issues in CAES layouts with liquid TES. Through this comprehensive investigation, the study provides valuable insights into enhancing the efficiency and sustainability of CAES systems. By constructing a volume–power–time conversion table, the research contributes to the advancement of CAES technology, facilitating more efficient energy storage and utilization, thereby addressing critical challenges in the field of energy storage.


Introduction
In recent years, the escalating demand for increased energy production has encountered challenges stemming from environmental constraints, particularly in limiting greenhouse gas (GHG) emissions and addressing the rising costs and geo-political market uncertainty associated with fossil fuels.While large-scale renewable power plants present a promising avenue for mitigating GHG emissions, their deployment introduces complexities such as grid congestion and voltage stability issues due to the inherent unpredictability of renewable sources.Furthermore, their intrinsic intermittent nature may lead to unavailability during periods of high demand.
When the production of renewable energy systems (RESs) varies, grid voltage and frequency are maintained within tolerance by balancing power.Regulation can be achieved through existing spinning reservoirs, voltage regulation, and load following, incurring negligible additional system costs up to a market penetration of 20% [1].Exceeding the 20% causes the incremental costs per kWh to become significant.The costs arise due to the necessity for supplemental generating capacity capable of compensating for power drops [2].Compressed air energy storage (CAES) systems emerge as a viable solution to attain the target generating capacity.
The fluctuations in generation patterns in wind parks create complexities in electrical grid management, requiring technological solutions to balance supply and demand.CAES systems provide dispatchable power to align with load and market prices, enhancing the capacity factor of a renewable power plant, and contributing to reduce infrastructure transmission costs by mitigating peak energy delivered by wind parks [3].
Energies 2024, 17, 1552 3 of 16 Through a two-step simulation based on an analytical model of air temperature and pressure and a finite element model for the cavern, Khaledi et al. [15] investigated the thermomechanical behavior of rock salt in caverns.The authors clarified that the cavern's internal pressure governs its stability, while high-temperature operations affect its longterm serviceability.The conclusion is that setting the minimum pressure in the cavern would reduce the mechanical stresses below the dilatancy threshold while maintaining the air temperature below a threshold would avoid the creep deformation of rock salt and material weakening.
Other works based on A-CAES modeling are briefly outlined.Zhou et al. [16] developed a model to demonstrate variable conditions and the feasibility of an A-CAES system.Tola et al. [17] modeled two A-CAES plant configurations using air and thermal oil.In the oil plant, structural issues arose due to increasing tank pressure, highlighting the importance of tank design and maintenance.Sciacovelli et al. [18] presented a dynamic model for an A-CAES plant with packed bed Thermal Energy Storage (TES), covering full and off-design performance, and highlighting the transient characteristics of thermal storage, cavern dynamics, and compression and expansion phases.
The previous literature underscores different aspects of CAES systems, encompassing the modeling of distinct plant configurations and analyzing thermomechanical characteristics.However, while individual studies have explored specific aspects of CAES systems, a comprehensive assessment of efficiency across different plant configurations is absent.
Therefore, an exhaustive model of CAES power facilities is proposed to determine the performance across different layouts, aiming to investigate key aspects such as efficiency, produced and absorbed energy, charging, and discharging time across a wide range of cavern volumes (from 500 to 200,000 m 3 ), and aiming to consider different heat exchange solutions, including solid or liquid TES, with or without external air heating, and with constant or variable cavern pressure.Moreover, the study focuses on developing algorithms for solid TES to enhance heat exchange during the charging and discharging phases.The simulation outputs allow for the building of a volume-power-time (VPT) conversion table, facilitating the correlation between energy production, efficiencies, and cavern volumes.Overall, the findings provide insights into optimizing CAES power plants for more efficient and sustainable energy storage, thereby facilitating the integration of renewable energy sources into the grid.
A constant-volume CAES is assessed, where a cavern of a specified capacity serves as a reservoir for supplying air to the expander phases, allowing for variable pressure operations.In this sense, the present model diverges from our earlier work [19], which emphasized the theoretically attainable efficiency in constant-pressure underwater CAES (UW-CAES) power plants, achieved through iterative adjustments of the heat exchanger effectiveness while maintaining a constant airflow rate between compression and expansion.
The numerical simulations are performed with the model presented in Section 2 developed in a MATLAB environment, based on several CAES plant layouts, such as the following: 1.
Solid TES with external air heating 3.
External air heating (no air cooling, no TES) 4.
High-temperature liquid TES with two aftercoolers 7.
Low-temperature liquid TES with two aftercoolers 8.
Liquid TES with external air heating 9.
Constant-pressure cavern and air cooling (no external air heating, no TES) Each CAES plant layout is evaluated across various cavern volumes within the range specified in Table 1 to establish the correlation between cavern volume and energy production.The numerical model is illustrated in Section 2, while the outputs of the simulations and conclusions are discussed in Sections 3 and 4, respectively.

System Description
Figure 1 presents the A-CAES plant layout of the model and includes single or multiple compressor phases, an underground cavity, air heaters, and expander phases.The fixedvolume cavern is modeled at a prescribed pressure range compliant with the lower limit imposed by the minimum working pressure of the turbomachinery plus the pressure drops due to the throttling valve and circuit.Intercoolers are added to reduce air temperature between two consecutive compression phases, while aftercoolers allow decreasing air temperature before the injection into the cavern downstream of the compressor.Solid and liquid TES represent optional components modeled to compare the power plant efficiency.The former is modeled with concrete cylinders supplied by tubes of a given diameter [20], ensuring high heat capacity allowed by its high thermal conductivity.The latter operates at higher temperatures (600-900 K) due to the utilization of nitrate salt, implying the installation of hot and cold tanks and a freezing point control to avoid the liquid salt freezing in the pipeline when the plant is not operating.Two external air heaters (heater 1 and heater 2) are modeled as combustion chambers or electric resistances driven by a renewable power source.The simulations are executed for the cavern volumes indicated in Table 1, with the compressor and expander's mass flow rates related to each cavern volume.The main parameters of the simulations are specified instead in Table 2.For a given cavern volume, after the initialization of the cavern pressure, the air mass stored in the cavern is determined in time based on the density of the injected air.For an accurate prediction of the cavern-wall temperature, a finite element assessment of the TES is executed to determine the hourly thermal power stored.Hence, the discharging time and extracted power are The existing CAES power plants (McIntosh [21], Huntorf [22,23]) are improved with a model implementing TES systems capable of storing the heat extracted from the hot Energies 2024, 17, 1552 5 of 16 compressed air during each compression phase.The stored heat reheats the cold air from the cavern, eliminating fossil fuel burners or external air heaters.The compressor module is modeled with two phases operating as inter-cooling and after-cooling systems, i.e., lowpressure TES (LP TES) and high-pressure TES (HP TES), respectively.The expander is modeled with two modules (LP exp and HP exp), with both TESs operating as preheaters for each turbine phase.When the cavern pressure reduces from the full-load operation, the plant configuration is modified.The HP expansion phase is bypassed so that HP TES and LP TES are supplied in series.While LP TES preheats the air after the LP compression, the HP TES increases the air temperature to the maximum cycle temperature.An expander high-pressure bypass is inserted to ensure the plant's operation at the low-pressure ratio.
The power delivered bypassing the HP phase is lower than the full-load operation.Therefore, the power plant can operate in full-load mode (operating HP and LP expansion phases) and partial-load mode (operating LP expansion phase).The throttling valves adapt the sliding cavern pressure to the compressor and expander's rated working pressures.
The simulations are executed for the cavern volumes indicated in Table 1, with the compressor and expander's mass flow rates related to each cavern volume.The main parameters of the simulations are specified instead in Table 2.For a given cavern volume, after the initialization of the cavern pressure, the air mass stored in the cavern is determined in time based on the density of the injected air.For an accurate prediction of the cavern-wall temperature, a finite element assessment of the TES is executed to determine the hourly thermal power stored.Hence, the discharging time and extracted power are determined, and so is the round-trip efficiency of the plant.The simulations are performed in full-load operation with the fluid supplied to both HP and LP turbine phases and in partial-load operation with the air bypassing the HP phase at a low cavern pressure.In a solid TES layout (without external air heating), the charging mass flow rate is set to 250 kg/s, ensuring 8-9 h of charging time for the maximum cavern volume.The discharging mass flow rates are computed by imposing the cavern discharging mass flow to a charging mass flow ratio equal to 2.2.The mass flow rate is set to 550 kg/s for the maximum cavern volume, allowing an expansion phase duration similar to that of the Huntorf power plant [22,23], i.e., 3-4 h.The other mass flow rates are scaled linearly with cavern volume for the reference volume of 200,000 m 3 .According to [27], if the mass flow rate ratio increases above 2.2 in solid TES, the air mass from the cavern would be withdrawn too quickly, implying a drop in the efficiency in the expansion phase and non-homogenous heat distribution in the concrete sections.
In a solid TES layout with air heating, the algorithm is the same but includes two external heaters to increase the air temperature from each corresponding TES, one per turbine phase.
The solution with external air heating without cooling (no TES), is used to supply hot compressed air directly into the cavern.With this layout, both air heaters are used.Conversely, the external air heating layout with air cooling (no TES) implies a drop in the air temperature using intercoolers and aftercoolers.
In a liquid TES layout (without external air heating), the heat exchanger heats the air from the cavern, which is expanded in the turbine phases.The variation in the heat exchanger size just slightly influences the round-trip efficiency.The chosen heat exchanger data (Table 2) allow the highest gain efficiency, however moderate.If two aftercoolers are introduced (case a), the algorithm computes the optimum liquid mass flow rate, maximizing its temperature.Alternatively, if the hot liquid temperature is reduced (case b) to 778 K to raise the hot liquid stored mass, the HP exchanger mass flow rate increases from 225 kg/s to 274 kg/s.
In a liquid TES layout with external air heating, an external heater is presented upstream of each expansion phase.Two heat exchangers are modeled upstream of the HP turbine phase because the outlet temperature from the HP turbine phase can be higher than the hot liquid temperature.
In the last layout, with air cooling and external air heating, the cavern pressure is maintained constant at 45 bar during charging and discharging.With intercooling and aftercooling, the air temperature diminished to 330 K to reduce the vapor content.
Throttling valves are required to allow the air pressure to match the compressor-and expander-rated pressure.

Compressor and Turbine Equations
The model incorporates temperature-dependent specific heat, viscosity, conductivity, and air density to ensure precise outputs.The compressor and expander's performances are estimated by assuming polytropic transformations given the compression and expansion ratios and the main cycle temperatures.The polytropic transformation is modeled according to (1), considering a polytropic index m and specific heat c.The index m serves to model an irreversible process through reversible and quasi-static heat transfer and energy transformation.The Real Gas Law is introduced in Equation (2), with RA being the gas constant.
The final compression temperature for a polytropic process is given in Equation ( 3), achieved by combining ( 1) and (2).
For an adiabatic compressive process, the polytropic efficiency η pol,c is defined as The exponent ν is achieved by solving Equation ( 4), and it is substituted in Equation ( 3), yielding the final compression temperature as a function of initial air temperature T in , compression ratio β c = p out /p in , polytropic efficiency η pol,c , and the isentropic exponent ε c .
Hence, the required enthalpy variation for a polytropic compression is given in Equation ( 6), with specific heat c p defined as the average value of those computed at the initial and final point, expressed in J/kgK.
Substituting Equation ( 3) into (6) yields the specific work required to perform a compression from p in to p out with a temperature gain equal to (T c,out − T c,in ) in Equation (7).
Concerning the turbine, the polytropic efficiency is defined as Equation ( 9) describes the air temperature at the end of the expansion phase as a function of the initial expansion temperature T t,in , the expander polytropic efficiency, and the isentropic exponent ε T .
The polytropic efficiency is calculated with the size parameter method [28] recalled in Equation (10), where V in denotes the inlet volumetric flow rate in m 3 /s, and ∆h is denotes the isentropic enthalpy drop in J/kg processed by the turbomachine.The size parameter SP correlates statistically η pol to the characteristic diameter in axial compressors of gas turbines near the optimum characteristic speed N S and represents the denominator of the characteristic diameter, i.e., D s = D/SP, where D is the mean diameter of the turbomachine.Polytropic efficiencies higher than 90% are achieved in gas turbines from 100 MW (β t = 15 ÷ 18) or aero-derivative gas turbines with β t ≈ 30 [28].SP = V in 0.5 /∆h is 0.25 (10) Equations ( 11) and ( 12) describe the polytropic efficiency for the compressor phase as a function of the SP, and the polytropic efficiency for the expander, respectively.The chosen parameters for cooled and uncooled phases are η pol,∞ = 0.89 and η pol,∞ = 0.925.

Solid TES
The solid thermal energy storage is modeled by solving the unsteady 1-D heat propagation equation through the radial direction of a solid storage media cylinder.The storage is designed to keep the air speed in the pipeline at 8 m/s in the most demanding operating condition, with oversized tubes to allow high performance during compression and discharge (Table 2).The thermal power exchanged between two neighboring elementary volumes is computed based on Equation (13), where the index i denotes the i-th elementary volume of solid media considered; λ concrete denotes the thermal conductivity of the solid elements in W/mK; A and r step denote the area of the solid storage cylinder in m 2 and the elementary step in the radial direction in m, respectively; and (T i − T i+1 ) denotes the discrete temperature variation in K.The heat exchanged between the hot/cold air and the first elementary volume is computed by Equation ( 14), where .Q 1 denotes the thermal power transferred at the first radial element, (T air − T 1 ) denotes the temperature variation between the air and the first elementary solid media of area A 1 , and h pipes denotes the heat transfer coefficient in W/m 2 K. . .
The computations are executed with a time step t step , achieving the discrete heat transfer ∆Q between two elements (15).The updated temperature of the i-th element T i (16) is derived from the energy conservation, where the mass is denoted by m solid , and the specific heat is denoted by c p,solid .Hence, the whole temperature profile is computed by iterating through the time and radial direction of the solid media storage.
The fluid-pipe interaction is modeled according to the heat convection in pipes as shown by equations from (17) to (19), where λ pipes denotes the pipe thermal conductivity, D pipes denotes the pipe diameter, µ denotes the dynamic viscosity, and c p,air represents the air-specific heat.The dimensionless parameters Nu, Re, and Pr represent Nusselt, Reynolds, and Prandtl's numbers, respectively.If the air is cooled, the exponent n equals 0.3.Diversely, if the air is heated by a solid medium, the exponent n equals 0.4.For solid storage media, the specific heat and density are assumed constant.

Liquid TES
While solid TES systems imply variations in the outlet air temperature T air,out , liquid TES systems allow it to be kept constant.The liquid thermal energy storage medium is modeled according to the experimental data about molten salt collected from [24,29,30].The specific heat and density are constant, while the heat exchangers that cool or heat air are modeled according to the experimental results synthesized in [25] and reported in Table 1.For simplicity, the thermal inertia of the heat exchanger is neglected.
According to the schematic of staggered bundles in [27], X t = X l = 1.At equilibrium, the heat extracted from the air is equivalent to the heat absorbed by the molten salt.Hence, steady heat exchange is modeled such that .m air × c p,air × ∆T air = .m salt × c p,salt × ∆T salt , with air mass flow rates indicated in Table 1.A drawback of liquid TES is the minimum molten salt temperature T salt,min of 550 K due to its freezing point, imposing a higher outlet air temperature T air,out of 560 K.The molten salt constraints are defined in Table 2. Once the liquid mass flow rate .m salt that satisfies the constraints is calculated, the heat transfer coefficient based on the Nusselt number correlated with staggered bundles in crossflow is derived according to Equation (20). The If 500 ≤ Re b ≤ 10 3 , C n denotes a correlation factor for the row number equal to 1 in this model.P rb and P rw represent the Prandtl's number for heat exchanger bulk and wall.These are estimated by Equations ( 24) and (25), respectively.
A simplified method to calculate the heat transfer coefficient was presented by Raju et al. [31] in the case of a lack of information about the correlation between the flow and heat transfer coefficient.The cavern wall exchange coefficient was approximated by fitting a two-parameter function that represents the global heat exchange factors until the numerical temperature behavior is close to the experimental temperature of the Huntorf power plant.
The fitted function ( 26) is used to estimate the cavern temperature in the considered range (Table 2), assuming the thermal exchange coefficient is independent of the cavern volume.Parameters a and b [31] denote natural and forced convection (air injection or withdrawal) factors.These are equal to 0.2356 kJ/m 3 and 0.0149 kJ/kgm 3 , respectively.

Pipeline Pressure Drop
The pipe head losses are computed according to Equation (27).Velocity V i and density ρ are determined at each section analyzed.Since the mass flow rate is constant, density variations of the airflow due to the continuous pressure drop generate a gain in the flow speed while the pipeline is assumed adiabatic.

Results
The first output of the simulation is represented by a finite element analysis of the solid TES as a function of the pipe length and radius.Figure 2 shows the temperature profile in the solid medium (concrete) after 8 h of the charging phase; meanwhile, Figure 3 shows the solid medium temperature profile after 3 h of the discharging phase, highlighting the remaining heat flux in the solid medium.This problem is due to the temperature difference between the charging and discharging phases and the solid TES piping configuration during partial-load operation (discharging phase).Since the two TESs are connected in series, the HP module experiences higher inlet air temperature with lower thermal exchange rates.
The second output of the simulation is achieved from the CAES configuration with liquid TES and without external air heating systems.The liquid medium (nitrite salt) has a freezing temperature close to 538 K [24], implying an inferior threshold for air temperature at 550 K due to freezing problems in the pipeline and heat exchangers.In a second postrefrigeration system, the storage medium (mineral oil) is added to further reduce the air temperature.This power plant layout experiences an icing problem in the last LP expander phase because the hot liquid stored during the compression phase cannot reheat the air to cover the entire expansion phase.
Energies 2024, 17, x FOR PEER REVIEW 10 of 17 heat transfer coefficient.The cavern wall exchange coefficient was approximated by fitting a two-parameter function that represents the global heat exchange factors until the numerical temperature behavior is close to the experimental temperature of the Huntorf power plant.The fitted function ( 26) is used to estimate the cavern temperature in the considered range (Table 2), assuming the thermal exchange coefficient is independent of the cavern volume.Parameters a and b [31] denote natural and forced convection (air injection or withdrawal) factors.These are equal to 0.2356 kJ/m 3 and 0.0149 kJ/kgm 3 , respectively.

Pipeline Pressure Drop
The pipe head losses are computed according to Equation (27).Velocity Vi and density ρ are determined at each section analyzed.Since the mass flow rate is constant, density variations of the airflow due to the continuous pressure drop generate a gain in the flow speed while the pipeline is assumed adiabatic.

Results
The first output of the simulation is represented by a finite element analysis of the solid TES as a function of the pipe length and radius.Figure 2 shows the temperature profile in the solid medium (concrete) after 8 h of the charging phase; meanwhile, Figure 3 shows the solid medium temperature profile after 3 h of the discharging phase, highlighting the remaining heat flux in the solid medium.This problem is due to the temperature difference between the charging and discharging phases and the solid TES piping configuration during partial-load operation (discharging phase).Since the two TESs are connected in series, the HP module experiences higher inlet air temperature with lower thermal exchange rates.The second output of the simulation is achieved from the CAES configuration with liquid TES and without external air heating systems.The liquid medium (nitrite salt) has a freezing temperature close to 538 K [24], implying an inferior threshold for air temperature at 550 K due to freezing problems in the pipeline and heat exchangers.In a second post-refrigeration system, the storage medium (mineral oil) is added to further reduce the air temperature.This power plant layout experiences an icing problem in the last LP expander phase because the hot liquid stored during the compression phase cannot reheat the air to cover the entire expansion phase.
As depicted in Figure 4, the hot liquid runs out after 4.2 h of expansion.A viable solution consists of decreasing the maximum hot liquid temperature to 738 K to store a higher mass of hot salt during compression, as illustrated in Figure 5.As depicted in Figure 4, the hot liquid runs out after 4.2 h of expansion.A viable solution consists of decreasing the maximum hot liquid temperature to 738 K to store a higher mass of hot salt during compression, as illustrated in Figure 5.
A generic air heating system to preheat the air, installed between the air preheater and the expander, represents the second viable solution.With this layout, the air temperature is under control.
The following results synthesize the simulations performed for each cavern volume indicated in Table 1.These are collected in the volume-power-time conversion tables (VPT), which are based on the following assumptions: the produced energy is almost independent of the expansion time; the power is constant during a generic phase (charging, discharging, or air heating); and for each phase duration, a gas turbine can change its mass flow rate according to the new phase duration.The first two assumptions are invalid for power plant layouts with solid TES without an external air heating system.Indeed, concrete does not provide constant cooling/heating air temperature.Hence, the absorbed/delivered power changes during the phases, so no VPT tables exist for this layout, as illustrated in Figure 6.The second output of the simulation is achieved from the CAES configuration with liquid TES and without external air heating systems.The liquid medium (nitrite salt) has a freezing temperature close to 538 K [24], implying an inferior threshold for air temperature at 550 K due to freezing problems in the pipeline and heat exchangers.In a second post-refrigeration system, the storage medium (mineral oil) is added to further reduce the air temperature.This power plant layout experiences an icing problem in the last LP expander phase because the hot liquid stored during the compression phase cannot reheat the air to cover the entire expansion phase.
As depicted in Figure 4, the hot liquid runs out after 4.2 h of expansion.A viable solution consists of decreasing the maximum hot liquid temperature to 738 K to store a higher mass of hot salt during compression, as illustrated in Figure 5.A generic air heating system to preheat the air, installed between the air preheater and the expander, represents the second viable solution.With this layout, the air temperature is under control.
The following results synthesize the simulations performed for each cavern volume indicated in Table 1.These are collected in the volume-power-time conversion tables (VPT), which are based on the following assumptions: the produced energy is almost independent of the expansion time; the power is constant during a generic phase (charging, discharging, or air heating); and for each phase duration, a gas turbine can change its mass flow rate according to the new phase duration.The first two assumptions are invalid for power plant layouts with solid TES without an external air heating system.Indeed, concrete does not provide constant cooling/heating air temperature.Hence, the absorbed/delivered power changes during the phases, so no VPT tables exist for this layout, as illustrated in Figure 6.
Figure 6 shows the HP and LP expander's power during the discharging phase and the HP phase bypass occurring 1.75 h after the start of the discharge phase.The chart also indicates the power reduction due to the progressive cooling of the solid TES.A generic air heating system to preheat the air, installed between the air preheater and the expander, represents the second viable solution.With this layout, the air temperature is under control.
The following results synthesize the simulations performed for each cavern volume indicated in Table 1.These are collected in the volume-power-time conversion tables (VPT), which are based on the following assumptions: the produced energy is almost independent of the expansion time; the power is constant during a generic phase (charging, discharging, or air heating); and for each phase duration, a gas turbine can change its mass flow rate according to the new phase duration.The first two assumptions are invalid for power plant layouts with solid TES without an external air heating system.Indeed, concrete does not provide constant cooling/heating air temperature.Hence, the absorbed/delivered power changes during the phases, so no VPT tables exist for this layout, as illustrated in Figure 6.
Figure 6 shows the HP and LP expander's power during the discharging phase and the HP phase bypass occurring 1.75 h after the start of the discharge phase.The chart also indicates the power reduction due to the progressive cooling of the solid TES.  Figure 6 shows the HP and LP expander's power during the discharging phase and the HP phase bypass occurring 1.75 h after the start of the discharge phase.The chart also indicates the power reduction due to the progressive cooling of the solid TES.

represents a VPT conversion table showing the available power during a full-load discharging phase (HP and LP expansion phases operating). A VPT is used by setting a cavern volume (upper row), choosing the duration of the full-load discharge phase
[hours] on the first column, and thus finding the corresponding CAES power for that phase.
The first two rows in Table 3 in bold indicate the simulated duration of the full-load expansion for each cavern volume.The other rows are scaled with the phase duration from the simulated data.These tables are also available for the compression phase, partial-load discharge, and external air heating.
Figure 7 collects the round-trip efficiencies as a function of the cavern volume, defined as the turbine work to compressor work ratio (28).The electrical works of compression and expansion (W t,el and W c,el ) are calculated taking into account the electromechanical losses of the motor/generator (η mot-gen in Table 2).The configuration without TES and air cooling but with external air heating demonstrates higher round-trip efficiency than its counterpart with air cooling.Indeed, the utilization of the hot air stored within the cavern is more efficient.However, it results in a lower mass of air stored, leading to a drop in electric work output, as illustrated in Figure 8.
High electric work output is achieved in configurations maintaining constant pressure and incorporating external air heating (triangle markers) and in configurations with air cooling and external air heating (circle markers).The two existing CAES plants are an example, wherein air cooled to approximately 50 °C, despite exhibiting low efficiency, yields increased energy production.Conversely, integrating a liquid TES to cool the air to 277 °C substantially enhances efficiency but decreases the stored air mass and the energy output.
In the partial-load operation of the liquid TES layout without external air heating, the temperature of the air stored in the cavern cannot increase above a threshold due to freezing issues in the heat exchanger and pipeline.If the plant includes two aftercoolers (case a), the air mass stored in the cavern increases due to the lower temperature of the supplied air.A higher cavern mass induces more compression work to achieve the design cavern pressure and a longer expansion phase.By reducing the hot liquid temperature in the same plant (case b), a higher mass flow rate of air is available for the expansion phase, even though the stored air mass is lower and the expansion is shorter, implying a lower energy extraction.
Efficiency depends slightly on volume, as shown by the constant pressure layout.Notably, layouts avoiding air cooling in favor of external air heating demonstrate higher efficiency than their counterparts without TES, leveraging the hot air stored within the The increasing volume within the cavern corresponds to enhanced turbine work and rated power because the higher the volume, the higher the mass stored in the cavern and the higher the mass flow rate processed by the turbine.
Among TES options, solid TES systems demonstrate higher performance than liquid TES and other options.A supplementary external air heating implies a rise in the expander energy produced.
The efficiency outputs of configurations without external heating exhibit a higher dependency on cavern volume, wherein larger volumes allow higher-rated power for expansion.Generally, these plants allow higher efficiency but lower electric work as they rely solely on heat stored during compression.Conversely, configurations incorporating external heating systems show reduced dependency on cavern volumes and lower efficiency, as they necessitate an additional power source.
In layouts with external air heating and air cooling, the power achieved at the expander is higher than the solution without cooling due to the higher stored air mass in the cavern.The increment in the air mass is due to the drop in cavern temperature induced by cooling.
The configuration without TES and air cooling but with external air heating demonstrates higher round-trip efficiency than its counterpart with air cooling.Indeed, the utilization of the hot air stored within the cavern is more efficient.However, it results in a lower mass of air stored, leading to a drop in electric work output, as illustrated in Figure 8.

Conclusions
This work presents a comprehensive model of a compressed air energy storage (CAES) system aimed at analyzing key performance parameters across a wide range of cavern volumes (from 500 to 200,000 m 3 ) and various heat exchange solutions, including solid or liquid thermal energy storage (TES), with or without external air heating, and constant or variable cavern pressure.
An optimization algorithm for solid TES parameters is developed to enhance heat exchange efficiency during the charging and discharging operations.However, despite the wide temperature range facilitated by solid TES, the air temperature fluctuations lead to power variability during the charging and discharging phases, exacerbated by the low thermal conductivity of the tested material (reinforced concrete), resulting in suboptimal heat extraction during the discharging phase.
The fragmentation of solid storage into multiple elements is recommended, with excess heat potentially distributed for other purposes, such as industrial or residential heating.Several CAES power plant layouts are evaluated to estimate their performances, including the energy-cavern volume relationship, efficiency, power, and energy.Simulation outputs enable the construction of volume-power-time (VPT) conversion tables, facilitating correlation between power plant parameters and cavern volume, charging, and discharging phase durations.
The liquid TES model highlights issues, such as hot liquid (molten salts) depletion during the discharging phase, potentially leading to icing problems in power plant layouts lacking external air heaters.Solutions include reducing the hot liquid temperature and installing additional aftercoolers with mineral oil to maintain stable air temperature, albeit at the expense of more components and a limited aftercooler outlet air temperature.
The power plant efficiency is influenced by the plant's layout, with simulations indicating a correlation between cavern volume, compressor-rated power, and cavern temperature.Higher cavern volumes result in higher compressor-rated power and efficiency, while cavern temperature affects efficiency by saving energy required to heat the air stored in the cavern.However, higher cavern temperatures also reduce the air mass storage capacity, ultimately impacting energy production.High electric work output is achieved in configurations maintaining constant pressure and incorporating external air heating (triangle markers) and in configurations with air cooling and external air heating (circle markers).The two existing CAES plants are an example, wherein air cooled to approximately 50 • C, despite exhibiting low efficiency, yields increased energy production.Conversely, integrating a liquid TES to cool the air to 277 • C substantially enhances efficiency but decreases the stored air mass and the energy output.
In the partial-load operation of the liquid TES layout without external air heating, the temperature of the air stored in the cavern cannot increase above a threshold due to freezing issues in the heat exchanger and pipeline.If the plant includes two aftercoolers (case a), the air mass stored in the cavern increases due to the lower temperature of the supplied air.A higher cavern mass induces more compression work to achieve the design cavern pressure and a longer expansion phase.By reducing the hot liquid temperature in the same plant (case b), a higher mass flow rate of air is available for the expansion phase, even though the stored air mass is lower and the expansion is shorter, implying a lower energy extraction.
Efficiency depends slightly on volume, as shown by the constant pressure layout.Notably, layouts avoiding air cooling in favor of external air heating demonstrate higher efficiency than their counterparts without TES, leveraging the hot air stored within the cavern for efficient operation, albeit with a reduced stored mass and a subsequently lower electric work output.

Conclusions
This work presents a comprehensive model of a compressed air energy storage (CAES) system aimed at analyzing key performance parameters across a wide range of cavern volumes (from 500 to 200,000 m 3 ) and various heat exchange solutions, including solid or liquid thermal energy storage (TES), with or without external air heating, and constant or variable cavern pressure.
An optimization algorithm for solid TES parameters is developed to enhance heat exchange efficiency during the charging and discharging operations.However, despite the wide temperature range facilitated by solid TES, the air temperature fluctuations lead to power variability during the charging and discharging phases, exacerbated by the low thermal conductivity of the tested (reinforced resulting in suboptimal heat extraction during the discharging phase. The fragmentation of solid storage into multiple elements is recommended, with excess heat potentially distributed for other purposes, such as industrial or residential heating.Several CAES power plant layouts are evaluated to estimate their performances, including the energy-cavern volume relationship, efficiency, power, and energy.Simulation outputs enable the construction of volume-power-time (VPT) conversion tables, facilitating correlation between power plant parameters and cavern volume, charging, and discharging phase durations.
The liquid TES model highlights issues, such as hot liquid (molten salts) depletion during the discharging phase, potentially leading to icing problems in power plant layouts lacking external air heaters.Solutions include reducing the hot liquid temperature and installing additional aftercoolers with mineral oil to maintain stable air temperature, albeit at the expense of more components and a limited aftercooler outlet air temperature.
The power plant efficiency is influenced by the plant's layout, with simulations indicating a correlation between cavern volume, compressor-rated power, and cavern temperature.Higher cavern volumes result in higher compressor-rated power and efficiency, while cavern temperature affects efficiency by saving energy required to heat the air stored in the cavern.However, higher cavern temperatures also reduce the air mass storage capacity, ultimately impacting energy production.
Layouts with external air heating generally yield higher electric work but lower efficiency.In contrast, layouts without external air heaters exhibit higher efficiency due to heat stored during the compressive phase.
Volume-power-time conversion tables and energy-volumes and efficiency-volumes maps are developed to establish relationships between energy, efficiencies, and cavern volumes, facilitating parameter setting for CAES power plants.These findings contribute to advancing CAES technology for more efficient energy storage and utilization, aiding the integration of renewable energy sources into the grid.

Energies 2024 , 17 Figure 1 .
Figure 1.General CAES plant layout with TES and external air heating in full-load operation (lefthand side) and partial-load (right-hand side).Solid line: active expander phases.Thin solid line: inactive expander phases.Dashed line: compressor phases.

Figure 1 .
Figure 1.General CAES plant layout with TES and external air heating in full-load operation (lefthand side) and partial-load (right-hand side).Solid line: active expander phases.Thin solid line: inactive expander phases.Dashed line: compressor phases.

Figure 2 .
Figure 2. Solid TES temperature profile after 8 h of charge, according to [27].Figure 2. Solid TES temperature profile after 8 h of charge, according to [27].

Figure 2 .
Figure 2. Solid TES temperature profile after 8 h of charge, according to [27].Figure 2. Solid TES temperature profile after 8 h of charge, according to [27].

Figure 5 .
Figure 5. LP turbine temperature during expansion phase with liquid TES with a lower hot liquid temperature.Inlet: solid line, outlet: dashed line.

Figure 6 .
Figure 6.Power delivered by the CAES system with solid TES and without external air heating.HP expander phase: solid line, LP expander phase: dashed line.

Figure 5 .
Figure 5. LP turbine temperature during expansion phase with liquid TES with a lower hot liquid temperature.Inlet: solid line, outlet: dashed line.

Figure 5 .
Figure 5. LP turbine temperature during expansion phase with liquid TES with a lower hot liquid temperature.Inlet: solid line, outlet: dashed line.

Figure 6 .
Figure 6.Power delivered by the CAES system with solid TES and without external air heating.HP expander phase: solid line, LP expander phase: dashed line.

Figure 6 .
Figure 6.Power delivered by the CAES system with solid TES and without external air heating.HP expander phase: solid line, LP expander phase: dashed line.

Figure 8 .
Figure 8. Produced energy chart of the analyzed CAES plant layouts, according to [27].

Figure 8 .
Figure 8. Produced energy chart of the analyzed CAES plant layouts, according to [27].

Table 1 .
Cavern volume range and mass flow rates.