Next Article in Journal
Aerodynamic Characteristics of a Single Airfoil for Vertical Axis Wind Turbine Blades and Performance Prediction of Wind Turbines
Next Article in Special Issue
Train Braking Time Variations Changing the Pressurized Air Temperature
Previous Article in Journal
Towards Reconstruction of Complex Flow Fields Using Unit Flows
Previous Article in Special Issue
Thermo-Environmental Performance of Four Different Shapes of Solar Greenhouse Dryer with Free Convection Operating Principle and No Load on Product
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TES-PD: A Fast and Reliable Numerical Model to Predict the Performance of Thermal Reservoir for Electricity Energy Storage Units

Department of Industrial Engineering, University of Padova, Via Venezia 1, 35131 Padova, Italy
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(7), 256; https://doi.org/10.3390/fluids6070256
Submission received: 26 April 2021 / Revised: 8 July 2021 / Accepted: 9 July 2021 / Published: 13 July 2021
(This article belongs to the Special Issue Advances in Thermo-Fluid Dynamics of Industrial Systems)

Abstract

:
The spread of renewable resources, such as wind and solar, is one of the main drivers to move from a fossil-based to a renewable-based power generation system. However, wind and solar production are difficult to predict; hence, to avoid a mismatch between electricity supply and demand, there is a need for energy storage units. To this end, new storage concepts have been proposed, and one of the most promising is to store electricity in the form of heat in a Thermal Energy Storage reservoir. However, in Thermal Energy Storage based systems, the critical component is the storage tank and, in particular, its mathematical model as this plays a crucial role in the storage unit performance estimation. Although the literature presents three modelling approaches, each of them differs in the considered parameters and in the method of modelling the fluid and the solid properties. Therefore, there is a need to clarify the model differences and the parameter influences on plant performance as well as to develop a more complete model. For this purpose, the present work first aim is to compare the models available in the literature to identify their strengths and weaknesses. Then, considering that the models’ comparison showed the importance of adopting temperature-dependent fluid and storage material properties to better predict the system performance, the authors developed a new and more detailed model, named TES-PD, which works with time and space variable fluid and solid properties. In addition, the authors included the tank heat losses and the solid effective thermal conductivity to improve the model accuracy. Based on the comparisons between the TES-PD model and the ones available in the literature, the proposal can better predict the first cycle charging time, as it avoids a 4% underestimation. This model also avoids overestimation of the delivery time, delivered energy, mean generated power and plant round-trip efficiency. Therefore, the results underline that a differential and time-accurate model, like the TES-PD, even if one-dimensional, allows a fast and effective prediction of the performance of both the tank and the storage plant. This is essential information for the preliminary design of innovative large-scale storage units operating with thermal storage.

1. Introduction

The ever-increasing energy demand and the need to replace approximately 40% of the currently in-operation power plants due to ageing means that, by 2040, 7200 GW of new capacity needs to be installed [1]. Since the power industry is responsible for approximately 25% of the global CO2 emissions [2], and considering that there is the need of knocking down the energy sector greenhouse gas emissions, the vast majority of the new capacity will be, on the one hand, of renewable origin (mainly wind and solar photovoltaic units) and, on the other hand, from low emissions fossil-based thermal units (e.g., combined-cycle units) and nuclear power plants.
However, despite the technology improvements and researchers’ efforts, nuclear and combined-cycle power plants lack in flexibility (see, e.g., [3,4,5,6,7,8,9]); therefore, the future power production systems will be made by Variable Renewable Energy Sources (VRESs) units and low flexibility thermal power plants, a generation mix that sets out the installation of Large-Scale Energy Storage Systems. Only spreading these large-scale energy storage installations can help balance the production of VRESs without encountering strong power fluctuations that can result in management and control problems, devices faults and, in the worst scenario, in local or even global blackout.
The International Energy Agency estimates that the United States, Europe, China and India require at least 310 GW of new storage capacity to guarantee the grid stability with a high VRES penetration [10]. In addition, the same international institute also underlined that available storage technologies, such as Pumped Hydro Storage (PHS), Compressed Air Energy Storage (CAES) and batteries, would not be able to cope with the need for storage capacity due to geographical constraints (PHS and CAES) or reduced lifespan (batteries) [11].
Therefore, there is the pressing need to make available storage units able to (i) boost VRES installations, (ii) be installed near both VRES and low flexibility thermoelectric units to reduce transmission losses and network congestion, (iii) ensure a grid safe operation and (iv) do not increase the grid overcapacity. The latter is a fundamental point because, at the time of writing, fossil-based thermal plants are used as backup units to compensate VRES power fluctuations.
Indeed, fossil capacity is not directly replaced by renewable plants, resulting in an installed overcapacity corresponding to a rapid decline of fossil units operating hours [12]. Moreover, additional overcapacity and a further reduction of fossil-fuelled plants working hours encounter in the case of extensive exploitation of large-scale energy storage. Thus, today’s better option to overcome these issues is to look beyond the traditional notion of storage systems and switch to the “Integrated Energy Storage System” (I-ESS) concept: storage units embeddable into existing fossil thermal power plants [13,14].
In a nutshell, devices, such as gas turbines, electric generators, step-up transformers and transmission lines, which are commonly mounted on existing thermal power units, can be used to build up the I-ESS plant, ensuring that the storage unit (i) does not add overcapacity, (ii) helps in the revamping of underutilised units, (iii) reduces soil consumption, (iv) boosts VRES spread and (v) guarantees safe network operation even with significant mismatches between electricity supply and demand.
For this purpose, the authors developed an I-ESS plant able to be integrated into underutilised thermal power plants and built near the VRES facilities or in an ad hoc site [13,14]. The plant stores electricity as sensible heat in a high-temperature artificial tank consisting of a packed bed and usually called “Thermal Energy Storage” (TES). The deployment of sensible TES guarantees greater cost-effectiveness compared to latent and thermochemical heat storage systems, as demonstrated by several works published in the literature (see, e.g., [15,16,17,18]).
Compared to PHS, which is the most widespread technology with over 96% of the total installed capacity [19], and CAES, such technology does not suffer from geographical constraints and does not require stable water (like PHS) or natural gas (as CAES) flow rate. A long-cycle life also characterises it compared to batteries with a higher simplicity than Pumped Thermal Energy Storage (PTES) [20,21,22,23,24,25,26,27,28,29]. The I-ESS plant proposed by Benato and Stoppato [13,14] is an open-cycle adopting air as working fluid in both charging and delivering mode.
In a nutshell, the charging scheme is made up of a high-temperature tank, a fan, an electric heater, an electric motor and a heat exchanger, while the delivery unit consists of a gas turbine in which the high-temperature tank replaces the combustion chamber. Compared to the available PTES layouts, the adoption of a unique hot tank instead of two (one at high and one at low temperature) and an open-cycle, which allows the use of air instead of argon, reduces the devices’ number, the scheme complexity, and the system management.
Despite a large number of published works and patents concerning TES-based storage units, the vast majority of them aim at the development and modelling of innovative plant configurations for storing electricity as sensible heat, while only a few are focused on tank behaviour investigation. The tank arranged as a packed bed is the core of the TES-based storage system, and its mathematical modelling plays a crucial role in the plant performance prediction as well as in the computation efforts and simulation time. In addition, the non-stationary behaviour of such storage remains a central issue also in its management. In particular, when the TES-based energy storage unit is coupled with VRES plants, it is crucial to accurately predict the charging and discharging times, the stored energy, and the energy that can be delivered, based on the state of charge. These parameters depend not only on the plant layout but also mainly on the physical nature of the storage material and the fluid evolution in the storage itself. Therefore, the availability of accurate storage mathematical models remains a pivotal point for predicting plant performance and its real-world implementation.
For these reasons, the storage mathematical model needs to consider the essential variables required to describe its behaviour accurately without significantly increasing the system complexity and computational time. These needs arise from two facts: (i) the tank model is inserted into the storage unit; therefore, the entire model has to provide accurate results in a reasonable time, (ii) computational Fluid Dynamics (CFD) analysis of the sole tank can provide accurate data, deriving from the solution of the Navier–Stokes system of equations using some closure logic for turbulent phenomena.
Such a strategy would take charge of the local variability of many time-related events and phenomena. However, to be truly effective and provide a complete description of the system phenomena, CFD-based approaches require a three-dimensional tank model: a design condition that clashes with the engineering needs of a reliable, accurate and, to some extend, fast model usable during the design phase.
Note that, due to an accurate and fast simulation of the entire plant, the charging and delivery times as well as the stored and deliverable energy can be estimated in light of storage requirements and packed bed constituting materials. Subsequently, using the CFD techniques, the storage tank can be further investigated to explore the local variability of a large number of time-related events and phenomena.
For these reasons, the work initially focuses on comparing the literature most adopted one-dimensional mathematical models concerning the influence of the single parameters on the tank and storage system performance. According to the authors’ best knowledge, three models have been used in the literature. Apparently, all of them can accurately describe the time-varying storage tank behaviour, and they have not previously been compared.
Thus, as a first point of novelty, the different tank models are inserted in the energy storage system architecture presented by Benato and Stoppato [13,14] and compared to predict the performance of such systems. As the models’ structures are quite different, our investigation was conducted with the logic of increasing modelling complexity to ensure its effectiveness. At first, the behaviour of the storage tank was analysed through the purely algebraic description proposed by Howell et al. [30]. This modelling strategy, due to its ease of deployment, is one of the top choices in the literature, but it is still the most limited.
The second modelling approach is based on the proposal made by McTigue et al. [31]. The model feels a pseudo-time-varying approach, for which the non-stationary of the storage system is treated in terms of average quantities, i.e., the mean temperature and average mass-flow rates. A partial variability of both the operating fluid and the storage solid’s thermal coefficients, albeit with a time interval averaged logic, is taken into account. The last modelling approach is the one proposed by Desrues et al. [21]. The model describes the reservoirs’ transient behaviour through a purely one-dimensional time marching partial differential equations approach.
The comparison among the three models allows us to underline the importance of (i) space and time variability of the thermophysical coefficients related to the operating fluid and the storage material, (ii) the overall heat loss coefficient and (iii) the solid effective thermal conductivity. Thus, starting from the formulation proposed by Desrues et al. [21], the authors developed a new tank model (named TES-PD model), which includes all these features.
Despite its one-dimensionality but due to its high level of completeness and ability to take charge of space and time variability of the fluid and storage material thermophysical properties, the TES-PD tank model is of real help in designing and managing such devices by predicting data that can be otherwise unobtainable; the latter can then be used as preliminary designing guidelines or input for a subsequent CFD simulations. In addition, the completeness of the proposed model makes it possible to investigate storage systems under innovative operating conditions, such as the one at low temperatures, a research field still partially unexplored.
The work is organized as follows: in Section 2 the TES-PD model is discussed in terms of its analytical description and numerical treatment. Section 3 presents the validation procedure of the models available in the literature in order to use them to highlight the peculiarities of the TES-PD model. Section 4 describes the plant scheme used to test the proposal, while the investigation outcomes and the comparison between the literature-available models, and the TES-PD one are summarised and discussed in Section 5. Finally, Section 6 states the conclusions.

