Vol. 5, No. 1, 2010

Download this article
Download this article For screen
For printing
Recent Issues
Volume 18, Issue 1
Volume 17, Issue 1
Volume 16, Issue 2
Volume 16, Issue 1
Volume 15, Issue 2
Volume 15, Issue 1
Volume 14, Issue 2
Volume 14, Issue 1
Volume 13, Issue 2
Volume 13, Issue 1
Volume 12, Issue 1
Volume 11, Issue 2
Volume 11, Issue 1
Volume 10, Issue 2
Volume 10, Issue 1
Volume 9, Issue 2
Volume 9, Issue 1
Volume 8, Issue 1
Volume 7, Issue 2
Volume 7, Issue 1
Volume 6, Issue 1
Volume 5, Issue 2
Volume 5, Issue 1
Volume 4, Issue 1
Volume 3, Issue 1
Volume 2, Issue 1
Volume 1, Issue 1
The Journal
About the Journal
Editorial Board
Submission Guidelines
Submission Form
Policies for Authors
Ethics Statement
ISSN: 2157-5452 (e-only)
ISSN: 1559-3940 (print)
Author Index
To Appear
Other MSP Journals
FETI and BDD preconditioners for Stokes–Mortar–Darcy Systems

Juan Galvis and Marcus Sarkis

Vol. 5 (2010), No. 1, 1–30

We consider the coupling across an interface of a fluid flow and a porous media flow. The differential equations involve Stokes equations in the fluid region, Darcy equations in the porous region, plus a coupling through an interface with Beaver–Joseph–Saffman transmission conditions. The discretization consists of P2/P1 triangular Taylor–Hood finite elements in the fluid region, the lowest order triangular Raviart–Thomas finite elements in the porous region, and the mortar piecewise constant Lagrange multipliers on the interface. We allow for nonmatching meshes across the interface. Due to the small values of the permeability parameter κ of the porous medium, the resulting discrete symmetric saddle point system is very ill conditioned. We design and analyze preconditioners based on the finite element by tearing and interconnecting (FETI) and balancing domain decomposition (BDD) methods and derive a condition number estimate of order C1(1 + (1κ)) for the preconditioned operator. In case the fluid discretization is finer than the porous side discretization, we derive a better estimate of order C2((κ + 1)(κ + (hp)2)) for the FETI preconditioner. Here hp is the mesh size of the porous side triangulation. The constants C1 and C2 are independent of the permeability κ, the fluid viscosity ν, and the mesh ratio across the interface. Numerical experiments confirm the sharpness of the theoretical estimates.

Stokes–Darcy coupling, mortar, balancing domain decomposition, FETI, saddle point problems, nonmatching grids, discontinuous coefficients, mortar elements
Mathematical Subject Classification 2000
Primary: 35Q30, 65N22, 65N30, 65N55, 76D07
Received: 22 November 2008
Revised: 17 November 2009
Accepted: 29 November 2009
Published: 27 December 2009
Juan Galvis
Department of Mathematics
Texas A&M University
College Station, TX 77843-3368
United States
Marcus Sarkis
Worcester Polytechnic Institute
Mathematical Sciences Department
100 Institute Road
Worcester, MA 01609
United States
Instituto Nacional de Matemática Pura e Aplicada
Estrada Dona Castorina 110
22460-320 Rio de Janeiro, RJ