Phase change problems with free convection: fixed grid numerical simulation

A numerical and experimental study of unsteady natural convection during freezing of water is presented. The mathematical model for the numerical simulations is based on the enthalpy-porosity method in vorticity-velocity formulation, equations are discretised on a fixed grid by means of a finite volume technique. A fully implicit method has been adopted for the mass and momentum equations. Experiments are performed for water in a differentially heated cube surrounded by air. The experimental data for natural convection with freezing in the cavity are collected to create a reference for comparison with numerical results. The method of simultaneous measurement of the flow and temperature fields using liquid crystal tracers is used. It allows us to collect transient data on the interface position, and the temperature and velocity fields. In order to improve the capability of the numerical method to predict experimental results, a conjugate heat transfer problem was solved, with finite thickness and internal heat conductivity of the non-isothermal walls. These results have been compared with the simulations obtained for the idealised case of perfectly adiabatic side walls, and with our experimental findings. Results obtained for the improved numerical model shown a very good agreement with the experimental data only for pure convection and initial time of freezing process. As time passes the discrepancies between numerical predictions and the experiment became more significant, suggesting a necessity for further improvements of the physical model used for freezing water.