2. The TES-PD Numerical Model

In this section, the TES-PD tank mathematical model used to predict the storage behaviour is presented. The model rearranges the most used tank models and improves their predictive capabilities with additional features. Those models, in order of complexity, are the Mumma & Marvin model [30], the Schumann & McTigue et al. model [31] and the Desrueset al. model [21]. For the sake of conciseness, such numerical strategies are reported in Appendix A and, in the following, they are addressed as “model m.1”, “model m.2” and “model m.3”, respectively.
The TES-PD model relies on the formulation proposed by Desrues et al. [21]. Despite the common starting point, the authors improved it by the integration in the TES-PD model of the entire set of parameters that affect the behaviour of the storage tank. These parameters are variously missing in the available literature contributions; in particular,
  • The space and time variability of all the thermophysical properties of both fluid and storage material.
  • The overall heat loss coefficient.
  • The solid material effective thermal conductivity.
Model m.3 was selected as the basis for the TES-PD model since it represents the most advanced description of the storage tank already available in the literature. However, model m.3 envisages constant solid and fluid properties computed at an arbitrary reference temperature. This limitation often leads to mispredicting the system dynamics and introducing uncertainties related to the reference state. Moreover, it can not fit the usage of the most advanced storage materials.
If adequately modified, model m.3 allows updating all the thermophysical properties in both space and time. Thus, in the TES-PD model, the authors embedded the update of both the solid and the fluid properties with a procedure similar to the one described by McTigue et al. [31] for model m.2. However, in contrast to model m.2 for which no specifications were reported as far the computation of the thermophysical coefficients, solid material properties were acquired from the National Institute of Standards and Technology database [32], while the fluid thermophysical properties were taken from the CoolProp library [33]. In this manner, the thermophysical properties variability was considered layer-by-layer, updating their values during the time-marching procedure as a function of the actual temperature.
Moreover, concerning the overall heat loss coefficient, this is accounted for only by model m.1 (see Appendix A.1 for the details). However, such a parameter is of fundamental importance for a storage tank because the charging and discharging periods can last from hours/days, making this a non-trivial contribution concerning its performance. Introducing this hypothesis, the TES-PD model reads as follows:
ρ f t + ( ρ f v f ) x = 0
( ρ f T f ) t + ( ρ f T f v f ) x = 1 ε α ( T s T f )
( ρ s c p , s T s ) t + x k s , e f f T s x = c p , f 1 ε α ( T f T s ) U i C u ( 1 ε ) ( T s T a m b )
The above-mentioned set of equations consists of the fluid’s mass conservation equation, the fluid energy conservation equation and the solid energy conservation equation. The model is closed by a constitutive expression for the pressure drop and the ideal gas equation of state:
p x = C f β 1 2 ρ f v f 2
p = ρ f r T f
where t is the time and x is the axial tank coordinate with L its length; ρ f , v f , T f and p are the fluid density, velocity, temperature and thermodynamic pressure, respectively; while ρ s and T s are the storage material density and temperature, respectively. k f and k s , e f f denote the fluid and the solid effective thermal conductivity, and c p , f , and c p , s are the fluid and the solid specific heat coefficient at constant pressure, respectively. ε denotes the void fraction, T a m b is the ambient temperature, and r is the specific gas constant.
The α and β coefficients generalise model m.3, assuming different values depending on the storage material geometry. α equals to 6 d 2 N u k f c p , f 1 for spheres packing and equals to 4 ε d h 2 N u k f c p , f 1 for microchannels. d denotes the particle equivalent diameter, while d h denotes their hydraulic diameter. As far as the β coefficient is concerned, the latter equals ( 1 ε ) ( d ε ) 1 for spheres packing, while it equals 4 d h 1 for microchannels. C u ( 1 ε ) 1 denotes the external surface to volume ratio, while U i is the overall heat loss coefficient.
C u values are 4 L 1 and 2 π A 1 A π 1 when the considered geometry is a microchannel structure or sphere packing, respectively. Finally, the hydraulic diameter, d h , together with the Nusselt number, N u , and the friction coefficient, C f , are computed according to Desrues et al. [21] as given in Appendix A.3. In the present work, the authors considered packed sphere geometries where stated.
Compared to model m.3, Equation (1a,b) remain formally unaltered, while Equation (1c) accounts for the time-variability of the thermophysical coefficients and embeds the heat losses through the external surface. Such improvements are of fundamental importance in detailing the system’s dynamics since truly layer-by-layer variability of the solid and fluid properties can be accounted for as well as the heat exchange between the tank and the external ambient temperature.
To solve the previously-mentioned set of partial differential equations, the authors adopted a first-order explicit Euler scheme for the temporal integration, while a first-order generalised upwind finite difference method was employed to compute the first spatial derivatives. A second-order centred finite difference approach was used for second derivatives treatment [34,35]. The cell-bound variability of the physical coefficients was accounted for through the method proposed by De Vanna et al. [36].
The fluid velocity was updated through a segregated approach to overcome the equation coupling. This set of discretization schemes is guided by the simplicity and computational efficiency. While using first-order methods, this choice proved to produce accurate results, which fits the available references in corresponding conditions. Note that the literature does not suggest methods to solve the set of equations; hence, the authors developed the above-mentioned approach and integrated it into the TES-PD model in order to provide both an updated TES model and a simple and computationally efficient resolution scheme.
Concerning the initial conditions for the partial differential equations problem, the entire set of tests was performed by the authors while adopting uniform ambient temperature and density, and the flow was considered initially at rest. The fluid properties were set according to the ideal gas law, while the values for the fluid mass flow rate, pressure, and temperature were enforced at the inlet location depending on the plant cycle in which the tank was installed.
Zero-gradient conditions, with first-order of accuracy, were set for the solid’s temperature distribution at the tank’s inlet and outlet. The strategy allows modelling the storage material as adiabatic, packing the heat losses in the dispersion coefficient only.
The modelling hypothesis introduced in the TES-PD model was fully validated. In particular, the proposed strategy produced similar heat losses to the one predicted by model m.1 (if used in corresponding conditions), the latter being the only model that accounts for this contribution. Several analyses were also performed looking at the role of non-constant thermophysical properties, tracing their contribution to the overall tank performance. More details of such investigations are given in Section 3 and Section 5.

3. Models Validation

Before using the TES-PD model to predict the storage tank dynamics and, subsequently, the I-ESS performance, it is fundamental to validate it. However, to the authors’ best knowledge, the literature lacks experimental data as well as numerical results obtained with detailed models like the proposal. For these reasons, the models m.1, m.2 and m.3, which partially enclose some of the TES-PD model features, can be used to demonstrate the accuracy and computational efficiency introduced by proposed model.
To this end, there is a need to implement those models and validate them in order to ensure the quality of the further comparisons. In particular, Section 3.1 describes the validation process of the model m.2 while Section 3.2 presents the procedure adopted to validate the model m.3. The validation process of model m.1 is not discussed in this work since it has been extensively analysed in [13,26], while the mathematical description of model m.1, m.2 and m.3 is summarised in Appendix A.

3.1. Validation of Model m.2

To validate the reservoir model according to the hypothesis of McTigue et al. [31], there is a need to consider their PTES plant arrangement [37]. The latter comprises one pair of turbomachines, two heat exchangers, and two packed beds: one at high pressure and temperature (named “hot tank”) and one at ambient pressure but low temperature (called “cold tank”).
During the charging phase, the system acts as a high-temperature heat pump cycle where the electricity is first converted into heat through the compressor and then stored as heat and cold in the hot and cold reservoirs, respectively. During the discharge, the plant acts as a closed Brayton cycle heat engine; thus, the system allows for extracting the stored energy and re-generating electricity. Both the hot and cold tanks are characterised by a height L of 5 m and a diameter D of 5 m; thus, L / D = 1 .
The void fraction, ε , is set equal to 0.33, while the particle diameter d is assumed to be 0.03 m. The working fluid is argon and the mass-flow-rate, m ˙ , is set equal to 12.5 kg·s−1. The pressure ratio between the two tanks is set equal to 10. Further details of the plant description and modelling assumptions can be found in [31].
Figure 1 shows the temperature trends inside the cold tank, considering a stabilized cyclic operation. The charging and discharging periods depend on a crucial parameter for cyclic operations: the utilization factor Π . This parameter, equal to 0.75 for both the phases, indicates the ratio between the charging and nameplate time (based on the thermal front speed estimation in ideal conditions).
Temperature trends refer to three different operating conditions during the charging process, considering constant or variable the solid storage material’s specific heat c p s . As given in Figure 1, the depicted trends are in good agreement with the reference for both constant and variable solid properties. In particular, the initial stages’ temperature trends show a maximum difference of 5%, while, at the end of the process, this decreases to almost 1%. The uncertainty for some settings (e.g., the update frequency or the component models) and material property values may affect the results, and the slight mismatches are to be ascribed to these factors.
To further validate the implemented model, another possible configuration for this PTES system emerges from the optimisation study performed by McTigue et al. [31]. In the reference, a more realistic scenario arises since the design implies the use of a polytropic efficiency equal to 0.9. The graphical extraction of the geometrical properties leads to a L / D ratio of 0.57 and 0.23 for the hot and cold reservoir, respectively, a void fraction ε of 0.35, and a particle diameter of 0.63 cm for both the tanks. The pressure ratio during the charging phase was maintained equal to 10.92, while during the delivery, it decreased to 9.105.
Figure 2 shows the round-trip efficiency and the energy density (i.e., the ratio between the restored energy and the total volume of the reservoirs) trends for the first 40 cycles, employing a utilisation factor of 0.25. The steady-state is achievable after many more cycles due to the lower utilisation factor Π . The reported trends for these two parameters consider a constant as well as a variable solid’s specific heat. The comparison involves the values reported by McTigue et al. [31], which were 0.59 for the efficiency and 65 MJ m−3 for the energy density.
Since McTigue et al. [31] sized the PTES system to produce an ideal energy density of 300 MJ m−3; the results in terms of energy density are in full agreement with the available data.
Hence, the model m.2 can be considered fully reproduced because the temperature profiles, the round-trip efficiency and the energy density trends were replicated concerning boththe constant and variable solid specific heat at constant pressure, c p , s .

3.2. Model m.3 Validation Process

