The error estimate of entropy-stable discontinuous Galerkin methods for hyperbolic conservation laws

Technische Samenvatting: Foutschattingen van Entropie-stabiele Discontinuous Galerkin-methoden voor Hyperbolische Conserveringswetten

Probleemstelling
Hyperbolische conserveringswetten beheersen fysische verschijnselen zoals gasdynamica, waarbij oplossingen vaak in eindige tijd schokken en contactdiscontinuïteiten ontwikkelen. Om de wel-gestelheid van deze problemen te waarborgen, zijn entropie-ongelijkheden vereist om de fysisch toelaatbare zwakke oplossing te selecteren. Hoewel Chen en Shu [J. Comput. Phys. 345 (2017)] een verenigd kader hebben vastgesteld voor de constructie van hoogwaardige entropie-stabiele discontinuous Galerkin (ESDG) methoden met behulp van specifieke numerieke kwadratuur en summation-by-parts (SBP) operatoren, ontbrak een rigoureuze convergentietheorie. Bestaande foutanalyses voor standaard DG-methoden steunen op de aanname dat de numerieke oplossing zich bevindt in een voorgeschreven ruimte van stuksgewijze polynomen. Echter, voor de nodale ESDG-methode komen de geproduceerde nodale waarden van het schema over het algemeen niet overeen met een stuksgewijze polynoom van graad $k$ in meerdere ruimtelijke dimensies. Bijgevolg zijn klassieke foutschattingen niet direct toepasbaar, wat een gat laat in de theoretische rechtvaardiging voor de convergentie van deze entropie-stabiele schema's.

Methodologie
Dit artikel stelt rigoureuze a priori foutschattingen vast voor semi-discrete ESDG-methoden op algemene ongestructureerde roosters. De analyse maakt gebruik van een eindige-verschillen-type consistentie–stabiliteit argument dat direct op nodaal niveau is geformuleerd, in plaats van te vertrouwen op standaard variationele argumenten.

De kernmethodologie omvat:

1. Polynoomreconstructie-hypothese: De auteurs nemen aan dat voor elk element de nodale vector een unieke polynoomreconstructie toelaat in een eindige-dimensionale ruimte $V_h(K)$ die polynomen van graad $k$ bevat. Dit maakt het gebruik van standaard inverse schattingen en benaderingstheorie mogelijk.
2. $L^\infty$ A Priori Begrenzing: De analyse neemt een begrenzing aan van de fout in de $L^\infty$-norm, die wordt gerechtvaardigd via een standaard continuïteitsargument (continuatiemethode).
3. Consistentie-Stabiliteitsframework:
   • Consistentie: Een lokale afkortingsfout wordt afgeleid door de exacte oplossing's nodale waarden in het semi-discrete schema in te vullen.
   • Stabiliteit: Een foutvergelijking wordt geformuleerd door het numerieke schema af te trekken van de evolutie van de exacte oplossing. De fout wordt gemeten in een kwadratuur-gebaseerde norm die equivalent is aan de gebroken $L^2$-norm op de reconstructieruimte.
   • Energie-schattingen: De foutvergelijking wordt vermenigvuldigd met de entropie-variabele fout en geïntegreerd. De volumetermen worden gecontroleerd met behulp van de skew-symmetrische SBP-structuur gecombineerd met entropie-symmetrisatie van de flux-Jacobianen. Interface-bijdragen worden beheerd via de entropie-stabiliteit van de numerieke flux. Niet-lineaire termen worden geschat met behulp van Taylor-expansies en mesh-geschaalde matrix-norm bounds.
   • Afsluiting: De resulterende differentiële ongelijkheid wordt gesloten met behulp van een Gronwall-argument.

Het framework wordt uitgebreid naar de Entropy-Stable Oscillation-Free DG (ESOFDG) methode, die extra dempingstermen bevat om spontane oscillaties te onderdrukken. De auteurs bewijzen dat deze dempingstermen de convergentieorde niet verslechteren.

Belangrijkste Bijdragen

• Rigoureuze Convergentietheorie: Het artikel levert de eerste rigoureuze a priori foutschatting voor de semi-discrete ESDG-methode geïntroduceerd door Chen en Shu [6], waarbij een $O(h^k)$ foutschatting in een kwadratuur-gebaseerde norm wordt bewezen voor gladde oplossingen van zowel scalaire conserveringswetten als systemen.
• Oplossing van de Polynoom-aanname Kloof: Door de polynoomreconstructie-hypothese te introduceren, overbruggen de auteurs de kloof tussen de nodale flux-differencing formulering (die de standaard stuksgewijze polynoom-aanname schendt) en klassieke foutanalyse technieken.
• Extensie naar ESOFDG: De analyse wordt succesvol uitgebreid naar de ESOFDG-methode, waarbij wordt aangetoond dat de aanvullende dempingstermen, ontworpen om oscillaties te beheersen, de theoretische convergentieorde behouden.
• Technische Lemma's: Het artikel stelt noodzakelijke technische resultaten vast met betrekking tot de schalingseigenschappen van SBP-operatoren, lokale afkortingsfouten en bounds op de dempingscoëfficiënten voor het ESOFDG-schema.

Resultaten

• Theoretische Bounds: Onder de aannames van polynoomreconstructie en een $L^\infty$ a priori begrenzing, bewijst het artikel dat de fout in de kwadratuur-gebaseerde norm schaalt als $O(h^k)$ voor polynoomgraad $k$. Dit resultaat geldt voor zowel scalaire vergelijkingen als systemen van conserveringswetten.
• Numerieke Experimenten: Numerieke tests op de Burgers-vergelijking (1D en 2D) en de 2D Euler-vergelijkingen (vortex evolutie) bevestigen de theoretische bevindingen.
  • De geobserveerde convergentiegraden komen over het algemeen overeen met of overtreffen zelfs licht de theoretische $O(h^k)$ voorspelling.
  • In sommige gevallen (bijv. 1D Burgers met $k=3$) benaderen de geobserveerde graden $O(h^{k+1})$, wat suggereert dat de theoretische bound mogelijk niet scherp is.
  • Het ESOFDG-schema vertoont convergentiegraden die ononderscheidbaar zijn van het standaard ESDG-schema, wat bevestigt dat de dempingstermen de nauwkeurigheid niet opofferen.

Betekenis en Claims
De auteurs stellen dat dit werk een kritieke kloof in de convergentietheorie voor entropie-stabiele DG-schema's dicht. Door een rigoureuze $O(h^k)$ foutschatting vast te stellen, biedt het artikel een theoretische fundering voor de betrouwbaarheid van deze hoogwaardige methoden bij het vastleggen van fysisch toelaatbare oplossingen.

Het artikel erkent bescheiden dat hoewel de numerieke experimenten suggereren dat de geobserveerde graden de theoretische bound met wel een halve orde kunnen overtreffen (wat potentieel duidt op superconvergentie of delicate cancellatie-mechanismen), de huidige analyse een solide, zij het mogelijk niet-scherpe, ondergrens biedt. De auteurs merken op dat de vraag of een scherpere $(k + 1/2)$-de orde schatting binnen het eindige elementen variationele kader kan worden afgeleid, een open vraag blijft. De primaire bijdrage is de vestiging van een rigoureuze convergentietheorie voor een klasse van schema's die voorheen zonder dergelijke rechtvaardiging bleven.