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.

DarcyLite: A Matlab Toolbox for Darcy Flow Computation.

A Fully Coupled Immersed Finite Element Method for Fluid Structure Interaction via the Deal.II Library

Deal.Ii Implementation Of A Weak Galerkin Finite Element Solver For Darcy Flow

A Parallel Fictitious Domain Method for the Interface-Resolved Simulation of Particle-Laden Flows and Its Application to the Turbulent Channel Flow

The Lowest-Order Weak Galerkin Finite Element Method for the Darcy Equation on Quadrilateral and Hybrid Meshes

Meshless Analysis of Darcy Flow with a Variational Multiscale Interpolating Element-Free Galerkin Method

Two-Dimensional Numerical Analysis of Non-Darcy Flow Using the Lattice Boltzmann Method: Pore-Scale Heterogeneous Effects

Multiphase Flow Modelling in Multiscale Porous Media: An Open-Sourced Micro-Continuum Approach

An efficient flux‐variable approximation scheme for Darcy's flow

A Discontinuous Galerkin Method for the Brinkman–Darcy-transport Problem

A Matlab Tutorial for Diffusion-Convection-Reaction Equations using DGFEM

Multiscale finite element method for Stokes-Darcy model

A Weak Galerkin Finite Element Method Based on H ( Div ) Virtual Element for Darcy Flow on Polytopal Meshes

Finite-volume Method with Lattice Boltzmann Flux Scheme for Incompressible Porous Media Flow at the Representative-Elementary-volume Scale

A divergence preserving cut finite element method for Darcy flow

Upscaled Lattice Boltzmann Method for Simulations of Flows in Heterogeneous Porous Media

A Mixed Discontinuous Galerkin Method for an Unsteady Incompressible Darcy Equation

On development of lattice-Boltzmann method based multi phase flow solver

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

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

DOLFIN: Automated Finite Element Computing