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.

Convergence Analysis of a DDFD Method for Flow Problems in Homogeneous Porous Media

Generalized Finite Difference Method for a Class of Multidimensional Space-Fractional Diffusion Equations

An Application of the Adomian's Decomposition Method to the Analysis of MHD Duct Flows

Convergence analysis of novel discontinuous Galerkin methods for a convection dominated problem

Block iterative methods for linear algebraic equation and domain decomposition method for incompressible viscous flow

Design and analysis of an exactly divergence-free hybridized discontinuous Galerkin method for incompressible flows on meshes with quadrilateral cells

Convergence of spectral discretization for the flow of diffeomorphisms

A two-grid characteristic block-centered finite difference method for Darcy-Forchheimer slightly compressible miscible displacement problem

A Reduced Basis Method for Darcy Flow Systems that Ensures Local Mass Conservation by Using Exact Discrete Complexes

Convergence analysis of numerical schemes for the Darcy-Forchheimer problem

Convergence of Dziuk's semidiscrete finite element method for mean curvature flow of closed surfaces with high-order finite elements

Exponentially-fitted Finite Elements for H(curl) and H(div) Convection-Diffusion Problems.

Structure-preserving schemes for drift-diffusion systems on general meshes: DDFV vs HFV

A divergence preserving cut finite element method for Darcy flow

A domain decomposed finite element method for solving Darcian velocity in heterogeneous porous media

A hybridizable discontinuous Galerkin method on unfitted meshes for single-phase Darcy flow in fractured porous media

An adaptive discontinuous Galerkin method for the Darcy system in fractured porous media

Hybrid GFD-RBF Method for Convection-Diffusion Problems

Superconvergence of a fully conservative finite difference method on non-uniform staggered grids for simulating wormhole propagation with the Darcy–Brinkman–Forchheimer framework

A hybrid-mixed finite element method for single-phase Darcy flow in fractured porous media