The hp-FEM applied to the Helmholtz equation with PML truncation does not suffer from the pollution effect

by   Jeffrey Galkowski, et al.

We consider approximation of the variable-coefficient Helmholtz equation in the exterior of a Dirichlet obstacle using perfectly-matched-layer (PML) truncation; it is well known that this approximation is exponentially accurate in the PML width and the scaling angle, and the approximation was recently proved to be exponentially accurate in the wavenumber k in [Galkowski, Lafontaine, Spence, 2021]. We show that the hp-FEM applied to this problem does not suffer from the pollution effect, in that there exist C_1,C_2>0 such that if hk/p≤ C_1 and p ≥ C_2 log k then the Galerkin solutions are quasioptimal (with constant independent of k), under the following two conditions (i) the solution operator of the original Helmholtz problem is polynomially bounded in k (which occurs for "most" k by [Lafontaine, Spence, Wunsch, 2021]), and (ii) either there is no obstacle and the coefficients are smooth or the obstacle is analytic and the coefficients are analytic in a neighbourhood of the obstacle and smooth elsewhere. This hp-FEM result is obtained via a decomposition of the PML solution into "high-" and "low-frequency" components, analogous to the decomposition for the original Helmholtz solution recently proved in [Galkowski, Lafontaine, Spence, Wunsch, 2022]. The decomposition is obtained using tools from semiclassical analysis (i.e., the PDE techniques specifically designed for studying Helmholtz problems with large k).


page 1

page 2

page 3

page 4


Decompositions of high-frequency Helmholtz solutions via functional calculus, and application to the finite element method

Over the last ten years, results from [Melenk-Sauter, 2010], [Melenk-Sau...

Wavenumber-explicit convergence of the hp-FEM for the full-space heterogeneous Helmholtz equation with smooth coefficients

A convergence theory for the hp-FEM applied to a variety of constant-coe...

Interaction with an obstacle in the 2d focusing nonlinear Schrödinger equation

We present a numerical study of solutions to the 2d focusing nonlinear S...

Perfectly-matched-layer truncation is exponentially accurate at high frequency

We consider a wide variety of scattering problems including scattering b...

Sharp preasymptotic error bounds for the Helmholtz h-FEM

In the analysis of the h-version of the finite-element method (FEM), wit...

A simple proof that the hp-FEM does not suffer from the pollution effect for the constant-coefficient full-space Helmholtz equation

In d dimensions, approximating an arbitrary function oscillating with fr...

Numerical solution of a one-dimensional nonlocal Helmholtz equation by Perfectly Matched Layers

We consider the computation of a nonlocal Helmholtz equation by using Pe...

Please sign up or login with your details

Forgot password? Click here to reset