$\phi$-FEM: a finite element method on domains defined by level-sets

TL;DR

This paper introduces $$-FEM, a finite element method for elliptic problems on level-set domains that avoids complex integration, achieves optimal convergence, and maintains good conditioning.

Contribution

The paper presents a novel fictitious domain finite element method that does not require non-standard integration, simplifying implementation and ensuring optimal convergence and conditioning.

Findings

Achieves optimal convergence in $H^1$ and $L^2$ norms.

Does not require non-standard numerical integration.

Maintains condition number comparable to standard FEM.

Abstract

We propose a new fictitious domain finite element method, well suited for elliptic problems posed in a domain given by a level-set function without requiring a mesh fitting the boundary. To impose the Dirichlet boundary conditions, we search the approximation to the solution as a product of a finite element function with the given level-set function, which also approximated by finite elements. Unlike other recent fictitious domain-type methods (XFEM, CutFEM), our approach does not need any non-standard numerical integration (on cut mesh elements or on the actual boundary). We consider the Poisson equation discretized with piecewise polynomial Lagrange finite elements of any order and prove the optimal convergence of our method in the $H^1$-norm. Moreover, the discrete problem is proven to be well conditioned, \textit{i.e.} the condition number of the associated finite element matrix is of the same order as that of a standard finite element method on a comparable conforming mesh. Numerical results confirm the optimal convergence in both $H^1$ and $L^2$ norms.