Embedded Trefftz DG method for steady Navier-Stokes flow. Part I: Oseen linearization

Technische Samenvatting: Embedded Trefftz DG-methode voor de stationaire Navier–Stokes-stroming (Deel I: Oseen-linearisatie)

Probleemstelling
Dit artikel behandelt de numerieke benadering van het stationaire Oseen-probleem, dat dient als de linearisatie van de stationaire incompressibele Navier–Stokes-vergelijkingen rond een voorgeschreven convectieveld $w$. De sturende vergelijkingen zoeken een snelheid $u$ en druk $p$ die voldoen aan:
$$-\nu\Delta u + (w \cdot \nabla)u + \nabla p = f, \quad \text{div } u = 0$$
in een begrensd domein $\Omega \subset \mathbb{R}^d$ ($d=2,3$), met homogene Dirichlet-randvoorwaarden. De primaire uitdaging die wordt aangepakt, is de ontwikkeling van een embedded Trefftz discontinuous Galerkin (DG) methode die een hoge benaderingskwaliteit behoudt met een gereduceerd aantal vrijheidsgraden (DoFs), terwijl de variabele coëfficiënten geïntroduceerd door de convectieterm $(w \cdot \nabla)u$ worden afgehandeld. In tegen tegenstelling tot het Stokes-probleem (waarbij $w=0$), bevat de Oseen-operator variabele coëfficiënten, waardoor standaard Trefftz-constructies met constante coëfficiënten ontoereikend zijn.

Methodologie
De auteurs stellen een embedded Trefftz-DG-benadering voor die de expliciete constructie van lokale Trefftz-bases vermijdt. In plaats daarvan vertrekt de methode vanuit een standaard polynoom DG-ruimte $X$ en legt zij de PDE lokaal op via projectie-restricties.

1. Lokale Operator en Trefftz-ruimte:
   Een lokale operator $K_{T,w}$ wordt gedefinieerd op elk element $T$, die de oplossingsruimte naar een duale testruimte $Q(T)'$ afbeeldt. De Trefftz-ruimte $T_w$ wordt gedefinieerd als de subruimte van $X$ waar de lokale residu zwak verdwijnt tegen $Q(T)$. Specifiek voldoen functies in $T_w$ aan de PDE alleen wanneer deze getest worden tegen lagere-orde polynomen (graad $k-2$ voor snelheid, $k-1$ voor druk).

2. Lokaal-Globale Splitsing en Complementaire Ruimten:
   Een cruciale theoretische stap is de constructie van een lokale complementaire ruimte $L(T)$ zodanig dat $X(T) = T_w(T) \oplus L(T)$.

   • Snelheidscomponent ($L_u$): Geconstrueerd met behulp van bubble-functies die verdwijnen op een sub-bal, verrijkt met specifieke polynoomtermen om facetfluxen te controleren, en beperkt om een nul-gemiddelde flux op $d$ facetten te garanderen.
   • Drukcomponent ($L_p$): Gedefinieerd als nul-gemiddelde polynomen van graad $k-1$.
   • Inverteerbaarheid: Het artikel bewijst dat de lokale operator $K_{T,w}$ begrensd inverteerbaar is op $L(T)$, mits het convectieveld $w$ aan een resolutievoorwaarde voldoet (specifiek een kleinheidsvoorwaarde op een discrete norm $|w|_{h,d}$ ten opzichte van de viscositeit). Deze inverteerbaarheid maakt het mogelijk om het lokale deel van de oplossing onafhankelijk op elk element te berekenen.

3. Gereduceerde Formulering:
   Door de hogere-orde drukcomponenten lokaal te elimineren (die via de lokale operator aan de snelheid gekoppeld zijn), leiden de auteurs een gereduceerde formulering af die uitsluitend wordt geponeerd in termen van de snelheid-onbekende. Deze formulering:

   • Verwijdert de saddle-point structuur op het lokale niveau.
   • Transfereert coerciviteit van de Stokes-setting naar de gereduceerde Oseen-Trefftz ruimten.
   • Behoudt slechts een globale koppeling voor de stuksgewijze constante drukcomponent om de discrete divergentie-restrictie af te dwingen.

