Convergence of Restricted Additive Schwarz with impedance transmission conditions for discretised Helmholtz problems

The Restricted Additive Schwarz method with impedance transmission conditions, also known as the Optimised Restricted Additive Schwarz (ORAS) method, is a simple overlapping one-level parallel domain decomposition method, which has been successfully used as an iterative solver and as a preconditioner for discretized Helmholtz boundary-value problems. In this paper, we give, for the first time, a convergence analysis for ORAS as an iterative solver – and also as a preconditioner – for nodal finite element Helmholtz systems of any polynomial order. The analysis starts by showing (for general domain decompositions) that ORAS as an unconventional finite element approximation of a classical parallel iterative Schwarz method, formulated at the PDE (non-discrete) level. This non-discrete Schwarz method was recently analysed in [Gong, Gander, Graham, Lafontaine, Spence, arXiv 2106.05218], and the present paper gives a corresponding discrete version of this analysis. In particular, for domain decompositions in strips in 2-d, we show that, when the mesh size is small enough, ORAS inherits the convergence properties of the Schwarz method, independent of polynomial order. The proof relies on characterising the ORAS iteration in terms of discrete ‘impedance-to-impedance maps’, which we prove (via a novel weighted finite-element error analysis) converge as h → 0 in the operator norm to their non-discrete counterparts.

[1]  Marcus Sarkis,et al.  Restricted Overlapping Balancing Domain Decomposition Methods and Restricted Coarse Problems for the Helmholtz Problem , 2007 .

[2]  Haijun Wu,et al.  Discontinuous Galerkin Methods for the Helmholtz Equation with Large Wave Number , 2009, SIAM J. Numer. Anal..

[3]  Zhiming Chen,et al.  A Source Transfer Domain Decomposition Method for Helmholtz Equations in Unbounded Domain , 2013, SIAM J. Numer. Anal..

[4]  T. Dupont,et al.  Polynomial approximation of functions in Sobolev spaces , 1980 .

[5]  Haijun Wu,et al.  hp-Discontinuous Galerkin methods for the Helmholtz equation with large wave number , 2008, Math. Comput..

[6]  Russell J. Hewett,et al.  L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation , 2019, J. Comput. Phys..

[7]  L. R. Scott,et al.  Finite element interpolation of nonsmooth functions satisfying boundary conditions , 1990 .

[8]  Tsuyoshi Murata,et al.  {m , 1934, ACML.

[9]  David Lafontaine,et al.  A sharp relative-error bound for the Helmholtz h-FEM at high frequency , 2019, Numerische Mathematik.

[10]  Bruno Després,et al.  Corners and stable optimized domain decomposition methods for the Helmholtz problem , 2020, Numerische Mathematik.

[11]  Shihua Gong,et al.  Domain decomposition preconditioners for high-order discretisations of the heterogeneous Helmholtz equation , 2020, ArXiv.

[12]  Martin J. Gander,et al.  A Class of Iterative Solvers for the Helmholtz Equation: Factorizations, Sweeping Preconditioners, Source Transfer, Single Layer Potentials, Polarized Traces, and Optimized Schwarz Methods , 2016, SIAM Rev..

[13]  Martin J. Gander,et al.  Optimized Multiplicative, Additive, and Restricted Additive Schwarz Preconditioning , 2007, SIAM J. Sci. Comput..

[14]  Martin J. Gander,et al.  Applying GMRES to the Helmholtz equation with shifted Laplacian preconditioning: what is the largest shift for which wavenumber-independent convergence is guaranteed? , 2015, Numerische Mathematik.

[15]  M. Gander,et al.  Why Restricted Additive Schwarz Converges Faster than Additive Schwarz , 2003 .

[16]  Martin J. Gander,et al.  A variational interpretation of Restricted Additive Schwarz with impedance transmission condition for the Helmholtz problem , 2021, ArXiv.

[17]  Pierre-Henri Tournier,et al.  A comparison of coarse spaces for Helmholtz problems in the high frequency regime , 2020, Comput. Math. Appl..

[18]  Lexing Ying,et al.  Sweeping Preconditioner for the Helmholtz Equation: Moving Perfectly Matched Layers , 2010, Multiscale Model. Simul..

[19]  Axel Modave,et al.  A non-overlapping domain decomposition method with high-order transmission conditions and cross-point treatment for Helmholtz problems , 2020, Computer Methods in Applied Mechanics and Engineering.

[20]  Haijun Wu,et al.  Pre-asymptotic error analysis of CIP-FEM and FEM for the Helmholtz equation with high wave number. Part I: linear version , 2014 .

[21]  Cornelis Vuik,et al.  A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems , 2005, SIAM J. Sci. Comput..

[22]  E. A. Spence,et al.  Domain Decomposition with Local Impedance Conditions for the Helmholtz Equation with Absorption , 2018, SIAM J. Numer. Anal..

[23]  W. McLean Strongly Elliptic Systems and Boundary Integral Equations , 2000 .

[24]  S. Nicaise,et al.  High-frequency behaviour of corner singularities in Helmholtz problems , 2018, ESAIM: Mathematical Modelling and Numerical Analysis.

[25]  S. Nicaise,et al.  Uniform a priori estimates for elliptic problems with impedance boundary conditions , 2020, Communications on Pure & Applied Analysis.

[26]  X. Claeys,et al.  Robust treatment of cross-points in optimized Schwarz methods , 2020, Numerische Mathematik.

[27]  Victorita Dolean,et al.  An introduction to domain decomposition methods - algorithms, theory, and parallel implementation , 2015 .

[28]  Haijun Wu,et al.  Preasymptotic Error Analysis of Higher Order FEM and CIP-FEM for Helmholtz Equation with High Wave Number , 2014, SIAM J. Numer. Anal..

[29]  Ivan G. Graham,et al.  Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption , 2017, Math. Comput..

[30]  Pierre Jolivet,et al.  Microwave Tomographic Imaging of Cerebrovascular Accidents by Using High-Performance Computing , 2016, Parallel Comput..

[31]  Stefan A. Sauter,et al.  A Refined Finite Element Convergence Theory for Highly Indefinite Helmholtz Problems , 2006, Computing.

[32]  A. H. Schatz,et al.  An observation concerning Ritz-Galerkin methods with indefinite bilinear forms , 1974 .

[33]  Bruno Després,et al.  A Domain Decomposition Method for the Helmholtz equation and related Optimal Control Problems , 1996 .

[34]  Martin J. Gander,et al.  Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation , 2021, ArXiv.