Asymptotic numerical hypocoercivity of the space-time discontinuous Galerkin method for the Kolmogorov equation

Technische Samenvatting: Asymptotische Numerieke Hypocoerciviteit van de Ruimtelijk-Temporele Discontinue Galerkin-methode voor de Kolmogorov-vergelijking

Probleemstelling
Het artikel behandelt de numerieke discretisatie van de klassieke Kolmogorov-vergelijking, een representatief model voor kinetische vergelijkingen die een degeneratieve diffusie vertonen. De vergelijking is gedefinieerd op een begrensd ruimtelijk domein $\Omega \subset \mathbb{R}^2$ en heeft de vorm $u_t + Lu = 0$, waarbij $L = -\partial_{xx} + x\partial_y$. Een kenmerkend aspect van deze operator is dat diffusie alleen plaatsvindt in de $x$-richting, terwijl de transportterm $x\partial_y$ in de $y$-richting werkt. Ondanks het gebrek aan expliciete diffusie in de $y$-richting, vertoont het systeem "hypocoerciviteit", een eigenschap waarbij de interactie tussen transport en degeneratieve dissipatie de afname over het gehele ruimtelijke domein voortplant, wat leidt tot exponentiële convergentie naar evenwicht.

De centrale uitdaging is om aan te tonen dat standaard numerieke schema's, specifiek de klassieke ruimtelijk-temporele discontinue Galerkin (dG) methode, deze hypocoercitieve structuur op het discrete niveau behouden. Hoewel hypocoerciviteit goed begrepen is in de continue setting (met name door het werk van Villani, Hérau, Nier, en anderen), is het behoud ervan in Galerkin-methoden grotendeels onverkend gebleven; eerdere numerieke resultaten waren grotendeels beperkt tot einddifferentie-, eindvolume- of niet-conforme eindelementmethoden die gespecialiseerde, niet-standaard formuleringen vereisen.

Methodologie
De auteurs analyseren de standaard ruimtelijk-temporele dG-methode toegepast op de Kolmogorov-vergelijking. De ruimtelijke discretisatie maakt gebruik van een upwind dG-schema voor de transportterm en een symmetrisch interior penalty dG-schema voor de diffusieve term. De temporele discretisatie maakt gebruik van een standaard dG-benadering met een polynoomgraad $q$.

De kern van de methodologie bestaat uit het construeren van een aangepaste stabiliteitsanalyse die de continue hypocoercitieve raamwerk nabootst. In de continue setting wordt hypocoerciviteit vaak bewezen door de inwendige product te modificeren met de toevoeging van gemengde afgeleiden of specifieke gradiënttermen, wat effectief een "sterker" norm creëert dat een spectrale kloof toelaat.

Om een vergelijkbaar resultaat op het discrete niveau te bereiken zonder de standaard bilineaire vorm van de dG-methode te wijzigen, gebruiken de auteurs een specifieke testfunctiestrategie. In plaats van de discrete vergelijking enkel te testen tegen de oplossing $U$, testen zij tegen een samengestelde functie $V = V_1 + V_2 + V_3$:

1. $V_1 = U$: De standaard energie-test.
2. $V_2 = \tau(xU_y + U_t)$: Een term die de transportoperator nabootst, geschaald door een parameter $\tau$ die afhankelijk is van de roostergrootte $h$ en de polynoomgraad $p$.
3. $V_3 = -\text{div}(A\nabla U)$: Een term die een gemodificeerde matrix $A$ betreft (afhankelijk van $h$ en $p$) en de divergentie van de gradiënt, ontworpen om controle te herstellen over de ontbrekende $y$-afgeleide.

Cruciaal is dat de auteurs observeren dat voor de discontinue polynoomruimten die in dG worden gebruikt, de term $V$ binnen dezelfde discrete ruimte als $U$ blijft. Dit maakt het mogelijk om $V$ als een toelaatbare testfunctie te gebruiken zonder de niet-standaard Petrov-Galerkin-termen te introduceren die nodig zijn in eerdere werken die continue elementen gebruikten.

Belangrijkste Bijdragen en Resultaten

1. Discrete Hypocoerciviteit voor Standaard dG: Het artikel stelt vast dat de klassieke ruimtelijk-temporele dG-methode een asymptotische numerieke hypocoerciviteit bezit. Door een eenzijdige inf-sup conditie in een "hypocoerciviteit-versterkte" norm te bewijzen, tonen de auteurs aan dat de numerieke oplossing voor grote tijden exponentieel naar evenwicht vervalt. Dit is het eerste dergelijk resultaat voor een standaard Galerkin-schema.

2. Versterkte Norm en Spectrale Kloof: De analyse introduceert een familie van normen die de volledige gradiënt van de numerieke oplossing bevatten. De auteurs bewijzen een roosterafhankelijke spectrale kloof (Poincaré-type ongelijkheid) in deze norm. De vervalsnelheid $\kappa$ hangt af van de discretisatieparameters en schaalt als $\kappa \sim h_{\min}^4 p^{-8}$. Hoewel de vervalsnelheid degradeert naarmate het rooster wordt verfijnd ($h \to 0$) of de polynoomgraad toeneemt ($p \to \infty$), zorgt de structuur voor stabiliteit en afname voor elke vaste discretisatie.

3. Inf-Sup Stabiliteit: Naast hypocoerciviteit bewijzen de auteurs dat de ruimtelijk-temporele dG-methode inf-sup stabiel is in de versterkte norm. Dit resultaat is sterker dan bestaande stabiliteitsresultaten voor ruimtelijk-temporele dG-methoden voor advectie-diffusieproblemen, omdat het controle biedt over de volledige gradiënt in plaats van enkel de advectierichting.

4. A Priori Foutbounds: Gebruikmakend van de vastgestelde inf-sup stabiliteit, leiden de auteurs optimale a priori foutbounds af voor voldoende regelmatige oplossingen in de versterkte ruimtelijk-temporele norm. De foutschatting is optimaal met betrekking tot zowel de ruimtelijke roostergrootte $h$ als de tijdstap $k$ (onder de voorwaarde $h \approx O(k)$).

5. Numerieke Validatie: Numerieke experimenten bevestigen de theoretische convergentiesnelheden en laten de exponentiële verval van de oplossing naar evenwicht zien. De experimenten tonen aan dat hoewel de theoretische vervalsnelheid afhangt van $h$ en $p$, de geobserveerde convergentie naar evenwicht robuust is.

Betekenis en Claims
De auteurs claimen dat dit werk de eerste demonstratie is van (asymptotische) numerieke hypocoerciviteit voor een bekende, klassieke Galerkin-methode. De betekenis ligt in het feit dat de methode de hypocoercitieve structuur behoudt zonder de noodzaak van de modificatie van de bilineaire vorm met extra Petrov-Galerkin-termen, die noodzakelijk waren in eerdere niet-conforme eindelementenbenaderingen.

Het artikel benadrukt dat de bewijstechniek generiek is. Door testfuncties te construeren die het effect nabootsen van de gemodificeerde inwendige producten die in de continue hypocoercitiviteitstheorie worden gebruikt, kan de aanpak potentieel worden uitgebreid naar andere kinetische type vergelijkingen en degeneratieve parabolische problemen. De resultaten bieden een rigoureuze fundering voor het gebruik van standaard hoogwaardige ruimtelijk-temporele dG-methoden voor langetermijnsimulaties van degeneratieve diffusieproblemen, waarbij wordt gewaarborgd dat de numerieke schema's de fysieke dissipatiemechanismen die inherent zijn aan het continue model niet kunstmatig onderdrukken.