Open Access

ARTICLE

A XFEM PHASE CHANGE MODEL WITH CONVECTION

Dave Martin^a,b,† , Hicham Chaouki^a,b, Jean-Loup Robert^a , Donald Ziegler^c , Mario Fafard^a,b

a Department of Civil and Water Engineering, Laval University, Quebec, QC, G1V 0A6, Canada
b NSERC/Alcoa Industrial Research Chair MACE3 and Aluminium Research Centre - REGAL, Laval University, Quebec, QC, G1V 0A6, Canada
c Alcoa Primary Metals, Alcoa Technical Center, 100 Technical Drive, Alcoa Center, PA, 15069-0001, USA

* Corresponding Author: Email:

Frontiers in Heat and Mass Transfer 2018, 10, 1-11. https://doi.org/10.5098/hmt.10.18

Abstract

A model for two dimensional solidification problems including convection was developed by coupling the Stefan problem with the Stokes problem. The extended finite element method (XFEM) was used to capture the strong discontinuity in velocity and pressure as well as the jump in heat flux at the phase change interface. The melting temperature and no-slip condition were imposed on the interface using a Lagrange multiplier and the penalization method, respectively. The resulting formulations were then coupled using a fixed point iteration algorithm. The model was able to reproduce the benchmark simulations while maintaining a sharp phase change interface.

---

Keywords

---