Analysis of a Helmholtz preconditioning problem motivated by uncertainty quantification

From MaRDI portal

Jump to:navigation, search

DOI10.1007/S10444-021-09889-0MaRDI QIDQ2230585zbMATH OpenOpenAlexFDO

Authors I. G. Graham, O. R. Pembery, E. A. Spence

Publication date 24 September 2021

Published in Advances in Computational Mathematics (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/2005.13390

zbMATH Keywords

exterior Dirichlet problem finite element discretisation heterogeneous Helmholtz equation

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Boundary value problems for second-order elliptic equations (35J25) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)

Abstract: This paper analyses the following question: let

m a t h b f 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 b l a c d o t (A_{j} a b l a 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

(m a t h b f A_{1})^{- 1} m a t h b f A_{2}

or

m a t h b f A_{2} (m a t h b f A_{1})^{- 1}

to converge in a

k

-independent number of iterations for arbitrarily large

k

? (In other words, for

m a t h b f A_{1}

to be a good left- or right-preconditioner for

m a t h b f 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.

Recommendations

Cites work

Cited in

(3)

Describes a project that uses

Uses Software

This page was built for publication: Analysis of a Helmholtz preconditioning problem motivated by uncertainty quantification

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2230585)

Retrieved from "https://portal.mardi4nfdi.de/w/index.php?title=Analysis_of_a_Helmholtz_preconditioning_problem_motivated_by_uncertainty_quantification&oldid=71622474"