3D macrosegregation simulation with anisotropic remeshing

The article presents a three-dimensional coupled numerical solution of momentum, mass, energy and solute conservation equations, for binary alloy solidification. The interdendritic flow in the mushy zone is assumed to obey the Darcy's law. Microsegregation is governed by the lever rule, assuming local equilibrium at phase interfaces. The resulting energy and solute advection–diffusion equations are solved using the Streamline-Upwind/Petrov–Galerkin (SUPG) finite element method. A SUPG-PSPG velocity-pressure formulation is applied for the momentum equation. The full algorithm was implemented in the 3D code THERCAST, together with an anisotropic remeshing method. Two applications have been considered: a small ingot of Pb-48wt%Sn alloy and a large steel ingot. The numerical results of these two cases are presented with the evolution of temperature, liquid velocity, and solute concentration fields during solidification. To cite this article: S. Gouttebroze et al., C. R. Mecanique 335 (2007).