An extended finite element method with algebraic constraints (XFEM-AC) for problems with weak discontinuities

Abstract

We present a new extended finite element method with algebraic constraints (XFEM-AC) for recovering weakly discontinuous solutions across internal element interfaces. If necessary, cut elements are further partitioned by a local secondary cut into body-fitting subelements. Each resulting subelement contributes an enrichment of the parent element. The enriched solutions are then tied using algebraic constraints, which enforce C0 continuity across both cuts. These constraints impose equivalence of the enriched and body-fitted finite element solutions, and are the key differentiating feature of the XFEM-AC. In so doing, a stable mixed formulation is obtained without having to explicitly construct a compatible Lagrange multiplier space and prove a formal inf-sup condition. Likewise, convergence of the XFEM-AC solution follows from its equivalence to the interface-fitted finite element solution. This relationship is further exploited to improve the numerical solution of the resulting XFEM-AC linear system. Examples are shown demonstrating the new approach for both steady-state and transient diffusion problems.