Hybrid Schwarz preconditioners for linear systems arising from hp-discontinuous Galerkin method

TL;DR

This paper introduces a hybrid Schwarz preconditioner for efficiently solving linear systems from hp-discontinuous Galerkin methods applied to elliptic problems, demonstrating improved convergence and adaptability to variable diffusion coefficients.

Contribution

The paper develops a novel two-level hybrid Schwarz preconditioner with spectral bounds and compares its performance to standard methods, including applications to nonlinear problems.

Findings

Spectral bound of the preconditioned operator is O((H/h)(p^2/q)).

Hybrid Schwarz preconditioner shows faster convergence than standard additive methods.

Effective for problems with varying diffusivity and domain decomposition respecting material interfaces.

Abstract

We deal with the numerical solution of linear elliptic problems with varying diffusion coefficient by the $hp$-discontinuous Galerkin method. We develop a two-level hybrid Schwarz preconditioner for the arising linear algebraic systems. The preconditioner is additive with respect to the local components and multiplicative with respect to the mesh levels. We derive the $hp$ spectral bound of the preconditioned operator in the form $O((H/h)(p^2/q))$, where $H$ and $h$ are the element sizes of the coarse and fine meshes, respectively, and $p$ and $q$ are the polynomial approximation degrees on the fine and coarse meshes. Further, we present a numerical study comparing the hybrid Schwarz preconditioner with the standard additive one from the point of view of the speed of convergence and also computational costs. Moreover, we investigate the convergence of both techniques with respect to the diffusivity variation and to the domain decomposition (non-)respecting the material interfaces. Finally, the combination with a $hp$-mesh adaptation for the solution of nonlinear problem demonstrates the potential of this approach.