High-Order Invariant-Domain Preserving Continuous Finite Elements via Graph-Poisson Convex Limiting

Technische Samenvatting: Hoogwaardige Invariantie-domein Preserverende Continue Elementen via Graph-Poisson Convex Limiting

Probleemstelling
Het artikel behandelt de numerieke benadering van nietlineaire scalaire behoudswetten van de vorm $\partial_t u + \nabla \cdot f(u) = 0$. Het primaire doel is het construeren van hoogwaardige continue eindelementmethoden (FEM) die de invariantie-domein eigenschap behouden. Deze eigenschap garandeert dat als de initiële data binnen een specifiek domein $A$ ligt, de entropie-oplossing voor alle tijd in $A$ blijft. Hoewel hoogwaardige Invariantie-domein Preserverende (IDP) schema's goed gevestigd zijn voor Finite Volume (FV) en Discontinuous Galerkin (DG) methoden, zijn vergelijkbare constructies voor continue Lagrange eindelementen minder ontwikkeld.

Bestaande hoogwaardige IDP-benaderingen voor continue FEM kampen met twee significante knelpunten:

1. Problemen met Mass Lumping: Standaard laagwaardige IDP-schema's vereisen een gelumpte massamatrix met niet-negatieve elementen om het Discrete Maximum Principe (DMP) te voldoen. Echter, voor hoogwaardige Lagrange-elementen (bijv. $P_k$ met $k \ge 4$ op simplexen) kunnen standaard gelumpte massa's negatief worden, wat de directe toepassing van standaard graph-viscositeit mechanismen verhindert.
2. Stencil Breedte en Diffusie: Het uitbreiden van edge-gebaseerde laagwaardige schema's naar hoogwaardige ruimtes resulteert in excessief brede graph-stencils, wat leidt tot overmatige numerieke diffusie.

Methodologie
De voorgestelde methode ontkoppelt de hoogwaardige en laagwaardige discretisaties om deze obstakels te overwinnen, door deze te verbinden via een nieuwe flux-reconstructie en limiting procedure.

1. Ruimte-separatie:

   • Hoogwaardig Schema: De hoogwaardige oplossing wordt berekend in de gewenste polynoomruimte ($P_k$) met behulp van een Residual Viscosity (RV) stabilisatiemethode, potentieel aangevuld met Continuous Interior Penalty (CIP) stabilisatie voor zelfs even de graden.
   • Laagwaardig Schema: In plaats van een laagwaardig schema op de hoogwaardige ruimte te construeren, bouwen de auteurs een optimaal eerst-orde IDP-schema op een fijnere $P_1$ submesh. De knopen van deze submesh komen exact overeen met de nodale punten van de hoogwaardige ruimte. Dit zorgt ervoor dat het laagwaardige schema opereert op een stuksgewijze lineaire graph met strikt positieve gelumpte massa's, waardoor het negatieve massa-probleem wordt vermeden.

2. Graph-Poisson Flux Reconstructie:
   Omdat de hoogwaardige en laagwaardige updates gedefinieerd zijn met betrekking tot verschillende massastructuren (de hoogwaardige consistente massa versus de fijne $P_1$ gelumpte massa), kan hun verschil niet direct worden gerepresenteerd als conserveerde edge-fluxen.

   • De auteurs definiëren een massadefect-vector die het verschil tussen de hoogwaardige en laagwaardige updates representeert.
   • Dit defect wordt gereconstrueerd als conserveerde, antisymmetrische graph-fluxen door een Graph-Poisson probleem (een discrete Neumann-probleem) op te lossen. Deze stap brengt het verschil tussen de twee schema's effectief in kaart op de fijnere $P_1$ graph.

3. Convex Limiting en Massa Herverdeling:
   Het limiting proces volgt een tweestapsstrategie geïnspireerd door Algebraic Flux Correction (AFC), maar aangepast voor continue elementen:

   • Stap 1 (Convex Limiting): De gereconstrueerde graph-fluxen worden begrensd met convexe limiting coëfficiënten. Deze stap dwingt lokale grenzen (het invariantie-domein) af op de oplossing, maar kan een massadefect introduceren ten opzichte van de hoogwaardige target omdat de limiting wordt uitgevoerd op de fijnere $P_1$ massastructuur.
   • Stap 2 (Capaciteit-gebaseerde Massa Herverdeling): Een secundaire stap herverdeelt het resterende massadefect. Deze maakt gebruik van de "admissible capacity" van elke knoop (de afstand tot de lokale grenzen) om de fijne $P_1$ massa van de hoogwaardige update te herstellen, terwijl de invariantie-domein grenzen strikt worden gerespecteerd. Deze stap steunt op de hypothese dat de hoogwaardige target-massa binnen het toegestane massainterval ligt dat door de grenzen wordt gedefinieerd.

