Analysis of a Helmholtz preconditioning problem motivated by uncertainty quantification

by   Ivan G. Graham, et al.

This paper analyses the following question: let A_j, j=1,2, be the Galerkin matrices corresponding to finite-element discretisations of the exterior Dirichlet problem for the heterogeneous Helmholtz equations ∇· (A_j ∇ u_j) + k^2 n_j u_j= -f. How small must A_1 -A_2_L^q and n_1 - n_2_L^q be (in terms of k-dependence) for GMRES applied to either (A_1)^-1A_2 or A_2(A_1)^-1 to converge in a k-independent number of iterations for arbitrarily large k? (In other words, for A_1 to be a good left- or right-preconditioner for A_2?). We prove results answering this question, give theoretical evidence for their sharpness, and give numerical experiments supporting the estimates. Our motivation for tackling this question comes from calculating quantities of interest for the Helmholtz equation with random coefficients A and n. Such a calculation may require the solution of many deterministic Helmholtz problems, each with different A and n, and the answer to the question above dictates to what extent a previously-calculated inverse of one of the Galerkin matrices can be used as a preconditioner for other Galerkin matrices.


page 1

page 2

page 3

page 4


Discontinuous Galerkin Finite Element Methods for 1D Rosenau Equation

In this paper, discontinuous Galerkin finite element methods are applied...

Curved Elements in Weak Galerkin Finite Element Methods

A mathematical analysis is established for the weak Galerkin finite elem...

Truncation Preconditioners for Stochastic Galerkin Finite Element Discretizations

Stochastic Galerkin finite element method (SGFEM) provides an efficient ...

Stochastic Discontinuous Galerkin Methods for Robust Deterministic Control of Convection Diffusion Equations with Uncertain Coefficients

We investigate a numerical behaviour of robust deterministic optimal con...

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

We consider one-level additive Schwarz domain decomposition precondition...

Please sign up or login with your details

Forgot password? Click here to reset