Open Collections

Convergence analysis of Padé approximations for Helmholtz problems with parametric/ stochastic wavenumber

Banff International Research Station for Mathematical Innovation and Discovery

The present work concerns the approximation of the solution map S associated to the parametric Helmholtz boundary value problem, i.e., the map which associates to each (real) wavenumber belonging to a given interval of interest the corresponding solution of the Helmholtz equation. We introduce a least squares rational Padé-type approximation technique applicable to any meromorphic Hilbert space-valued univariate map, and we prove the uniform convergence of the Padé approximation error on any compact subset of the interval of interest that excludes any pole. This general result is then applied to the Helmholtz solution map S, which is proven to be meromorphic in â , with a pole of order one in every (single or multiple) eigenvalue of the Laplace operator with the considered boundary conditions. Numerical tests are provided that confirm the theoretical upper bound on the Padé approximation error for the Helmholtz solution map. The Padé-type approximation can then be used to compute the probability distribution of quantities of interest associated to the Helmholtz problem with random wavenumber.

Mathematics; Probability theory and stochastic processes; Numerical analysis; Industrial mathematics / modelling and simulation

video/mp4

Author affiliation: Ecole Polytechnique Fédéral de Lausanne

2019-03-10

Attribution-NonCommercial-NoDerivatives 4.0 International

http://hdl.handle.net/2429/68575

Unreviewed

http://creativecommons.org/licenses/by-nc-nd/4.0/