The jump filter in the discontinuous Galerkin method for hyperbolic conservation laws

Technische Samenvatting: De Jump Filter in de Discontinuous Galerkin-methode voor Hyperbolische Behoudswetten

Probleemstelling
Hogere-orde lineaire numerieke schema's, in het bijzonder de Discontinuous Galerkin (DG)-methode, worden veelvuldig gebruikt voor het oplossen van hyperbolische behoudswetten vanwege hun hoge nauwkeurigheid en flexibiliteit met ongestructureerde roosters. Echter, bij het simuleren van oplossingen die discontinuïteiten bevatten (zoals schokgolven), genereren deze schema's vaak ongewenste spurious oscillaties (het Gibbs-fenomeen). Deze oscillaties kunnen leiden tot divergentie van de solver of de creatie van kunstmatige structuren. Hoewel er diverse shock-capturing technieken bestaan—zoals Total Variation Bounded (TVB) limiters, WENO-limiters en kunstmatige viscositeitsmethoden—lijden veel hiervan aan nadelen zoals hoge computationele kosten, moeilijkheden bij implementatie op ongestructureerde of gebogen roosters, de noodzaak van karakteristieke decompositie, of excessieve numerieke dissipatie die de nauwkeurigheid in gladde regio's degradeert.

Methodologie
De auteurs stellen een nieuwe jump filter voor die ontworpen is om binnen het DG-raamwerk te opereren om deze oscillaties te mitigeren, terwijl de hoge-orde nauwkeurigheid en compactheid van de methode behouden blijven.

• Lokale Viscositeitsformulering: De kern van de methode is de introductie van een lokale viscositeitsterm gebaseerd op jump-informatie bij celgrenzen. In tegen tegenstelling tot globale spectrale viscositeit, is deze aanpak gelokaliseerd. De viscositeitscoëfficiënt, $\sigma_j$, wordt adaptief bepaald door de grootte van de jumps in de oplossing en de afgeleiden over celgrenzen.
  • In gladde regio's waar jumps minimaal zijn, is de toegevoegde viscositeit verwaarloosbaar, waardoor de hoge-orde nauwkeurigheid behouden blijft.
  • In regio's met sterke discontinuïteiten neemt de viscositeit toe om de hoog-orde polynomiale modi die verantwoordelijk zijn voor oscillaties te dempen.
• Spectrale Filtering via Tijd-Splitsing: Om de stijfheid en computationele kosten te vermijden die geassocieerd worden met het impliciet oplossen van de diffusieterm (wat de tijdstap zou beperken), maken de auteurs gebruik van een tijd-splitsingsstrategie. Na elke stap van het standaard expliciete Runge-Kutta DG-schema wordt de jump filter toegepast als een post-processing stap. Dit houdt het oplossen van een eenvoudige Gewone Differentiaalvergelijking (ODE) voor de modale coëfficiënten van de polynomiale oplossing in:
  $$\frac{d}{dt} u^l_j + \sigma^l_j u^l_j = 0$$
  De oplossing wordt bijgewerkt door de coëfficiënten te vermenigvuldigen met een exponentiële vervalfactor, $e^{-\sigma^l_j \Delta t}$. Deze operatie is expliciet, computationeel goedkoop en vereist geen karakteristieke decompositie.
• Hybride Limiter Integratie: De jump filter is ook geïntegreerd in een hybride limiter-framework. In deze aanpak worden "problematische" cellen geïdentificeerd met behulp van een jump-indicator. Voor deze cellen wordt de jump filter gebruikt om een lager-orde polynomiale oplossing te berekenen, die vervolgens wordt gemengd met de hoog-orde DG-oplossing. Dit heeft als doel verdere reductie van numerieke dissipatie vergeleken met het gebruik van enkel de filter.
• Toepasbaarheid: De methode is geformuleerd voor zowel scalaire als systeem-behoudswetten (specifiek de Euler-vergelijkingen) in één en twee dimensies. Het opereert direct in de fysieke (conserveerde) variabelenruimte, waardoor de complexiteit van karakteristieke decompositie wordt vermeden.

Belangrijkste Bijdragen

1. Nieuwe Jump Filter: Introductie van een lokale viscositeitsfilter die gebruikmaakt van jump-informatie bij celgrenzen om hoog-orde modi adaptief te dempen.
2. Behoud van DG-eigenschappen: Theoretische analyse bevestigt dat het voorgestelde schema conservatie, $L^2$-stabiliteit en optimale convergentiesnelheden (hoge-orde nauwkeurigheid) voor gladde oplossingen handhaaft.
3. Computationele Efficiëntie: De methode is ontworpen voor lage computationele kosten. Het vereist geen karakteristieke decompositie, opereert volledig in de fysieke ruimte en behoudt de compactheid van de DG-methode (lokale data-communicatie), wat het zeer geschikt maakt voor parallel computing en GPU-omgevingen.
4. Hybride Integratie: Demonstratie van de naadloze integratie van de filter in een hybride limiter-framework, wat effectief oscillaties in sterke schokken onderdrukt terwijl numerieke dissipatie in gladde regio's wordt geminimaliseerd.

Numerieke Resultaten
Het artikel valideert de methode via uitgebreide numerieke experimenten op de 1D en 2D Euler-vergelijkingen, variërend van:

• Gladde Stroming: Nauwkeurigheidstests (bijv. Burgers-vergelijking, isentropische vortex) bevestigen dat de filter de convergentieorde van de standaard DG-methode op fijne roosters niet degradeert.
• Sterke Schokken: Tests met de Lax-probleem, blast waves en double Mach reflection laten zien dat de filter effectief spurious oscillaties controleert en scherpe schokprofielen vastlegt zonder de overshoots die typisch zijn voor ongefilterde hoog-orde schema's.
• Complexe Interacties: De methode presteert goed in scenario's met schok-vortex interacties, schok-bubbel interacties en Kelvin-Helmholtz instabiliteit, waarbij hoog-orde schema's met de filter complexe stromingsstructuren (zoals slip lines en turbulente menglagen) beter oplossen dan lager-orde tegenhangers.
• Robuustheid: Het schema vertoont robuustheid in hoog-Mach getal astrofysische jet-simulaties en behoudt positiviteit wanneer het wordt gecombineerd met een positivity-preserving limiter.

Betekenis en Claims
De auteurs beweren dat de jump filter een robuust, nauwkeurig en efficiënt alternatief biedt voor traditionele limiters en kunstmatige viscositeitsmethoden. De primaire betekenis ligt in de eenvoud en computationele efficiëntie:

• Het vermijdt de dure karakteristieke decompositie die vereist is door veel andere hoog-orde limiters.
• Het behoudt de compactheid van de DG-methode, wat efficiënte parallellisatie vergemakkelijkt.
• Het biedt een naadloze manier om de noodzakelijke dissipatie alleen daar te introduceren waar nodig, waarbij een balans wordt gevonden tussen shock-capturing vermogen en het behoud van hoge-orde nauwkeurigheid in gladde regio's.

Het artikel concludeert dat deze aanpak effectief het probleem aanpakt van het simuleren van hyperbolische behoudswetten met discontinue oplossingen, en een praktisch hulpmiddel biedt voor hoog-resolutie shock-capturing. Toekomstig werk richt zich op de implementatie op ongestructureerde roosters en steady-state problemen.