New SFB 1313 publication, published in "Computational Geosciences". The work has been developed within the SFB 1313 research project I-01 and A-02.
Authors
- Wietse M. Boon (Politecnico di Milano, former University of Stuttgart, research project I-01)
- Dennis Gläser (University of Stuttgart, research project I-03)
- Rainer Helmig (University of Stuttgart, research project A-02)
- Kilian Weishaupt (University of Stuttgart, research project A-02)
- Ivan Yotov (University of Pittsburgh, SFB1313 External Partner, University Stuttgart Visiting Mercator Fellow, research project A-02)
Abstract
A discretization method with non-matching grids is proposed for the coupled Stokes-Darcy problem that uses a mortar variable at the interface to couple the marker and cell (MAC) method in the Stokes domain with the Raviart-Thomas mixed finite element pair in the Darcy domain. Due to this choice, the method conserves linear momentum and mass locally in the Stokes domain and exhibits local mass conservation in the Darcy domain. The MAC scheme is reformulated as a mixed finite element method on a staggered grid, which allows for the proposed scheme to be analyzed as a mortar mixed finite element method. We show that the discrete system is well-posed and derive a priori error estimates that indicate first order convergence in all variables. The system can be reduced to an interface problem concerning only the mortar variables, leading to a non-overlapping domain decomposition method. Numerical examples are presented to illustrate the theoretical results and the applicability of the method.