Numerical Study of Phase Change of PCM in Spherical Cavities

Abstract

This work presents a numerical study of the phase change process of PCM (Phase Change Materials) stored in spherical cavities. The numerical model is two-dimensional and it is composed by the equations of conservation of mass, momentum, energy and volumetric fraction, which are modeled using the enthalpy-porosity technique. The computational mesh is tetrahedral, with refinements on regions that have large thermic and fluid dynamic gradients. The numeric model was validated with result from literature. It was studied the melting process of PCM RT35, RT 55 and RT 82 in spherical cavity with constant wall temperature. Four diameters of spheres D were used (40, 60, 80 and 100 mm) and three temperature differences ΔT (10, 20 and 30 oC) between the wall temperature and the melting temperature of the PCM. Liquid fraction results from the 36 cases studied are presented. It was observed that the time required to reach a certain liquid fraction increases with the diameter and reduces with the increment of ΔT, being possible to predict the fusion time by knowing the characteristic length of the sphere. The largest percentage reduction of the fusion time was obtained with ΔT = 10 oC – 20 oC for all the D considered. The shortest fusion time was obtained with the largest ΔT combined with the smallest D. It is possible to see the dependence of the liquid fraction results in relation with the PCM properties and the its independence in relation its melting temperature, since all the PCM studied presented equal fusion time for the same ΔT and D.