The geometric error is less than the pollution error when solving the high-frequency Helmholtz equation with high-order FEM on curved domains

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Chaumont-Frelet, Théophile, Spence, Euan A.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915214438432768
author Chaumont-Frelet, Théophile
Spence, Euan A.
author_facet Chaumont-Frelet, Théophile
Spence, Euan A.
contents We consider the $h$-version of the finite-element method, where accuracy is increased by decreasing the meshwidth $h$ while keeping the polynomial degree $p$ constant, applied to the Helmholtz equation. Although the question "how quickly must $h$ decrease as the wavenumber $k$ increases to maintain accuracy?" has been studied intensively since the 1990s, none of the existing rigorous wavenumber-explicit analyses take into account the approximation of the geometry. In this paper we prove that for nontrapping problems solved using straight elements the geometric error is order $kh$, which is then less than the pollution error $k(kh)^{2p}$ when $k$ is large; this fact is then illustrated in numerical experiments. More generally, we prove that, even for problems with strong trapping, using degree four (in 2-d) or degree five (in 3-d) polynomials and isoparametric elements ensures that the geometric error is smaller than the pollution error for most large wavenumbers.
format Preprint
id arxiv_https___arxiv_org_abs_2401_16413
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle The geometric error is less than the pollution error when solving the high-frequency Helmholtz equation with high-order FEM on curved domains
Chaumont-Frelet, Théophile
Spence, Euan A.
Numerical Analysis
We consider the $h$-version of the finite-element method, where accuracy is increased by decreasing the meshwidth $h$ while keeping the polynomial degree $p$ constant, applied to the Helmholtz equation. Although the question "how quickly must $h$ decrease as the wavenumber $k$ increases to maintain accuracy?" has been studied intensively since the 1990s, none of the existing rigorous wavenumber-explicit analyses take into account the approximation of the geometry. In this paper we prove that for nontrapping problems solved using straight elements the geometric error is order $kh$, which is then less than the pollution error $k(kh)^{2p}$ when $k$ is large; this fact is then illustrated in numerical experiments. More generally, we prove that, even for problems with strong trapping, using degree four (in 2-d) or degree five (in 3-d) polynomials and isoparametric elements ensures that the geometric error is smaller than the pollution error for most large wavenumbers.
title The geometric error is less than the pollution error when solving the high-frequency Helmholtz equation with high-order FEM on curved domains
topic Numerical Analysis
url https://arxiv.org/abs/2401.16413