Simulations were carried out using Engineering Equation Solver (EES) software, which enables advanced numerical calculations and optimization. The analysis showed that the economizer significantly increased the efficiency of the system, reducing the need for energy for pressurization and cooling production, and enabling a reduction in the area of the solar collector. Simulation results showed a 20% increase in the coefficient of performance of the system with the economizer.
Simulations showed that the cooling system operated more efficiently at lower heating medium temperatures, which is beneficial for longer absorption chiller operation cycles. Supporting cooling with solar energy was able to reduce the total cooling energy demand by 33% to 43%. In addition, better results were obtained at lower operating temperatures for the solar panels. During the heating season, the solar system was able to cover more than 90% of the heating demand, provided appropriate energy-saving strategies were used. The solar collectors operated efficiently at the lowest possible discharge temperature, and during the transitional season the building could operate without an external energy source, using excess heat to produce chilled water, which was stored for later.
The tank model did not allow any node to have a temperature higher than the node above it, as this could result in disruption of thermal stratification. In the event of such a situation, adiabatic mixing of the two nodes and equalization of temperatures occurred. The changing temperatures of the heating medium significantly affect the efficiency of the system and the absorption processes, which is crucial for optimizing the operation of thermal energy storage.
3.1. Using CFD to Simulate Thermal Stratification in Thermal Energy Storage Facilities
The growing interest in thermal energy storage technologies, especially in the context of Concentrated Solar Power (CSP) systems, underscores the importance of effective management and optimization of these systems. This technology is gaining popularity due to its ability to integrate with thermal energy storage (TES - Thermal Energy Storage) systems, which significantly improve the reliability and flexibility of energy supply.
According to the authors of the review [
24], thermal storage is a key component in solar energy systems. Such thermal stores create layers of working fluid due to temperature differences, where the warmest particles rise to the top and cooler particles sink to the bottom. Such an arrangement lowers the return temperature of the medium to the solar collector, which increases the efficiency of the system. Research indicates that thermal stratification in thermal storage contributes to better utilization of stored energy and facilitates its further storage. The stratification reduces the rate at which liquids are mixed, which reduces the time it takes to charge the storage compared to storages that lack distinct layers of liquids at different temperatures. A key factor affecting the effectiveness of stratification is the proportionality factor. According to Al-Marafi’s research, increasing the height-to-diameter ratio of a storage tank beyond a value of 4 brings no further benefit. Thermal stratification can also be improved by adding obstacles to the inside of the storage, reducing the flow rate of the medium, and moving the inlet and outlet to the edge of the tank [
25,
26]. Confirmation of one of the above statements can be found in
Figure 9, showing the temperature distribution in the thermal storage for two variants of the location of the inlet and outlet of the working medium.
From the figure above, it can be seen that there is more pronounced thermal stratification when the spigots are closer to the horizontal surfaces of the tank. In this case, layers with different temperatures are better separated, indicating less mixing of heat between layers. As the distance increases, the temperature gradient is more blurred, suggesting more intense mixing and thus greater thermal stratification, which can lead to less efficient use of energy in storage.
In the context of the thermal stratification in question, low-cost thermal stores that use thermocline technology are an interesting alternative. A thermocline is a thin layer that separates liquids of different temperatures, thus preserving thermal stratification - a process in which temperature differences in different layers of liquids are maintained, without excessive mixing of hot and cold water. However, poorly designed systems can lead to intense mixing of hot and cold liquids, which adversely affects temperature stratification and ultimately reduces the thermal performance of the storage facility. The thickness of the thermocline determines the level of mixing of liquids of different temperatures - a thin and stable thermocline reduces heat transfer between hot and cold water during tank loading and unloading, improving stratification.
Factors negatively affecting stratification in TES systems [
27] primarily include natural convection with the environment, mixing due to the kinetic force of the fluid entering the tank, thermal diffusion and thermal conductivity in the water, and the interaction of walls, pipes and other internal components. The shape factor of the tank (Aspect Ratio) also plays an important role. To improve the performance of thermal storage, special attention is paid to the type of heat transfer fluid (HTF - Heat Transfer Fluid) and the stratification process. Physical and geometrical parameters such as inlet velocity, supply temperature, tank aspect ratio and the presence of stratification baffles that optimize temperature distribution and minimize mixing of layers with different temperatures are also important factors. CFD simulations provide a better understanding and optimization of thermal stratification processes, which increases the efficiency and reliability of TES systems in CSP technology.
In response to the growing demand for efficient heat storage materials, the authors of the text [
28] presented a novel approach. They used air as the working medium and used a porous bed as the charge inside the tank. CFD simulations using the k-ε model investigated the effect of heat conduction coefficient and Capacity Ratio on the charging efficiency of the thermal storage. Various values of Reynolds number between 8300 and 50000, heat conduction coefficient (3.5-1062) and heat capacity of the cartridge (1483-7415) were analyzed. Two simulation scenarios included air turbulence in the porous structure and the thermal behavior of the system.
Figure 10 shows the distribution of turbulent kinetic energy in the inlet region of the tank. Hot air injection generates turbulent flow, especially in the zone where the medium contacts the porous material. Analysis of the turbulence kinetic energy (KV) for different values of the Reynolds number shows that the KV rate increases significantly as it increases, penetrating deeper into the porous medium. The highest KV values occur near the air inlet (r = 0 to r = 0.25 m), and as the Re number increases, the area of clear flow is filled with high KV values.
The results showed that lower thermal conductivity values reduced heat loss, resulting in higher temperatures inside the tank after longer charging times. Heat capacity had a significant effect on transition moments, and higher Reynolds numbers, which imply stronger turbulence in the flow, reduced the benefit of using materials with low thermal conductivity. To compensate for this, it was necessary to increase the thermal conductivity of the materials.
A paper [
29] described a two-dimensional numerical model for analyzing the process of energy storage and release in a thermocline of molten salt in a porous medium. The model of a cylindrical tank, considered the process of heat storage by introducing high-temperature molten salt into the tank through an upper valve while removing lower-temperature salt through a valve at the bottom. The process of recovering the stored heat was reversed: cold molten salt flowed into the bottom of the tank, while hot salt flowed out the top.
Simulations showed that during storage charging, the thermocline layer stabilizes and grows over time, although the rate of growth decreases. The authors noted that although stable thermocline layers are formed in a tank with porous material, its heat storage efficiency is lower compared to tanks with pure molten salt. Ultimately, choosing materials with high heat capacity and optimizing the structure of porous fillers can improve the performance and cost-effectiveness of the storage.
Another important case described in the article [
30], was the study of simultaneous charging and energy extraction from the thermal storage. The authors developed a model of a thermocline reservoir using molten salt as a medium and various filling materials: quartzite rock, slag stones and alumina ceramics. The effects of charging and discharging intensities on the thickness of the thermocline layer and stored power were analyzed. The tank was subjected to stable charging or intermittent charging while discharging continuously. At
Figure 6 shows the simulation results in the form of temperature distribution in tanks with different media and the process of thermocline formation.
Simulations showed that a tank filled with pure molten salt provided better thermal stratification compared to tanks with additional materials. The thickness of the thermocline layer increased with increasing flow rate, and the use of filler promoted uniform layer distribution along the radius of the tank. The greatest thickness of the thermocline was achieved with alumina ceramics. For the same flow rate, the temperature gradient was smaller for slagstone than for quartzite rock, suggesting that slagstone performs better as a heat storage material. In addition, continuous charging proved to be more efficient than intermittent charging at higher flow rates.
Much research has focused on the thermal processes occurring inside heat storage tanks. Researchers from the Department of Science and Engineering at Queensland University of Technology in Australia [
31] analyzed the theoretical performance of such tanks to identify key factors affecting their thermal efficiency in heating and cooling applications. Using CFD simulations, they studied five different fully insulated storage tank geometries with water as the heat transfer medium. The goal was to determine how the inlet velocity of the water, the ratio of tank dimensions and the temperature difference between the tank water and the charging water affect water mixing and thermocline layer formation. The tanks differed in their aspect ratio, that is, the ratio of height to diameter (TES 1: AR=3.8, TES 2: AR=2.75, TES 3: AR=2.38, TES 4: AR=1.33, TES 5: AR=1).
Figure 11 clearly shows the favorable effect of a high aspect ratio on the process of thermocline formation. In addition, it was noted that thermal stratification improved significantly as the temperature difference between the inlet medium and the initial temperature in the tank increased, along with lower inlet velocities. Observations showed that reducing the temperature difference between tank water and inlet water from 80°C to 10°C led to an increase in mixing by as much as 303%. Reducing the dimension ratio from 3.8 to 1.0 increased mixing by 143%. A similar trend was also observed when the inlet water velocity was increased. It is interesting to note that in the tank with the highest AR ratio (TES-1), it was possible to increase the inlet velocity without a clear effect on the shape of the thermocline after it was formed.
In an article [
32], a research team from the University of Nantes in France studied the problem of uneven flow distribution in a single-tank solar salt heat storage system. Their proposal to solve this problem was to use stratification baffles with openings to improve the uniformity of temperature distribution in the tank. Using CFD, they conducted 2D simulations of the temperature distribution and fluid velocity during storage loading and unloading.
At the beginning of the study, they focused on the charging process and conducted three key CFD simulations for different configurations of tank geometry. The first simulation involved a tank without a top baffle (step -1), the second with two uniform baffles with openings (step 0), and the third involved an optimized baffle geometry whose distances and widths were adjusted based on an iterative optimization process (
Figure 12). The optimization algorithm in MATLAB, implemented in six iterations (step 6), analyzed the data from the first two simulations and calculated changes in geometric parameters, which were then used to run new simulations in Ansys FLUENT. This yielded an optimized baffle configuration that provided the highest heat storage efficiency.
When analyzing the cyclic operation of the storage, including both loading and unloading, it was found that the geometry of the lower baffle has a significant impact on the efficiency of the unloading process. Therefore, using the optimized upper baffle from previous studies, an additional optimization of the lower baffle was carried out. The results showed that the optimized lower baffle improves the efficiency of both charging and discharging, leading to better flow management and less heat loss in the tank.
Various design parameters such as Reynolds number, aspect ratio, cone angle, number of holes in the baffle and baffle porosity were also analyzed. It was found that optimizing the two baffles significantly improves the efficiency of heat charging and discharging in thermal reservoirs, and their appropriate combination yields the best results. A higher Reynolds number negatively affected temperature stratification, suggesting the need to control the fluid flow velocity at the inlet. In addition, a positive cone angle (θ=+15°) improved control of warm liquid flow compared to a negative angle (θ=-15°). The larger number of holes, higher aspect ratio and lower porosity of the baffles improved stratification, leading to more efficient heat storage.
Another example of the use of computational fluid dynamics in studying the degree of thermal stratification of thermal storage tanks is published [
33]. A team of researchers from Shandong University in China developed three models of three-dimensional cylindrical storage tanks that differ in the number of stratification surfaces with an opening. Numerical simulations were used to study the effects of the number of obstructions and fluid inlet velocity on thermal stratification (
Figure 13).
It turned out that increasing the number of obstacles significantly improved stratification in the tank. Specifically, at an initial velocity of 0.3 m/s, the tank with two obstacles was more effective in reducing the mixing of hot and cold water compared to the others. Thermocline thickness measurements were 0.87 m, 0.23 m and 0.29 m for one, two and three baffles, respectively. As the inlet velocity increased, mixing increased, especially in tanks with one barrier. Tank number 3 showed the least sensitivity to changes in velocity, and the thinnest thermocline layer was observed at 0.20-0.25 m/s.
A study by a team of researchers from the CanmetENERGY institute in Canada [
34] aimed to investigate the effects of design and operational parameters during charging on the flow, thermal stratification and performance of the heat tank in solar energy systems. Two key design parameters were considered in the analysis: the aspect ratio (the ratio of the tank’s height to its diameter) and the position of the inlet and outlet ports. The operating parameters were mass flow rate, hot water temperature at the inlet, and initial water temperature. The study used CFD simulations of a three-dimensional insulated cylindrical tank with water as the medium.
The results showed that a high aspect ratio (AR = 3.5) provided the best thermal stratification. The higher the coefficient, the greater the temperature gradient between layers, which promoted better stratification. Decreasing the coefficient (AR = 2) led to worse stratification and lower temperatures. Changing the height of the tank at a constant diameter had no significant effect on stratification.
The location of the inlet/outlet port proved equally crucial. The best results were obtained with the inlet located near the top of the tank, where the thermal layer developed the fastest (
Figure 14).
Ports placed lower promoted more complex water mixing. The high hot water flow rate affected more intensive mixing, but stable thermal stratification was maintained. The optimal flow rate was 0.15 kg/s.
In the context of the inlet temperature, regardless of the difference between the initial temperature and the inlet water temperature, thermal stratification was stable and water mixing was minimal.
An article by researchers from the Department of Engineering and Energy in India [
35] analyzed the performance of single-phase tank heat storage systems using different inlet configurations and working materials. They analyzed water and solar salt (a mixture of NaNO3 and KNO3) in different tank variants - both with and without a distributor. The study was based on CFD simulations, taking into account the location of the hot medium inlet and the presence of the distributor (
Figure 15). The distributor in a heat storage system is designed to evenly distribute the hot medium in the tank, which translates into improved thermal stratification and increases the thermal efficiency of the system.
The results showed that the use of a distributor significantly improved thermal efficiency regardless of the working material. In the case of water, the distributor significantly reduced the thickness of the thermocline layer, which promoted better thermal stratification. In contrast, solar salt showed lower efficiency compared to water under the same conditions. The article emphasizes that the distributor improves thermal efficiency, especially in systems with water, while for systems with solar salt, further research is needed to improve the results.
Other researchers at the same institute in India investigated a thermal energy storage system using molten salts as the heat transfer medium and quartzite stones as the filler. In a paper [
36], they analyzed the effect of stone size on the performance of the thermal storage (TES) charging and discharging process using a three-dimensional model under laminar and turbulent flow conditions, verified by experimental results. The study included analysis of thermodynamic and geometric parameters, and the CFD-Taguchi method was used to optimize the system.
The results showed that charging and discharging efficiency decreased with time, and the optimal Reynolds number range (100-300) allowed better thermal stratification in the tank. Turbulent flow reduced system efficiency by 17.14% compared to laminar flow. Inlet temperature mainly affected heat loss through the outer walls, while smaller quartzite stones improved storage efficiency, although stones that were too fine or too large reduced porosity, which negatively affected charging and discharging processes.
The most significant parameters affecting performance were porosity and aspect ratio. The optimal settings, determined by the CFD-Taguchi method, were an aspect ratio of 0.25, Reynolds number of 10, porosity of 0.8 and fill size of 0.01 m.
3.2. Using CFD to Simulate Heat Loss in Thermal Energy Storage Systems
One of the main challenges in thermal energy storage systems negatively affecting their performance is heat loss. In order to reduce these losses and better understand the thermal processes occurring in thermal storage, computational fluid dynamics methods are also widely used.
Researchers at the University of Stuttgart have focused their research on insulating materials for thermal energy stores operating at high temperatures of up to 2000°C [
37]. Their research focuses on multilayer insulation, in which materials with better temperature stability are placed closer to the heat source.
Numerous simulations have been carried out to study heat loss and temperature propagation as a function of insulation materials, insulation layer thickness and gaseous atmosphere. The thermal storage (TES) was modeled as a truncated cone with a height of 150 mm and diameters of 100 mm and 150 mm, maintaining a constant temperature of 2300 K.
The authors of the publication found that increasing insulation thickness can lead to an increase in heat loss. As a result, they suggest that the use of vacuum in insulation can significantly reduce thermal conductivity, which can almost halve the heat loss rate with less insulation thickness. They expressed the need for further research into vacuum materials and the costs associated with maintaining the vacuum, as well as alternative insulation techniques such as multilayer foil insulation.
Analysis of all tested variants (air, argon and vacuum layer between insulation layers) showed that the heat loss coefficient decreases with increasing total insulation thickness (
Figure 16). However, increasing the thickness of the insulation leads to a decrease in effectiveness in reducing heat loss. This is due to the increasing external heat-transfer surface area with greater total insulation thickness. Heat loss reduction efficiency refers to how effectively an insulation reduces heat loss compared to its thickness. When insulation thickness increases, the amount of heat lost initially decreases noticeably. However, after a certain point, further increases in insulation thickness bring less and less benefit in reducing heat loss.
The article [
38] presents research on thermal energy storage systems in the context of insulation efficiency and heat loss in two different types of systems: rock cavern and above-ground (
Figure 17), which focused on comparing heat loss and heat transfer efficiency in these two systems.
The study used a TES storage model using loosely stacked stones that stores heat up to 685°C. Numerical simulations were carried out, which included analysis of heat transfer in the stone packs and evaluation of ambient heat loss. The simulations were carried out over 10 years to determine the effects of various parameters on heat storage efficiency.
The results showed that cave TES systems have lower heat loss (69.2% lower) and better long-term stability compared to above-ground systems. In the first five days, there was no significant difference in the rate of heat loss between the two technologies (185.2-243.5 kW for cave TES and 188.0-241.3 kW for aboveground TES). However, in the last five days, the rate of heat loss in the cave system was significantly lower (162.48-110.23 kW) than in the aboveground system (190.3-242.0 kW).
Interestingly, the authors simulated the predicted temperature distribution on the rocks surrounding the cave storage after one year and 10 years of operation (
Figure 18).
Cave systems have shown less sensitivity to changes in insulation performance in long-term use, making them a more efficient solution compared to above-ground systems.
A Turkish scientist from Amasya University [
39] studied the effect of using glass wool insulation for a cylindrical tank (truncated cone insulation) and a rectangular tank (pyramid-shaped insulation) (
Figure 19) on improving thermal stratification in hot water tanks. The results of the temperature distribution and velocity of the medium in the tank, obtained by simulation, are shown in
Figure 20.
The effects of three parameters were analyzed: tank size ratio, diameter ratio, and thickness of bottom and top insulation. According to the study, the highest loss of exergy efficiency was only 0.034% for truncated cone insulation and 0.02% for truncated pyramid insulation, indicating minimal loss of exergy when using these geometries. In addition, the results showed that lowering the ratios of tank dimensions, bottom and top insulation diameters, and insulation thickness significantly improves thermal stratification, leading to a significant increase in the temperature difference between the bottom and top of the tank. The maximum increase in thermal stratification was recorded at certain values of these parameters, confirming the effectiveness of truncated cone and pyramid-shaped insulation in improving the efficiency of hot water tanks.
Computational fluid dynamics (CFD) methods are commonly used to analyze heat loss in thermal energy storage (TES) systems. These simulations allow detailed modeling of heat transfer processes, taking into account the influence of various parameters such as tank geometry, insulation materials and ambient conditions. With CFD, it is possible to evaluate insulation efficiency, predict temperature distribution and optimize thermal storage designs, minimizing energy losses. This allows better design of TES systems with higher efficiency and long-term thermal stability.