Multiphase modeling of melting : solidification with high density variations using XFEM

Authors: Martin, Dave
Advisor: Fafard, MarioRobert, Jean-Loup
Abstract: The modelling of the cryolite bath, used in the smelting of aluminum, offers multiple challenges, particularly the presence of discontinuities in the solution and a difference in density between the solid and liquid phases. To over come these challenges, several novel elements were developed in this thesis. First of all, the phase change problem, commonly named the Stefan problem, was solved in two dimensions using the extended finite element method. A specially designed Lagrange multiplier formulation, using an enriched Lagrange multiplier solution, was implemented to impose the melting temperature on the interface. The interface velocity is determined by the jump in the heat flux across the interface and was calculated using the Lagrange multiplier values. Secondly, convection was included by solving the Stokes equations in the liquid phase using the extended finite element method as well. Thirdly, the density change between solid and liquid phases, usually neglected in the literature, was taken into account by the addition of a non-zero velocity boundary condition at the solid-liquid interface to maintain mass conservation in the system. Benchmark analytical and numerical problems were solved to validated the various components of the model and the coupled system of equations. The solutions to the numerical problems were compared to the solutions obtained using Comsol’s moving mesh algorithm. Theses comparisons show that the extended finite element method correctly solves the phase change problem with non-constant densities.
Document Type: Thèse de doctorat
Issue Date: 2016
Open Access Date: 24 April 2018
Grantor: Université Laval
Collection:Thèses et mémoires

Files in this item:
32754.pdf2.02 MBAdobe PDFView/Open
All documents in CorpusUL are protected by Copyright Act of Canada.