An Interface Green's Function Framework for Complete Discrete $W^{1,\infty}$ Analysis of Discontinuous Galerkin Methods

Technische Samenvatting: Een Interface Green's Functie Framework voor Volledige Discrete $W^{1,\infty}$ Analyse van Discontinuous Galerkin Methoden

Probleemstelling
Puntgewijze foutanalyse van Discontinuous Galerkin (DG) methoden is een belangrijk onderzoeksgebied geweest, waarbij de bestaande literatuur optimale foutschattingen heeft vastgesteld in de $L^\infty$-norm en de gebroken $W^{1,\infty}$-seminorm. Echter, er blijft een kritieke theoretische kloof bestaan: de gebroken $W^{1,\infty}$-seminorm, gedefinieerd als $\|\nabla(u - u_h)\|_{L^\infty_h(\Omega)}$, slaagt er niet in om discontinuïteiten over mesh-interfaces te controleren. Bijgevolg bieden bestaande analyses geen directe controle over interface-sprongen, die inherent zijn aan discontinue benaderingen en essentieel zijn voor stabiliteit, conservering en post-processing. Bovendien is het op convexe polyhedrale domeinen, waar zwakke discrete maximumprincipes momenteel niet beschikbaar zijn voor DG-methoden, het afleiden van optimale schattingen zonder suboptimale logaritmische factoren (bijv. $|\ln h|^\beta)$ een hardnekkige uitdaging gebleken, met name voor polynoomgraden groter dan één.

Methodologie
De auteurs stellen een nieuw interface Green's functie framework voor om een volledige discrete $W^{1,\infty}$ analyse te bereiken. De kern van de methodologische innovatie wijkt af van klassieke Green's functie-analyses, die doorgaans een enkele discrete deltafunctie gebruiken om de oplossing of haar gradiënt te schatten. In plaats daarvan introduceert dit framework interface Green's functies gegenereerd door het verschil tussen aangrenzende discrete deltafuncties.

1. Definitie van Interface Green's Functies:

   • Voor interne mesh-interfaces definieert het framework een Green's functie $\Gamma_{if}$ die voldoet aan $-\Delta \Gamma_{if} = \delta_{h, K_+^*} - \delta_{h, K_-^*}$, waarbij $\delta_{h, K}$ een gladde discrete deltafunctie is ondersteund op element $K$. Deze bronterm representeert van nature de sprong van discontinue benaderingen over de interface.
   • Voor randvlakken (boundary faces) wordt een vergelijkbare Green's functie $\Gamma_{bf}$ gedefinieerd met behulp van een enkele discrete deltafunctie.

2. Representatieformules:
   De sprongen van de numerieke oplossingen (en hun gradiënten) worden gerepresenteerd met behruik van deze interface Green's functies en hun DG-benaderingen. Dit reduceert de schatting van interface-hoeveelheden tot de benaderingsfout van de interface Green's functies zelf.

3. Cruciaal Analytisch Onderdeel: Lokale Energie-schattingen:
   Het bewijs steunt op een nieuwe lokale energie-schatting voor deze interface Green's functies, vastgesteld via dyadische decompositie van het domein.

   • De auteurs maken gebruik van een Taylor-expansie argument met betrekking tot de bronvariabele.
   • Cruciaal is dat het verschil tussen naburige discrete deltafuncties ($\delta_{h, K_+^*} - \delta_{h, K_-^*}$) een annulerend effect creëert op het niveau van de singuliere bron.
   • Deze annulering elimineert de dominante singuliere bijdrage, wat leidt tot scherpere lokale schattingen (specifiek met een extra factor $h$) vergeleken met klassieke schattingen.
   • Deze verbetering maakt het mogelijk om interface-schattingen af te leiden zonder het introduceren van logaritmische factoren die doorgaans voortvloeien uit de directe toepassing van discrete trace-inegaliteiten gecombineerd met bestaande maximum-norm schattingen.

Belangrijkste Bijdragen en Resultaten
Het paper vestigt een "volledige" discrete $W^{1,\infty}$ theorie, gedefinieerd door de norm:
$$|||v|||_{W^{1,\infty}_h(\Omega)} = \max_{K \in \mathcal{T}_h} \|\nabla v\|_{L^\infty(K)} + \max_{F \in \mathcal{E}_h} h_F^{-1} \|[v]\|_{L^\infty(F)}$$
Deze norm meet simultaan de gebroken gradiënt en de geschaalde interface-sprongen.

Het framework wordt toegepast op de volgende methoden op convexe polyhedrale domeinen voor de Poisson-vergelijking:

• Symmetric Interior Penalty DG (IPDG) methode.
• Hybrid Mixed DG methode.
• Twee representatieve Hybridizable DG (HDG) methoden (HDG-1 en HDG-2).

Belangrijkste Theoretische Resultaten:

• Optimale Jump-schattingen: De auteurs bewijzen optimale maximum-norm schattingen voor de interface-sprongen $[u - u_h]$ en numerieke fluxen, waarbij de logaritmische factor $\beta$ die voorheen voorkwam in schattingen op convexe polyhedrale domeinen wordt verwijderd.
• Volledige Foutgrenzen: Ze leiden optimale foutschattingen voor de gebroken gradiënten en interface-sprongen simultaan af, waarmee de puntgewijze analyse voor deze methoden wordt voltooid.
• Extensie naar Vierde-orde Problemen: Het framework wordt uitgebreid naar de $C^0$ interior penalty methode voor de biharmonische vergelijking. De auteurs leiden nieuwe maximum-norm schattingen af voor de sprongen van de normale afgeleide, een hoeveelheid die voorheen onbehandeld bleef in de puntgewijze analyse voor deze methode.

Betekenis en Claims
Het paper claimt dat de primaire bijdrage methodologisch is. Door de introductie van interface Green's functies gebaseerd op het verschil van discrete deltafuncties, overwinnen de auteurs de fundamentele barrière van het gebrek aan zwakke discrete maximumprincipes voor DG-methoden op convexe polyhedrale domeinen.

• Verwijdering van Logaritmische Factoren: Het framework elimineert succesvol de suboptimale logaritmische factoren die overgenomen zijn van klassieke trace-inegaliteit benaderingen, wat optimale schattingen biedt voor polynoomgraden groter dan één.
• Generaliseerbaarheid: De constructie is afhankelijk van de lokale discontinue structuur van de benadering in plaats van de specifieke numerieke formulering. De auteurs stellen dat de methodologie voldoende algemeen is om toepasbaar te zijn op een brede klasse van discontinue en niet-conforme eindige element methoden.
• Volledigheid: Het werk biedt de eerste volledige discrete $W^{1,\infty}$ analyse die simultaan de gebroken gradiënten en interface-sprongen controleert voor de gespecificeerde DG-methoden, waarmee een kritieke kloof in het theoretische begrip van puntgewijze benaderingen voor discontinue methoden wordt opgevuld.