To validate the bed model according to the model m.3 description, it was necessary to build up the PTES system proposed by Desrues et al. [21]. The latter is similar to the McTigue et al. [31] one, with some changes in the components’ management and cycle design parameters.
The numerical setup, in terms of the grid and time-step, accounts for the following choices: two hundred layers were used for the spatial discretisation allowing for a spatial resolution of Δ x = 0.05 m. Such resolution was found to be sufficient in guaranteeing grid independence, and the reader is directed to Appendix B for further details. Temporal stability is recovered setting the time step, Δ t , equal to 0.01 s. The storage material is considered made of microchannels; thus, the α and β coefficients were selected, consequently. The global cycle update time (i.e., all the inputs and outputs update of each component of the system) was set equal to 120 s.
The process’ steady state was reached after approximately ten cycles, considering a cyclic operation of charge–discharge. The pseudo-stationary results were compared to those reported in the reference (graphically extracted). The plant was managed by controlling the temperature at the outlet of the tanks; thus, the utilisation factor was not defined.
Figure 3 shows the trends of stored and delivered energy. Stored energy represents the whole energy absorbed by the system during the charging phase, while delivered energy is the net energy output during discharge. Another trend shown is the efficiency over cycles (i.e., the ratio between restored and stored energy), which is in good agreement with the reference values.
The mismatch in energy between the authors’ implementation and the reference [21] depends on the not explicitly reported choices: for example, the storage material used or the other component model. Thus, the efficiency trend is more representative of the numerical model adopted since it is more independent of such particular choices.
Figure 4 also illustrates thermal fronts inside the hot and cold tanks. The reported results refer to the 10th cycle at the end of both the charge and discharge phases. As clearly depicted, even these trends show a good match concerning the reference solution, particularly where the temperature gradient is high. Thus, the authors convey that the model m.3 implementation can be considered validated since the maximum temperature difference occurs at the tank limits (mostly influenced by other components’ simulation choices), while the internal area gradient shows a less than 1% difference.
Considering that the TES-PD model is based on model m.3, the authors introduced their model in the plant scheme proposed by Desrues et al. [21], and, by avoiding the use of variable thermophysical properties and heat losses through the surface of the tank, they verified that the obtained results were exactly the ones obtained with model m.3. Therefore, it is realistic to consider the TES-PD model validated under model m.3 conditions. In any case, to support this claim, further investigations will be presented in Section 5 with regard to the selected test case.

4. The Selected Test Case

To carry out an objective comparison among numerical reservoir models, there is a need to include each of them in a plant. For this purpose, the I-ESS configuration developed by Benato and Stoppato [13] has been selected. The plant model is linked to the National Institute of Standards and Technology [32] database and CoolProp [33] library to acquire the storage material and the working fluid thermophysical properties. More information on plant component modelling and validation procedure were reported by Benato and Stoppato [13] and Benato and Stoppato [14], while, in the following, a brief description of the plant scheme and control strategy are given.
The plant layout used to store electricity in the form of thermal energy is depicted in Figure 5a, while Figure 5b illustrates the discharge configuration that allows converting the thermal energy into electricity.
The charging layout consists of a fan (FAN), a heat exchanger (HX), an electric heater (EH) and a high-temperature storage tank as shown in Figure 5a. The plant acts as an open-cycle where the air is used as the working fluid. The fan sucks the air from the ambient atmosphere and forces it to follow the path from 1 to 6. The electricity generated in surplus, e.g., by the VRES plants, is converted into heat through the electric heater.
In practice, the device uses electrical energy to increase the air temperature from T3 to T4. The latter is equal to the maximum cycle temperature, Tmax, and it is a design parameter. Compared to other TES-based layouts where the electricity conversion into heat is done by the compressor, the electric heater allows maintaining T4 constant and equal to the maximum one independently to the value assumed by T3.
After the conversion, the electricity is stored as thermal energy in the hot storage, which made by a packed bed, is initially at ambient temperature and, during the charge, is slowly heated up by the hot air entering it at a temperature equal to Tmax. Since the plant is based on an open-cycle, a heat exchanger is inserted to reduce the heat losses. The device, if needed, should be bypassed and installed only in storage units characterised by medium-to-low size to contain its dimensions and, subsequently, the costs.
Note that the charge stops when the mismatch between the maximum temperature, Tmax, and the air temperature at the tank’s outlet, T5, is equal to a pre-computed value, named the charging tolerance ( t o l ). The charging tolerance (which can be equal to, e.g., 10, 100 or 300 K) widely affects the performance of the plant in terms of the charging and delivering times ( t C H G and t D I S ), stored and delivered energy ( E C H G and E D I S ) and round-trip efficiency ( η r t ). Therefore, based on the service requested to the I-ESS plant, the charging tolerance needs to be carefully selected and, in the case of switching the operation mode, reassessed.
Figure 5b depicts the delivery arrangement and shows the route follows by the air: from 11 to 16. The system is composed of a compressor, an expander, an electric motor/generator (MG), a heat exchanger (HX) and a high-temperature tank. The generation unit acts as a modified gas turbine in which the combustion chamber has been replaced with the high-temperature storage tank.
As previously stated, the idea of configuring the I-ESS plant using gas turbine components stems from the fact that gas turbines operating in open-cycle are not able to compete in the actual electricity market due to high operating costs and low efficiency. However, in this scenario, their utilisation hours dropped rather steeply. Thus, gas turbines conversion in storage units can be helpful to revamp such plants as well as to redeploy industrial sites.
The compressor sucks the air at ambient conditions and increases its pressure and temperature. Then, the fluid passes through the heat exchanger and enters the hot tank, where it is heated up. At the beginning of the delivery, the air leaves the hot tank at a temperature approximately equal to T 14 = T m a x t o l . At this point, the air is expanded in the turbine and injected into the heat exchanger. The expansion process generates mechanical power. Part of it is used to drive the compressor while the remaining one is converted into electricity through the electric generator.
After the expansion, the air is cooled down into the heat exchanger and, then it is released into the environment. As for the charging scheme, the heat exchanger should be bypassed if needed and can be installed only in storage units characterised by medium-to-small size to contain its dimensions and, subsequently, costs. As is evident, the heat exchanger and the tank shown in Figure 5a are the same components of Figure 5b.
As for the charge, the delivery phase can also be stopped at different points depending on the service requested to the I-ESS plant. The state of discharge influences the time requested for the second to i-th charge, the energy that can be stored during the second to i-th charge, the discharge time and the round-trip efficiency. In addition, during the discharge, the generated power ( P G E N ) by the I-ESS plant decreases with the reduction of the tank temperature. Nevertheless, to guarantee an almost constant generation during the delivery phase, the plant control system stops the discharge when the mismatch between the generated power and the design one is higher than 5%.
The tank’s size has to correspond to the amount of energy that needs to be stored at the required temperature. Therefore, as suggested by Singh et al. [38], the size of the packed bed needs to be fixed so that, during the charge, the bed absorbs the maximum amount of energy made available by the heat transfer fluid flowing through the bed.
In this way, at the beginning of the charge, the fluid leaves the storage at a temperature approximately equal to the initial one. Conversely, at the end of the charge, the bed can achieve a temperature near the one of the fluid entering the storage. Note that, during the charge, the heat transfer fluid flows from the top to the bottom, while, during delivery, the flow is reversed. The tank is cylindrically shaped with an upper plenum, a packed bed storage and a lower plenum as depicted in Figure 5, while it is arranged vertically to avoid buoyancy-driven instabilities of the thermal fronts.
We selected an existing tank placed in an underutilized facility that can accommodate a packed bed with a volume of 500 m3 and a height of 8.30 m. We assumed that the tank temperature at the beginning of the first charge is equal to ambient temperature, 25 °C. The packed bed void fraction was assumed equal to 0.4 because the tank is filled with spheres of aluminium oxide (0.05 m in diameter) with a density and a mean specific heat equal to 3990 kg m−3 and 1150 J kg−1 K−1, respectively [39].
The authors selected the above-mentioned packed bed specifications following the indications given by Ataer [40] and Singh et al. [41]. As suggested by Desrues [42], for this geometrical arrangement, the solid material’s effective thermal conductivity can be assumed up to 1 W m−1 K−1, while the overall heat loss coefficient U i can be set equal to 0.7 W m−2 K−1.
Considering again the idea of converting underutilized thermal units into a storage system, the authors, with the help of an Italian gas turbine manufacturer, selected a 3 MW gas turbine unit installed in an Italian industrial site as a device to be used for the delivery phase. Its operating parameters included a mass-flow rate of 16 kg s−1, a pressure ratio of 9.5 and a maximum temperature at the turbine inlet of 1123 K. Therefore, the maximum temperature of the cycle, Tmax, was fixed equal to 1123 K.

5. Results and Discussion

Before proceeding with the model comparisons, it is crucial to remark that each model is not different from the others concerning the structure only; however, there are even different modelling choices. Model m.1 considers the heat losses, while model m.3 considers the solid thermal conductivity contribution but it is not reciprocal. Model m.2 does not consider any of these effects. In addition to these aspects, model m.1 and m.3 required constant properties for both solid and fluid, i.e., c p , s , c p , f , k f and μ f while model m.2 and model TES-PD introduce the property variability.
Intending to understand the influence on the plant performance of both the heat loss coefficient (considered in m.1) and the effective thermal conductivity (considered in m.3 and TES-PD), the authors first compare the models overriding these contributions. However, considering that models m.2 and TES-PD required variable properties and that model TES-PD becomes equal to m.3 when the constant properties are accounted for and the heat loss coefficient and the effective conductivity neglected, the authors ran models m.1, m.2 and m.3 examining different temperatures at which the fluid and the solid material properties were computed.
Note that the authors imposed the same constant property behaviour even for model m.2, aiming to compare the results. After that, the obtained findings can also be compared with the outcomes retrieved from model m.2 and TES-PD model run with variable properties of the fluid and the solid. The outcomes of this analysis are given in Section 5.1.
In Section 5.2, instead, we investigated the role of the charging tolerance in the tank performance. The dynamic performance of a TES-based storage plant was affected not only by the fluid and solid properties but also by the heat losses and effective conductivity and charging tolerance. The lower the mismatch between the maximum temperature and the air temperature at the tank’s outlet was, the higher the charging time and the stored energy were.
Therefore, it is interesting to compare the models’ behaviour also with different charging tolerance. To overcome issues related to fluid and solid properties evaluations, the comparison among the models was made considering constant properties also in the case of model m.2 and neglecting the heat loss coefficient and the effective conductivity. Under these conditions, model TES-PD falls back to the results obtained from model m.3.
Finally, in Section 5.3, the models are compared considering their original formulation to underline the pros and cons of a more complete formulation. Regarding the storage performance affected by the solid material conductivity and given the lack of experimental data, the authors performed a sensitivity analysis based on literature data and available correlations.
As far as the simulations setup, for each model, a specific layer height was based on the solution sensibility considerations, ensuring less than a 2% difference in the observed results, doubling the number of layers. Thus, the chosen values of grid refinement were Δ x equal to 0.1 m (83 layers, [13]) and 0.0166 m (500 layers, [37]) for model m.1 and m.2, respectively. As in the previous validation of model m.3, 200 layers (i.e., Δ x = 0.0415 m) were used for the integration of both model m.3 and the TES-PD model.
Comparisons were made analysing the plant performance parameters: charging and delivering times, stored and delivered energy, generated power and round-trip efficiency evaluated for the first cycle, with the first stable cycle and the average cycle calculated by averaging 50 stable cycles.

