In addition, [PERSON] and [PERSON] [15] have shown that not only the absolute value but also the shape of the Darcy term function influences the phase change behavior. Recently, [PERSON] et al. [14] have stated that the influence of the Darcy constant reduces for finer meshes and for an isothermal phase change. Beside the switch-off techniques, the fixed grid methods depend on a method that takes into account the latent heat. Overviews of these so-called enthalpy methods can be found in literature [16; 17; 18]. [PERSON] et al. [19] analyzed five common enthalpy methods and concluded that for a broad range of parameters, four of the five tested methods gave almost identical results, only the applied apparent heat capacity method performed poorly. In addition, the effect of two different formulations of the convective term in the energy equation was studied. It was found that a formulation that only convects sensible enthalpy is more resilient compared to a formulation were all enthalpy is convected. For cement mortar with PCM, a purely conductive analysis was performed by [PERSON] et al. [20]. Three methods were tested and the best results were achieved with an enthalpy method that models the cement mortar PCM composite as a binary mixture. The tested apparent heat capacity method gave unphysical results. Other examples of tested numerical setups and settings are the tolerance for solving the energy equation [21] and different discretization schemes for the convective term of the energy equation [22]. In general, both the tolerance within the energy equation and the discretization of the convective term show a large impact; however, for the tolerance, it was not observed for all solvers tested. Moreover, [PERSON] et al. [13] compared the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) with the Pressure-Implicit with Splitting of Operators (PISO) scheme for solving the Navier-Stokes equation and the PRESure STaggering Option (PRESTO!) with the body-force-weighted approach for interpolating the pressure. There were no differences found in the results achieved with the SIMPLE or the PISO scheme, but the PRESTOI approach showed a better agreement with the experiments than the body-force-weighted approach. Lastly, it should be noted that independency studies for the mesh and the time step are commonly reported within papers on LTESS simulation, and a few authors compare two-dimensional (2D) with three-dimensional (3D) simulations showing that distinctive differences in the flow patterns may appear between both approaches [22; 23]. Analyses regarding the material properties mostly only focused on one or a few specific aspects. From analyses with simple analytical models, [PERSON] et al. [24] drew some basic conclusions regarding the effects of uncertain material data. The effect of differential scanning calorimetry (DSC) results that were achieved with different heating rates were investigated numerically [25; 26], and the influence of the width of the melting range on the numerical results was studied by several authors [9; 10; 15; 21; 27; 28]. Interestingly, it is reported that a larger melting range may either lead to a faster phase change [9; 21] or a slower one [9]. [PERSON] et al. [22] compared simulation results achieved with temperature dependent material properties to results that were obtained with constant material properties. It seemed that next to the density change between solid and liquid octadecane, the variation of the thermal conductivity of the liquid PCM has the largest influence on the melting rate. [PERSON] et al. [29] checked the impact of assuming identical material properties in the solid and the liquid state compared to differing values for some of the material properties. The authors claim that the model with different material properties in the solid and liquid phase has better potential to reduce the discrepancy between simulation and experiment than the model with constant material properties.