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

摘要

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 C~0 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.