Abstract:In this work, we develop algebraic solvers for linear systems arising from the discretization of second-order elliptic problems by saddle-point mixed finite element methods of arbitrary polynomial degree $p \ge 0$. We present a multigrid and a two-level domain decomposition approach in two or three space dimensions, which are steered by their respective a posteriori estimators of the algebraic error. First, we extend the results of [A. Miraçi, J. Papež, and M. Vohralík, SIAM J. Sci. Comput. 43 (2021), S117--S145] to the mixed finite element setting. Extending the multigrid procedure itself is rather natural. To obtain analogous theoretical results, however, a multilevel stable decomposition of the velocity space is needed. In two space dimensions, we can treat the velocity space as the curl of a stream-function space, for which the previous results apply. In three space dimensions, we design a novel stable decomposition by combining a one-level high-order local stable decomposition of [Chaumont-Frelet and Vohralík, SIAM J. Numer. Anal. 61 (2023), 1783--1818] and a multilevel lowest-order stable decomposition of [Hiptmair, Wu, and Zheng, Numer. Math. Theory Methods Appl. 5 (2012), 297--332]. This allows us to prove that our multigrid solver contracts the algebraic error at each iteration and, simultaneously, that the associated a posteriori estimator is efficient. A $p$-robust contraction is shown in two space dimensions. Next, we use this multilevel methodology to define a two-level domain decomposition method where the subdomains consist of overlapping patches of coarse-level elements sharing a common coarse-level vertex. We again establish a ($p$-robust) contraction of the solver and efficiency of the a posteriori estimator. Numerical results presented both for the multigrid approach and the domain decomposition method confirm the theoretical findings.

Robust domain decomposition methods for high-contrast multiscale problems on irregular domains with virtual element discretizations

Scalable solvers for complex electromagnetics problems

Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption

A Multiresolution Domain Decomposition Preconditioner for the MoM Solution of Multiscale Complex Structures

Space-time balancing domain decomposition

Hybrid Localized Spectral Decomposition for multiscale problems

A robust and adaptive GenEO-type domain decomposition preconditioner for $\mathbf{H}(\mathbf{curl})$ problems in three-dimensional general topologies

Robust Preconditioners for DG-Discretizations with Arbitrary Polynomial Degrees

A robust two-level overlapping preconditioner for Darcy flow in high-contrast media

Coarse Spaces Based on Higher-Order Interpolation for Schwarz Preconditioners for Helmholtz Problems

Substructuring preconditioners with novel interface solvers for general elliptic-type equations in three dimensions

On Edge Multiscale Space based Hybrid Schwarz Preconditioner for Helmholtz Problems with Large Wavenumbers

Analysis and preconditioning of a probabilistic domain decomposition algorithm for elliptic boundary value problems

Two-level hybrid Schwarz Preconditioners for The Helmholtz Equation with high wave number

A-posteriori-steered $p$-robust multigrid and domain decomposition methods with optimal step-sizes for mixed finite element discretizations of elliptic problems

An explicit substructuring method for overlapping domain decomposition based on stochastic calculus

A modified domain decomposition spectral collocation method for parabolic partial differential equations

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

A Robust Algebraic Domain Decomposition Preconditioner for Sparse Normal Equations

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry