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.

Pressure and convection robust bounds for continuous interior penalty divergence-free finite element methods for the incompressible Navier-Stokes equations

Pressure and convection robust bounds for continuous interior penalty divergence-free finite element methods for the incompressible Navier–Stokes equations

Pressure and convection robust Finite Elements for Magnetohydrodynamics

Pressure robust SUPG-stabilized finite elements for the unsteady Navier-Stokes equation

A pressure-robust mixed finite element method for the coupled Stokes–Darcy problem

Penalty method for the Navier-Stokes-Fourier system with Dirichlet boundary conditions: convergence and error estimates

A Posteriori Error Analysis for Pressure-Robust HDG Methods for the Stationary Incompressible Navier–Stokes Equations

A nonconforming pressure-robust finite element method for the Stokes equations on anisotropic meshes

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

Optimal Error Estimates of the Penalty Finite Element Method for the Unsteady Navier–Stokes Equations with Nonsmooth Initial Data

Analysis and computation of a pressure-robust method for the rotation form of the stationary incompressible Navier-Stokes equations by using high-order finite elements

Numerical analysis of a porous natural convection system with vorticity and viscous dissipation

Divergence-free cut finite element methods for Stokes flow

A divergence preserving cut finite element method for Darcy flow

Pressure-induced locking in mixed methods for time-dependent (Navier-)Stokes equations

A pressure-robust HHO method for the solution of the incompressible Navier-Stokes equations on general meshes

A divergence-free finite element method for the Stokes problem with boundary correction

Convergence and error estimates of a penalization finite volume method for the compressible Navier-Stokes system

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

A tangential and penalty-free finite element method for the surface Stokes problem