CFD simulation of solid/liquid phase change in commercial PCMs using the slPCMlib library

Modelling the solidification and melting behavior of phase change materials (PCMs) can present several challenges, including non-isothermal phase change, multi-step transitions, complex enthalpy-temperature relations and convective heat transfer. To tackle these challenges, this study proposes a novel numerical method based on real experimental data for modelling complex phase change behavior of commercial PCMs in Ansys Fluent CFD software. The phase transition of PCMs is simulated with the apparent heat capacity (AHC) method, which is implemented as a User-Defined Function (UDF). Heat convection is simulated with the Boussinesq approximation. The method is validated by comparing to the Ansys Fluent built-in Solidification and Melting (S&M) model for the case of piecewise constant heat capacity. The results show that AHC method provides excellent agreement with S&M model for a piecewise constant heat capacity model. In addition, the AHC method solves convergence problems even for narrow phase change temperature ranges when a smooth heat capacity function is assumed. Finally, the AHC model is also tested on a set of commercial PCMs exhibiting complex phase change behavior. Temperature dependent specific heat capacities are determined from heat capacity data provided by PCM manufacturing companies and are imported as cubic Hermite spline functions from the solid-liquid phase change material library slPCMlib available at https://slpcmlib.ait.ac.at/. From the presented case study and tested PCMs it can be concluded that the AHC method is numerically robust when used with sufficiently smooth heat capacity functions, and it can be recommended for the analysis of heat transfer with PCMs showing a complex phase change behavior.