Analysis of optimal preconditioners for CutFEM
arXiv:2202.09069 · doi:10.1002/nla.2486
Abstract
In this paper we consider a class of unfitted finite element methods for scalar elliptic problems. These so-called CutFEM methods use standard finite element spaces on a fixed unfitted triangulation combined with the Nitsche technique and a ghost penalty stabilization. As a model problem we consider the application of such a method to the Poisson interface problem. We introduce and analyze a new class of preconditioners that is based on a subspace decomposition approach. The unfitted finite element space is split into two subspaces, where one subspace is the standard finite element space associated to the background mesh and the second subspace is spanned by all cut basis functions corresponding to nodes on the cut elements. We will show that this splitting is stable, uniformly in the discretization parameter and in the location of the interface in the triangulation. Based on this we introduce an efficient preconditioner that is uniformly spectrally equivalent to the stiffness matrix. Using a similar splitting, it is shown that the same preconditioning approach can also be applied to a fictitious domain CutFEM discretization of the Poisson equation. Results of numerical experiments are included that illustrate optimality of such preconditioners for the Poisson interface problem and the Poisson fictitious domain problem.
arXiv commentnew version: added 2D experiments, including results for quadratic finite elements. arXiv admin note: text overlap with arXiv:2107.01182
References in corpus (4)
- High order unfitted finite element methods on level set domains using isoparametric mappings
- Preconditioning immersed isogeometric finite element methods with application to flow problems
- Scalable multigrid methods for immersed finite element methods and immersed isogeometric analysis
- Robust and scalable domain decomposition solvers for unfitted finite element methods