Free-surface Stokes problem: stability estimates and time-step improvements

Technisches Resümee: Free-Surface-Stokes-Problem: Stabilitätsschätzungen und Verbesserungen des Zeitschritts

Problemstellung
Die genaue Simulation hochviskoser freier Oberflächenflüsse, wie etwa Eisstromdynamik, Mantelkonvektion oder Lavaströme, erfordert die Lösung der gekoppelten Stokes/Free-Surface-Gleichungen. Diese Probleme sind durch ein bewegliches Gebiet charakterisiert, in dem die Fluidgeschwindigkeit und der Druck mit der Entwicklung der Höhe der freien Oberfläche gekoppelt sind. Während Standarddiskretisierungen unter Verwendung von Finitelementen im Raum und Euler-Verfahren in der Zeit üblich sind, leiden diese oft unter schweren Stabilitätsbeschränkungen. Insbesondere explizite und semi-implizite Euler-Schemata erfordern häufig prohibitiv kleine Zeitschritte, um stabil zu bleiben, was zu hohen Rechenkosten führt. Bestehende Stabilisierungstechniken, wie der Free-Surface Stabilization Algorithm (FSSA), verbessern zwar die Zeitschrittgrößen, bieten jedoch keine a priori theoretischen Schranken für die Stabilität und garantieren weder unbedingte Stabilität noch die Volumenkonstanz.

Methodik
Die Autoren präsentieren eine rigorose theoretische Analyse der Stabilitäts- und Erhaltungseigenschaften des gekoppelten Stokes/Free-Surface-Systems sowohl für Newtonsche als auch für nicht-Newtonsche (Power-Law-Viskosität) Fluide. Die Analyse erfolgt auf zwei Ebenen:

1. Kontinuierliche Ebene: Die Autoren leiten eine totale Energiebilanz her, die die Energie der Stokes-Impulsgleichungen mit der Energie der freien Oberflächenhöhe koppelt. Diese Ableitung beruht auf der Etablierung spezifischer Beziehungen zwischen dem Advektionsterm der freien Oberfläche und Oberflächenintegralen und verknüpft diese anschließend mit dem Gravitationsterm in den Stokes-Gleichungen.
2. Diskrete Ebene: Die Studie untersucht voll diskrete Schemata unter Verwendung von Finitelementen im Raum kombiniert mit expliziten, semi-impliziten und impliziten Euler-Zeitschrittverfahren.
   • Regularisierung: Um das mangelnde inhärente Regularitätsniveau der freien Oberflächenhöhenfunktion (welches für die Wohlunterschiedlichkeit des Stokes-Problems erforderlich ist) zu adressieren, führen die Autoren einen hochgradigen Edge-Stabilisierungsterm, $J(s^{n+1}, w)$, ein, der in der Zeit implizit behandelt wird.
   • Vorschlag zur Stabilisierung: Motiviert durch die Identifizierung „spurer“ (scheinbarer) Terme in den Energieschätzungen der expliziten Euler-Methode, schlagen die Autoren einen neuen Stabilisierungsterm vor, der zur schwachen Form der Stokes-Gleichung hinzugefügt wird. Dieser Term bestraft die Normalkomponente der Geschwindigkeit an der freien Oberfläche und ist darauf ausgelegt, die destabilisierenden Terme, die aus der expliziten Behandlung der freien Oberfläche entstehen, exakt zu kompensieren.

Wesentliche Beiträge

• Totale Energiebilanz: Das Paper leitet eine neuartige totale Energiebilanz für das gekoppelte Stokes/Free-Surface-System her, welche die Impulsenergie mit der Oberflächenhöhenenergie koppelt. Dies wurde als grundlegender Baustein für das Verständnis der Stabilität identifiziert, der zuvor in der Literatur fehlte.
• Theoretische Stabilitätsbeweise:
  • Impliziter Euler: Es wurde bewiesen, dass dieser bezüglich der Zeitschrittgröße unbedingt stabil ist (unter der Annahme der Lipschitz-Stetigkeit der Oberflächenhöhe) und das Domänenvolumen konserviert.
  • Expliziter Euler (unstabilisiert): Es wurde bewiesen, dass dieser bezüglich der totalen Energienorm für jede Zeitschrittgröße instabil ist, aufgrund spurer Terme im Zusammenhang mit der normalen Oberflächengeschwindigkeit.
  • Semi-impliziter Euler (FSSA): Es wurde gezeigt, dass dieser nur bedingt stabil ist und – entscheidend – das Domänenvolumen nicht konserviert.
  • Vorgeschlagener stabilisierter expliziter Euler: Die Autoren beweisen, dass ihr neuer Stabilisierungsterm den expliziten Euler-Verfahren unbedingte Stabilität bezüglich der Zeitschrittgröße verleiht (unter der Lipschitz-Stetigkeitsannahme) und gleichzeitig das Domänenvolumen bewahrt.
