Stability analysis of Arbitrary-Lagrangian-Eulerian ADER-DG methods on classical and degenerate spacetime geometries

Technische Samenvatting: Stabiliteitsanalyse van Arbitrary-Lagrangian-Eulerian ADER-DG-methoden

Probleemstelling
Het Arbitrary-Lagrangian-Eulerian (ALE) raamwerk wordt breed toegepast voor het oplossen van hyperbolische partiële differentiaalvergelijkingen (PDE's) op bewegende roosters, met name wanneer topologische veranderingen optreden. In deze context verbinden Direct ALE-methoden roosters op verschillende tijdstappen via ruimtetijd-controlevolumes. Recente ontwikkelingen door Gaburro et al. hebben "gedegenereerde" ruimtetijd-elementen (sliver-elementen) geïntroduceerd om topologische veranderingen af te handelen waarbij cellen een ruimtelijke grootte van nul kunnen hebben aan het begin of einde van een tijdstap, maar wel een niet-nul ruimtetijdvolume bezitten. Hoewel ADER-DG (Arbitrary high-order DERivatives Discontinuous Galerkin) methoden effectief zijn voor deze problemen, ontbrak tot nu toe een rigoureuze theoretische stabiliteitsanalyse voor deze schema's—specifiek op gedegenereerde ruimtetijd-geometrieën. Bestaande literatuur vertrouwt vaak op empirische CFL (Courant-Friedrichs-Lewy) grenzen of beperkte stabiliteitsstudies voor lage polynomiale graden, wat een gat laat in het begrip van de stabiliteitsrestricties voor hoog-orde expliciete en impliciete ALE ADER-DG methoden op zowel klassieke als gedegenereerde geometrieën.

Methodologie
De auteurs maken gebruik van een von Neumann-stabiliteitsanalyse om de stabiliteit van expliciete en impliciete ALE ADER-DG-methoden te onderzoeken. De studie richt zich op de lineaire advectievergelijking als een benadering voor hyperbolische systemen.

1. Formulering: Het artikel stelt eerst het wiskundige raamwerk vast voor zowel expliciete als impliciete ALE ADER-DG-schema's.
   • Expliciet Schema: Maakt gebruik van een predictor-corrector-benadering. De predictor-stap construeert een lokale ruimtetijd-polynomiale benadering binnen elk controlevolume (of sliver-element) met behulp van een fixed-point iteratie. De corrector-stap werkt de oplossing bij met bewegende basisfuncties en numerieke fluxen (Rusanov-type) op interfaces.
   • Impliciet Schema: Lost de globale ruimtetijd-polynomial simultaan op over alle controlevolumes met behulp van een Newton-GMRES-iteratie, waarbij de PDE direct over de ruimtetijd-controlevolumes wordt geïntegreerd.
2. Modellering van Gedegenereerde Geometrie: Om gedegenereerde geometrieën te analyseren, introduceren de auteurs een 1D-surrogaatsetting waarin standaard interfaces tussen controlevolumes worden vervangen door "sliver-elementen". Deze elementen bootsen de gat-achtige slivers na die worden gebruikt bij 2D/3D topologische veranderingen, waarbij ze een ruimtelijke breedte van nul hebben bij $t_n$ en $t_{n+1}$, maar een niet-nul ruimtetijdvolume bezitten.
3. Stabiliteitsanalyse:
   • Klassieke Geometrieën: De auteurs leiden amplificatiematrices af voor de expliciete en impliciete schema's. Ze berekenen de spectrale radius $\rho$ van deze matrices over een bereik van CFL-getallen en fasehoeken.
   • Gedegenereerde Geometrieën: De analyse wordt uitgebreid naar sliver-elementen. Het domein wordt gemodelleerd als een periodieke blok die een standaard controlevolume en een sliver-element bevat. Een amplificatiematrix die afhankelijk is van het CFL-getal, het sliver-breedteparameter $\delta$ en de fasehoek wordt geconstrueerd.
   • Numerieke Verificatie: Theoretische stabiliteitsgrenzen worden gevalideerd door discrete amplificatiefactoren te berekenen over fijne roosters van CFL-waarden en polynomiale graden ($N=1$ tot $9$). Consistentieordes worden ook numeriek geverifieerd.

Belangrijkste Bijdragen en Resultaten

• Verfijning van Expliciete Stabiliteitsgrenzen: Voor klassieke geometrieën bevestigt de studie dat voor lage polynomiale graden ($N \le 3$), de stabiliteitslimieten overeenkomen met veelgebruikte empirische CFL-grenzen. Echter, voor hogere ordes ($N \ge 4$), demonstreert het artikel dat veelgebruikte empirische CFL-waarden in de praktijk de strikte von Neumann-stabiliteitscondities schenden (waarbij $\rho > 1$). De auteurs identificeren aanzienlijk lagere, rigoureuze CFL-limieten die vereist zijn voor strikte stabiliteit in hoog-orde expliciete schema's.
• Rooster-snelheidsrestricties: De analyse biedt een theoretische karakterisering van het toegestane bereik van rooster-snelheden voor een vast doel-CFL, waarbij wordt getoond hoe roosterbeweging de stabiliteitsrestricties beïnvloedt.
• Impliciete Onvoorwaardelijke Stabiliteit: De von Neumann-analyse en een theoretisch bewijs (Theorem 1) bevestigen dat de impliciete ALE ADER-DG-methode onvoorwaardelijk stabiel is voor de lineaire advectievergelijking op klassieke geometrieën, ongeacht de grootte van de tijdstap.
• Stabiliteit op Gedegenereerde Geometrieën:
  • Expliciet Geval: De introductie van sliver-elementen verslechtert de stabiliteit niet. De discrete amplificatiefactor voor de gedegenereerde setting wordt gevonden te zijn kleiner dan of gelijk aan die van de klassieke setting. Bijgevolg zijn dezelfde CFL-grenzen die van toepassing zijn op klassieke geometrieën ook geldig voor gedegenereerde geometrieën. De auteurs merken op dat de impliciete behandeling van de sliver-predictorstap binnen het expliciete globale schema de acceptabele CFL-limiet licht kan verhogen.
  • Impliciet Geval: De impliciete methode blijft onvoorwaardelijk stabiel, zelfs in aanwezigheid van sliver-elementen, mits de sliver-breedteparameter $\delta$ voldoet aan een specifieke geometrische restrictie ten opzichte van de tijdstap en de roostergrootte (Theorem 2).
• Consistentie: Numerieke experimenten bevestigen dat zowel expliciete als impliciete schema's hun verwachte orde van consistentie ($N+1$) behouden op zowel klassieke als gedegenereerde geometrieën.

Betekenis en Claims
Het artikel claimt het vullen van een kritiek theoretisch gat door de eerste rigoureuze von Neumann-stabiliteitsanalyse te bieden voor ALE ADER-DG-methoden op gedegenereerde ruimtetijd-geometrieën. De primaire betekenis ligt in het valideren van het gebruik van gedegenereerde elementen (slivers) voor het verbinden van bewegende roosters met topologische veranderingen. De resultaten tonen aan dat het gebruik van deze elementen met een grootte van nul geen extra stabiliteitsrestricties oplegt of de toegestane tijdstap vermindert in vergelijking met klassieke geometrieën. Dit resultaat biedt een theoretische fundering voor de praktische toepassing van directe ALE-methoden in complexe scenario's met topologische veranderingen en plaveit de weg voor de ontwikkeling van nieuwe ruimtetijd cut-cell-gebaseerde methoden. De auteurs benadrukken dat hoewel empirische CFL-waarden in de praktijk vaak worden gebruikt (gemitigeerd door viscositeit en limiters), de afgeleide rigoureuze grenzen essentieel zijn voor een volledige theoretische onderbouwing van het gedrag van de methode.