The melting of pure gallium in a rectangular cavity has been numerically investigated using the enthalpy-porosity approach for modeling combined convection-diffusion phase change. The major advantage of this technique is that it allows a fixed-grid solution of the coupled momentum and energy equations to be undertaken without resorting to variable transformations. In this work, a two-dimensional dynamic model is used and the influence of laminar natural-convection flow on the melting process is considered. Excellent agreement exists between the numerical predictions and experimental results available in the literature. The enthalpy-porosity approach has been found to converge rapidly, and is capable of producing accurate results for both the position and morphology of the melt front at different times with relatively modest computational requirements. These results may be taken to be a sound validation of this technique for modeling isothermal phase changes in metallurgical systems.
Bibliographical noteFunding Information:
Computations were carried out using an AT&T 3B2/400 microcomputer and a Cray 2/1 supercomputer. The AT&T 3B2/400 was pan of a computer donation by AT&T lnformation Systems, Inc., to the Civil and Mineral Engineering Department of the University of Minnesota. Time on the Cray 2/ 1 was provided by a resource grant from the Minnesota Supercomputer Institute. The support by both AT&T lnformation Systems. Inc.. and the Minnesota Supercomputer Institute is gratefully acknowledged.