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

Technische Zusammenfassung: Stabilitätsanalyse von ALE-ADER-DG-Verfahren mit beliebigen Lagrange-Euler-Strukturen

Problemstellung
Das Arbitrary-Lagrangian-Eulerian (ALE)-Framework wird weit verbreitet eingesetzt, um hyperbolische partielle Differentialgleichungen (PDEs) auf bewegten Gittern zu lösen, insbesondere wenn Topologieänderungen auftreten. In diesem Zusammenhang verbinden Direct-ALE-Methoden Gitter zu verschiedenen Zeitschritten über Raumzeit-Kontrollvolumina. Jüngste Entwicklungen von Gaburro et al. führten „degenerierte“ Raumzeit-Elemente (Sliver-Elemente) ein, um Topologieänderungen zu handhaben, bei denen Zellen zu Beginn oder am Ende eines Zeitschritts eine räumliche Größe von Null aufweisen, aber ein nicht-verschwindendes Raumzeitvolumen besitzen. Während ADER-DG-Verfahren (Arbitrary high-order DERivatives Discontinuous Galerkin) für diese Probleme effektiv sind, fehlte bisher eine rigorose theoretische Stabilitätsanalyse für diese Schemata – speziell auf degenerierten Raumzeit-Geometrien. Die bestehende Literatur stützt sich oft auf empirische CFL-Grenzwerte (Courant-Friedrichs-Lewy) oder begrenzte Stabilitätsstudien für niedrige Polynomgrade, was eine Lücke im Verständnis der Stabilitätsbeschränkungen für hochordentliche explizite und implizite ALE-ADER-DG-Methoden sowohl auf klassischen als auch auf degenerierten Geometrien hinterlässt.

Methodik
Die Autoren verwenden eine von-Neumann-Stabilitätsanalyse, um die Stabilität von expliziten und impliziten ALE-ADER-DG-Verfahren zu untersuchen. Die Studie konzentriert sich auf die lineare Advektionsgleichung als Stellvertreter für hyperbolische Systeme.

1. Formulierung: Das Paper etabliert zunächst den mathematischen Rahmen für sowohl explizite als auch implizite ALE-ADER-DG-Schemata.
   • Explizites Schema: Verwendet einen Prädiktor-Korrektor-Ansatz. Der Prädiktor-Schritt konstruiert eine lokale Raumzeit-Polynomapproximation innerhalb jedes Kontrollvolumens (oder Sliver-Elements) mittels einer Fixpunktiteration. Der Korrektor-Schritt aktualisiert die Lösung unter Verwendung beweglicher Basisfunktionen und numerischer Flüsse (Rusanov-Typ) an den Grenzflächen.
   • Implizites Schema: Berechnet das globale Raumzeit-Polynom simultan über alle Kontrollvolumina hinweg mittels einer Newton-GMRES-Iteration, indem die PDE direkt über die Raumzeit-Kontrollvolumina integriert wird.
2. Modellierung degenerierter Geometrien: Um degenerierte Geometrien zu analysieren, führen die Autoren ein 1D-Surrogat-Setting ein, bei dem Standard-Grenzflächen zwischen Kontrollvolumina durch „Sliver-Elemente“ ersetzt werden. Diese Elemente ahmen die lochähnlichen Slivers nach, die in 2D/3D-Topologieänderungen verwendet werden, und besitzen eine räumliche Breite von Null bei $t_n$ und $t_{n+1}$, aber ein nicht-verschwindendes Raumzeitvolumen.
3. Stabilitätsanalyse:
   • Klassische Geometrien: Die Autoren leiten Verstärkungsmatrizen für die expliziten und impliziten Schemata ab. Sie berechnen den Spektralradius $\rho$ dieser Matrizen über einen Bereich von CFL-Zahlen und Phasenwinkeln.
   • Degenerierte Geometrien: Die Analyse wird durch die Einbeziehung von Sliver-Elementen auf die degenerierten Geometrien ausgeweitet. Das Gebiet wird als periodischer Block modelliert, der ein Standard-Kontrollvolumen und ein Sliver-Element enthält. Eine Verstärkungsmatrix, die vom CFL-Wert, dem Sliver-Breitenparameter $\delta$ und dem Phasenwinkel abhängt, wird konstruiert.
   • Numerische Verifizierung: Theoretische Stabilitätsgrenzen werden durch die Berechnung diskreter Verstärkungsfaktoren über feine Gitter von CFL-Werten und Polynomgraden ($N=1$ bis $9$) validiert. Die Konsistenzordnungen werden ebenfalls numerisch verifiziert.

