Modeling phase change materials (PCMs) has been a topic of research interest in the past, carried out experimentally and by means of computational fluid dynamics (CFD). The implemented solidification and melting (SM) model in Ansys Fluent-based on the enthalpy-porosity formulation is widely used in the literature. To the authors’ knowledge, few publications apply the apparent heat capacity (AHC) method in Ansys Fluent and even fewer have discussed both. The SM approach applies a linear relationship of the liquid fraction between solidus and liquidus temperature although it is known that the phase transition follows a non-linear behavior, which can be captured using the AHC method as a curve shape and location of the specific heat capacity containing information about the nature of phase transition behavior. Important factors in modeling are the temperature dependent thermophysical material properties density, viscosity, and thermal conductivity. They are often considered constant in the respective phase (solid or liquid) with a (linear) transition over the melting range. Temperature-dependent density is taken into account by using the Boussinesq approximation to model convective heat transfer. SM and AHC are compared to the analytical solution of the two-phase Stefan problem. As this does not include gravity and thus natural convection behavior, an additional comparison to two different PCMs, one from literature and a second data set gained in a new experiment is provided. The present work helps to evaluate the differences between the SM and AHC approach and to decide which is better suited for intended studies.

Monitoring of the state of charge of the thermal energy storage component in solar thermal systems for space heating and/or cooling in residential buildings is a key element from the overall system control strategy point of view. According to the literature, there is not a unique method for determining the state of charge of a thermal energy storage system that could generally be applied in any system. This contribution firstly provides a classification of the state-of-the-art of available techniques for the determination of the state of charge, and secondly, it presents an experimental analysis of different methods based on established sensor technologies, namely temperature, mass flow rates, and pressure measurements, tested using a lab-scale heat exchanger filled with a commercial phase change material for cooling applications. The results indicate that, depending on the expected accuracy and available instrumentation, each of the methods studied here can be used in the present application, the deviations between the methods generally being below 20%. This study concludes that a proper combination of two or more of these methods would be the ideal strategy to obtain a more reliable and accurate estimation of the state of charge of the latent heat thermal energy storage.

The area-proportional baseline method generates phase fraction–temperature curves from heat capacity data of phase change materials. The curves describe the continuous conversion from solid to liquid over an extended temperature range. They are consistent with the apparent heat capacity and enthalpy modeling approach for the numerical solution of heat transfer problems. However, the curves are non-smooth, discrete signals. They are affected by noise in the heat capacity data and should not be used as input to continuous simulation models. This contribution proposes an alternative method based on spline approximation for the generation of consistent and smooth phase fraction–temperature, apparent heat capacity–temperature and enthalpy–temperature curves. Applications are presented for two commercial paraffins from Rubitherm GmbH considering heat capacity data from Differential Scanning Calorimetry and 3-layer-calorimetry. Apparent heat capacity models are validated for melting experiments using a compact heat exchanger. The best fitting models and the most efficient numerical solutions are obtained for heat capacity data from 3-layer-calorimetry using the proposed spline approximation method. Because of these promising results, the method is applied to melting data of all 44 Rubitherm paraffins. The computer code of the corresponding phase transition models is provided in the Supplementary Information.