Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes

Technische Samenvatting: Embedded Trefftz DG-methode voor reactie-diffusieproblemen op anisotrope roosters

Probleemstelling
Het artikel behandelt de numerieke oplossing van singular geperturbeerde reactie-diffusieproblemen van de vorm $-\varepsilon^2 \Delta u + u = f$ in een begrensd domein $\Omega \subset \mathbb{R}^2$, met Dirichlet-randvoorwaarden. De focus ligt op het regime waar $0 < \varepsilon \ll 1$, een scenario waarin oplossingen scherpe grensverlagen (boundary layers) en sterk anisotrope kenmerken vertonen. Standaard isotrope roosters zijn vaak inefficiënt voor dergelijke problemen, wat de noodzaak voor anisotrope roosterverfijning (kwadrillaat-elementen uitgerekt in specifieke richtingen om lokale kenmerken te resolveren met minder vrijheidsgraden) met zich meebrengt. Hoewel Discontinuous Galerkin (DG) methoden flexibel zijn voor hogere-orde discretisaties op algemene roosters, resulteren standaard tensorproduct DG-formuleringen in grote globale algebraïsche systemen. Trefftz-methoden, die lokale benaderingsruimten gebruiken die behoren tot de kern van de differentiaaloperator, bieden een reductie in vrijheidsgraden (DoF); echter is de constructie van geschikte Trefftz-ruimten voor algemene elliptische PDE's niet triviaal, in het bijzonder op anisotrope, gebogen kwadrillaat-elementen.

Methodologie
De auteurs stellen een Embedded Trefftz Discontinuous Galerkin (DG) methode voor die op maat is gemaakt voor anisotrope kwadrillaat-roosters. De methode wordt geconstrueerd via de volgende stappen:

1. Rooster en Ruimten: Het domein wordt verdeeld in curvilineaire kwadrillaat-elementen. Elk element $K$ wordt vanuit een eenheidvierkant getransformeerd via een affine transformatie (voor het afhandelen van anisotropie) en een $C^1$-diffeomorfisme (voor het afhandelen van kromming). De discrete ruimte bestaat uit tensorproduct-polynomen in de hulpcoördinaten.
2. Embedded Trefftz-constructie: In plaats van de globale Trefftz-ruimte expliciet te construeren, gebruikt de methode een embedding-strategie. Een lokale testruimte $L_h(K)$ wordt op elk element gedefinieerd. Deze ruimte bestaat uit tensorproduct-polynomen die nul zijn op de twee langere zijden van het pullback-element.
3. Lokale Operator en Randvoorwaarden: Een lokale reactie-diffusie-operator $A_K$ wordt gedefinieerd. De embedded Trefftz-ruimte $T_h(K)$ is de subruimte van de standaard DG-ruimte waar de residu van de PDE, getest tegen functies in $L_h(K)$, elementgewijs nul is: $(A_K v_h, q_h)_K = 0$ voor alle $q_h \in L_h(K)$.
4. Globale Formulering: De globale methode zoekt een oplossing $u_h$ die aan twee voorwaarden voldoet:
   • Lokale restricties: $(A_K u_h, v_L)_K = \ell_K(v_L)$ voor alle lokale testfuncties $v_L \in L_h(K)$.
   • Globale consistentie: De standaard interior penalty DG-bilineaire vorm $a_h(u_h, v_T) = \ell_h(v_T)$ geldt voor alle globale testfuncties $v_T$ in de Trefftz-ruimte.
5. Implementatie: In de praktijk wordt de Trefftz-ruimte niet expliciet samengesteld. In plaats daarvan maakt de methode gebruik van de Singular Value Decomposition (SVD) van de lokale operatormatrices om een embedding te construeren, waarbij alleen de globale DoF's geassocieerd met de Trefftz-ruimte worden behouden.

Belangrijkste Bijdragen
Het artikel levert de volgende specifieke technische bijdragen:

• Ontwerp van Lokale Testruimten: De primaire innovatie is het ontwerp van een lokale testruimte die compatibel is met tensorproduct-polynoombases. In tegenstelling tot eerdere werken die gebruikmaken van simplex-bases, maakt deze aanpak systematische constructie op kwadrillaat- en hexaëder-roosters mogelijk.
• Stabiliteit op Anisotrope Elementen: De auteurs bewijzen stabiliteit en quasi-optimaliteit op anisotrope, mogelijk gebogen, kwadrillaat-elementen. Een cruciaal aspect is de specifieke keuze van de testruimte (die nul is op de langere zijden), wat de afleiding van anisotrope inverse ongelijkheden en Poincaré–Steklov ongelijkheden mogelijk maakt die essentieel zijn voor stabiliteit.
• Foutanalyse: Het artikel leidt anisotrope a priori foutschattingen af in de energienorm. Deze schattingen vangen het gedrag van de oplossing in de werkelijke richtingen van anisotropie op, wat scherpere grenzen biedt dan isotrope schattingen.
• Dimensiereductie: De methode reduceert de dimensie van de lokale benaderingsruimte van $(p+1)^2$ (standaard tensorproduct) naar $2p+2$, wat het aantal globaal gekoppelde vrijheidsgraden aanzienlijk verlaagt vergeleken met standaard DG en Hybridizable DG (HDG) methoden, met name in 3D.

Resultaten
Numerieke experimenten valideren de theoretische bevindingen:

• $h$-Convergentie: Op vierkante en circulaire domeinen met grensverlagen bereikt de embedded Trefftz-methode optimale convergentiesnelheden in de energienorm. Het vereist aanzienlijk minder vrijheidsgraden om dezelfde nauwkeurigheid te bereiken als standaard DG-methoden.
• Gebogen Domeinen: De methode presteert robuust op gebogen roosters, inclusief gevallen met gemengde driehoekige en kwadrillaat-elementen (waarbij de driehoekige elementen via standaard embedding-technieken worden behandeld).
• $hp$-Convergentie: Op een negen-elementenrooster aangepast aan grensverlagen vertoont de methode exponentiële convergentie. De auteurs merken echter op dat voor hoge polynoomgraden ($p \gtrsim 5$) en grote grensverlagen de standaard embedding instabiel kan worden. Dit wordt gemitigeerd door de testruimte in hoekelementen aan te passen om de Trefftz-restrictie te versoepelen, wat de stabiliteit en exponentiële convergentie herstelt.
• Anisotrope Diffusie: De methode wordt getest op een zuiver anisotrope diffusieprobleem (zonder reactieterm), wat buiten het strikte theoretische kader valt. De resultaten laten zien dat de embedded Trefftz-methode consistent lagere fouten oplevert dan standaard DG onder zowel $h$- als $p$-verfijningsstrategieën.

Betekenis en Claims
Het artikel claimt een algemeen kader te bieden voor Trefftz-type DG-methoden die effectief blijven, zelfs wanneer de kern van de PDE-operator niet expliciet bekend of triviaal is. Door het succesvol uitbreiden van de embedded Trefftz-aanpak naar anisotrope kwadrillaat-roosters, opent dit werk de deur naar hogere-orde discretisaties op hexaëder-roosters voor gesingular de geperturbeerde problemen. De auteurs benadrukken dat de methode een competitief alternatief biedt voor HDG, waarbij het de kleinste globale systemen in drie dimensies oplevert terwijl de benaderingseigenschappen van hogere-orde discretisaties behouden blijven. De analyse is bescheiden wat betreft $p$-robuustheid, waarbij wordt erkend dat stabiliteit voor zeer hoge polynoomgraden op grove roosters met grote lagen specifieke aanpassingen aan de testruimte vereist kan hebben, zoals aangetoond in de numerieke experimenten.