A posteriori error estimates and adaptive mesh refinement for the Stokes-Brinkman problem

TL;DR

This paper develops residual-based a posteriori error estimates for the Stokes-Brinkman equations, enabling adaptive mesh refinement to efficiently simulate flow in heterogeneous porous media, demonstrated through numerical experiments in 2D and 3D.

Contribution

It introduces a novel a posteriori error estimation method for the Stokes-Brinkman problem and applies it to adaptive mesh refinement, improving simulation accuracy and efficiency.

Findings

Adaptive mesh refinement improves solution accuracy.

Error estimates effectively guide mesh refinement.

Numerical experiments confirm method effectiveness in 2D and 3D.

Abstract

The Stokes-Brinkman equations model flow in heterogeneous porous media by combining the Stokes and Darcy models of flow into a single system of equations. With suitable parameters, the equations can model either flow without detailed knowledge of the interface between the two regions. Thus, the Stokes-Brinkman equations provide an alternative to coupled Darcy-Stokes models. After a brief review of the Stokes-Brinkman problem and its discretization using Taylor-Hood finite elements, we present a residual-based a posteriori error estimate and use it to drive an adaptive mesh refinement process. We compare several strategies for the mesh refinement, and demonstrate its effectiveness by numerical experiments in both 2D and 3D.