Table of Contents
Fetching ...

Monolithic two-level Schwarz preconditioner for Biot's consolidation model in two space dimensions

Stefan Meggendorfer, Guido Kanschat, Johannes Kraus

TL;DR

The paper develops a monolithic two-level overlapping Schwarz preconditioner for the quasi-static Biot consolidation model discretized with mass-conserving $H^{\text{div}}$-conforming methods in two space dimensions. By transforming the saddle-point system into an equivalent singularly perturbed SPD problem and establishing stable space decompositions, the authors prove uniform convergence of both additive and multiplicative Schwarz iterations under mild mesh-size relations and bounded permeability. Theoretical results are complemented by 2D numerical experiments demonstrating parameter-robust performance across a range of Lamé, permeability, and storage coefficients, with a coarse-space solve that scales favorably. The work advances robust preconditioning for poroelastic systems and lays groundwork for extensions to 3D and high-frequency/contrast settings.

Abstract

This paper addresses the construction and analysis of a class of domain decomposition methods for the iterative solution of the quasi-static Biot problem in three-field formulation. The considered discrete model arises from time discretization by the implicit Euler method and space discretization by a family of strongly mass-conserving methods exploiting $H^{div}$-conforming approximations of the solid displacement and fluid flux fields. For the resulting saddle-point problem, we construct monolithic overlapping domain decomposition (DD) methods whose analysis relies on a transformation into an equivalent symmetric positive definite system and on stable decompositions of the involved finite element spaces under proper problem-dependent norms. Numerical results on two-dimensional test problems are in accordance with the provided theoretical uniform convergence estimates for the two-level multiplicative Schwarz method.

Monolithic two-level Schwarz preconditioner for Biot's consolidation model in two space dimensions

TL;DR

The paper develops a monolithic two-level overlapping Schwarz preconditioner for the quasi-static Biot consolidation model discretized with mass-conserving -conforming methods in two space dimensions. By transforming the saddle-point system into an equivalent singularly perturbed SPD problem and establishing stable space decompositions, the authors prove uniform convergence of both additive and multiplicative Schwarz iterations under mild mesh-size relations and bounded permeability. Theoretical results are complemented by 2D numerical experiments demonstrating parameter-robust performance across a range of Lamé, permeability, and storage coefficients, with a coarse-space solve that scales favorably. The work advances robust preconditioning for poroelastic systems and lays groundwork for extensions to 3D and high-frequency/contrast settings.

Abstract

This paper addresses the construction and analysis of a class of domain decomposition methods for the iterative solution of the quasi-static Biot problem in three-field formulation. The considered discrete model arises from time discretization by the implicit Euler method and space discretization by a family of strongly mass-conserving methods exploiting -conforming approximations of the solid displacement and fluid flux fields. For the resulting saddle-point problem, we construct monolithic overlapping domain decomposition (DD) methods whose analysis relies on a transformation into an equivalent symmetric positive definite system and on stable decompositions of the involved finite element spaces under proper problem-dependent norms. Numerical results on two-dimensional test problems are in accordance with the provided theoretical uniform convergence estimates for the two-level multiplicative Schwarz method.
Paper Structure (12 sections, 12 theorems, 119 equations, 4 tables)

This paper contains 12 sections, 12 theorems, 119 equations, 4 tables.

Key Result

Theorem 4.1

Assume that $H \le c h$ for some positive constant $c$ and that the bounds eq:normequivalence_k on $\hat{\bm K}$ hold. Then, the additive Schwarz method eq:additiveschwarz defines a uniform preconditioner for the system eq:discreteproblem. Consequently, the GMRES method converges uniformly, i.e., it

Theorems & Definitions (26)

  • Theorem 4.1
  • Theorem 4.2
  • Remark 4.3
  • Remark 4.4
  • Lemma 4.5
  • proof
  • Proposition 4.6
  • Lemma 4.7
  • proof
  • Lemma 4.8
  • ...and 16 more