A HHO formulation for variable density incompressible flows where the density is purely advected

Diese Arbeit stellt ein Hybrid High-Order (HHO)-Verfahren für inkompressible Strömungen mit variabler Dichte vor, das eine exakte Volumenerhaltung und reine Dichte-Advektion gewährleistet und sich durch Druckrobustheit, hohe zeitliche Genauigkeit sowie effiziente Lösungsstrategien für die Simulation von Mehrphasenströmungen und Rayleigh-Taylor-Instabilitäten auszeichnet.

Das Wesentliche

Titel: Wie man flüssige Mischungen mit einem digitalen „Schneidbrett" simuliert

Stellen Sie sich vor, Sie gießen Honig und Wasser in ein Glas. Das Wasser ist leicht, der Honig schwer. Wenn Sie das Glas schütteln, vermischen sie sich nicht wirklich, aber sie wirbeln durcheinander. Die schwerere Flüssigkeit will nach unten, die leichtere nach oben. Das ist eine Rayleigh-Taylor-Instabilität – ein physikalisches Phänomen, das auch in Sternen oder bei Explosionen passiert.

Die Herausforderung für Computer ist: Wie berechnet man das genau, ohne dass die Simulation „verrücktspielt" oder zu viel Rechenzeit braucht?

Hier kommt die neue Methode vor, die Lorenzo Botti und Francesco Carlo Massa entwickelt haben. Sie nennen sie HHO-ESDIRK. Klingt kompliziert? Lassen Sie uns das mit einfachen Bildern erklären.

1. Das Problem: Der „undichte" Eimer

Bei herkömmlichen Computer-Simulationen von Flüssigkeiten passiert oft ein Problem: Der Computer vergisst manchmal, dass Wasser (oder Honig) nicht komprimierbar ist. Es ist, als würde man Wasser in einen Eimer füllen, der kleine Löcher hat. Irgendwann ist das Wasser weg oder der Eimer platzt, weil der Druck zu hoch wird.

Außerdem: Wenn man zwei Flüssigkeiten mischt, darf die Dichte (wie schwer sie ist) nicht plötzlich negativ werden oder unendlich groß. Das wäre physikalisch unmöglich.

2. Die Lösung: Ein intelligentes „Schneidbrett" (HHO)

Die Autoren haben eine neue Art entwickelt, den Raum in kleine Stücke zu zerlegen, ähnlich wie man ein Brot in Scheiben schneidet. Aber statt nur die Scheiben zu betrachten, schauen sie sich auch die Kanten zwischen den Scheiben an.

• Die Hybrid-Methode (HHO): Stellen Sie sich vor, Sie haben ein Puzzle. Bei alten Methoden musste jedes Puzzleteil perfekt passen, um das Bild zu zeigen. Bei dieser neuen Methode (HHO) können die Teile etwas „schief" sein, aber die Kanten zwischen ihnen passen perfekt zusammen.
• Der Vorteil: Das erlaubt dem Computer, viel komplexere Formen zu berechnen, ohne dass das Bild (die Flüssigkeit) verzerrt wird. Es ist wie ein flexibles Gummiband, das sich an jede Form anpasst, aber trotzdem straff bleibt.

3. Die Magie: Der „perfekte" Druck

Ein großes Problem bei Flüssigkeiten ist der Druck. Wenn der Computer den Druck falsch berechnet, fließt die Flüssigkeit in die falsche Richtung.

• Druck-Robustheit: Die neue Methode ist wie ein erfahrener Koch, der den Druck genau spürt. Selbst wenn die Rechnung des Drucks nicht 100 % perfekt ist, fließt die Flüssigkeit trotzdem genau dorthin, wo sie soll. Der Druck stört die Bewegung der Flüssigkeit nicht.
• Volumen-Erhaltung: Die Methode garantiert, dass das Volumen der Flüssigkeit zu 100 % erhalten bleibt. Es ist, als würde man Wasser in einem undurchlässigen Behälter bewegen: Nicht ein Tropfen geht verloren.