4. Stabiliteitsanalyse:
   De analyse vestigt:

   • Continuïteit van de bilineaire vormen.
   • Coerciviteit van de gereduceerde vorm op de Trefftz-subruimte.
   • Inf-sup stabiliteit voor het volledige gekoppelde systeem in standaard DG-normen.
     Deze resultaten steunen op het abstracte lokaal-globale kader en de specifieke constructie van $L(T)$, wat welgesteldheid en quasi-optimaliteit waarborgt.

Belangrijkste Bijdragen
Het artikel breidt de embedded Trefftz-DG-methodologie uit van het Stokes-probleem naar de stationaire Oseen-vergelijkingen, wat de lineaire fundering biedt voor de daaropvolgende niet-lineaire Navier–Stokes-analyse (Deel II). De specifieelijke bijdragen zijn:

• Constructie van Lokale Complementaire Ruimten: In tegenstelling tot de Stokes-casus, die vertrouwde op constante-coëfficiënt Taylor-polynoomargumenten, construeert dit werk $L(T)$ om de variabele-coëfficiënt lokale operator af te handelen die voortvloeit uit de convectieterm.
• Behandeling van Inhomogene Rechterzijden: De methode incorporeert inhomogene termen direct in de lokale probleemformulering in plaats van een correctie-argument te gebruiken, wat natuurlijk past binnen het abstracte lokaal-globale splitsingskader.
• Gereduceerde Formulering en Stabiliteit: De auteurs leiden een gereduceerde snelheid-alleen formulering af en bewijzen gereduceerde inf-sup stabiliteit voor de volledige Oseen-operator. Dit elimineert de saddle-point structuur lokaal en maakt de transfer van coerciviteit mogelijk onder een geschikte kleinheidsvoorwaarde op het convectieveld.

Resultaten

• Theoretisch: Het artikel bewijst volledige stabiliteit en quasi-optimaliteit in standaard DG-normen. Het stelt vast dat de methode welgesteld is en dat de fout begrensd is door de beste benaderingsfout in de discrete ruimte, geschaald door constanten die afhangen van de resolutie van het convectieveld.
• Numeriek:
  • Oseen–Kovasznay Test: Numerieke experimenten op de voorgeschreven Oseen–Kovasznay-stroming laten zien dat de Trefftz-DG methode de verwachte asymptotische convergentiegraden bereikt ($O(h^{k+1})$ voor snelheid, $O(h^k)$ voor druk).
  • Efficiëntie: Hoewel de standaard DG-methode iets kleinere foutconstanten oplevert bij een vaste polynoomgraad, bereikt de Trefftz-DG methode een gegeven foutniveau met aanzienlijk minder vrijheidsgraden, met name voor hogere polynoomordes.
  • Computationele Kosten: Tijdmetingen geven aan dat de Trefftz-DG methode consistent sneller is dan de standaard DG-methode, waarbij het voordeel toeneemt op fijnere roosters en voor hogere orden.
  • Beperkingen: De methode vertoont gevoeligheid in het convectie-gedomineerde regime (lage viscositeit), waarbij fouten eerder verslechteren dan in de standaard DG. Dit wordt toegeschreven aan de resolutie en conditionering van de lokale embedded Trefftz-constructie.

Betekenis
Het artikel claimt de lineaire analytische grondslag te bieden die noodzakelijk is voor de niet-lineaire Navier–Stokes-theorie in Deel II. Door het succesvol uitbreiden van de embedded Trefftz-benadering naar de Oseen-vergelijkingen met variabele coëfficiënten, demonstreren de auteurs dat Trefftz-technieken kunnen worden toegepast op inhomogene, variabele-coëfficiënt operatoren zonder expliciete basisconstructie. De gereduceerde formulering wordt benadrukt als een cruciale stap voor het analyseren van het niet-lineaire probleem, aangezien het de lokale druk ontkoppelt van het globale systeem, wat de stabiliteitsanalyse voor iteratieve solvers vereenvoudigt. Het werk positioneert het Oseen-probleem als een brug tussen klassieke DG-discretisaties en Trefftz-type modelreductie, wat een pad biedt naar high-order, low-DoF simulaties voor viskeuze incompressibele stromingen.