5.1. Influence on Fluid and Solid Properties of the Reference Temperature

Dealing with constant properties leads to assuming a temperature corresponding to the properties’ reference value. Admissible reference temperatures can be the ambient temperature, the maximum temperature of the cycle or a mean value between them. In this analysis, due to the absence of literature specifications on reference temperature selection, the authors tested three values to underline the differences introduced by this assumption.
Table 1 lists the results obtained when the fluid and storage material properties c p , s , c p , f , k f and μ f were computed with a reference temperature equal to T a m b = 25 C, T m a x = 850 C and T m e a n = 437.5 C. We also list the outcomes of simulations computed using properties variable in space and time for the models that allow these features (m.2 and TES-PD).
T a m b columns list the system performance parameters when setting the reference temperature for properties evaluation to the ambient temperature. The maximum difference involves the first cycle’s round-trip efficiency, whose relative difference is about 7% considering different models. Then, in a stabilised cycle, the discrepancies become more relevant for discharging times and energy. T m a x columns report the system performance parameters when the reference temperature for properties evaluation was set to 850 C.
The maximum differences drop from 7% to 4% concerning the first charge, compared to the T a m b case. In a stabilised cycle, the differences decrease to almost 3.5%, which is a minor mismatch. When assuming T = T m e a n = 437.5 C as the reference temperature, the relative differences are almost comparable to the T = T m a x reference temperature case.
Observing these simulation outcomes, it is evident that the differences among the models at a specific reference temperature are negligible. Model m.3 always predicted a higher efficiency, even if the differences were generally variable but, still, less than 9%. This concerns only the differences among the models. More meaningful is the evaluation of how the choice of the reference temperature affects the results.
In this case, since models m.2 and TES-PD can deal with variable properties, it is also interesting to add these results to the discussion (columns T v a r ). Note that, assuming variable properties means not only updating values during the time-marching procedure but also layer by layer. The solid and fluid’s specific heat and the air’s viscosity and conductivity were evaluated on the single layer’s local temperature basis.
From the obtained results, it is clear that the selection of the reference temperature is a fundamental step because adopting T a m b instead of T m e a n or T m a x caused an underestimation of the first and average charge stored energy ( E C H G ) of 37% and 42%, respectively. Similarly, the energy deliverable during the discharges ( E D I S ) was underestimated by a value higher than 47.5%.
Underling these underestimations in the energy storage capacity are essential since it predominantly affects the design volume of the storage reservoir. Independently to the adopted model, selecting T a m b as reference temperature instead of T m e a n or T m a x can lead to the selection of a storage volume 50% smaller than the one requested from, e.g., the VRES plant. In addition to the above-mentioned design issue, the reference temperature selection also affects the time required to charge, t C H G and discharge the storage, t D I S .
In fact, the adoption of T a m b as a reference reduced the first and the subsequent charges times of 33% and 37%, respectively, while the discharge time was underestimated by 47.7%. These underestimations of the charging and the delivery times are the source of wrong evaluations of storage plant market availability and response time: parameters that need to be accurately estimated to avoid storage plant inability of compensating the mismatch between, e.g., the VRES supply and network demand, which can be the cause of devices fault or, in the worst scenario, grid blackout.
Analysing the columns of Table 1, which refer to T v a r , and comparing them with the outcomes of the simulations derived by the selection of T m e a n or T m a x as the reference temperature, it is clear that the mismatches in the charge and discharge times and energy were lower than 3%. Therefore, it appears that adopting a constant temperature equal to T m e a n or T m a x for properties computations did not introduce relevant discrepancies.
However, since the choice of any reference temperature is purely arbitrary and, considering that the simulation time with variable properties increases by 10% compared to constant properties at T m e a n or T m a x , the use of variable properties can lead to a genuinely consistent model especially in view of including in it parameters, such as heat losses and solid thermal conductivity, as done in the TES-PD model.
Figure 6 shows the temperature trends for these tests, adopting model m.3. Since differences between models in this configuration are not too marked, the figure focuses on a single model results’ comparison with different reference temperature choices. The most significant difference concerns the results derived using the environmental temperature as the reference.
The charge and discharge processes last almost 30% less than in the other cases, thus, involving less energy. The temperature profiles inside the tank at the end of the first stable cycle’s charge and discharge were not greatly affected.
A minor difference concerns the last layer temperature during the first charge, which increased more gradually when the properties were variable. These considerations lead to say that, for these specific conditions, i.e., materials and system behaviour, T m a x or T m e a n are good choices as the reference temperature. However, this can not be generalized. Adopting variable properties, a possibility offered by model TES-PD, is a solution that allows avoiding wrong temperature choices, although it is 10% more computationally expensive.

5.2. Influence of the Charging Tolerance

The temperature used as the reference for fluid and solid properties computations is not the unique parameter that directly affects the performance of TES-based energy storage systems. In fact, the time requested for charging the storage tank and the energy stored in it is directly linked to the charging tolerance, t o l , namely the mismatch between the maximum temperature and the air temperature at the tank’s outlet. The lower the tolerance is, the higher the charging time is.
As for the properties’ reference temperature, the charging tolerance is also a design parameter whose selection depends exclusively on how the plant is operated, which, in turn, relies on the network or VRES plant needs. For this purpose, we investigated the difference in the system dynamics caused by different the charging tolerance and adopting a t o l equal to 10, 100 and 200 K. To better detect model mismatches concerning only the charging tolerance, simulations were performed using model m.1, m.2 and m.3, computing the fluid and solid properties at T m e a n and neglecting the heat loss coefficient and the solid effective thermal conductivity. Table 2 summarises the simulation outcomes.
Even if the maximum relative difference between the characteristic parameters slightly increased (from a 4% to 7% difference in the discharging times and energy during a stabilized operation) when the charging tolerance was 200 K, the most noticeable differences concerned the system’s round-trip efficiency, η r t . In particular, the higher the tolerance is, the higher the round-trip efficiency is. In contrast, the lower the tolerance is, the higher the charging and delivery times and the charging and delivery energy are.
Apart from these aspects, comparing the simulations results for a given charging tolerance revealed minor mismatches among them: usually, less than 6% and 8% in the first and first stable cycle parameters, respectively. Therefore, the model comparison does not underline a significant mismatch in the plant performance prediction if the same charging tolerance is set.
However, looking deeper into the results listed in Table 2 reveals that the charging tolerance influenced the stabilisation of the process. When t o l equalled 10 K, the performance parameters at the first stabilised cycle were identical to the ones obtained from the average cycle calculated by averaging 50 stabilised cycles; this was no longer valid when setting the charging tolerance to 200 K. The different temperature profile at the end of the charge caused a delay in reaching the stabilisation. The delay effect was more significantly remarkable when the tolerance was higher. In fact, when t o l = 10, the temperature was almost T m a x in the whole tank at the charge’s end, making the second discharge process similar to the first one.
To underline this critical aspect, Figure 7 shows the temperature trends over time, adopting the m.3 model with different charging tolerances. Even if the discharge process ends when the generated power is 5% lower than the maximum, the charging tolerance reduction affects this phase. The temperature distributions inside the tank (Figure 7b) show the lower energy availability when setting the tolerance to 200 K. The data in Table 2 state that when t o l = 10 K, the average cycle’s charging energy was about 190 MWh, compared to around 130 MWh when t o l = 200 K.
The same happens to the discharge energy but less markedly since the efficiency slightly increases. In conclusion, the choice of a particular tolerance in the charging phase did not have an exclusive impact on the charging time; however, the whole system’s dynamics were conditioned. This aspect is fundamental because such a tank is hooked up to a highly variable energy source, as the system’s history heavily affects its overall performance.

5.3. Models’ Original Formulation

As previously stated, in this work, we aim to compare the available packed bed models, understand the influence of the constituting parameters, and then develop a new mathematical structure to take into account the most impacting parameters while preserving the reliability and computation efficiency. However, to evaluate the impact on plant performance of the temperature at which the fluid and the solid properties are computed (Section 5.1) and of the charging tolerance (Section 5.2), the models need to work under the same operating conditions.
Thus, the authors set to zero the heat loss term ( U i ) considered in model m.1, the effective thermal conductivity ( k s , e f f ) deemed in model m.3 and the property variability through layers in model m.2. In this manner, the analysis revealed that, for models working with variables computed at a constant temperature, a good compromise was to consider T m e a n (the average temperature between the ambient and the maximum ones) despite the best choice is a model considering layer-by-layer fluid and solid properties’ variability (like model m.2 and TES-PD).
This analysis is crucial in light of a model comparison based on their original formulations. In fact, to carry out this analysis, the properties for model m.1 and m.3 were computed considering T m e a n while m.2 and TES-PD assumed variable properties for each layer, according to the local temperature. The heat loss coefficient U i was set equal to 0.7 Wm−2K−1, while k s , e f f was assumed equal to 1 Wm−1K−1 [30,42]. Regarding the charging tolerance, the authors selected 10 K because they assumed to fully charge the hot storage. Table 3 lists the main outcomes of the performed investigation.
Compared to values reported in Table 1 (column related to m.1 adopting T m e a n ), the consideration of the heat loss coefficient increases the first and average cycle charging times of 6.2 and 7%, while the first and average cycle discharging times were reduced to 8%: a fact that guarantees increasing the energy stored in the tank by 2.2% (first cycle) and 1% (average cycle), but that reduces both the discharge energy and round-trip efficiency by approximately 10%. The plant’s mean generated power resulted as reduced by approximately 1.2%.
Thus, the heat loss coefficient is a parameter that needs to be taken into account in the tank model to correctly estimate its dynamic behaviour and, subsequently, the storage plant performance. In a nutshell, neglecting U i implies a 10% overestimation of the plant deliverable energy, a fact that can be the source of storage plant management issues.
Compared to results listed in Table 1 (column related to m.3 adopting T m e a n ), a model formulation in which the solid conductivity was set equal to 1 Wm−1K−1 as suggested by Desrues [42], did not introduce considerable differences. For the first and the average cycle, the charge and discharge times as well as the charged and delivered energy, the gaps were less than 1%. Therefore, a model formulation in which the solid conductivity was considered constant did not drastically affect the prediction of plant performance.
More interesting is the comparison among the computations performed in Table 1 (column related to TES-PD adopting T v a r ) and the ones given in Table 3. They reveal that neglecting heat losses and solid effective thermal conductivity largely influenced the plant performance prediction. The first cycle charging time resulted as 4.3% lower while the delivering one was reduced by 6.1%.
The stored and delivered energy were 1.6 and 7.7% lower while the plant round-trip efficiency resulted as 9.4% lower in the first charge–discharge cycle and 8.3% considering an average cycle. Thus, the original formulation of the TES-PD model guarantees considering the account property variability, heat losses, and solid conductivity without compromising the simulation time, which was only 12% higher than for the other models.
For the sake of completeness, in Figure 8, we present the non-dimensional temperature trends in the hot tank during the first and the average charge, the average discharge and the temperature distributions inside the tank at the end of the average cycle’s charge (CHG) and discharge (DIS) for the four models.
The trends did not substantially differ among each other, although the completeness of the TES-PD model guaranteed a better estimation of the overall performance. However, in previous simulations, the thermal conductivity value was assumed constant and equal to a value suggested by Desrues [42] despite that solid thermal conductivity can lead to very different results if it is properly modelled. In particular, in contrast to other models, the TES-PD allows the implementation of correlations in which the solid thermal conductivity is variable according to a temperature-based law (see, e.g., [43]).
However, in this case, due to the physical process behaviour, the thermal conductivity can assume a value up to 7, a fact that imposes a drastic reduction of both the simulation time-steps and layer thickness, parameters, which, in turn, affects the computation time, which can be to three orders of magnitude higher. For these reasons and considering the literature lack of accurate data for solid thermal conductivity in the simulated conditions, the authors performed a parametric analysis (see, Table 4), which demonstrated low discrepancies among k s , e f f equal to 0, 1 or 2 W m−1 K−1.
In contrast, the computational time increased by 20% from 1 to 2 W m−1 K−1. Therefore, despite the model’s ability to take into account the effective conductivity variability according to a temperature-based law, the lack of objective data and the vast increase in computational time led the authors to choose, for the design phase, a constant value equal to 1 W m−1 K−1. In future TES-PD updated versions and CFD investigations, the effect of this parameter needs to be further investigated to sharpen the model’s ability to predict storage dynamic behaviour.

