A divergence preserving cut finite element method for Darcy flow

05/24/2022
by   Thomas Frachon, et al.
0

We study cut finite element discretizations of a Darcy interface problem based on the mixed finite element pairs RT_0×P_0, BDM_1×P_0, and RT_1×P_1. 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, pollutes the computed velocity field so the divergence-free property of the considered elements is lost. 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. Numerical experiments indicate that with the new stabilization terms the unfitted finite element discretization, for the given element pairs, results in 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 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 the computational mesh.

READ FULL TEXT

page 10

page 12

page 13

page 14

research
04/27/2023

Divergence-free cut finite element methods for Stokes flow

We develop two unfitted finite element methods for the Stokes equations ...
research
07/11/2022

Cut finite element method for divergence free approximation of incompressible flow: optimal error estimates and pressure independence

In this note we design a cut finite element method for a low order diver...
research
07/26/2022

On the convergence of a low order Lagrange finite element approach for natural convection problems

The purpose of this article is to study the convergence of a low order f...
research
09/13/2019

Divergence-free tangential finite element methods for incompressible flows on surfaces

In this work we consider the numerical solution of incompressible flows ...
research
06/22/2023

Analysis of divergence-preserving unfitted finite element methods for the mixed Poisson problem

In this paper we present a new H(div)-conforming unfitted finite element...
research
12/18/2021

Variational multiscale modeling with discretely divergence-free subscales: Non-divergence-conforming discretizations

A recent paper [J. A. Evans, D. Kamensky, Y. Bazilevs, "Variational mult...
research
09/09/2020

A mixed finite element method on polytopal mesh

In this paper, we introduce new stable mixed finite elements of any orde...

Please sign up or login with your details

Forgot password? Click here to reset