A divergence preserving cut finite element method for Darcy flow

Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Numerical methods for partial differential equations, boundary value problems (65N99)

Abstract: We study cut finite element discretizations of a Darcy interface problem based on the mixed finite element pairs

e x t b f {R T}_{0} i m e s Q_{0}

,

e x t b f {B D M}_{1} i m e s Q_{0}

, and

e x t b f {R T}_{1} i m e s Q_{1}

. Here

Q_{k}

is the space of discontinuous polynomial functions of degree k,

e x t b f {R T}_{k}

is the Raviart-Thomas space, and

e x t b f {B D M}_{k}

is the Brezzi-Douglas-Marini 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 resulting linear system matrix, destroys the divergence-free property of the considered element pairs. Therefore, we propose two corrections to the standard stabilization strategy: using macro-elements and new stabilization terms for the pressure. By decomposing the computational mesh into macro-elements and applying ghost penalty terms only on interior edges of macro-elements, stabilization is active only where needed. By modifying the standard stabilization terms for the pressure we recover the optimal approximation of the divergence without losing control of the condition number of the linear system matrix. We derive a priori error estimates for the proposed unfitted finite element discretization based on

e x t b f {R T}_{k} i m e s Q_{k}

,

k g e q 0

. Numerical experiments indicate that with the new method 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.

This page was built for publication: A divergence preserving cut finite element method for Darcy flow

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6507842)