6. Conclusions

The environmental impact of power generation units is in boosting wind and solar photovoltaic plant production because they can guarantee green electricity generation. However, such plants are characterised by a variable and intermittent nature, a feature that results in considerable power fluctuations, which, in turn, can provoke grid control issues and device faults. To overcome these problems and guarantee a rapid transition from a fossil-based generation mix to a renewable one, installing large-scale energy storage units is required.
For this purpose, we developed a system able to compete with pumped hydro energy storage, compressed air energy storage and batteries. However, as the proposed storage system is based on the thermal energy storage concept, there is a need to predict the storage performance correctly. This is not an easy task considering that the literature reported three tank models, each one different from the others in terms of how the fluid and storage material properties are computed and the consideration or not of the overall heat loss coefficient and the effective solid conductivity.
Therefore, to highlight the strengths and weaknesses of the models available in the literature, the authors built them and, after a validation process, performed a comparison. The analysis revealed the influence of the temperature at which the solid and fluid properties were computed as well as the need to consider the thermal losses and the solid thermal conductivity.
Starting from these outcomes, the authors developed a new and more complete model for the tank, namely the TES-PD model. The model accounts for the thermophysical properties’ local time and space variability, and it embeds the thermal losses and the solid thermal conductivity. Features that avoid any uncertainty concerning the reference temperature definition, and considering the solid and fluid properties variability layer-by-layer, allowed for a better performance prediction of both the storage and the plant in which it works.
In addition, the TES-PD model can take charge of real properties, paving the way for the study of innovative fluids, materials and configurations regarding thermal storage issues if these parameters are available from accurate measurements, dataset or CFD analysis. Finally, it is possible to claim that the performed analysis on the TES-PD model revealed that such differential and time-accurate description of the storage tank, even if one-dimensional, allowed a fast and effective prediction of the performance of both the tank and the storage plant; this is essential information for the preliminary design of innovative large-scale storage units operating with thermal storage.

Author Contributions

Conceptualization, A.B. and F.D.V.; Formal analysis, A.B., F.D.V., E.G., A.S. and G.C.; Methodology, A.B., F.D.V., E.G., A.S. and G.C.; Validation, A.B., F.D.V., E.G., A.S. and G.C.; Writing—original draft, A.B. and F.D.V.; Writing—review and editing, A.B., F.D.V., E.G., A.S. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CAESCompressed Air Energy Storage
CFDComputational Fluid Dynamics
CHGCharge
DISDischarge
EHElectric Heater
GENGenerated
HXHeat Exchanger
I-ESSIntegrated Energy Storage System
MGMotor/Generator
PHSPumped Hydro Storage
PTESPumped Thermal Energy Storage
TESThermal Energy Storage
VRESVariable Renewable Energy Sources

Appendix A. Detailed Description of Most Used Numerical Models for Energy Storage Tanks

In this section, the most used numerical models for storage tanks are described with a high detail level. The logic of an increased complexity governs the description. Despite the models are already available in literature, only a few details were reported regarding their numerical implementation; thus, the authors present them in this section to highlight their peculiarities and the method that they developed for their mathematical resolution.

Appendix A.1. Model m.1

The model m.1, presented for the first time in 1982 by Howell et al. [30], represents one of the most straightforward schemes for packed bed numerical description as demonstrated by its large use in the literature (see, e.g., [13,26,38,41,44,45]). The tank model consists of discretising the bed in N number of elements over the height of L; thus, Δ x = L / N represents the spacing within the bed layers. Inside any layers, the temperature and the pressure are assumed constant and uniform throughout the single bed element (the cross-sectional area is A). With this hypothesis, the pressure drop in the tank, Δ p , can be computed as:
Δ p = C f L G 2 d ρ f
where C f is the friction coefficient, G = m ˙ / A is the specific mass-flow rate, m ˙ is the mass-flow rate, d is the particle equivalent diameter and ρ f is the fluid density. The fluid temperature distribution along the bed is computed as:
T f , i t = T s , i t + ( T f , i 1 t T s , i t ) Φ 1
where
Φ 1 = h v A L N ( m ˙ c p , f ) = N T U N
The temperature of the bed, instead, holds as follows:
T s , i t + Δ t = T s , i t + [ Φ 2 ( T f , i 1 t T f , i t ) Φ 3 ( T s , i t T a m b ) ] Δ t
where
Φ 2 = N m ˙ c p , f ρ s A L c p , s ( 1 ε ) Φ 3 = U i Δ A i m ˙ c p , f Φ 2
Here, T f , i is the temperature of the i-th fluid layer while T s , i is the temperature of the i-th bed’s layer. The density of the solid material constituting the packed bed is denoted as ρ s , while c p , s and ε are the solid specific heat at constant pressure and the void fraction of the bed, respectively. The time-step, Δ t , is evaluated as Δ t = ( ρ f L A ε ) / ( N m ˙ ) . N T U is the Number of Transfer Units, U i is the overall heat loss coefficient, and Δ A i is the external area of a layer. The volumetric heat transfer coefficient, h v , depends on the thermal conductivity of the fluid, k f ; the hydraulic diameter, d; and the Nusselt number, N u , according to the following expression:
h v = N u k f d 2
To compute such parameters, the correlations reported by Singh et al. [46] are employed; thus, the Nusselt number is evaluated as:
N u = 0.437 ( R e ) 0.75 ( Ψ ) 3.35 ( ε ) 1.62 e x p [ 29.03 ( log Ψ ) 2 ]
Note that R e = ( G d ) / μ f is the Reynolds number, μ f is the dynamical fluid viscosity and Ψ = A s / A e is particle sphericity (i.e., the ratio between the surface area of a spherical particle of the same volume and the particle surface area). Finally, always according to Singh et al. [46], the friction coefficient is computed as:
C f = 4.466 ( R e ) 0.2 ( Ψ ) 0.696 ( ε ) 2.945 e x p [ 11.85 ( log Ψ ) 2 ]
The above set of correlations is useful since they provide an algebraic description of storage tank and they are valid for packing material with a geometry different from the spherical one, only through the sphericity definition.

Appendix A.2. Model m.2

Model m.1 gives considerable success in the literature due to its simplicity. However, this fact is limiting in more complex applications or whenever time-accurate details are required. The model m.2 follows the suggestion of McTigue et al. [31]. The model characterises the heat exchanged and the pressure losses within the reservoirs based on the Schumann method [47]. Thus, the flow is still assumed as one-dimensional, and the exchanged heat is limited by surface effects. The model is detailed by [37], and here a brief description follows.
According to Schumann [31,47], the fluid and the solid temperature, T f and T s , are assumed to obey the following set of differential equations:
T f x = T s T f
T s t = T f T s τ
where and τ are, respectively, the lengths and the time scale defined as:
= 1 S t ( 1 ε ) S v
τ = ρ s c p , s S t c p , f G S v
Pressure losses are accounted with the following expression:
d p d x = S v ( 1 ε ) C f G 2 2 ε 3 ρ f
where C f is the pressure loss coefficient, S t = h / ( G c p , f ) denotes the Stanton number and S v is the surface to volume ratio of the solid particles. Other symbols follow the definitions provided in Appendix A.1. Based on these hypotheses, White et al. [37] proposed an improved version of Schumann’s model by introducing variable properties for the solid and the fluid, including momentum and mass continuity equations. Accounting for these contributions, it is necessary to consider the equations governing the heat and mass exchanges that occur in an infinitesimal control volume along with the tank. The latter is expressed as:
δ Q ˙ = A δ x ( 1 ε ) S v h ( T f T s )
where h is the superficial heat transfer coefficient. Then, mass continuity, momentum and energy (for both fluid and solid) equations in a differential formulation can be expressed as:
ε ρ f t + G x = 0
G t + x G 2 ρ f + ε 2 ( p ρ f g x ) = ε ( 1 ε ) τ S v
ε t ( ρ f e f ) + x ( G h f ) = ( 1 ε ) S v h ( T s T f )
ρ s c p s T s t = S v h ( T f T s )
The particle surface shear stress is represented by τ , p is the pressure, g is the gravity acceleration, x is the axial coordinate, e f is the internal fluid energy per unit mass, and h f is the fluid enthalpy per unit mass. Equation (A14a–d) can be solved in a fully coupled way employing a time-marching algorithm combined with a spatial discretisation method. However, as suggested by White et al. [37], to reduce computation effort, the model can be cast in a semi-analytical version, based on the fact that the main terms of the Equation (A14a–d) can be simplified according to the Schumann model (i.e., Equations (A8) and (A9)). To evaluate the time-dependent terms contribution, a correction factor F can be defined and expressed as:
F = ε G c p f p t ρ f c p f T f t
The latter expression, based on Equation (A14a,c), can be averaged and added to the fluid energy equation. Thus, a semi-analytical model is obtained by integrating the Schumann model and accounting for the factor F. The resulting discrete scheme reads as follows:
T f , i t + Δ t = ( T s ¯ + F ¯ ) 1 e x p Δ x + T f , i 1 t + Δ t e x p Δ x
T s , i t + Δ t = T f ¯ 1 e x p Δ t τ + T s , i t e x p Δ t τ
where the mean temperature of both the solid, T s ¯ , and the fluid, T f ¯ , are defined as:
T s ¯ = T s , i 1 t + Δ t + T s , i t + Δ t 2
T f ¯ = T f , i t + T f , i t + Δ t 2
The Equation (A14a) is implemented in the form:
G i t + Δ t = G i 1 t + Δ t + ε Δ x Δ t ( ρ f , i t ρ f , i t + Δ t )
As final remarks, White et al. [37] suggest adopting a Δ x and a Δ t less than 20 % of and τ , respectively; assumptions that are sufficient to obtain the grid independence of the results. The correlations used to evaluate the Stanton number, S t , as well as the particle and modified Reynolds numbers are reported in [37,48]. Fluid and solid properties are updated during the time-marching procedure using the data available in databases. Since the available literature does not provide any details concerning properties acquisition, in the present work, the authors proposed to acquire the storage material properties from the National Institute of Standards and Technology database [32]. The fluid properties are acquired from CoolProp [33] library.

