Dave Martin, Hicham Chaouki, Jean-Loup Robert, Mario Fafard, Donald Ziegler

Frontiers in Heat and Mass Transfer (FHMT) 7 - 31 (2016)

The two dimensional phase change problem was solved using the extended finite element method with a Lagrange formulation to apply the interface boundary condition. The Lagrange multiplier space is identical to the solution space and does not require stabilization. The solid-liquid interface velocity is determined by the jump in heat flux across the interface. Two methods to calculate the jump are used and compared. The first is based on an averaged temperature gradient near the interface. The second uses the Lagrange multiplier values to evaluate the jump. The Lagrange multiplier based approach was shown to be more robust and precise.