Abstract:This paper presents a pressure-robust enriched Galerkin (EG) method for the Brinkman equations with minimal degrees of freedom based on EG velocity and pressure spaces. The velocity space consists of linear Lagrange polynomials enriched by a discontinuous, piecewise linear, and mean-zero vector function per element, while piecewise constant functions approximate the pressure. We derive, analyze, and compare two EG methods in this paper: standard and robust methods. The standard method requires a mesh size to be less than a viscous parameter to produce stable and accurate velocity solutions, which is impractical in the Darcy regime. Therefore, we propose the pressure-robust method by utilizing a velocity reconstruction operator and replacing EG velocity functions with a reconstructed velocity. The robust method yields error estimates independent of a pressure term and shows uniform performance from the Stokes to Darcy regimes, preserving minimal degrees of freedom. We prove well-posedness and error estimates for both the standard and robust EG methods. We finally confirm theoretical results through numerical experiments with two- and three-dimensional examples and compare the methods' performance to support the need for the robust method.

What problem does this paper attempt to address?

This paper attempts to solve several key problems in the numerical solution of the Brinkman equation, especially the problem of poor performance of the standard method when the viscosity parameter ν is close to 0 (i.e., the Darcy region). Specifically: 1. **Limitations of the Standard Method**: When dealing with the Brinkman equation, the standard mixed finite element method will produce sub - optimal convergence orders or even not converge in the Darcy region (ν → 0) if the mesh size does not satisfy the condition of \(h < \sqrt{\nu}\). This is because the standard method requires the velocity field to satisfy H(div) - consistency in the Darcy region, and commonly used Stokes elements such as MINI and Taylor - Hood elements cannot meet this requirement. 2. **The Need for Pressure Robustness**: In the Darcy region, the velocity error of the standard method is related to the pressure term, which leads to instability and inaccuracy of the numerical solution. Therefore, it is necessary to develop a pressure - robust method so that the velocity error estimate is independent of the pressure term and has uniform performance throughout the range from Stokes to Darcy. 3. **The Requirement for the Minimum Degrees of Freedom**: In order to improve computational efficiency, the new method needs to use as few degrees of freedom as possible while maintaining pressure robustness. This requires that the method is not only mathematically stable but also efficient in practical applications. For this purpose, the paper proposes a pressure - robust Enriched Galerkin (EG) method. By introducing a velocity reconstruction operator, the velocity function in the EG velocity space is mapped to the Brezzi - Douglas - Marini (BDM) space, thereby ensuring that the velocity function is H(div) - consistent in the Darcy term. This method not only eliminates the strict limitation on the mesh size but also ensures stability and accuracy in different regions. In summary, the main objective of the paper is to develop a new numerical method that can provide stable, accurate and efficient numerical solutions when dealing with the Brinkman equation, whether in the Stokes region or the Darcy region.

A uniform and pressure-robust enriched Galerkin method for the Brinkman equations

A Uniform and Pressure-Robust Enriched Galerkin Method for the Brinkman Equations

Pressure-robust enriched Galerkin methods for the Stokes equations

Robust globally divergence-free weak Galerkin methods for stationary incompressible convective Brinkman-Forchheimer equations

Robust globally divergence-free weak Galerkin methods for unsteady incompressible convective Brinkman-Forchheimer equations

A Stable Numerical Algorithm for the Brinkman Equations by Weak Galerkin Finite Element Methods

A New Weak Galerkin Finite Element Scheme for the Brinkman Model

A mixed discontinuous Galerkin method with symmetric stress for Brinkman problem based on the velocity–pseudostress formulation

A conforming discontinuous Galerkin finite element method for Brinkman equations

Stabilized equal lower-order finite element methods for simulating Brinkman equations in porous media

The pressure-robust weak Galerkin finite element method for Stokes-Darcy problem

A weak Galerkin pseudostress-based mixed finite element method on polygonal meshes: application to the Brinkman problem appearing in porous media

A priori and a posteriori error estimates of a really pressure-robust virtual element method for the incompressible Brinkman problem

Robust Approximation of Generalized Biot-Brinkman Problems

A pressure-residual augmented GLS stabilized method for a type of Stokes equations with nonstandard boundary conditions

A pressure-robust virtual element method for the Stokes problem

Enriched Galerkin methods for two-phase flow in porous media with capillary pressure

A locally mass-conservative enriched Petrov-Galerkin method without penalty for the Darcy flow in porous media

A pressure-robust mixed finite element method for the coupled Stokes–Darcy problem

A pressure-robust virtual element method for the Navier-Stokes problem on polygonal mesh

Coupling nonconforming and enriched Galerkin methods for robust discretization and fast solvers of poroelasticity problems