4. Die Zeitreise: Der „Schritt-für-Schritt"-Tanz (ESDIRK)

Flüssigkeiten bewegen sich schnell. Um das zu simulieren, muss der Computer viele kleine Zeitschritte machen.

• Der Tanz: Die Autoren nutzen eine spezielle Tanzmethode (ESDIRK), bei der der Computer nicht nur einen Schritt macht, sondern mehrere kleine „Probe-Schritte" innerhalb eines großen Schrittes macht, um sicherzugehen, dass er nicht stolpert.
• Warum wichtig? Das ist besonders wichtig, wenn sich die Flüssigkeit sehr schnell bewegt oder wenn die Dichte stark schwankt (wie bei Honig und Wasser).

5. Der Clou: Das „Geheimnis" der Rechenleistung

Normalerweise braucht eine hochgenaue Simulation einen riesigen Computer, weil er so viele Daten speichern muss.

• Statistische Kondensation: Die Autoren haben einen Trick angewendet. Sie sagen dem Computer: „Berechne nur die Kanten des Puzzles genau, den Rest kann ich mir im Kopf vorstellen."
• Das Ergebnis: Der Computer braucht viel weniger Speicherplatz und ist viel schneller, obwohl die Ergebnisse genauso genau sind wie bei den großen, langsamen Methoden. Es ist wie beim Packen eines Koffers: Man packt nur das Nötigste ein und lässt den Rest zu Hause, aber man hat trotzdem alles, was man braucht.

6. Der Test: Der Honig-Wasser-Tanz

Um zu beweisen, dass ihre Methode funktioniert, haben sie den Rayleigh-Taylor-Effekt simuliert:

• Szenario: Schwere Flüssigkeit oben, leichte unten. Die schwere will runter, die leichte hoch. Es entstehen schöne, spiralförmige Wirbel.
• Ergebnis: Ihre Methode konnte diese Wirbel auch bei sehr hohen Geschwindigkeiten (hohe Reynolds-Zahlen) und bei sehr unterschiedlichen Dichten (wie Honig und Wasser) genau abbilden.
• Besonderheit: Selbst bei sehr groben „Scheiben" (wenige Rechenpunkte) lieferte die Methode genaue Ergebnisse, solange sie die „Kanten" (die Hybrid-Technik) richtig nutzte.

Fazit

Diese neue Methode ist wie ein super-intelligenter, sparsamer und genauer Koch, der komplexe Flüssigkeitsmischungen simulieren kann. Sie ist schnell, braucht wenig Speicherplatz und macht keine Fehler beim Volumen oder der Dichte.

Das ist ein großer Schritt für die Zukunft, um zum Beispiel zu verstehen, wie sich Öl und Wasser in Pipelines verhalten, wie sich Wolken bilden oder wie sich Explosionen in der Atmosphäre ausbreiten – alles mit weniger Rechenleistung und mehr Genauigkeit als je zuvor.

Technische Zusammenfassung

Hier ist eine detaillierte technische Zusammenfassung des Papers „A HHO formulation for variable density incompressible flows where the density is purely advected" von Botti und Massa auf Deutsch.

1. Problemstellung

Das Paper adressiert die numerische Simulation von inkompressiblen Strömungen mit variabler Dichte, insbesondere für Mischungen nichtmischbarer (immiscible) Fluide. Die zugrunde liegenden Gleichungen sind die Navier-Stokes-Gleichungen, bei denen die Dichte $\rho$ als rein advektive Größe behandelt wird (d.h. sie wird nur durch das Geschwindigkeitsfeld transportiert, ohne Diffusionsterme).

Die Hauptherausforderungen bei solchen Problemen sind:

