Trefftz methods with evanescent plane waves

Titel: Trefftz-methode met evanescente vlakke golven

Auteurs: Andrea Moiola, Nicola Galante, Emile Parolin
Context: WAVES 2026, Montreal, Canada

---

1. Probleemstelling

Het paper adresseert de fundamentele numerieke instabiliteit die optreedt bij klassieke Trefftz-methoden voor het oplossen van de Helmholtz-vergelijking ($\Delta u + \kappa^2 u = 0$).

• Klassieke aanpak: Traditionele Trefftz-methoden benaderen oplossingen met behulp van propagatieve vlakke golven (PPW, Propagative Plane Waves) van de vorm $e^{i\kappa \mathbf{d}\cdot\mathbf{x}}$, waarbij $\mathbf{d}$ een reële eenheidsvector is.
• Het probleem: Hoewel deze methoden snelle convergentie beloven, lijden ze onder ernstige numerieke instabiliteiten. Dit wordt veroorzaakt door de bijna-lineaire afhankelijkheid van de basisfuncties, wat leidt tot slecht geconditioneerde lineaire systemen.
• Fundamentele beperking: Het paper citeert een stelling (gebaseerd op eerdere werken) die aantoont dat voor bepaalde Helmholtz-oplossingen (met name "evanescente modi"), een nauwkeurige benadering met PPW's alleen mogelijk is met coëfficiënten die zo groot zijn dat ze de machineprecisie overschrijden. Dit maakt een stabiele representatie in computeralgoritmen onmogelijk, ongeacht hoe geavanceerde lineaire algebra-technieken (zoals preconditionering) worden toegepast.

2. Methodologie

De auteurs stellen een nieuwe aanpak voor die de basisfuncties uitbreidt van puur propagatieve golven naar evanescente vlakke golven (EPW, Evanescent Plane Waves).

A. Evanescente Vlakke Golven (EPW)

In plaats van alleen reële propagatierichtingen, worden complexe propagatievectoren $\mathbf{d} \in \mathbb{C}^n$ toegestaan, met de voorwaarde $\mathbf{d} \cdot \mathbf{d} = 1$.

• Een EPW wordt gedefinieerd door:
  • Een propagatierichting $\mathbf{p}$ (reëel).
  • Een evanescencerichting $\mathbf{e}$ (reëel, orthogonaal op $\mathbf{p}$).
  • Een evanescentiestrakte $\eta \geq 0$.
• De functie heeft de vorm: $e^{i\kappa \zeta \mathbf{p}\cdot\mathbf{x}} e^{-\kappa \eta \mathbf{e}\cdot\mathbf{x}}$.
• Deze golven oscilleren in richting $\mathbf{p}$ met een schijnbare golfgetal $\kappa\zeta$ en vervallen exponentieel in richting $\mathbf{e}$.
• Theoretische onderbouwing: Volgens de "Herglotz-transformatie" kan elke Helmholtz-oplossing worden voorgesteld als een continue superpositie van EPW's met beperkte coëfficiënten. Dit is niet het geval voor PPW's, wat de stabiliteit van EPW's theoretisch garandeert.

B. Selectie van de Basis (Sampling Recipe)

Om een discrete ruimte te construeren binnen een Ultraweak Variational Formulation (UWVF), stellen de auteurs een specifieke strategie voor om de parameters van de EPW's te kiezen:

• 2D Strategie: Gegeven een budget van vrijheidsgraden (DOF) $P$, wordt een "Fourier-truncatieniveau" $L = P/4$ gedefinieerd.
• De parameters worden gesampled uit een uniforme verdeling over een domein dat de hoek van propagatie, de richting van evanescence en een schaalparameter voor de evanescentiestrakte omvat.
• Adaptiviteit: Als het DOF-budget laag is, worden de basisfuncties voornamelijk PPW's ($\eta=0$). Naarmate het budget groeit, neemt het aandeel van golven met een hogere evanescentiestrakte ($\eta > 0$) toe. Dit zorgt ervoor dat de methode zich aanpast aan de complexiteit van de oplossing.
• Oversampling en Regularisatie: Om numerieke instabiliteiten ondanks de keuze van de basis te voorkomen, wordt oversampling (10% extra basisfuncties) gebruikt. Daarnaast wordt de matrix $D$ (uit het UWVF-systeem) geregulariseerd via een afgeknotte SVD (Singular Value Decomposition). Kleine singuliere waarden worden genegeerd om de omkeerbaarheid van het systeem te garanderen zonder de nauwkeurigheid te verliezen.

3. Belangrijkste Bijdragen

1. Identificatie van de oorzaak: Het paper verduidelijkt dat de instabiliteit van PPW-methoden fundamenteel is en niet opgelost kan worden door alleen lineaire algebra-technieken; de basisfuncties zelf moeten worden gewijzigd.
2. EPW-Selectiealgoritme: Een eenvoudige, reproduceerbare "recept" (recipe) voor het selecteren van een discrete ruimte van evanescente vlakke golven in 2D (en een schets voor 3D) die de kwaliteit van de benadering maximaliseert.
3. Integratie met UWVF: De toepassing van deze EPW-basis binnen de Ultraweak Variational Formulation, inclusief de noodzakelijke regularisatiestappen voor een stabiele oplossing.
4. Theoretisch bewijs: Verwijzing naar bestaande theorema's die aantonen dat EPW's een stabiele representatie mogelijk maken voor een breder scala aan Helmholtz-oplossingen dan PPW's.

4. Resultaten (Numerieke Experimenten)

De auteurs vergelijken de prestaties van PPW's en EPW's in drie scenario's:

• Puntbron (Fundamentele oplossing):
  • Bij een moderate DOF-grootte convergeren beide methoden snel.
  • Bij toenemende DOF stagneert de PPW-convergentie door numerieke instabiliteit (fouten blijven hoog).
  • De EPW-methode blijft convergeren en bereikt fouten die $10^6$ keer kleiner zijn dan die van de PPW-methode.
• Invloed van het golfgetal ($\kappa$):
  • Wanneer de grootte van de benaderingsruimte lineair toeneemt met $\kappa$, verbetert de nauwkeurigheid van de EPW-methode naarmate $\kappa$ toeneemt.
  • De PPW-methode faalt bij hogere waarden van $\kappa$.
• Verstrooiingsprobleem (Scattering):
  • Een niet-convexe verstrooier met hoekpunten (waar singulariteiten optreden) wordt getest.
  • Ook hier stagneert de PPW-methode, terwijl de EPW-methode een aanzienlijk lagere fout levert en de complexe golffronten nauwkeurig kan oplossen.

5. Betekenis en Conclusie

Dit paper biedt een oplossing voor een decennialang bestaand probleem in de numerieke golfmechanica. Door de beperking tot propagatieve golven los te laten en evanescente golven te introduceren, wordt de numerieke stabiliteit van Trefftz-methoden drastisch verbeterd.

De voorgestelde methode maakt het mogelijk om:

• Hogere nauwkeurigheid te bereiken bij hoge golfgetallen.
• Complexere geometrieën en singulariteiten te behandelen zonder dat de oplossing instabiel wordt.
• Een robuust algoritme te bieden dat minder afhankelijk is van de precieze keuze van de mesh of de lineaire algebra-preconditioners.

De resultaten suggereren dat EPW's de nieuwe standaard moeten worden voor Trefftz-discontinuous Galerkin-methoden, vooral in toepassingen waar hoge precisie en stabiliteit cruciaal zijn.