Modelling convection in solidification processes using stabilized finite element techniques

Solidification of dendritic alloys is modelled using stabilized finite element techniques to study convection and macrosegregation driven by buoyancy and shrinkage. The adopted governing macroscopic conservation equations of momentum, energy and species transport are derived from their microscopic counterparts using the volume-averaging method. A single domain model is considered with a fixed numerical grid and without boundary conditions applied explicitly on the freezing front. The mushy zone is modelled here as a porous medium with either an isotropic or an anisotropic permeability. The stabilized finite-element scheme, previously developed by authors for modelling flows with phase change, is extended here to include effects of shrinkage, density changes and anisotropic permeability during solidification. The fluid flow scheme developed includes streamline-upwind/Petrov–Galerkin (SUPG), pressure stabilizing/Petrov–Galerkin, Darcy stabilizing/Petrov–Galerkin and other stabilizing terms arising from changes in density in the mushy zone. For the energy and species equations a classical SUPG-based finite element method is employed with minor modifications. The developed algorithms are first tested for a reference problem involving solidification of lead–tin alloy where the mushy zone is characterized by an isotropic permeability. Convergence studies are performed to validate the simulation results. Solidification of the same alloy in the absence of shrinkage is studied to observe differences in macrosegregation. Vertical solidification of a lead–tin alloy, where the mushy zone is characterized by an anisotropic permeability, is then simulated. The main aim here is to study convection and demonstrate formation of freckles and channels due to macrosegregation. The ability of stabilized finite element methods to model a wide variety of solidification problems with varying underlying phenomena in two and three dimensions is demonstrated through these examples. Copyright © 2005 John Wiley & Sons, Ltd.

[1]  S. Felicelli,et al.  Numerical model for dendritic solidification of binary alloys , 1993 .

[2]  Frank P. Incropera,et al.  Analysis of the effect of shrinkage on macrosegregation in alloy solidification , 1995 .

[3]  Qing-chun Li,et al.  NUMERICAL METHOD FOR SOLUTION OF STRONGLY COUPLED BINARY ALLOY SOLIDIFICATION PROBLEMS , 1991 .

[4]  S. Felicelli,et al.  Three-dimensional simulations of freckles in binary alloys , 1998 .

[5]  V. Voller,et al.  The modelling of heat, mass and solute transport in solidification systems , 1989 .

[6]  Qingchun Li,et al.  GRAVITY- AND SOLIDIFICATION-SHRINKAGE-INDUCED LIQUID FLOW IN A HORIZONTALLY SOLIDIFIED ALLOY INGOT , 1991 .

[7]  Frank P. Incropera,et al.  A continuum model for momentum, heat and species transport in binary solid-liquid phase change systems. II: Application to solidification in a rectangular cavity , 1987 .

[8]  Nicholas Zabaras,et al.  Variational multiscale stabilized FEM formulations for transport equations: stochastic advection- , 2004 .

[9]  Thomas J. R. Hughes,et al.  A stabilized mixed finite element method for Darcy flow , 2002 .

[10]  Nicholas Zabaras,et al.  A stabilized volume‐averaging finite element method for flow in porous media and binary alloy solidification processes , 2004 .

[11]  F. Incropera,et al.  A continuum model for momentum, heat and species transport in binary solid-liquid phase change systems—I. Model formulation , 1987 .

[12]  Tayfun E. Tezduyar,et al.  Finite element stabilization parameters computed from element matrices and vectors , 2000 .

[13]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[14]  J. Heinrich,et al.  Numerical simulation of incompressible flow driven by density variations during phase change , 1999 .

[15]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[16]  S. Felicelli,et al.  Finite element analysis of directional solidification of multicomponent alloys , 1998 .