Amplitude-based generalized plane waves: New quasi-trefftz functions for scalar equations in two dimensions

Research output: Contribution to journalArticlepeer-review

Abstract

Generalized plane waves (GPWs) were introduced to take advantage of Trefftz methods for problems modeled by variable coefficient equations. Despite the fact that GPWs do not satisfy the Trefftz property, i.e., they are not exact solutions to the governing equation, they instead satisfy a quasi-Trefftz property: They are only approximate solutions. They lead to high-order numerical methods, and this quasi-Trefftz property is critical for their numerical analysis. The present work introduces a new family of GPWs: amplitude-based. The motivation lies in the poor behavior of the phase-based GPW approximations in the preasymptotic regime, which will be tamed by avoiding high-degree polynomials within an exponential. The new ansatz introduces higher-order terms in the amplitude rather than in the phase of a plane wave as was initially proposed. The new functions' construction and the study of their approximation properties are guided by the road map proposed in [L.-M. Imbert-Gérard and G. Sylvand, Numer. Math., to appear]. For the sake of clarity, the first focus is on the two-dimensional Helmholtz equation with spatially varying wave number. The extension to a range of operators allowing for anisotropy in the first- and second-order terms follows. Numerical simulations illustrate the theoretical study of the new quasi-Trefftz functions.

Original languageEnglish (US)
Pages (from-to)1663-1686
Number of pages24
JournalSIAM Journal on Numerical Analysis
Volume59
Issue number3
DOIs
StatePublished - 2021
Externally publishedYes

Keywords

  • Best approximation properties
  • Generalized plane waves
  • Quasi-Trefftz methods

ASJC Scopus subject areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Amplitude-based generalized plane waves: New quasi-trefftz functions for scalar equations in two dimensions'. Together they form a unique fingerprint.

Cite this