Mathematical Analysis of Phase Change Thermal Energy Storage System and Effect of Stefan’s Number on TESS Performance

Mathematical Analysis of Phase Change Thermal Energy Storage System and Effect of Stefan’s Number on TESS Performance

Mayank SrivastavaM.K. Sinha 

Mechanical Engineering Department, National Institute of Technology, Jamshedpur, India

Corresponding Author Email:
22 August 2018
13 November 2018
31 December 2018
| Citation



Thermal energy storage system with phase change materials has attracted a great attention because of its important role in latent heat energy conservation. TESS gives a high thermal storage density with a wide range of temperature. This paper considers the analytical solution of outward melting/solidification of phase change materials in thermal energy storage system. Due to its non-linear behaviour, it is complicated to have exact solution of melting process and predict the behavior of interface movement, heat transfer rate. Heat balance integral method is applied to solve one dimensional outward melting problem in cylindrical geometry. Interface location, heat transfer rate and heat transfer with time is obtained for the geometry. Matlab code has been written to obtained the results.


conduction, encapsulation, HBIM, interface position, melting, phase change materials, Stefan’s number

1. Introduction
2. Mathematical Modelling
3. Numerical Results and Discussions
4. Conclusions

[1] Zalba B, Ma Marın J, Cabeza LF, Mehling H. (2003). Review on thermal energy storage with phase change: materials, heat transfer analysis and applications. Applied Thermal Engineering 23: 251–283.

[2] Farid MM, Khudhair AM, Razack SAK, Al-Hallaj S. (2004). A review on phase change energy storage: Materials and applications. Energy Conversion and Management 45: 1597–1615.

[3] Agyenim F, Hewitt N, Eames P, Smyth M. (2010). A review of materials, heat transfer and phase change problem formulation for latent heat thermal energy storage systems (LHTESS). Renewable and Sustainable Energy Reviews 14: 615–628.

[4] Bell GE. (1978). A refinement of the heat balance integral method applied to a melting problem. Int. J. Heat Mass Transfer 21: 1357-1362.

[5] Poots G. (1962). On the application of integral methods to the solution of the problem involving the solidification of liquids initially at fusion temperature. Int. J. Heat Mass Transfer 15: 525-531.

[6] Prasad A. (1979). Melting of solid bodies due to convective heating with the removal of melt. American Institute of Aeronautics and Astronautics, Inc. 16: 445-448.

[7] Caldwell J, Chiu CK. (2000). Numerical solution of one phase Stefan problems bt the heat balance integral method, Past-I – cylindrical and spherical geometries. Communications in Numerical Methods in Engineering 16: 569-583.

[8] Kumar A, Prasad A, Upadhaya SN. (1987). Spherical Phase change energy storage with constant temperature heat injection. Journal of Energy Resources Technology 109: 101-104.

[9] Ren HS. (2007). Application of the heat balance integral to an inverse Stefan problem. Int. J. of Thermal Sciences 46: 118-127.

[10] Prasad A, Singh SP. (1994). Conduction controlled phase change energy storage with radiative heat addition. Transactions of the ASME 116: 218-223.

[11] Sharma RK, Ganesan P, Sahu JN, Metselaar HSC, Mahila TMI. (2014). Numerical study for enhancement of solidification of phase change materials using trapezoidal cavity. Int. J. of Powder Technology 268: 38-47.

[12] Morgan K. (1981). A numerical analysis of freezing and melting with convection. Computer Methods in Applied Mechanics and Engineering 28: 275-284.

[13] Wood AS. (2001). A new look at the heat balance integral method. Applied Mathematical Modelling 25: 818-824.

[14] Myers TG, Mitchell SL, Muchatibaya G, Myers MY. (2007). Acubic heat balance integral method for one dimensional melting of a finite thickness layer. Int. J. Heat Mass Transfer 50: 5305-5317.

[15] Gupta RS. (1981). Dhirendra Kumar, Variable time step methods for one-dimensional Stefan problem with mixed boundary condition. Int. J. Heat Mass Transfer 24: 251-259.

[16] Caldwell J, Chan CCh. (2000). Spherical solidification by the enthalpy method and the heat balance integral method. Applied Mathematical Modelling 24: 45-53.

[17] Trp A. (2005). An Experimental and Numerical investigation of heat transfer during technical grade paraffin melting and solidification in a shell and tube latent thermal energy storage unit. Int. J. of Solar Energy 79: 648-660.

[18] Thomas DG, Sajith BC, Gopi S. (2016). Performance analysis of a latent heat thermal energy storage system for solar energy applications. Inernational Conferrence on Emerging Trends in Engineering, Science and Technology 24: 469-476.

[19] Lamberg P, Lehtiniemi R, Anna-Maria H. (2004). Numerical and experimental investigation of melting and freezing processes in phase change material storage. Int. J. of Thermal Sciences 43: 277-287.

[20] Savovic S, Caldwell J. (2009). Numerical solution of Stefan problem with time-dependent boundary conditions by variable space grid method. Int. J. of Thermal Science 13: 165-174.

[21] Du YP, Ding YL. (2016). The value of thermal radiation in assessing the charge/discharge rate of high‐grade thermal energy storage using encapsulated phase change materials (PCMs). International Journal of Energy Research 40(9): 1235-1244.

[22] Mantelli H, Braga WF. (2011). Temperature profiles for diffusion problem precise solutions using heat balance integral method. Journal of Thermophysics and Heat Transfer 25(3): 443-449.

[23] Srivastava M, Sinha MK. (2018). Computtional analysis of encapsulated thermal energy phase change storage system. International Journal of Mechanical Engineering and Technology 9(5): 662-668.

[24] Goodman TR. (1958). The heat balance integral and its application to problems involving a change of phase, Trans. ASME Journal of Heat Transfer 80: 335-342.

[25] Garg R, Thakur H, Tripathi B. (2018). Numerical simulation of two-dimensional fluid flow problem using truly meshfree method. Mathematical Modelling of Engineering Problems 5(4): 357-364.