A hybrid extended finite element/level set method for modeling phase transformations

A hybrid numerical method for modelling the evolution of sharp phase interfaces on xed grids is presented. We focus attention on two-dimensional solidi cation problems, where the temperature eld evolves according to classical heat conduction in two subdomains separated by a moving freezing front. The enrichment strategies of the eXtended Finite Element Method (X-FEM) are employed to represent the jump in the temperature gradient that governs the velocity of the phase boundary. A new approach with the X-FEM is suggested for this class of problems whereby the partition of unity is constructed with C1( ) polynomials and enriched with a C0( ) function. This approach leads to jumps in temperature gradient occurring only at the phase boundary, and is shown to signi cantly improve estimates for the front velocity. Temporal derivatives of the temperature eld in the vicinity of the phase front are obtained with a projection that employs discontinuous enrichment. In conjunction with a ner nite di erence grid, the Level Set method is used to represent the evolution of the phase interface. An iterative procedure is adopted to satisfy the constraints on the temperature eld on the phase boundary. The robustness and utility of the method is demonstrated with several benchmark problems of phase transformation