A discontinuous Galerkin method for elliptic-hyperbolic equations

Titel: Een Discontinuous Galerkin-methode voor elliptisch-hyperbolische vergelijkingen

Auteurs: Chiara Perinati, Lise-Marie Imbert-Gérard, Andrea Moiola, Paul Stocker.

1. Probleemstelling

Het artikel richt zich op de numerieke oplossing van een klasse van lineaire partiële differentiaalvergelijkingen (PDV's) van gemengde orde. Specifiek wordt de Frankl-operator onderzocht:
$$Lu := K u_{xx} + u_{yy} = f$$
waarbij de coëfficiënt $K = K(y)$ van teken verandert binnen het rekenkundige domein $\Omega$.

• Elliptisch gebied: Waar $y > 0$ (en $K(y) > 0$).
• Hyperbolisch gebied: Waar $y < 0$ (en $K(y) < 0$).
• Parabolische kromme: De lijn $y=0$ waar het type van de vergelijking overgaat.

Een klassiek voorbeeld is de Tricomi-vergelijking ($K(y)=y$). Het probleem wordt geformuleerd als een randwaardeprobleem (het Tricomi-probleem) met Dirichlet-randvoorwaarden op een deel van de rand ($\Gamma_0$ in het elliptische gebied en een karakteristieke kromme $\Gamma_1$ in het hyperbolische gebied). Dit wordt beschouwd als een "open" probleem, wat in tegenstelling tot "gesloten" problemen (waarvoor randvoorwaarden op de volledige rand nodig zijn) goed gesteld is.

2. Methodologie

De auteurs stellen een Discontinuous Galerkin (DG) methode voor die direct de tweede-orde vergelijking benadert, in plaats van deze eerst om te vormen naar een stelsel van eerste-orde vergelijkingen.

Kernaspecten van de formulering:

• Morawetz Multiplier Techniek: De variatiele formulering is gebaseerd op de "abc-methode" van Friedrichs. De vergelijking wordt vermenigvuldigd met een geschikte testfunctie (de Morawetz multiplier) van de vorm $Mv = b v_x + c v_y$, waarbij $b$ en $c$ specifieke functies zijn die voldoen aan voorwaarden voor stabiliteit en coerciviteit.
• DG Formule: De methode gebruikt element-wise polynomen (of Trefftz-achtige ruimten) die niet noodzakelijk continu zijn over de elementgrenzen.
  • Stabilisatie: Om coerciviteit te garanderen, worden straff termen toegevoegd die de sprongen (jumps) van de oplossing en haar eerste afgeleiden over de facetten straffen.
  • Randvoorwaarden: De randvoorwaarden worden zwak opgelegd via de multiplier-techniek en de DG-fluxen.
• Ruimten: De methode wordt geanalyseerd voor drie soorten discrete ruimten:
  1. Standaard stuksgewijze polynomen ($P_p$).
  2. Quasi-Trefftz ruimten: Ruimten opgespannen door polynomen die lokaal de homogene vergelijking benaderen (via Taylor-ontwikkeling).
  3. Embedded Trefftz ruimten: Ruimten waarbij de Trefftz-eigenschap zwak wordt afgedwongen via projectie binnen de standaard polynoomruimte.

3. Belangrijkste Bijdragen

1. Well-posedness (Goed gesteldheid):

   • Er wordt bewezen dat het discrete probleem goed gesteld is.
   • Coerciviteit wordt bewezen in een energie-norm, afhankelijk van de mesh, met behulp van de Morawetz multiplier.
   • Continuïteit wordt bewezen in een sterkere "residu"-norm.
   • De methode is stabiel voor algemene polygonale meshes en willekeurige polynoomgraden, mits de straffparameters ($\gamma_1, \gamma_2, \gamma_3$) voldoende groot zijn.

2. Foutanalyse en Convergentie:

   • Er worden a priori foutbounds afgeleid in de energie-norm.
   • Voor standaard polynomen worden hp-foutschattingen verkregen. De convergentie in de energie-norm is van de orde $O(h^{p-1})$ (suboptimaal ten opzichte van de theoretische $O(h^p)$ voor standaard elliptische problemen, wat typisch is voor gemengde type problemen).
   • Voor Quasi-Trefftz ruimten worden $h$-foutschattingen afgeleid die profiteren van de hoge nauwkeurigheid van deze ruimten voor gladde oplossingen.

3. Efficiëntie van Trefftz-methoden:

   • De auteurs tonen aan dat Quasi-Trefftz en Embedded Trefftz ruimten een aanzienlijke reductie in het aantal vrijheidsgraden mogelijk maken (orde $O(p)$ in plaats van $O(p^2)$ voor standaard polynomen) terwijl ze vergelijkbare nauwkeurigheid behouden.

4. Numerieke Resultaten

De theorie wordt gevalideerd met numerieke experimenten voor de Tricomi-vergelijking ($K(y)=y$).

• Convergentie: De numerieke resultaten bevestigen de theoretische convergentieorden.
  • In de energie-norm wordt een convergentie van $O(h^{p-1})$ waargenomen voor alle ruimten.
  • In de $L^2$-norm wordt vaak een hogere orde waargenomen (bijv. $O(h^p)$ of beter).
• Vergelijking Ruimten:
  • Standaard polynomen leveren de beste absolute nauwkeurigheid voor een gegeven $h$ en $p$, maar vereisen meer vrijheidsgraden.
  • Quasi-Trefftz en Embedded Trefftz bereiken vergelijkbare nauwkeurigheid met veel minder vrijheidsgraden. Ze tonen zelfs exponentiële convergentie ($O(e^{-AN_{dofs}})$) wanneer de polynoomgraad $p$ wordt verhoogd, terwijl standaard polynomen slechts wortel-exponentiële convergentie tonen.
• Robuustheid: De methode is robuust ten opzichte van de keuze van de straffparameters. Hoewel er een theoretische ondergrens is voor stabiliteit, werkt de methode goed over een breed bereik van waarden.

5. Betekenis en Conclusie

Dit werk biedt een solide theoretische en numerieke basis voor het oplossen van elliptisch-hyperbolische problemen met DG-methoden.

• Innovatie: De integratie van Morawetz multipliers direct in de DG-formulering voor tweede-orde vergelijkingen is een belangrijke stap, omdat dit de stabiliteit garandeert zonder de noodzaak van $C^1$-continuïteit (wat vaak nodig is bij conformerende elementen voor dit type probleem).
• Toepassing: De methode is relevant voor toepassingen in transsonische gasstromen, waar de overgang van subsonisch naar supersonisch stromen precies dit type gemengde vergelijking vereist.
• Efficiëntie: De succesvolle toepassing van Trefftz-achtige methoden (Quasi en Embedded) op dit complexe probleem toont aan dat het mogelijk is om de rekenkosten drastisch te verlagen zonder in te boeten aan nauwkeurigheid, wat cruciaal is voor complexe simulaties in de natuurkunde en techniek.

Samenvattend presenteert dit artikel een geavanceerde, stabiele en efficiënte numerieke methode voor een klassiek maar uitdagend probleem in de wiskundige fysica.