Eliminating the pollution effect in Helmholtz problems by local subscale correction

Abstract

We introduce a new Petrov-Galerkin multiscale method for the numerical approximation of the Helmholtz equation with large wave number κ \kappa in bounded domains in R d \mathbb {R}^d . The discrete trial and test spaces are generated from standard mesh-based finite elements by local subscale correction in the spirit of numerical homogenization. The precomputation of the correction involves the solution of coercive cell problems on localized subdomains of size ℓ H \ell H , H H being the mesh size and ℓ \ell being the oversampling parameter. If the mesh size and the oversampling parameter are such that H κ H\kappa and log ⁡ ( κ ) / ℓ \log (\kappa )/\ell fall below some generic constants and if the cell problems are solved sufficiently accurately on some finer scale of discretization, then the method is stable and its error is proportional to H H . Pollution effects are eliminated in this regime.