A residual-based a posteriori error estimator for a fully-mixed formulation of the Stokes–Darcy coupled problem

Abstract

In this paper we develop an a posteriori error analysis of a new fully mixed finite element method for the coupling of fluid flow with porous media flow in 2D. Flows are governed by the Stokes and Darcy equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We consider dual-mixed formulations in both media, which yields the pseudostress and the velocity in the fluid, together with the velocity and the pressure in the porous medium, and the traces of the porous media pressure and the fluid velocity on the interface, as the resulting unknowns. The set of feasible finite element subspaces includes Raviart–Thomas elements of lowest order and piecewise constants for the velocities and pressures, respectively, in both domains, together with continuous piecewise linear elements for the traces. We derive a reliable and efficient residual-based a posteriori error estimator for the coupled problem. The proof of reliability makes use of the global inf–sup condition, Helmholtz decompositions in both media, and local approximation properties of the Clément interpolant and Raviart–Thomas operator. On the other hand, inverse inequalities, the localization technique based on element-bubble and edge-bubble functions, and known results from previous works, are the main tools for proving the efficiency of the estimator. Finally, some numerical results confirming the theoretical properties of this estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities of the solution, are reported.