A uniformly well-conditioned, unfitted Nitsche method for interface problems

A finite element method for elliptic partial differential equations that allows for discontinuities along an interface not aligned with the mesh is presented. The solution on each side of the interface is separately expanded in standard continuous, piecewise-linear functions, and jump conditions at the interface are weakly enforced using a variant of Nitsche's method. In our method, the solutions on each side of the interface are extended to the entire domain which results in a fixed number of unknowns independent of the location of the interface. A stabilization procedure is included to ensure well-defined extensions. We prove that the method provides optimal convergence order in the energy and the L (2) norms and a condition number of the system matrix that is independent of the position of the interface relative to the mesh. Numerical experiments confirm the theoretical results and demonstrate optimal convergence order also for the pointwise errors.