Appendix A.3. Model m.3

Model m.3 follows the description proposed by Desrues et al. [21], which consists of a purely one-dimensional time-marching approach. The model fits a simple geometry consisting of a regular square array of square sub-channels. Assuming constant solid properties as stated by Desrues et al. [21], the equations used to describe the transient phenomenon for the specific geometry are:
ρ f t + ( ρ f v f ) x = 0
( ρ f T f ) t + ( ρ f T f v f ) x = 4 N u k f c p , f ( d h ) 2 ( T s T f )
ρ s c p , s T s t + k s 2 T s x 2 = ε 1 ε 4 N u k f ( d h ) 2 ( T f T s )
The model is closed by the pressure drop constitutive equation and the ideal gas equation of state
p x = C f 4 d h 1 2 ρ f v f 2
p = ρ f r T f
Since no detailed descriptions are reported in the literature concerning the numerical treatment of this model, the authors integrated it with the same numerical strategy described for the TES-PD model. Thus, a first-order generalised upwind scheme is used to treat the first-order derivatives, while a second-order central method is employed for the second derivative contributions. Temporal contributions are accounted through a first-order explicit Euler method. In addition, the pressure-velocity coupling is demanded to a segregated approach that consists of iterating up a certain convergence tolerance among every time-step. As mentioned in the TES-PD model description, such a numerical treatment is found a good compromise between accuracy and computational cost. Heat transfer and pressure drop coefficients evaluation require the Nusselt number and the friction factor correlations, the latter being related to the employed regenerator. The interested reader can refer to Desrues [42] for more details, while in the following, the adopted correlations are reported concerning spheres packing geometry.
  • Microchannels structure. The friction factor is defined as:
    C f = R e 1 · 24 ( 1 1.3553 α + 1.9467 α 2 1.7012 α 3 + 0.9564 α 4 0.2537 α 5 )
    where α = L d e p t h / L w i d t h defines the aspect ratio of the channels. The Nusselt number is defined as
    N u = 7.541 ( 1 2.610 α + 4.970 α 2 5.119 α 3 + 2.702 α 4 0.548 α 5 )
    and it is used to compute the superficial heat transfer coefficient and the effective thermal conductivity:
    k s , e f f = ( 1 ε ) k s
  • Packing of spheres The hydraulic diameter, d h , is given as:
    d h = ε 1 ε d
    while the hydraulic Reynolds, R e h , which is linked to the superficial one, can be defined as:
    R e h = 1 1 ε R e
    and it is used to determine the friction factor based on the following equation:
    C f = 2 ( A / R e h + B ) , A = 150 , B = 1.75 ,
    which represents the Ergun equation [49] for spheres packing, and the Wakao & Kaguei [50] correlation for Nusselt number
    N u = 2 + 1.1 ( P r ) 1 / 3 ( R e ) 3 / 5
    where P r = ( μ f c p , f ) / k f is the Prandtl number and R e = ( G d ) / μ f is the Reynolds number.

Appendix B. Convergence Analysis and Grid Sensitivity

The section aims at performing a grid sensitivity analysis applied to the TES-PD model. In particular, a parametrical study is performed to tune the number of points required for the grid independence of the results. The analysis outcomes are depicted in Figure A1. In particular, Figure A1 shows the temperature trends inside the tank at the end of the first charge and discharge cycles, adopting an increasing number of subdivisions, N. Table A1 lists the main performance parameters. The reader can notice how the computational Time of Completion increases rapidly with a finer discretisation. The authors convey in tolerating a 2% variation of all the performance parameters, counting that further refinement of the grid leads to almost doubling calculation times, which affects the final use and the scope of a one-dimensional model. For these reasons, a number of layers equal to 200 are found a good compromise between numerical accuracy and computational efficiency. The same number of layers is used in the model m.3 basing on the formally equivalent numerical structure compared to TES-PD.
Figure A1. TES-PD model behaviour with increasingly refined grid level.
Figure A1. TES-PD model behaviour with increasingly refined grid level.
Fluids 06 00256 g0a1
Table A1. Performance parameters as function of the grid resolution.
Table A1. Performance parameters as function of the grid resolution.
Parameter[·] N = 20 N = 50 N = 100 N = 200 N = 300
t C H G [h]39.8234.4032.2230.9830.55
t D I S [h]11.9314.1515.2215.8716.10
E C H G [MWh]374.76354.38346.06341.32339.65
E D I S [MWh]30.7936.5039.2340.9041.49
η r t [%]8.2210.3011.3411.9812.22
T O C C H G [min]2.253.395.459.6413.97
T O C D I S [min]0.621.292.314.386.51

