Domain decomposition preconditioners for high-order discretisations of the heterogeneous Helmholtz equation

by   Shihua Gong, et al.

We consider one-level additive Schwarz domain decomposition preconditioners for the Helmholtz equation with variable coefficients (modelling wave propagation in heterogeneous media), subject to boundary conditions that include wave scattering problems. Absorption is included as a parameter in the problem. This problem is discretised using H^1-conforming nodal finite elements of fixed local degree p on meshes with diameter h = h(k), chosen so that the error remains bounded with increasing k. The action of the one-level preconditioner consists of the parallel solution of problems on subdomains (which can be of general geometry), each equipped with an impedance boundary condition. We prove rigorous estimates on the norm and field of values of the left- or right-preconditioned matrix that show explicitly how the absorption, the heterogeneity in the coefficients and the dependence on the degree enter the estimates. These estimates prove rigorously that, with enough absorption and for k large enough, GMRES is guaranteed to converge in a number of iterations that is independent of k,p, and the coefficients. The theoretical threshold for k to be large enough depends on p and on the local variation of coefficients in subdomains (and not globally). Extensive numerical experiments are given for both the absorptive and the propagative cases; in the latter case we investigate examples both when the coefficients are nontrapping and when they are trapping. These experiments (i) support our theory in terms of dependence on polynomial degree and the coefficients; (ii) support the sharpness of our field of values estimates in terms of the level of absorption required.


Higher-order finite element methods for the nonlinear Helmholtz equation

In this work, we analyze the finite element method with arbitrary but fi...

Wavenumber-explicit convergence of the hp-FEM for the full-space heterogeneous Helmholtz equation with smooth coefficients

A convergence theory for the hp-FEM applied to a variety of constant-coe...

Polynomial bounds for the solutions of parametric transmission problems on smooth, bounded domains

We consider a family (P_ω)_ω∈Ω of elliptic second order differential ope...

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...

Analysis of a Helmholtz preconditioning problem motivated by uncertainty quantification

This paper analyses the following question: let A_j, j=1,2, be the Galer...

Optimal artificial boundary conditions based on second-order correctors for three dimensional random elliptic media

We are interested in numerical algorithms for computing the electrical f...

Wave propagation in high-contrast media: periodic and beyond

This work is concerned with the classical wave equation with a high-cont...

Please sign up or login with your details

Forgot password? Click here to reset