Modified enthalpy method for the simulation of melting and solidification

Abstract

Enthalpy method is commonly used in the simulation of melting and solidification owing to its ease of implementation. It however has a few shortcomings. When it is used to simulate melting/solidification on a coarse grid, the temperature time history of a point close to the interface shows waviness. While simulating melting with natural convection, in order to impose no-slip and impermeability boundary conditions, momentum sink terms are used with some arbitrary constants called mushy zone constants. The values of these are very large and have no physical basis. Further, the chosen values affect the predictions and hence have to be tuned for satisfactory comparison with experimental data. To overcome these deficiencies, a new cell splitting method under the framework of the enthalpy method has been proposed. This method does not produce waviness nor requires mushy zone constants for simulating melting with natural convection. The method is then demonstrated for a simple one-dimensional melting problem and the results are compared with analytical solutions. The method is then demonstrated to work in two-dimensions and comparisons are shown with analytical solutions for problems with planar and curvilinear interfaces. To further benchmark the present method, simulations are performed for melting in a rectangular cavity with natural convection in the liquid melt. The solid-liquid interface obtained is compared satisfactorily with the experimental results available in literature. © 2013 Indian Academy of Sciences.