Numerical investigation of the influence of mushy zone parameter Amush on heat transfer characteristics in vertically and horizontally oriented thermal energy storage systems

The effect of the value used for the mushy zone parameter (A(mush)) on predicted heat transfer and melting characteristics of a phase change material (PCM) Lauric acid, in both vertical and horizontal enclosures was studied. There is a lack of clarity regarding which value of this parameter should be used for accurate simulations of phase change heat transfer, addressing this will aid in accurate simulation and design of systems for LHTES (Latent heat thermal energy storage). The numerical analysis undertaken used a commercial CFD code ANSYS FLUENT 18.2 and the enthalpy-porosity formulation. The range of mushy zone parameter used was from 10(5) to 10(7). The predicted locations of the melt front were compared to published experimental data available in the literature. The simulations provided quantitative information about the amount of energy stored and the melt fraction and providing improved understanding of the heat transfer process. Comparison between predictions using different values of A(mush), and experimental data showed that correct selection of the value of A(mus)h to be used in the momentum equations is an important parameter for accurate modelling of LHTES and has a significant influence on the solid-liquid interface shape and progression. The study reveals that increasing the value of A(mush) leads to a decrease in fluid velocity, decreasing convection and the rate of heat transfer, therefore, proper selection of the mushy zone parameter is necessary to accurately simulate LHTES systems and provide a better understanding of the phase change behaviour and heat transfer characteristics.