The Stefan problem regarding the formation of several liquid-solid interfaces produced by the oscillations of the ambient temperature around the melting point of a phase change material has been addressed by several authors. Numerical and semi-analytical methods have been used to find the thermal response of a phase change material under these type of boundary conditions. However, volume changes produced by the moving fronts and their effects on the thermal performance of phase change materials have not been addressed. In this work, volume changes are incorporated through an additional equation of motion for the thickness of the system. The thickness of the phase change material becomes a dynamic variable of motion by imposing total mass conservation. The modified equation of motion for each interface is obtained by coupling mass conservation with a local energy-mass balance at each front. The dynamics of liquid-solid interface configurations is analyzed in the transient and steady periodic regimes. Finite element and heat balance integral methods are used to verify the consistency of the solutions to the proposed model. The heat balance integral method is modified and adapted to find approximate solutions when two fronts collide, and the temperature profiles are not smooth. Volumetric corrections to the sensible and latent heat released (absorbed) are introduced during front formation, annihilation, and in the presence of two fronts. Finally, the thermal energy released by the interior surface is estimated through the proposed model and compared with the solutions obtained through models proposed by other authors.
Keywords: phase change material; several front formation; thermal performance.