References

  1. International Energy Agency. World Energy Outlook 2016. 2016. Available online: https://www.iea.org/reports/world-energy-outlook-2016 (accessed on 12 July 2021).
  2. International Energy Agency. IEA Statistics, 2019. CO2 Emissions from Fuel Combustion: Highlights. 2019. Available online: https://iea.blob.core.windows.net/assets/eb3b2e8d-28e0-47fd-a8ba-160f7ed42bc3/CO2_Emissions_from_Fuel_Combustion_2019_Highlights.pdf (accessed on 12 July 2021).
  3. Jenkins, J.D.; Zhou, Z.; Ponciroli, R.; Vilim, R.; Ganda, F.; de Sisternes, F.; Botterud, A. The benefits of nuclear flexibility in power system operations with renewable energy. Appl. Energy 2018, 222, 872–884. [Google Scholar] [CrossRef]
  4. Cany, C.; Mansilla, C.; Mathonnière, G.; Da Costa, P. Nuclear power supply: Going against the misconceptions. Evidence of nuclear flexibility from the French experience. Energy 2018, 151, 289–296. [Google Scholar] [CrossRef]
  5. Hentschel, J.; Spliethoff, H. A parametric approach for the valuation of power plant flexibility options. Energy Rep. 2016, 2, 40–47. [Google Scholar] [CrossRef] [Green Version]
  6. Stoppato, A.; Benato, A.; Mirandola, A. Assessment of stresses and residual life of plant components in view of life-time extension of power plants. Proc. ECOS 2012, 2012, 26–29. [Google Scholar]
  7. Buschmeier, M.; Kleinwächter, T.; Feldmüller, A.; Köhne, P.; Energy, S. Improving Flexibility of the Combined Cycle Power Plant Hamm Uentrop to Cover the Operational Profiles of the Future. POWER-GEN Europe. 2014. Available online: https://assets.siemens-energy.com/siemens/assets/api/uuid:1e178af1d8e18630f18de72eb4d65fde3f99dcd5/improving-flex-ccpp-hamm-pspg-b10229-00-7600-en-lr-final.pdf (accessed on 12 July 2021).
  8. Therkorn, D.; Gassner, M.; Lonneux, V.; Zhang, M.; Bernero, S. CCPP operational flexibility extension below 30% load using reheat burner switch-off concept. In ASME Turbo Expo 2015: Turbine Technical Conference and Exposition; American Society of Mechanical Engineers Digital Collection: New York, NY, USA, 2015. [Google Scholar]
  9. Encabo Caceres, I.; Mocholí Montañés, R.; Nord, L.O. Flexible Operation of Combined Cycle Gas Turbine Power Plants with Supplementary Firing. Appl. Energy 2020, 278, 115630. [Google Scholar] [CrossRef]
  10. International Energy Agency. Technology Roadmap: Energy Storage. Available online: https://www.iea.org/reports/technology-roadmap-energy-storage (accessed on 12 July 2021).
  11. Benato, A.; Stoppato, A. Pumped Thermal Electricity Storage: A technology overview. Therm. Sci. Eng. Prog. 2018, 6, 301–315. [Google Scholar] [CrossRef]
  12. Vinnemeier, P.; Wirsum, M.; Malpiece, D.; Bove, R. Integration of heat pumps into thermal plants for creation of large-scale electricity storage capacities. Appl. Energy 2016, 184, 506–522. [Google Scholar] [CrossRef] [Green Version]
  13. Benato, A.; Stoppato, A. Energy and cost analysis of an Air Cycle used as prime mover of a Thermal Electricity Storage. J. Energy Storage 2018, 17, 29–46. [Google Scholar] [CrossRef]
  14. Benato, A.; Stoppato, A. Integrated Thermal Electricity Storage System: Energetic and cost performance. Energy Convers. Manag. 2019, 197, 111833. [Google Scholar] [CrossRef]
  15. Aneke, M.; Wang, M. Energy storage technologies and real life applications—A state of the art review. Appl. Energy 2016, 179, 350–377. [Google Scholar] [CrossRef] [Green Version]
  16. Benato, A.; Stoppato, A. Energy and Cost Analysis of a New Packed Bed Pumped Thermal Electricity Storage Unit. J. Energy Resour. Technol. 2017, 140, 020904. [Google Scholar] [CrossRef]
  17. Gür, T.M. Review of electrical energy storage technologies, materials and systems: Challenges and prospects for large-scale grid storage. Energy Environ. Sci. 2018, 11, 2696–2767. [Google Scholar] [CrossRef]
  18. Koohi-Fayegh, S.; Rosen, M.A. A review of energy storage types, applications and recent developments. J. Energy Storage 2020, 27, 101047. [Google Scholar] [CrossRef]
  19. International Renewable Energy Agency (IRENA). Electricity Storage and Renewables: Costs and Markets to 2030; Technical Report; IRENA: Abu Dhabi, United Arab Emirates, 2017. [Google Scholar]
  20. Divya, K.C.; Østergaard, J. Battery energy storage technology for power systems—An overview. Electr. Power Syst. Res. 2009, 79, 511–520. [Google Scholar] [CrossRef]
  21. Desrues, T.; Ruer, J.; Marty, P.; Fourmigué, J.F. A thermal energy storage process for large scale electric applications. Appl. Therm. Eng. 2010, 30, 425–432. [Google Scholar] [CrossRef] [Green Version]
  22. Fuchs, G.; Lunz, B.; Leuthold, M.; Sauer, D.U. Technology Overview on Electricity Storage: Overview on the Potential and on the Deployment Perspectives of Electricity Storage Technologies; Institute for Power Electronics and Electrical Drives (ISEA), RWTH Aachen University: Aachen, Germany, 2012. [Google Scholar]
  23. White, A.; Parks, G.; Markides, C.N. Thermodynamic analysis of pumped thermal electricity storage. Appl. Therm. Eng. 2013, 53, 291–298. [Google Scholar] [CrossRef] [Green Version]
  24. Luo, X.; Wang, J.; Dooner, M.; Clarke, J. Overview of current development in electrical energy storage technologies and the application potential in power system operation. Appl. Energy 2015, 137, 511–536. [Google Scholar] [CrossRef] [Green Version]
  25. Amirante, R.; Cassone, E.; Distaso, E.; Tamburrano, P. Overview on recent developments in energy storage: Mechanical, electrochemical and hydrogen technologies. Energy Convers. Manag. 2017, 132, 372–387. [Google Scholar] [CrossRef]
  26. Benato, A. Performance and cost evaluation of an innovative Pumped Thermal Electricity Storage power system. Energy 2017, 138, 419–436. [Google Scholar] [CrossRef]
  27. McTigue, J.; White, A. A Comparison of Radial-flow and Axial-flow Packed Beds for Thermal Energy Storage. Energy Procedia 2017, 105, 4192–4197. [Google Scholar] [CrossRef]
  28. McTigue, J.D.; Markides, C.N.; White, A.J. Performance response of packed-bed thermal storage to cycle duration perturbations. J. Energy Storage 2018, 19, 379–392. [Google Scholar] [CrossRef]
  29. Alva, G.; Lin, Y.; Fang, G. An overview of thermal energy storage systems. Energy 2018, 144, 341–378. [Google Scholar] [CrossRef]
  30. Howell, J.R.; Bannerot, R.B.; Vliet, G.C. Solar Thermal Energy Systems: Analysis and Design; Mcgraw-Hill College: New York, NY, USA, 1982. [Google Scholar]
  31. McTigue, J.D.; White, A.J.; Markides, C.N. Parametric studies and optimisation of pumped thermal electricity storage. Appl. Energy 2015, 137, 800–811. [Google Scholar] [CrossRef] [Green Version]
  32. NIST. Nist-Janaf Thermochemical Tables. Available online: https://janaf.nist.gov (accessed on 22 December 2020).
  33. Bell, I.H.; Wronski, J.; Quoilin, S.; Lemort, V. Pure and pseudo-pure fluid thermophysical property evaluation and the open-source thermophysical property library CoolProp. Ind. Eng. Chem. Res. 2014, 53, 2498–2508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dahlquist, G.; Björck, Å. Numerical Methods in Scientific Computing, Volume I; SIAM Publications Online; University of Auckland: Auckland, New Zealand, 2008. [Google Scholar]
  35. Moin, P. Fundamentals of Engineering Numerical Analysis; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  36. De Vanna, F.; Benato, A.; Picano, F.; Benini, E. High order conservative formulation of viscous terms for variable viscosity flows. Acta Mech. 2021, 232, 2115–2133. [Google Scholar] [CrossRef]
  37. White, A.; McTigue, J.; Markides, C. Wave propagation and thermodynamic losses in packed-bed thermal reservoirs for energy storage. Appl. Energy 2014, 130, 648–657. [Google Scholar] [CrossRef] [Green Version]
  38. Singh, R.; Saini, R.; Saini, J. Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size—Part I. Open Fuels Energy Sci. J. 2014. [Google Scholar] [CrossRef] [Green Version]
  39. Poling, B.E.; Thomson, G.H.; Friend, D.G.; Rowley, R.L.; Wilding, W.V. Perry’s Chemical Engineers’ Handbook; McGraw Hill: New York, NY, USA, 2008. [Google Scholar]
  40. Ataer, O.E. Storage of thermal energy. In Energy Storage Systems, Encyclopedia of Life Support Systems; Developed under the Auspices of the UNESCO; Eolss Publishers: Oxford, UK, 2006; pp. 1–20. [Google Scholar]
  41. Singh, H.; Saini, R.P.; Saini, J.S. A review on packed bed solar energy storage systems. Renew. Sustain. Energy Rev. 2010, 14, 1059–1069. [Google Scholar] [CrossRef]
  42. Desrues, T. Stockage Massif d’éLectricité Sous Forme Thermique. Ph.D. Thesis, Univeristé de Grenoble, Grenoble, France, 2011. [Google Scholar]
  43. Esence, T.; Bruch, A.; Molina, S.; Stutz, B.; Fourmigue, J.F. A review on experience feedback and numerical modeling of packed-bed thermal energy storage systems. Sol. Energy 2017, 153, 628–654. [Google Scholar] [CrossRef]
  44. Singh, R.; Saini, R.; Saini, J. Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size—Part II. Open Fuels Energy Sci. J. 2014. [Google Scholar] [CrossRef] [Green Version]
  45. Singh, R.; Saini, R.; Saini, J. Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size—Part III. Open Fuels Energy Sci. J. 2008, 1, 102–106. [Google Scholar] [CrossRef] [Green Version]
  46. Singh, R.; Saini, R.P.; Saini, J.S. Nusselt number and friction factor correlations for packed bed solar energy storage system having large sized elements of different shapes. Sol. Energy 2006, 80, 760–771. [Google Scholar] [CrossRef]
  47. Schumann, T. Heat transfer: A liquid flowing through a porous prism. J. Frankl. Inst. 1929, 208, 405–416. [Google Scholar] [CrossRef]
  48. White, A.J. Loss analysis of thermal reservoirs for electrical energy storage schemes. Appl. Energy 2011, 88, 4150–4159. [Google Scholar] [CrossRef] [Green Version]
  49. Hicks, R. Pressure drop in packed beds of spheres. Ind. Eng. Chem. Fundam. 1970, 9, 500–502. [Google Scholar] [CrossRef]
  50. Wakao, N.; Kaguei, S.; Funazkri, T. Effect of fluid dispersion coefficients on particle-to-fluid heat transfer coefficients in packed beds: Correlation of Nusselt numbers. Chem. Eng. Sci. 1979, 34, 325–336. [Google Scholar] [CrossRef]