• Erhaltung des Volumens: Da die Strömung inkompressibel ist, muss die Divergenzfreiheit des Geschwindigkeitsfeldes exakt eingehalten werden, um eine physikalisch korrekte reine Advektion der Dichte zu gewährleisten.
• Druckrobustheit: Fehler in der Druckapproximation dürfen nicht die Geschwindigkeits- oder Dichtefehler beeinflussen, insbesondere bei rotationsfreien Körperkräften.
• Stabilität: In konvektionsdominierten Regimen (hohe Reynolds-Zahlen) und bei scharfen Dichtegrenzen (z. B. Rayleigh-Taylor-Instabilität) müssen numerische Oszillationen kontrolliert werden, ohne die Genauigkeit zu stark zu beeinträchtigen.
• Effizienz: Die Behandlung hochauflösender, zeitabhängiger Probleme erfordert effiziente Löser und geringen Speicherbedarf.

2. Methodik

Die Autoren entwickeln eine Hybrid High-Order (HHO) Formulierung, gekoppelt mit Explicit Singly Diagonal Implicit Runge-Kutta (ESDIRK) Methoden für die Zeitdiskretisierung.

Räumliche Diskretisierung (HHO)

• Hybride Freiheitsgrade: Die Diskretisierung verwendet Polynome sowohl auf den Zellen (Mesh-Elemente) als auch auf den Rändern (Skeleton) des Gitters.
• Variablen: Geschwindigkeit ($u$), Druck ($p$) und Dichte ($\rho$) werden alle hybrid diskretisiert.
• $H(\text{div})$-Konformität: Durch die Wahl eines hybriden Druckraums und geeigneter lokaler Polynomräume wird eine $H(\text{div})$-konforme Geschwindigkeitsapproximation rekonstruiert. Dies ist entscheidend, da es eine exakte Divergenzfreiheit ($\nabla \cdot u = 0$) punktweise auf jedem Element ermöglicht.
• Konvektionsterme: Die Advektionsterme für Impuls und Masse werden mit schief-symmetrischen Formulierungen behandelt, ergänzt durch Upwind-Stabilisierungsterme, um Robustheit bei konvektionsdominierten Strömungen zu gewährleisten.
• Randbedingungen: Schwache Auferlegung von Dirichlet- und Neumann-Randbedingungen.
• Varianten: Es werden zwei Varianten vorgestellt:
  • HHO-πk: Projektion der advectierten Geschwindigkeit auf den Raum $P_k$ (beweisbar stabil).
  • HHO-πk+1: Projektion auf $P_{k+1}$ (verbesserte Genauigkeit in konvektionsdominierten Fällen, aber nicht schief-symmetrisch).

Zeitdiskretisierung (ESDIRK)

• Da das Problem als Differential-Algebraische Gleichung (DAE) vom Index 2 formuliert ist (aufgrund der Inkompressibilitätsbedingung), werden ESDIRK-Verfahren verwendet, die speziell für DAEs entwickelt wurden.
• Dies ermöglicht eine implizite, hochgenaue Zeitintegration mit hoher Stabilität.

Algebraische Lösung und Effizienz

• Statische Kondensation: Elementinterne Freiheitsgrade werden lokal eliminiert, sodass nur die Randfreiheitsgrade global gekoppelt sind. Dies reduziert die Dimension der globalen Matrix erheblich, insbesondere bei hohen Polynomgraden $k$.
• Löser: Ein $p$-Multigrid-Vorkonditionierer wird in Kombination mit einem FGMRES-Löser eingesetzt, um die Konvergenz zu beschleunigen.

3. Wichtige Beiträge

1. Exakte Volumenerhaltung: Die Methode erzwingt die Volumenerhaltung zellweise bis auf Maschinengenauigkeit. Dies garantiert, dass die Dichte rein advectiert wird, was für die Simulation von nichtmischbaren Fluiden essenziell ist.
2. Druckrobustheit: Die Formulierung ist so konstruiert, dass irrotationale Körperkräfte nur den Druck beeinflussen, nicht jedoch die Geschwindigkeit oder Dichte. Dies wurde theoretisch hergeleitet und numerisch bestätigt.
3. Stabilitätsanalyse: Es wird eine Stabilitätsanalyse für den diskreten Fall durchgeführt, die zeigt, dass die Diskretisierung dissipativ ist und die Energieerhaltungseigenschaften erfüllt (unter homogenen Randbedingungen).
4. Effiziente Hochordnungs-Simulation: Durch die Kombination von HHO, statischer Kondensation und $p$-Multigrid wird ein hoher Rechenaufwand bei hohen Polynomgraden vermieden.