• Volumenkonstanz: Das Paper liefert theoretische Beweise und numerische Verifizierungen, dass das vorgeschlagene stabilisierte explizite Schema das Domänenvolumen konserviert, während das durch FSSA stabilisierte semi-implizite Schema dies nicht tut.

Ergebnisse
Numerische Experimente wurden unter Verwendung der Open-Source-Bibliothek FEniCS sowohl für Newtonsche als auch für nicht-Newtonsche Fluidfälle durchgeführt, einschließlich einer komplexen Geometrie basierend auf Daten des Grönland-Eisschildes.

• Stabilität: Das unstabilisierte explizite Euler-Schema war für alle getesteten Zeitschritte instabil. Die FSSA-stabilisierten Schemata (sowohl semi-implizit als auch explizit) waren bedingt stabil und versagten bei größeren Zeitschritten. Im Gegensatz dazu blieb das vorgeschlagene EE-Stabilized-Schema sowohl in den Newtonschen als auch in den nicht-Newtonschen Fällen für beliebig große Zeitschritte stabil.
• Volumenkonstanz: Das EE-Stabilized-Schema konservierte das Domänenvolumen bis auf Maschinengenauigkeit. Das SIE-FSSA-Schema konnte das Volumen nicht erhalten, was die theoretischen Befunde bestätigt.
• Genauigkeit: Konvergenztests zeigten, dass die vorgeschlagenen Stabilisierungsterme konsistent mit der diskretisierten Problemstellung sind. Das EE-Stabilized-Schema behielt eine erste Ordnung der Genauigkeit für die Oberflächenhöhe und die Geschwindigkeit unter gleichzeitiger Netz- und Zeitschrittverfeinerung bei, wobei die Fehler vergleichbar mit oder kleiner als die anderer stabilisierter Schemata waren.
• Oszillationen: Der implizite Edge-Stabilisierungsterm unterdrückte hochfrequente Oszillationen in der freien Oberflächenhöhe effektiv, ohne zusätzliche Zeitschrittbeschränkungen einzuführen.

Bedeutung und Behauptungen
Das Paper behauptet, dass das vorgeschlagene stabilisierte explizite Euler-Schema eine robuste und effiziente Alternative zu rechenintensiven impliziten Lösern darstellt. Durch die theoretische Garantie der unbedingten Stabilität bezüglich der Zeitschrittgröße und die Sicherstellung der Volumenkonstanz ermöglicht die Methode stabile Simulationen von Free-Surface-Flüssen selbst bei komplexen Geometrien und Quelltermen. Die Autoren betonen, dass ihre Stabilisierung die erste ist, die theoretisch beliebige Zeitschrittgrößen für die explizite Euler-Diskretisierung dieses gekoppelten Problems garantiert. Der Ansatz zeichnet sich dadurch aus, dass er symmetrisch (vorteilhaft für iterative Löser), einfach zu implementieren und frei von abstimmbaren Parametern ist. Die Autoren kommen zu dem Schluss, dass, obwohl ihre Analyse die Lipschitz-Stetigkeit der Oberflächenhöhe voraussetzt, die abgeleiteten Stabilitätsschätzungen in der $L^2$-Norm eine notwendige Grundlage für die Wohlunterschiedlichkeit darstellen, sobald eine angemessene Regularisierung angewendet wird. Zukünftige Arbeiten sollen die Untersuchung von Regularisierungstermen für die Wohlunterschiedlichkeit sowie die Ausweitung dieser Stabilitätsanalyse auf Navier-Stokes/Free-Surface-Probleme und höhere Zeitschrittintegratoren untersuchen.