Figure 1. Temperature distribution inside the cold reservoir during the charging process, at different charging stages (5%, 50% and 95%). Our are compared with the results reported in [31,37]. In black empty dots, temperature profiles are reported considering constant the storage material’s specific heat (the mean value is selected). In black filled dots, temperature profiles concern the case of a variable storage material’s specific heat: each layer has a different value for this parameter, according to the actual temperature.
Figure 1. Temperature distribution inside the cold reservoir during the charging process, at different charging stages (5%, 50% and 95%). Our are compared with the results reported in [31,37]. In black empty dots, temperature profiles are reported considering constant the storage material’s specific heat (the mean value is selected). In black filled dots, temperature profiles concern the case of a variable storage material’s specific heat: each layer has a different value for this parameter, according to the actual temperature.
Fluids 06 00256 g001
Figure 2. The round-trip efficiency and energy density for the first 40 cycles of the PTES arrangement proposed by McTigue et al. [31]. (a) The round-trip efficiency obtained by the authors is compared with the one reported in the reference. (b) Energy density, defined as the ratio between the energy re-obtained during the discharging phase and the tanks’ total volume. The latter agrees with that reported in the reference since the energy density is equal to 300 MJ m−3 in design conditions.
Figure 2. The round-trip efficiency and energy density for the first 40 cycles of the PTES arrangement proposed by McTigue et al. [31]. (a) The round-trip efficiency obtained by the authors is compared with the one reported in the reference. (b) Energy density, defined as the ratio between the energy re-obtained during the discharging phase and the tanks’ total volume. The latter agrees with that reported in the reference since the energy density is equal to 300 MJ m−3 in design conditions.
Fluids 06 00256 g002
Figure 3. Energy and efficiency during the first 10 charge–discharge cycles for Desrues [21] and the present case. (a) In black are reported the authors’ stored energy (empty dots) during the charging process and the delivered energy (filled dots) during the discharging process. In red are shown the results of the reference. (b) Dots in black indicate the efficiency in the authors’ case, defined as the ratio between restored and stored energy. These trends show good accuracy in reproducing the process efficiency, compared with those described in the reference (red empty squares).
Figure 3. Energy and efficiency during the first 10 charge–discharge cycles for Desrues [21] and the present case. (a) In black are reported the authors’ stored energy (empty dots) during the charging process and the delivered energy (filled dots) during the discharging process. In red are shown the results of the reference. (b) Dots in black indicate the efficiency in the authors’ case, defined as the ratio between restored and stored energy. These trends show good accuracy in reproducing the process efficiency, compared with those described in the reference (red empty squares).
Fluids 06 00256 g003
Figure 4. Temperature profiles inside the tanks at the end of the two cycle phases, for a stabilized cycle. (a) Temperature trend inside the hot tank. (b) Temperature trend inside the cold reservoir.
Figure 4. Temperature profiles inside the tanks at the end of the two cycle phases, for a stabilized cycle. (a) Temperature trend inside the hot tank. (b) Temperature trend inside the cold reservoir.
Fluids 06 00256 g004
Figure 5. Schemes of the plant used as a test case: (a) Charging layout, (b) Discharge arrangement. The working fluid adopted in both charge and discharge is air. The system components are: a fan (FAN), an electric heater (EH), a heat exchangers (HX), an electric motor/generator (MG), a compressor, an expander and the storage tank.
Figure 5. Schemes of the plant used as a test case: (a) Charging layout, (b) Discharge arrangement. The working fluid adopted in both charge and discharge is air. The system components are: a fan (FAN), an electric heater (EH), a heat exchangers (HX), an electric motor/generator (MG), a compressor, an expander and the storage tank.
Fluids 06 00256 g005
Figure 6. Comparisons of the non-dimensional curves of temperature over time. Model m.3 assumes properties (cps, cpf , kf and μf) at different reference temperatures, while model TES-PD adopts properties computed at each layer with its temperature. The reference temperatures are marked with Tamb, Tmax, Tmean and Tvar. (a) First charge temperature trends over time, (b) first stable cycle charge temperature trends, (c) first stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the first stable cycle cycle’s charge (CHG) and discharge (DIS).
Figure 6. Comparisons of the non-dimensional curves of temperature over time. Model m.3 assumes properties (cps, cpf , kf and μf) at different reference temperatures, while model TES-PD adopts properties computed at each layer with its temperature. The reference temperatures are marked with Tamb, Tmax, Tmean and Tvar. (a) First charge temperature trends over time, (b) first stable cycle charge temperature trends, (c) first stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the first stable cycle cycle’s charge (CHG) and discharge (DIS).
Fluids 06 00256 g006
Figure 7. Comparisons of the model m.3 non-dimensional curves of temperature over time when the charging tolerance are different. The charging tolerance’s values were 10, 100 and 200, which are temperature differences between Tmax and the air temperature at the tank’s outlet. (a) First charge temperature trends over time, (b) first stable cycle charge temperature trends, (c) first stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the first stable cycle’s charge (CHG) and discharge (DIS).
Figure 7. Comparisons of the model m.3 non-dimensional curves of temperature over time when the charging tolerance are different. The charging tolerance’s values were 10, 100 and 200, which are temperature differences between Tmax and the air temperature at the tank’s outlet. (a) First charge temperature trends over time, (b) first stable cycle charge temperature trends, (c) first stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the first stable cycle’s charge (CHG) and discharge (DIS).
Fluids 06 00256 g007
Figure 8. Non-dimensional temperature trends over time. Comparison among the models in their original formulation. The charging tolerance was tol = 10, and the reference temperature was Tmean. (a) First charge temperature trends over time, (b) average cycle charge temperature trends, (c) average stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the average cycle’s during charge (CHG) and discharge (DIS).
Figure 8. Non-dimensional temperature trends over time. Comparison among the models in their original formulation. The charging tolerance was tol = 10, and the reference temperature was Tmean. (a) First charge temperature trends over time, (b) average cycle charge temperature trends, (c) average stable cycle discharge temperature trends and (d) temperature distributions inside the tank at the end of the average cycle’s during charge (CHG) and discharge (DIS).
Fluids 06 00256 g008
Table 1. Comparison among models when the fluid and solid properties ( c p , s , c p , f , k f and μ f ) are computed at a fixed reference temperature ( T a m b , T m a x and T m e a n ) and at a temperature, T v a r , reliant on the layer and time (only for m.2 and TES-PD). In order to compare the models, the heat loss coefficient U i and the solid effective thermal conductivity, k s , e f f , were both set equal to 0, while the charging tolerance was set equal to 10 K. The plant performance parameters are listed for the first cycle, the first stable cycle and average one.
Table 1. Comparison among models when the fluid and solid properties ( c p , s , c p , f , k f and μ f ) are computed at a fixed reference temperature ( T a m b , T m a x and T m e a n ) and at a temperature, T v a r , reliant on the layer and time (only for m.2 and TES-PD). In order to compare the models, the heat loss coefficient U i and the solid effective thermal conductivity, k s , e f f , were both set equal to 0, while the charging tolerance was set equal to 10 K. The plant performance parameters are listed for the first cycle, the first stable cycle and average one.
T amb T max T mean T var
m.1m.2m.3m.1m.2m.3m.1m.2m.3m.2TES-PD
First cycle
t C H G [h]22.2323.0222.0530.2530.2529.5230.2730.6229.7230.1329.42
t D I S [h]11.4510.8511.4816.5716.3716.9516.4216.0016.6716.3516.93
E C H G [MWh]252.44254.84251.34348.50347.48345.22347.72348.14345.14336.60334.95
E D I S [MWh]29.9428.3530.0243.3342.8044.3442.9341.8343.5942.7544.29
P G E N [MW]2.612.612.612.612.612.622.612.612.622.612.62
η r t [%]11.8611.1211.9512.4312.3212.8412.3512.0212.6312.7013.22
First stable cycle
t C H G [h]18.2018.3718.1225.1325.0324.8725.0825.0824.8725.0224.87
t D I S [h]11.4310.8311.4716.5516.3516.9316.4015.9816.6516.3316.92
E C H G [MWh]137.45134.00137.28193.38191.59194.80192.34189.51193.09190.76194.31
E D I S [MWh]29.8928.3029.9843.2842.7544.2942.8841.7843.5542.7144.25
P G E N [MW]2.612.612.612.612.612.622.612.612.622.612.62
η r t [%]21.7521.1221.8422.3822.3122.7422.3022.0522.5522.3922.77
Average cycle
t C H G [h]18.2018.3718.1225.1325.0324.8725.0825.0824.8725.0224.87
t D I S [h]11.4310.8311.4716.5516.3516.9316.4015.9816.6516.3316.92
E C H G [MWh]137.45134.00137.28193.38191.59194.80192.34189.51193.09190.76194.31
E D I S [MWh]29.8928.3029.9843.2842.7544.2942.8841.7843.5542.7144.25
P G E N [MW]2.612.612.612.612.612.622.612.612.622.612.62
η r t [%]21.7521.1221.8422.3822.3122.7422.3022.0522.5522.3922.77
Table 2. Comparison between models when assuming different charging tolerances. The charging tolerance, which is the temperature difference between the cycle’s maximum temperature and the air temperature at the tank’s outlet, was set equal to 10, 100 and 200 K.
Table 2. Comparison between models when assuming different charging tolerances. The charging tolerance, which is the temperature difference between the cycle’s maximum temperature and the air temperature at the tank’s outlet, was set equal to 10, 100 and 200 K.
tol 10 K tol 100 K tol 200 K
m.1m.2m.3m.1m.2m.3m.1m.2m.3
First cycle
t C H G [h]30.2730.6229.7226.1726.3025.9024.4824.5224.32
t D I S [h]16.4216.0016.6715.7815.3516.0714.9314.4815.28
E C H G [MWh]347.72348.14345.14331.31330.87329.86322.06321.07321.16
E D I S [MWh]42.9341.8343.5941.2540.1142.0039.0137.8239.94
P G E N [MW]2.612.612.622.612.612.612.612.612.61
η r t [%]12.3512.0212.6312.4512.1212.7312.1111.7812.43
First stable cycle
t C H G [h]25.0825.0824.8718.8718.5519.0014.6314.1714.95
t D I S [h]16.4015.9816.6515.1814.7215.5012.7712.2513.17
E C H G [MWh]192.34189.51193.09163.14159.01165.52132.87127.83136.43
E D I S [MWh]42.8841.7843.5539.6438.4140.4833.2731.9034.33
P G E N [MW]2.612.612.622.612.612.612.612.602.61
η r t [%]22.3022.0522.5524.3024.1624.4625.0424.9625.16
Average cycle
t C H G [h]25.0825.0824.8718.8718.5519.0014.6414.1814.96
t D I S [h]16.4015.9816.6515.1814.7215.5012.7712.2513.17
E C H G [MWh]192.34189.51193.09163.14158.95165.52132.96127.92136.55
E D I S [MWh]42.8841.7843.5539.6438.4340.4833.2831.9134.34
P G E N [MW]2.612.612.622.612.612.612.612.602.61
η r t [%]22.3022.0522.5524.3024.1724.4625.0324.9525.15
Table 3. Test of all models in their original formulation. In this case, the properties for m.1 and m.3 were assumed constant (evaluated at T m e a n ), while m.2 and TES-PD assumed variable properties for each layer, according to the local temperature. The heat loss coefficient U i was set equal to 0.7 Wm−2K−1, while k s , e f f was 1 Wm−1K−1. The charging tolerance was fixed equal to 10 K.
Table 3. Test of all models in their original formulation. In this case, the properties for m.1 and m.3 were assumed constant (evaluated at T m e a n ), while m.2 and TES-PD assumed variable properties for each layer, according to the local temperature. The heat loss coefficient U i was set equal to 0.7 Wm−2K−1, while k s , e f f was 1 Wm−1K−1. The charging tolerance was fixed equal to 10 K.
m.1m.2m.3TES-PD
First cycle
t C H G [h]32.2830.1329.5230.73
t D I S [h]15.2016.3516.7515.95
E C H G [MWh]355.77336.60344.41340.43
E D I S [MWh]39.1542.7543.8141.11
P G E N [MW]2.582.612.622.58
η r t [%]11.0012.7012.7212.08
Average cycle
t C H G [h]26.9825.0224.7826.15
t D I S [h]15.2016.3316.7315.93
E C H G [MWh]194.36190.76193.34195.35
E D I S [MWh]39.1542.7143.7741.07
P G E N [MW]2.582.612.622.58
η r t [%]20.1422.3922.6421.02
Table 4. Simulations for different choices of k s , e f f , adopting model TES-PD with variable properties and U i =  0.7 Wm−2K−1. The k s , e f f values are 0, 1 and 2 Wm−1K−1.
Table 4. Simulations for different choices of k s , e f f , adopting model TES-PD with variable properties and U i =  0.7 Wm−2K−1. The k s , e f f values are 0, 1 and 2 Wm−1K−1.
k s , eff = 0 k s , eff = 1 k s , eff = 2
First cycle
t C H G [h]32.2231.9831.73
t D I S [h]15.2215.2815.35
E C H G [MWh]346.06345.23344.35
E D I S [MWh]39.2339.4039.58
P G E N [MW]2.582.582.58
η r t [%]11.3411.4111.49
Average cycle
t C H G [h]26.8326.7226.60
t D I S [h]15.2015.2715.35
E C H G [MWh]193.12193.14193.28
E D I S [MWh]39.1939.3639.57
P G E N [MW]2.582.582.58
η r t [%]20.2920.3820.48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Benato, A.; De Vanna, F.; Gallo, E.; Stoppato, A.; Cavazzini, G. TES-PD: A Fast and Reliable Numerical Model to Predict the Performance of Thermal Reservoir for Electricity Energy Storage Units. Fluids 2021, 6, 256. https://doi.org/10.3390/fluids6070256

AMA Style

Benato A, De Vanna F, Gallo E, Stoppato A, Cavazzini G. TES-PD: A Fast and Reliable Numerical Model to Predict the Performance of Thermal Reservoir for Electricity Energy Storage Units. Fluids. 2021; 6(7):256. https://doi.org/10.3390/fluids6070256

Chicago/Turabian Style

Benato, Alberto, Francesco De Vanna, Ennio Gallo, Anna Stoppato, and Giovanna Cavazzini. 2021. "TES-PD: A Fast and Reliable Numerical Model to Predict the Performance of Thermal Reservoir for Electricity Energy Storage Units" Fluids 6, no. 7: 256. https://doi.org/10.3390/fluids6070256

APA Style

Benato, A., De Vanna, F., Gallo, E., Stoppato, A., & Cavazzini, G. (2021). TES-PD: A Fast and Reliable Numerical Model to Predict the Performance of Thermal Reservoir for Electricity Energy Storage Units. Fluids, 6(7), 256. https://doi.org/10.3390/fluids6070256

Article Metrics

Back to TopTop