1. Introduction
Storage of energy is crucial to overcome the mismatch between variations in energy demand and renewable generation. Renewable sources of energy are intermittent, as they can be affected by weather conditions or other factors. Therefore, storing excess energy during times of low demand and using it later when the energy demand is higher is important to mitigate the high demand for these intermittent resources. One effective solution to overcome this challenge is to use energy storage systems [
1].
Pumped Thermal Energy Storage (PTES) is a promising technology for large-scale energy storage. Compared to other thermal energy storage methods, PTES offers high round-trip efficiency (RTE), high capacity, a lifespan of up to 30 years, a short response time [
2,
3,
4], and a fast start-up time [
5,
6]. PTES systems are not only environmentally friendly but also have a smaller carbon footprint compared with other novel large-scale commercialised energy storage technologies, such as Compressed Air Energy Storage (CAES) and Pumped Hydro Energy Storage Systems (PHES) [
7,
8].
During the charging process of PTES systems, electrical energy is converted into thermal energy. This is performed using a heat pump that moves heat from a low-temperature reservoir to a high-temperature reservoir. Later, during the discharging process, the thermal energy stored in the reservoirs is used to power a heat engine. The heat engine then converts the thermal energy back into electrical energy. The conventional layout of PTES is shown in
Figure 1 [
1].
There are various TES systems; nonetheless, the packed bed is one of the preferred designs with some limitations based on the type of application [
9,
10,
11,
12]. The packed bed storage system (PBSS) is a compact structure that offers a high heat storage capacity, a large surface area, an efficient energy transfer process, and very efficient storage.
The PBSS is an encouraging TES technology that helps enhance renewable energy generation systems and reduce energy storage costs [
13,
14]. The PBSS is an efficient technology that can be used for a wide range of temperatures, making it suitable for both high- and low-temperature applications [
14,
15].
The PBSS system consists of three main components: the mechanism of thermal energy transfer, the storage medium, and the containment system. The bed itself, see
Figure 2, is made up of storage materials, such as rocks, PCM, ores, or concrete, depending on factors such as economic considerations, intended use or application, and energy storage capacity [
16].
During both the charging and discharging phases, it is important to maintain contact between the heat transfer fluid (HTF) and the packing elements that circulate the bed voids to add or remove heat to the solid material. This increases the effective area of surfaces that are in contact, which in turn help to improve efficiency [
17]. Additionally, by adding improved stratification to the PBSS, the collector efficiency can be increased [
18].
High-temperature thermal energy storage is becoming increasingly important as a key component in large-scale applications. Packed bed storage represents an economically viable large-scale energy storage solution.
A wide range of applications, including chemical and drying processes, reactor designs, and thermal storage, have been studied in the literature concerning heat transfer within packed beds.
The first analysis of packed bed temperature transients dates back to 1926 and was undertaken by Anzelius [
19]. This later developed into the well-known Schumann model in 1929 [
20], which is essentially a one-dimensional, ‘two-phase’ model allowing for temperature differences between the gas and solid.
Various numerical schemes based on the Schumann approach have since been developed and validated by measurements [
21]. Schumann’s model [
20,
22,
23], which models fluid flow through a packed bed of solid spheres as a flow through porous media, is a frequently used basic two-phase model. Sanderson et al. [
24] and Regin et al. [
25] use this model to study the effect of various parameters, such as particle diameter on axial dispersion and pressure drop. The size of particles has a significant impact on the heat transfer area, which in turn affects the temperature of the fluid as it moves through the packed bed. The main factor that is affected by particle diameter is the permeability, which equals [
26] the following:
according to the Ergun model for non-Darcy flow,
k is the permeability,
dp is the particle diameter, and
ɛ is the porosity.
Based on the equation, if the particle diameter increases, the permeability will increase while the porosity remains constant—permeability is the property of porous materials that describes the ease with which a substance, such as a fluid or gas, can flow through the interconnected void spaces within the material. The explanation for this relation is that larger particles create more open structures between them, which creates greater void spaces. This facilitates easier fluid movement through the porous medium. Essentially, permeability is a measure of how easily fluids can move through a porous material.
The relationship between particle size and efficiency can be understood as follows: heat transfer within the packed bed depends heavily on the surface area of the particles. Smaller particles have a greater surface area compared to larger ones, which leads to increased interaction between the fluid and the porous material. As a result, heat is transferred more effectively from the fluid to the solid, resulting in a higher heat transfer. Therefore, smaller particles tend to have higher thermal efficiency in applications related to heat storage.
Recently, Navier–Stokes-based 3D optimisation studies were performed by Bataineh et al. [
27]. From the above literature, it is observed that there is no particular theory for the precise size selection of TES systems.
Singh et al. [
16] provided a brief review of analytical models and experimental work. However, experiments can be expensive and require a numerical simulation approach to predict the complex transient of storage systems. Recent reviews of analytical or numerical models are given by Gracia et al. [
28] and Elouali et al. [
29].
The storage of thermal energy can be divided into three categories, with each method designed for different purposes. Sensible Heat Storage (SHS) is a widely used technique for storing and releasing thermal energy by altering the temperature of a medium. Latent Heat Storage (LHS) involves the energy gained or lost during the transformation of a substance from solid to liquid, or vice versa. This method is commonly used in low-temperature applications, and the materials used in this category are known as Phase Change Materials (PCMs) Lastly, Thermochemical Energy Storage is a method that uses reversible chemical reactions to store and release thermal energy. This method has a high energy density and can be used in various applications.
Several packed bed designs are available in the literature, varying in geometry, storage medium, and heat transfer fluid (HTF). Typically, the HTF travels axially through the store, which is cylindrical. The filler is usually a solid substance such as ceramics, alumina, or crushed rock. Packed beds of this type are classified as SHS because the thermal change in the filler is responsible for storing energy. As an alternative, PCMs that are encapsulated can also be used to create LHS systems [
30], but this can increase the cost and complexity of the system.
HTFs can exist in the form of gases, such as air [
18,
31] or argon [
28], or liquids, such as thermal oils [
29] or molten salts [
25]. A packed bed storage facility operates in two different ways, depending on whether it is being charged or discharged. When heat is added to the bed, the flow usually moves in one direction, and when heat is removed, the flow moves in the opposite direction.
This paper applies numerical modelling using COMSOL Multiphysics software to investigate the performance of a thermal energy storage-packed bed system. A parametric study is conducted to determine the optimum design parameter values, such as the particle diameter of the porous media, void fraction, storage material, and aspect ratio, in order to enhance the hot storage packed bed’s overall performance.
3. Results and Discussion
This paper presents the development of a hot storage-packed bed model derived from White’s model. White’s model used argon as an HTF, and Fe3O4 was used as the packed storage material, with a particle diameter of 4 mm, a mass flow rate of 13.7 kg/s, and a void fraction of 0.4. The model has a height and width of 4.62 m.
Changing the heat transfer fluid to air is the first step in the model. Next, the particle diameter is examined to determine the optimal value, with a range chosen based on a literature review while maintaining other parameters unchanged. Following the selection of the optimal particle diameter, this paper discusses selecting the optimum porosity and storage material. The optimum particle diameter is then used alongside different packed storage materials and material porosities while maintaining a constant mass flow rate and aspect ratio. The aspect ratio is then tested under constant volume using the optimum particle diameter, material porosity, and storage material. The optimisation results are presented in the sections that follow.
3.1. Optimum Particle Diameters
In this paper, solid particle diameters of 4, 20, 40, 80, 100, 120, and 160 mm were used to investigate the effect of particle size on the performance of the hot storage tank during the charging process. The capacity factor and the total energy stored calculation methods results are shown in
Figure 7 and
Figure 8, respectively.
According to
Figure 7, the capacity factor increases as the particle diameter and the amount of energy stored over the eight hours decrease.
Figure 8 shows that the total energy stored increases as the particle diameter size decreases.
The temperature profile of different particle diameters along the packed bed after eight hours of charging is presented in
Figure 9.
Based on
Figure 9, the temperature of particles with the smallest diameter increases more rapidly than the temperature of particles with the largest diameter. This is due to the fact that the smallest particle diameter has the largest surface area and the greatest number of particles, which enhances the heat transfer performance of the smaller particle sizes.
Accordingly,
Figure 10a illustrates the heat transfer efficiency along the
Z-axis (packed bed height) after 8 h of charging time; similar results have been provided after 2, 4, and 6 h of charging. It is clear that the heat transfer efficiency profile pattern follows the temperature profile in
Figure 9 and agrees with the conclusions obtained in
Figure 7 and
Figure 8, which show that the smallest particle diameter is associated with the highest performance.
Figure 10b shows the heat transfer efficiency of each particle size after 8 h of charging.
Based on these results, a conclusion is drawn that the smallest particle diameter provides a high energy storage and a more efficient process, which is consistent with the literature review viewpoint.
3.2. Material Selection and Porosity
Storage materials play a crucial role in converting heat energy, as their physical and thermal properties largely determine their effectiveness. To be considered suitable for use, storage materials must possess a high energy density, high heat capacity, high thermal conductivity, long-term cycle stability, good mechanical stability, and low carbon footprint, and also be composed of sustainable, non-toxic materials that are compatible with the intended operating conditions.
In this study, five solid Sensible Heat Storage materials that are suitable for TES systems operating at temperatures above 500 °C have been selected. The thermophysical properties of the selected solid packing materials are shown in
Table 6 [
6,
38,
39].
Porosity refers to the measure of void spaces in a material, and it affects the way fluids flow through it. In fluid dynamics, the impact of porosity on pressure drop is significant, especially in the flow through porous media. The Darcy–Weisbach equation is commonly used to describe pressure drop in porous media. This study selected six different porosity values from 0.1 to 0.6 to find the optimum value for the hot storage design.
Figure 11 presents the relation between the hot tank capacity factor against the material porosity for the different materials considered.
Figure 11 shows the charging efficiency among five solid sensible storage materials after 8 h of charging. The capacity factor of Magnesia is comparatively low. For instance, at a porosity level of 0.2, the capacity factor of Magnesia is approximately 7%, which signifies it still has the capacity to store up to 14 times as much energy. This also means that if the full capacity is not required, the container could be made smaller. When compared to other materials, like NaCl, reinforced concrete, and silica, Magnesia has a higher charging factor of around 13–14% after 8 h of charging. Although cast iron performs slightly better than Magnesia, it shows poorer performance in
Figure 12 with respect to the total energy stored.
Figure 12 shows that Magnesia performs well. It stores slightly more than NaCl and slightly lower than reinforced concrete and silica, and they perform much better than cast iron. Although cast iron was seen to perform slightly better than Magnesia in
Figure 11, and its poor performance in terms of energy stores, as seen in
Figure 12, suggests that it is not an ideal material for the backed bed.
Accordingly, Magnesia performs the second best according to both efficiency and energy stored, and is only slightly behind the silica and reinforced concrete in
Figure 12 and the cast iron in
Figure 11. In terms of the porosity, it is clear that the total energy stored increases by decreasing the porosity, suggesting that both 0.1 and 0.2 would be appropriate values. However,
Figure 13 shows that the pressure drop through the storage container starts to increase at the lower porosity values, so a porosity of 0.2 has been selected.
After analysing
Figure 11,
Figure 12 and
Figure 13, it has been concluded that porosity affects the three primary design factors: pressure drop, heat transfer, and energy storage.
3.3. Optimum Aspect Ratio
When designing thermal energy storage tanks for high-temperature applications, it is important to consider the aspect ratio. This ratio affects the efficiency of heat transfer within the storage medium. By optimising the aspect ratio, heat can be distributed effectively throughout the storage material, which is crucial for charging and discharging processes in high-temperature applications.
The aspect ratio in tanks refers to the ratio of the tank’s height to its diameter or width. It is a geometric parameter that characterises the shape of the storage tank and is denoted by an aspect ratio = H/D.
Another numerical investigation has been conducted to explore the influence of studying different aspect ratios on the performance of a hot storage tank.
Four different aspect ratios have been tested and conducted in this paper, where the volume of the container was kept fixed, as shown in
Table 7.
The aspect ratio plays a critical role in determining the temperature distribution within the tank. A higher aspect ratio results in a more even temperature distribution, while a lower aspect ratio may cause significant temperature variations from the top to the bottom of the tank. This is because a tank with a higher aspect ratio experiences a smaller temperature gradient along its height, leading to a more homogeneous temperature distribution.
Table 8 shows the effect of changing the aspect ratio on the heat transfer, total energy stored, and system efficiency.
According to the results shown in
Table 8, the highest aspect ratio shows a high-capacity factor, high total energy stored, and high heat transfer efficiency with a slight difference between the amount of energy stored and the system efficiency. The results in
Table 8 show that the highest aspect ratio shows a high-capacity factor, high total energy stored, and high heat transfer efficiency. However, the differences are relatively small, and other features should also be considered, such as the increased pressure drop across longer storage devices.
3.4. Temperature Profile
Figure 14 shows the temperature distribution within the storage tank with an aspect ratio of 1 and with the parameters selected based on the optimised parameters found in the preceding analysis. Magnesia is used as the packed bed material with a particle diameter of 4 mm and a porosity of 0.2. The figure shows the evolution of the temperature distribution within the storage system at selected times over the first four days of charging and indicates how the temperature profile within the packed bed evolves over the charging process.