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.

A posteriori error estimate for the non-stationary concentration equation coupled with the Darcy system discretized by the Raviart–Thomas finite element

A posteriori error analysis of a semi‐augmented finite element method for double‐diffusive natural convection in porous media

A posteriori error estimates for the Brinkman-Darcy-Forchheimer problem

A posteriori error analysis of a Banach spaces-based fully mixed FEM for double-diffusive convection in a fluid-saturated porous medium

A new optimal error analysis of a mixed finite element method for advection–diffusion–reaction Brinkman flow

A posteriori error estimates for a dual finite element method for singularly perturbed reaction–diffusion problems

A posteriori error analysis and adaptive non-intrusive numerical schemes for systems of random conservation laws

A posteriori error analysis of a mixed FEM for the coupled Brinkman-Forchheimer/Darcy problem

Analysis of a semi-augmented mixed finite element method for double-diffusive natural convection in porous media

An augmented mixed FEM for the convective Brinkman-Forchheimer problem: a priori and a posteriori error analysis

Finite Element Approximation of Data-Driven Problems in Conductivity

An adaptive stabilized finite element method for the Stokes–Darcy coupled problem

Newton iterative method based on finite element discretization for the stationary Darcy–Brinkman equations

A divergence preserving cut finite element method for Darcy flow

A Posteriori Error Bounds and Global Error Control for Approximation of Ordinary Differential Equations

A posteriori error estimates for stabilised mixed finite element methods for a nonlinear elliptic problem

A-posteriori error estimates for systems of hyperbolic conservation laws

Guaranteed a-posteriori error estimation for finite element solutions of nonstationary heat conduction problems based on their elliptic reconstructions

A posteriori error estimates for nonconforming discretizations of singularly perturbed biharmonic operators

A posteriori error analysis of a space-time hybridizable discontinuous Galerkin method for the advection-diffusion problem