Galkowski, J; Spence, EA; (2023) Does the Helmholtz boundary element method suffer from the pollution effect? SIAM Review , 65 (3) pp. 806-828. 10.1137/22M1474199.

Preview

Text 22m1474199.pdf Download (525kB) | Preview

Abstract

In dimensions, accurately approximating an arbitrary function oscillating with frequency requires degrees of freedom. A numerical method for solving the Helmholtz equation (with wavenumber and in dimensions) suffers from the pollution effect if, as , the total number of degrees of freedom needed to maintain accuracy grows faster than this natural threshold (i.e., faster than for domain-based formulations, such as finite element methods, and for boundary-based formulations, such as boundary element methods). It is well known that the -version of the finite element method (FEM) (where accuracy is increased by decreasing the meshwidth and keeping the polynomial degree fixed) suffers from the pollution effect, and research over the last 30 years has resulted in a near-complete rigorous understanding of how quickly the number of degrees of freedom must grow with to maintain accuracy (and how this depends on both and properties of the scatterer). In contrast to the -FEM, at least empirically, the -version of the boundary element method (BEM) does not suffer from the pollution effect (recall that in the boundary element method the scattering problem is reformulated as an integral equation on the boundary of the scatterer, with this integral equation then solved numerically using a finite element--type approximation space). However, the current best results in the literature on how quickly the number of degrees of freedom for the -BEM must grow with to maintain accuracy fall short of proving this. In this paper, we prove that the -version of the Galerkin method applied to the standard second-kind boundary integral equations for solving the Helmholtz exterior Dirichlet problem does not suffer from the pollution effect when the obstacle is nontrapping (i.e., does not trap geometric-optic rays). While the proof of this result relies on information about the large- behavior of Helmholtz solution operators, we show in an appendix how the result can be proved using only Fourier series and asymptotics of Hankel and Bessel functions when the obstacle is a 2-d ball. It is well known that the h-version of the finite element method (FEM) (where accuracy is increased by decreasing the meshwidth h and keeping the polynomial degree p fixed) suffers from the pollution effect, and research over the last ∼ 30 years has resulted in a near-complete rigorous understanding of how quickly the number of degrees of freedom must grow with k (and how this depends on both p and properties of the scatterer). In contrast to the h-FEM, at least empirically, the h-version of the boundary element method (BEM) does not suffer from the pollution effect (recall that in the boundary element method the scattering problem is reformulated as an integral equation on the boundary of the scatterer, with this integral equation then solved numerically using a finite-element-type approximation space). However, the current best results in the literature on how quickly the number of degrees of freedom for the h-BEM must grow with k fall short of proving this. In this paper, we prove that the h-version of the Galerkin method applied to the standard second-kind boundary integral equations for solving the Helmholtz exterior Dirichlet problem does not suffer from the pollution effect when the obstacle is nontrapping (i.e., does not trap geometric-optic rays).

Download activity - last month

Export as JSONCSVXML

Download activity - last 12 months

Export as CSVXMLJSON

Downloads by country - last 12 months

Export as CSVXMLJSON

Archive Staff Only

View Item