Development and validation of a stabilized immersed boundary CFD model for freezing and melting with natural convection

The efficiency of the stabilization is first demonstrated by studying the Stefan problem. The two approaches to impose stasis are then compared using 2D test cases before they are both used to study melting in a rectangular (2D) and prismatic (3D) cavity. Significant differences are observed in the flow profiles and the solid-liquid interface position between the 2D and the 3D simulations.