Abstract:We study cut finite element discretizations of a Darcy interface problem based on the mixed finite element pairs $\textbf{RT}_k\times Q_k$, $k\geq 0$. Here $Q_k$ is the space of discontinuous polynomial functions of degree less or equal to $k$ and $\textbf{RT}$ is the Raviart-Thomas space. We show that the standard ghost penalty stabilization, often added in the weak forms of cut finite element methods for stability and control of the condition number of the linear system matrix, destroys the divergence-free property of the considered element pairs. Therefore, we propose new stabilization terms for the pressure and show that we recover the optimal approximation of the divergence without losing control of the condition number of the linear system matrix. We prove that the method with the new stabilization term has pointwise divergence-free approximations of solenoidal velocity fields. We derive a priori error estimates for the proposed unfitted finite element discretization based on $\textbf{RT}_k\times Q_k$, $k\geq 0$. In addition, by decomposing the mesh into macro-elements and applying ghost penalty terms only on interior edges of macro-elements, stabilization is applied very restrictively and only where needed. Numerical experiments with element pairs $\textbf{RT}_0\times Q_0$, $\textbf{RT}_1\times Q_1$, and $\textbf{BDM}_1\times Q_0$ (where $\textbf{BDM}$ is the Brezzi-Douglas-Marini space) indicate that we have 1) optimal rates of convergence of the approximate velocity and pressure; 2) well-posed linear systems where the condition number of the system matrix scales as it does for fitted finite element discretizations; 3) optimal rates of convergence of the approximate divergence with pointwise divergence-free approximations of solenoidal velocity fields. All three properties hold independently of how the interface is positioned relative to the computational mesh.

On the convergence of 2D P4+ triangular and 3D P6+ tetrahedral divergence‐free finite elements

A nonconforming P3 and discontinuous P2 mixed finite element on tetrahedral grids

A nonconforming P2 and discontinuous P1 mixed finite element on tetrahedral grids

Nonconforming Finite Element Methods of Order Two and Order Three for the Stokes Flow in Three Dimensions.

A P 2 H-Div-Nonconforming-H-Curl Finite Element for the Stokes Equations on Triangular Meshes

A second-order convergent nonconforming polynomial stokes element on quadrilateral meshes

Divergence-free cut finite element methods for Stokes flow

A locally calculable P3‐pressure in a decoupled method for incompressible Stokes equations

Nonconforming finite element Stokes complexes in three dimensions

Piecewise Divergence-Free Nonconforming Virtual Elements for Stokes Problem in Any Dimensions

Third order Maximum-Principle-Satisfying Direct discontinuous Galerkin methods for time dependent convection diffusion equations on unstructured triangle mesh

A General Degree Divergence-Free Finite Element Method for the Two-Dimensional Stokes Problem on Smooth Domains

A divergence-free reconstruction of the nonconforming virtual element method for the Stokes problem

A locally calculable $P^3$-pressure in a decoupled method for incompressible Stokes equations

Mixed finite element projection methods for the unsteady Stokes equations

General degree divergence-free finite element methods for the Stokes problem on smooth domains

A stabilizer-free pressure-robust finite element method for the Stokes equations

-conforming HDG methods for the stress-velocity formulation of the Stokes equations and the Navier–Stokes equations

A divergence preserving cut finite element method for Darcy flow

Superconvergence of Discontinuous Galerkin and Local Discontinuous Galerkin Schemes for Linear Hyperbolic and Convection-Diffusion Equations in One Space Dimension