Convergence analysis of Padé approximations for Helmholtz frequency response problems^★

¹ Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien, Austria
² CSQI – Calcul Scientifique et Quantification de l’Incertitude, MATHICSE, École Polytechnique Fédérale de Lausanne, Station 8, CH-1015 Lausanne, Switzerland
³ Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien, Austria

^* Corresponding author: francesca.bonizzoni@univie.ac.at

Received: 15 July 2016
Accepted: 28 September 2017

Abstract

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.