Wesentliche Beiträge und Ergebnisse

• Verfeinerung der expliziten Stabilitätsgrenzen: Für klassische Geometrien bestätigt die Studie, dass für niedrige Polynomgrade ($N \le 3$) die Stabilitätsgrenzen mit den weit verbreiteten empirischen CFL-Grenzwerten übereinstimmen. Für höhere Ordnungen ($N \ge 4$) zeigt das Paper jedoch auf, dass häufig verwendete empirische CFL-Werte tatsächlich die strengen von-Neumann-Stabilitätsbedingungen verletzen (wo $\rho > 1$). Die Autoren identifizieren signifikant niedrigere, rigorose CFL-Limits, die für eine strikte Stabilität in hochordentlichen expliziten Schemata erforderlich sind.
• Gittergeschwindigkeitsbeschränkungen: Die Analyse liefert eine theoretische Charakterisierung des zulässigen Bereichs der Gittergeschwindigkeiten für einen festen Ziel-CFL und zeigt auf, wie die Netzbewegung die Stabilitätsbeschränkungen beeinflusst.
• Implizite unbedingte Stabilität: Die von-Neumann-Analyse und ein theoretischer Beweis (Theorem 1) bestätigen, dass das implizite ALE-ADER-DG-Verfahren für die lineare Advektionsgleichung auf klassischen Geometrien unabhängig von der Schrittweite unbedingt stabil ist.
• Stabilität auf degenerierten Geometrien:
  • Expliziter Fall: Die Einführung von Sliver-Elementen verschlechtert die Stabilität nicht. Der diskrete Verstärkungsfaktor für die degenerierte Einstellung wurde als kleiner oder gleich dem der klassischen Einstellung befunden. Folglich sind dieselben für klassische Geometrien geltenden CFL-Grenzwerte auch für die degenerierten Fälle gültig. Die Autoren merken an, dass die implizite Behandlung des Sliver-Prädiktor-Schritts innerhalb des expliziten globalen Schemas das akzeptable CFL-Limit leicht erhöhen kann.
  • Impliziter Fall: Die implizite Methode bleibt selbst in Gegenwart von Sliver-Elementen unbedingt stabil, sofern der Sliver-Breitenparameter $\delta$ eine spezifische geometrische Beschränkung relativ zum Zeitschritt und der Gittergröße erfüllt (Theorem 2).
• Konsistenz: Numerische Experimente bestätigen, dass sowohl explizite als auch implizite Schemata ihre erwartete Konsistenzordnung ($N+1$) sowohl auf klassischen als als auch auf degenerierten Geometrien beibehalten.

Bedeutung und Ansprüche
Das Paper beansprucht, eine kritische theoretische Lücke zu schließen, indem es die erste rigorose von-Neumann-Stabilitätsanalyse für ALE-ADER-DG-Methoden auf degenerierten Raumzeit-Geometrien liefert. Die primäre Bedeutung liegt in der Validierung der Verwendung degenerierter Elemente (Slivers) zur Verbindung bewegter Gitter bei Topologieänderungen. Die Ergebnisse zeigen, dass die Verwendung dieser räumlich null-großen Elemente keine zusätzlichen Stabilitätsbeschränkungen auferlegt oder die zulässigen Zeitschritte im Vergleich zu klassischen Geometrien reduziert. Dieser Befund bietet eine theoretische Grundlage für die praktische Anwendung von Direct-ALE-Methoden in komplexen Szenarien mit Topologieänderungen und ebnet den Weg für die Entwicklung neuer Methoden auf Basis von Raumzeit-Cut-Cell-Verfahren. Die Autoren betonen, dass zwar in der Praxis oft empirische CFL-Werte verwendet werden (abgemildert durch Viskosität und Limiter), die abgeleiteten rigorosen Grenzen jedoch für ein vollständiges theoretisches Verständnis des Verhaltens der Methode essenziell sind.