A least‐squares front‐tracking finite element method analysis of phase change with natural convection

This paper focuses on the numerical modelling of phase-change processes with natural convection. In particular, two-dimensional solidification and melting problems are studied for pure metals using an energy preserving deforming finite element model. The transient Navier–Stokes equations for incompressible fluid flow are solved simultaneously with the transient heat flow equations and the Stefan condition. A least-squares variational finite element method formulation is implemented for both the heat flow and fluid flow equations. The Boussinesq approximation is used to generate the bulk fluid motion in the melt. The mesh motion and mesh generation schemes are performed dynamically using a transfinite mapping. The consistent penalty method is used for modelling incompressibility. The effect of natural convection on the solid/liquid interface motion, the solidification rate and the temperature gradients is found to be important. The proposed method does not possess some of the false diffusion problems associated with the standard Galerkin formulations and it is shown to produce accurate numerical solutions for convection dominated phase-change problems.

[1]  Keith Miller,et al.  Moving Finite Elements. I , 1981 .

[2]  Roland W. Lewis,et al.  Finite element solution of non‐linear heat conduction problems with special reference to phase change , 1974 .

[3]  Kevin O'Neill,et al.  Moving boundary‐moving mesh analysis of phase change using finite elements with transfinite mappings , 1986 .

[4]  J. P. Gupta,et al.  Thermal and fluid flow effects during solidification in a rectangular enclosure , 1982 .

[5]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[6]  V. Voller,et al.  A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems , 1987 .

[7]  Mario A. Storti,et al.  Making curved interfaces straight in phase‐change problems , 1987 .

[8]  Nicholas Zabaras,et al.  Moving and deforming finite‐element simulation of two‐dimensional Stefan problems , 1990 .

[9]  Thomas J. R. Hughes,et al.  Unconditionally stable algorithms for nonlinear heat conduction , 1977 .

[10]  R. L. Sani,et al.  Consistent vs. reduced integration penalty methods for incompressible media using several old and new elements , 1982 .

[11]  Boris Rubinsky,et al.  Numerical Computation Using Finite Elements for the Moving Interface in Heat Transfer Problems with Phase Transformation , 1983 .

[12]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[13]  Y. S. Touloukian.  Thermophysical properties of matter , 1970 .

[14]  Jonathan A. Dantzig,et al.  MODELLING LIQUID-SOLID PHASE CHANGES WITH MELT CONVECTION , 1989 .

[15]  Juan C. Heinrich,et al.  Petrov-Galerkin method for multidimensional time-dependent convective-diffusion equations , 1987 .

[16]  Ephraim M Sparrow,et al.  Analysis of melting in the presence of natural convection in the melt region , 1977 .

[17]  Vaughan R Voller,et al.  Finite difference solutions of solidification phase change problems: Transformed versus fixed grids , 1990 .

[18]  J. Heinrich,et al.  Heat Transfer Rates in Natural Convection at High Rayleigh Numbers in Rectangular Enclosures: a Numerical Study , 1982 .

[19]  A. Segal,et al.  Comparison of finite element techniques for solidification problems , 1986 .

[20]  Chie Gau,et al.  Melting and Solidification of a Pure Metal on a Vertical Wall , 1986 .

[21]  Vaughan R Voller,et al.  ENTHALPY-POROSITY TECHNIQUE FOR MODELING CONVECTION-DIFFUSION PHASE CHANGE: APPLICATION TO THE MELTING OF A PURE METAL , 1988 .

[22]  J. P. Gupta,et al.  TWO-DIMENSIONAL SOLIDIFICATION WITH NATURAL CONVECTION IN THE MELT AND CONVECTIVE AND RADIATIVE BOUNDARY CONDITIONS , 1981 .

[23]  Daniel R. Lynch,et al.  Continuously deforming finite elements for the solution of parabolic problems, with and without phase change , 1981 .

[24]  William G. Gray,et al.  Finite element simulation of flow in deforming regions , 1980 .

[25]  O. Zienkiewicz,et al.  Natural Convection in a Square Enclosure by a Finite-Element, Penalty Function Method Using Primitive Fluid Variables , 1978 .

[26]  V. Sastri,et al.  Efficient numerical method for two-dimensional phase change problems , 1984 .

[27]  Richard H. Gallagher,et al.  A general two‐dimensional, graphical finite element preprocessor utilizing discrete transfinite mappings , 1981 .

[28]  K. O'Neill,et al.  Computation of Porous Media Natural Convection Flow and Phase Change , 1984 .

[29]  M. Cross,et al.  Accurate solutions of moving boundary problems using the enthalpy method , 1981 .

[30]  Boris Rubinsky,et al.  A finite element method for the study of solidification processes in the presence of natural convection , 1986 .