Belangrijkste Bijdragen

• Ruimte Ontkoppeling: Het artikel introduceert een framework waarbij het laagwaardige IDP-schema wordt gebouwd op een fijnere $P_1$ submesh geïnduceerd door hoogwaardige knopen. Dit omzeilt het negatieve gelumpte massa probleem inherent aan hoogwaardige Lagrange-bases en vermijdt de excessieve diffusie van brede stencils.
• Graph-Poisson Reconstructie: Een nieuwe methode wordt voorgesteld om het verschil tussen hoogwaardige en laagwaardige updates te representeren als conserveerde graph-fluxen door een graph-Laplace probleem op te lossen. Dit vervangt de directe algebraïsche subtractie die in traditionele FCT/AFC methoden wordt gebruikt, en accommodeert de mismatch in massamatrices.
• Tweestaps Limiting: De combinatie van convex limiting op graph-fluxen gevolgd door een capaciteit-gebaseerde massa herverdelingsstap zorgt zowel voor het behoud van invariantie-domeinen als voor globale conservatie met betrekking tot de fijnere $P_1$ massa.

Numerieke Resultaten
De methode is getest op verschillende benchmarks met behulp van de FEniCSx library op ongestructureerde meshes en polynoomgraden tot $P_4$:

• Smooth Transport: De methode vertoonde optimale convergentiecijfers in $L_1$, $L_2$, en $L_\infty$ normen voor gladde problemen, wat bevestigt dat de limiter inactief blijft in gladde regimes.
• Three-Body Rotation: De methode slaagde erin de vorm van discontinuïteiten en gladde extrema te behouden zonder spuriaanse oscillaties. Het gebruik van ontspannen grenzen (mesh-afhankelijk) bleek noodzakelijk om clipping van lokale extrema te voorkomen en de hoogwaardige nauwkeurigheid te handhaven.
• Burgers' Vergelijking: De methode convergeerde naar de correcte entropie-oplossing voor een probleem met een sonisch punt. Vergelijkingen lieten zien dat het gebruik van Residual Viscosity (RV) stabilisatie cruciaal was; methoden die enkel CIP-stabilisatie gebruikten, faalden om te convergeren naar de entropie-oplossing voor hogere polynoomgraden.
• KPP Probleem: In een uitdagende benchmark met composiet schokgolven convergeerde de gelimiteerde RV-methode naar de entropie-oplossing, terwijl de gelimiteerde CIP-alleen methode convergeerde naar een niet-entropische zwakke oplossing.

Betekenis en Claims
Het artikel claimt een robuuste, hoogwaardige, invariantie-domein preserverende continue eindelementmethode te bieden die de specifieke beperkingen van het toepassen van standaard IDP-technieken op hoogwaardige Lagrange-elementen overwint. Door de laagwaardige en hoogwaardige ruimtes te scheiden en een Graph-Poisson reconstructie te gebruiken, bereikt de methode:

1. Robuustheid: Het handelt discontinuïteiten en sonische punten af zonder de fysieke grenzen te schenden.
2. Nauwkeurigheid: Het behoudt hoogwaardige nauwkeurigheid in gladde regio's en convergeert naar de entropie-oplossing voor nietlineaire problemen.
3. Generaliteit: De aanpak is toepasbaar op willekeurige polynoomgraden en ongestructureerde meshes, waarbij de restricties van negatieve gelumpte massa's in standaard hoogwaardige bases worden omzeild.

De auteurs benadrukken dat de Graph-Poisson flux reconstructie een onderscheidend kenmerk is, waardoor de methode kan functioneren als een AFC-stijl correctie waarbij antidiffusieve fluxen worden gereconstrueerd via een PDE-gebaseerde (graph) oplossing in plaats van directe algebraïsche subtractie.