4. Numerische Ergebnisse

Die Methode wurde an drei Testfällen validiert:

• Kovasznay-Strömung (Stationär):

  • Diente zur Validierung der räumlichen Konvergenzraten.
  • Ergebnisse zeigen die erwarteten Konvergenzraten ($k+2$ für Geschwindigkeit, $k+1$ für Druck in $L^2$-Norm).
  • Die Dichte bleibt konstant (wie erwartet), und die Divergenzfreiheit wird bis auf Maschinengenauigkeit eingehalten.
  • HHO-πk+1 zeigt bei $k=0$ eine bessere Genauigkeit als HHO-πk.

• Guermond-Quartapelle-Lösung (Instationär):

  • Validierung der zeitlichen und räumlichen Konvergenzraten mit einer analytischen manufactured solution.
  • Bestätigung der Druckrobustheit: Geschwindigkeitsfehler sind unabhängig von der Druckgenauigkeit.
  • Zeitliche Konvergenz wurde mit 3-stufigen (2. Ordnung) und 6-stufigen (4. Ordnung) ESDIRK-Schemata getestet. Die Ergebnisse stimmen mit der theoretischen Ordnung überein.

• Rayleigh-Taylor-Instabilität (Realistisches Szenario):

  • Simulation bei niedrigen ($At=0.5$) und hohen ($At=0.75$) Atwood-Zahlen sowie verschiedenen Reynolds-Zahlen ($Re=1000, 5000$).
  • Niedrige Atwood-Zahl: Hohe Genauigkeit ($k=6$) auf groben Gittern lieferte Ergebnisse, die mit feinen Gittern ($k=1$) vergleichbar waren, was den Rechenaufwand drastisch reduzierte. Bei sehr hohen Reynolds-Zahlen ($Re=5000$) traten bei hohen Polynomgraden Oszillationen an scharfen Grenzschichten auf, was auf den Bedarf an Oszillationskontrollstrategien hinweist.
  • Hohe Atwood-Zahl: Bei $k=0$ (lineare Elemente) war die Methode schrankenbewahrend (bounds-preserving), d.h. die Dichte blieb im physikalisch sinnvollen Bereich ($1 \le \rho \le 7$), was Stabilität bei hohen Dichteverhältnissen garantiert.

5. Bedeutung und Fazit

Das Paper stellt einen bedeutenden Fortschritt in der numerischen Strömungsmechanik dar, insbesondere für die Simulation von Mehrphasenströmungen mit variabler Dichte.

• Physikalische Treue: Die Kombination aus exakter Inkompressibilität und reiner Dichte-Advektion verhindert physikalisch inkorrekte Artefakte (wie künstliche Diffusion der Dichte).
• Skalierbarkeit: Die Nutzung von statischer Kondensation und $p$-Multigrid macht hochordentliche Methoden auch für komplexe, dreidimensionale Probleme praktikabel.
• Robustheit: Die Druckrobustheit und die Fähigkeit, Dichtegrenzen bei niedrigen Ordnungen einzuhalten, machen die Methode für Anwendungen mit extremen Dichteverhältnissen (z. B. Wasser-Luft-Strömungen) geeignet.

Die Autoren schließen, dass die entwickelte ESDIRK-HHO-Formulierung ein leistungsfähiges Werkzeug für die Simulation von Mischungen nichtmischbarer Fluide ist und planen zukünftige Erweiterungen auf das Cahn-Hilliard-Navier-Stokes-System (CHNS).