On the implementation of the eXtended Finite Element Method (XFEM) for interface problems
Provides a practical implementation guide for researchers using deal.II to solve interface problems with XFEM, but the approach is incremental.
The paper presents an implementation of the XFEM in the deal.II library for solving interface problems with unfitted meshes, demonstrating optimal convergence for strong and weak discontinuities.
The eXtended Finite Element Method (XFEM) is used to solve interface problems with an unfitted mesh. We present an implementation of the XFEM in the FEM-library deal.II. The main parts of the implementation are (i) the appropriate quadrature rule; (ii) the shape functions for the extended part of the finite element formulation; (iii) the boundary and interface conditions. We show how to handle the XFEM formulation providing a code that demonstrates the solution of two exemplary interface problems for a strong and a weak discontinuity respectively. In the weak discontinuity case, the loss of conformity due to the blending effect and its remedy are discussed. Furthermore, the optimal convergence of the presented unfitted method is numerically verified.