Numerical simulation of processes in the latent-heat thermal energy storage tank

Apr 23, 2021, 2:50 PM
1224 (Lomonosov st. 9)


Lomonosov st. 9

ITMO University
oral Energy Saving and Energy Efficiency Energy Saving and Energy Efficiency


Рафаэль Файзуллин (ITMO University)


In the near future, a significant increase in the number of air conditioning systems and an increase in their total electricity consumption are predicted [1]. The use of latent-heat thermal energy storage (LHTES) in these systems will reduce the peak load on equipment. This will enable to establish systems of lower power and reduce energy costs in the presence of a double tariff rate [2].

For the creation and effective operation of the LHTES, it is necessary to calculate the key parameters of the system, like the charging-discharging time and the power of the storage. Because of the high nonlinearity of thermophysical parameters of phase change materials (PCM), such a calculation presents a complex computational problem for CFD programs [3]. An alternative is to calculate using simplified mathematical models. This calculation significantly requires fewer computing resources, time and does not require the purchase of expensive licenses.

This article proposes a computer model for numerical calculations of heat fluxes in packed bed of capsules with PCM. A mathematical model of heat transfer in the capsule is based on finding a zero-dimensional solution of the Stefan problem with the influence of convective flows arising in the liquid phase. In Stefan's approach, the location of the mobile interface between the solid and liquid phases is unknown. To account for convective heat transfer in the liquid phase, the effective coefficient of thermal conductivity is calculated [4].

To describe the conditions of heat transfer between the coolant and the wall of the capsule, an experimental dependence was used [5]. The calculation is reduced to obtaining the temperature of the coolant after passing one layer of the packed bed. The resulting temperature represents an input parameter for calculating the subsequent layer. This operation is rehearsed until the calculation is made for all layers of the packed bed.

The numerical calculation was carried out in the mathematical software Scilab. The results of calculating the temperature of the coolant after passing through the LHTES according to the proposed model are in good agreement with the experimental data for LHTES with spherical capsules filled with paraffin.


This work was supported by the ITMO University under the project “Improving the efficiency of energy systems by using thermal energy storage batteries” (620150).


[1] IEA, World Energy Outlook 2019, November 2019, ISBN 978-92-64-97300-8
[2] Mingli Li, Qi Cao, Hong Pan, Xingyu Wang, Zhibin Lin. Effect of melting point
on thermodynamics of thin PCM reinforced residential frame walls in different climate zones//Applied Thermal Engineering, Volume 188, 2021, 116615
[3] Zilong Wang, Hua Zhang, Binlin Dou, Guanhua Zhang, Weidong Wu. Influence of inlet structure on thermal stratification in a heat storage tank with PCMs: CFD and experimental study//Applied Thermal Engineering, Volume 162, 2019, 114151
[4] Zakharova V.Y., Kuznetsov P.A., Baranenko A.V., Faizullin R.O. Calculation of cold accumulators with phase change materials//IOP Conference Series: Materials Science and Engineering, Volume 826, 2020, 012040
[5] Yankaus R.I., Survila V. Yu. Influence of wall permeability on heat transfer of spherical packing // VANT. Series: Atomic-hydrogen energy and technology, 1985, Volume 1 (20), p. 25-26

Position of speaker graduate student
Publication IOP Conference Series: Earth and Environmental Science
Affiliation of speaker ITMO University

Primary author


Dr Alexander Baranenko (ITMO University) Victoria Zakharova (ITMO University)

Presentation materials