The extended finite element method (XFEM) for solidification problems

An enriched finite element method for the multi-dimensional Stefan problems is presented. In this method the standard finite element basis is enriched with a discontinuity in the derivative of the temperature normal to the interface. The approximation can then represent the phase interface and the associated discontinuity in the temperature gradient within an element. The phase interface can be evolved without re-meshing or the use of artificial heat capacity techniques. The interface is described by a level set function that is updated by a stabilized finite element scheme. Several examples are solved by the proposed method to demonstrate the accuracy and robustness of the method. Copyright (C) 2001 John Wiley Sons, Ltd.