Adaptive FEM for Helmholtz Equation with Large Wavenumber

• Published in: Journal of scientific computing 2023, Vol.94 (1), p.21, Article 21
• Main Authors: Duan, Songyao, Wu, Haijun
• Format: Article
• Language: English
• Subjects: Adaptive algorithms; Algorithms; Approximation; Computational Mathematics and Numerical Analysis; Errors; Estimates; Exact solutions; Finite element analysis; Finite element method; Helmholtz equations; Lower bounds; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematics; Mathematics and Statistics; Partial differential equations; Theoretical; Wavelengths
• ISSN: 0885-7474; 1573-7691
• DOI: 10.1007/s10915-022-02074-5

A posteriori upper and lower bounds are derived for the linear finite element method (FEM) for the Helmholtz equation with large wavenumber. It is proved rigorously that the standard residual type error estimator seriously underestimates the true error of the FE solution for the mesh size h in the preasymptotic regime, which is first observed by Babuška et al. (Int J Numer Methods Eng 40:3443–3462, 1997) for a one dimensional problem. By establishing an equivalence relationship between the error estimators for the FE solution and the corresponding elliptic projection of the exact solution, an adaptive algorithm is proposed and its convergence and quasi-optimality are proved under the condition that k 3 h 0 1 + α is sufficiently small, where k is the wavenumber, h 0 is the initial mesh size, and 1 2 < α ≤ 1 is a regularity constant depending on the maximum reentrant angle of the domain. Numerical tests are given to verify the theoretical findings and to show that the adaptive continuous interior penalty finite element method (CIP-FEM) with appropriately selected penalty parameters can greatly reduce the pollution error and hence the residual type error estimator for this CIP-FEM is reliable and efficient even in the preasymptotic regime.