next up previous contents
Next: Summary of the algorithm Up: Discretization with the X-FEM Previous: The Time Projection   Contents

The X-FEM Approximation

The standard finite element approximation to the temperature field takes the form:

(48)

where are nodal shape functions and are the corresponding coefficients. The sum is taken over all nodes in the mesh. When the shape functions are piecewise linear, the above construction has the gradient

(49)

which can only be discontinuous across element boundaries. This is, essentially, one of the motivations for using a moving mesh technique, where the element and interface boundaries always coincide.

In the X-FEM, we extend the above using the construction

(50)

where is the enrichment function, are enriched degrees of freedom and is a subset of the nodes in the mesh. In the case of phase transformations, we take to consist of all nodes whose shape function supports contain the liquid-solid interface:

(51)

where is the location of the freezing front. This construction has a gradient

(52)

which can be discontinuous wherever is discontinuous. This represents the basic idea for this work: the motion of the interface is modeled through a change in the enrichment function .


next up previous contents
Next: Summary of the algorithm Up: Discretization with the X-FEM Previous: The Time Projection   Contents
Renaud Merle 2000-08-26