Convergence of parallel overlapping domain decomposition methods for the Helmholtz equation

by   Shihua Gong, et al.

We analyse parallel overlapping Schwarz domain decomposition methods for the Helmholtz equation, where the subdomain problems satisfy first-order absorbing (impedance) transmission conditions, and exchange of information between subdomains is achieved using a partition of unity. We provide a novel analysis of this method at the PDE level (without discretization). First, we formulate the method as a fixed point iteration, and show (in dimensions 1,2,3) that it is well-defined in a tensor product of appropriate local function spaces, each with L^2 impedance boundary data. Given this, we then obtain a bound on the norm of the fixed point operator in terms of the local norms of certain impedance-to-impedance maps arising from local interactions between subdomains. These bounds provide conditions under which (some power of) the fixed point operator is a contraction. In 2-d, for rectangular domains and strip-wise domain decompositions (with each subdomain only overlapping its immediate neighbours), we present two techniques for verifying the assumptions on the impedance-to-impedance maps which ensure power contractivity of the fixed point operator. The first is through semiclassical analysis, which gives rigorous estimates valid as the frequency tends to infinity. These results verify the required assumptions for sufficiently large overlap. For more realistic domain decompositions, we directly compute the norms of the impedance-to-impedance maps by solving certain canonical (local) eigenvalue problems. We give numerical experiments that illustrate the theory. These also show that the iterative method remains convergent and/or provides a good preconditioner in cases not covered by the theory, including for general domain decompositions, such as those obtained via automatic graph-partitioning software.


Sharp bounds on Helmholtz impedance-to-impedance maps and application to overlapping domain decomposition

We prove sharp bounds on certain impedance-to-impedance maps (and their ...

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

The Restricted Additive Schwarz method with impedance transmission condi...

A Fixed Point Theorem for Iterative Random Contraction Operators over Banach Spaces

Consider a contraction operator T over a Banach space X with a fixed po...

Analysis of the SORAS domain decomposition preconditioner for non-self-adjoint or indefinite problems

We analyze the convergence of the one-level overlapping domain decomposi...

H(curl^2) conforming element for Maxwell's transmission eigenvalue problem using fixed-point approach

Using newly developed H(curl^2) conforming elements, we solve the Maxwel...

Robust treatment of cross points in Optimized Schwarz Methods

In the field of Domain Decomposition (DD), Optimized Schwarz Method (OSM...

On the Iteration Complexity of Hypergradient Computation

We study a general class of bilevel problems, consisting in the minimiza...

Please sign up or login with your details

Forgot password? Click here to reset