A Fully Implicit Enthalpy-Porosity Model for Phase-Change

This article proposes a new formulation for a phase change model based on the enthalpy-porosity idea. A general one-energy equation model (1EEM) is extended to deal with the melting and solidification of pure substances and alloys. Before melting and after solidification, solid material is seen as a porous media with low porosity and very small permeability. During phase change, thermal equilibrium in the mushy zone is assumed. Viscous and form drag in the volume-averaged momentum equation are reduced as the temperature rises above the melting point. In the energy equation, latent heat is treated implicitly in the accumulation term instead of explicitly as in most works in the literature. Liquid fraction for the entire field is updated after a new temperature field is calculated. Thermophysical properties are updated with the new liquid fraction field. Governing equations are discretized according to the control-volume method. Algebraic equation sets are relaxed with the Simple Method. Inner iterations make use of the Strong Implicit Procedure. Preliminary results indicate good agreement with the literature for pure substances.