A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

Dieser Artikel stellt eine hybride Volume-of-Fluid-Phasenfeld-Methode mit adaptiver Gitterverfeinerung für direkte numerische Simulationen von löslichen Tensid-beladenen Strömungen vor, die die Kopplung zwischen Volumen- und Grenzflächentransport präzise erfasst, um zu zeigen, wie Marangoni-Spannungen die Aufstiegsdynamik von Blasen in dreidimensionalen Geometrien erheblich verändern.

Das Wesentliche

Das große Ganze: Eine neue Art, „seifiges" Wasser zu simulieren

Stellen Sie sich vor, Sie beobachten eine Blase, die durch ein Glas Wasser aufsteigt. Wenn das Wasser rein ist, schießt die Blase gerade nach oben, wackelt ein wenig, bewegt sich aber schnell. Wenn Sie jedoch Seife (Tenside) zum Wasser hinzufügen, verhält sich die Blase anders. Sie könnte langsamer werden, stärker wackeln oder sogar ihren Weg ändern.

Dies geschieht, weil Seifenmoleküle es lieben, an der Oberfläche der Blase zu haften. Wenn sich die Blase bewegt, werden diese Moleküle herumgeschoben und erzeugen eine ungleichmäßige „Spannung" auf der Haut der Blase. Diese ungleichmäßige Spannung wirkt wie unsichtbare Hände, die die Blase in verschiedene Richtungen drücken und ihre Bewegung verändern.

Das Problem:
Die Simulation davon am Computer ist unglaublich schwierig. Es ist wie der Versuch, eine Seifenblase mit einer Kamera zu filmen, die zwei widersprüchliche Einstellungen hat:

1. Die scharfe Kamera: Muss den Rand der Blase als hauchdünne Linie sehen (um Druck und Form zu berechnen).
2. Die unscharfe Kamera: Muss sehen, wie sich die Seifenmoleküle sanft über diesen Rand ausbreiten (um zu berechnen, wie sich die Seife bewegt).

Die meisten Computermethoden zwingen Sie, sich für eine Kameraeinstellung zu entscheiden, was die Simulation entweder physikalisch ungenau oder für komplexe 3D-Formen rechnerisch unmöglich macht.

Die Lösung:
Die Autoren dieses Papiers entwickelten eine Hybridmethode. Stellen Sie sich dies als eine „Split-Screen"-Simulation vor, die gleichzeitig das Beste aus beiden Welten nutzt:

• Der scharfe Rand (Volume-of-Fluid): Sie verwenden eine Methode, die den Rand der Blase scharf hält und die Flüssigkeitsmenge perfekt erhält (wie eine hochauflösende Umrissskizze).
• Die sanfte Seife (Phase-Field): Sie verwenden eine zweite, „verschwommene" Schicht, die als glatte Autobahn für die Seifenmoleküle dient. Dies ermöglicht es der Seife, sich natürlich zwischen Wasser und Blasenoberfläche zu bewegen, ohne stecken zu bleiben oder verloren zu gehen.

Wie es funktioniert: Die „Verkehrsleiter"-Analogie

Um dies zu ermöglichen, schufen die Autoren ein digitales Verkehrssystem für die Seifenmoleküle:

1. Die Autobahn (Die Grenzfläche): Die Oberfläche der Blase ist eine belebte Autobahn. Die Seifenmoleküle sind Autos.
2. Die Auffahrten und Abfahrten (Adsorption/Desorption):
   • Adsorption: Seifenmoleküle aus dem Wasser (dem Volumen) wollen auf die Autobahn (die Blasenoberfläche) springen.
   • Desorption: Seifenmoleküle werden müde und springen von der Autobahn zurück ins Wasser.
   • Die neue Methode berechnet genau, wie viele Autos zu jedem einzelnen Moment auf- oder abfahren, sodass keine Autos verschwinden oder aus dem Nichts auftauchen.
3. Der Stau (Marangoni-Spannung): Wenn zu viele Seifenautos an einer Stelle auf der Blase zusammenlaufen, wird dieser Punkt „klebrig" (hohe Spannung). Die Blasenhaut versucht, sich von diesem klebrigen Punkt wegzuziehen, was eine Kraft erzeugt, die die Blase verlangsamt oder zum Wackeln bringt. Die Simulation erfasst diesen Tauziehen-Verlauf perfekt.

Was sie testeten (Die „Fahrprüfungen")

Bevor sie ihr neues Auto auf die Autobahn ließen, führten sie es durch eine Fahrschule mit drei spezifischen Tests:

1. Der Dehnungstest (Erweiternde Kugel): Sie bliesen eine mit Seife bedeckte Blase auf. Sie prüften, ob sich die Seife gleichmäßig ausbreitete, während die Blase größer wurde. Die Simulation stimmte mit der Mathematik perfekt überein.
2. Der Rotationstest (Drehende Blase): Sie drehten eine Blase mit Seife darauf. Sie prüften, ob sich die Seife korrekt im Kreis bewegte, ohne zu lecken. Auch hier war die Simulation punktgenau.
3. Der Austauschtest (Flache Wand): Sie beobachteten, wie sich Seife vom Wasser zu einer flachen Wand und zurück bewegte. Sie testeten drei Szenarien:
   • Nur Aufspringen: Bleibt die Seife haften? Ja.
   • Nur Abspringen: Verlässt die Seife die Stelle? Ja.
   • Beides: Findet sie ein Gleichgewicht? Ja.

Das Hauptereignis: Die aufsteigende Blase

Schließlich ließen sie ihre neue Methode eine Blase in einem 3D-Wassertank aufsteigen simulieren.

• Die saubere Blase: Sie stieg relativ schnell und gerade auf.
• Die „unlösliche" Seifenblase: Die Seife war an der Oberfläche festgeklebt und konnte nicht weg. Dies erzeugte einen starken „Stau" am hinteren Ende der Blase, der sie erheblich verlangsammte.
• Die „lösliche" Seifenblase (Das echte Ding): Hier glänzt die neue Methode. Die Seife konnte auf die Blase springen und von ihr abspringen.
  • Wenn die Seife leicht absprang (hohe „Desorption"), verhielt sich die Blase fast wie eine saubere.
  • Wenn die Seife leicht aufsprang (hohe „Adsorption"), verhielt sich die Blase wie die Version mit festgeklebter Seife.
  • In der Mitte zeigte die Blase einen komplexen Tanz: Sie verlangsamte sich, änderte ihren Kurs und hinterließ beim Aufsteigen eine „Spur" aus Seife im Wasser.

Warum das wichtig ist (laut dem Papier)

Die Autoren behaupten, diese Methode sei robust, skalierbar und genau.

• Robust: Sie stürzt nicht ab, wenn die Blase seltsame Formen annimmt oder zerbricht.
• Skalierbar: Sie kann auf Supercomputern laufen, um riesige, komplexe 3D-Simulationen effizient zu bewältigen.
• Genau: Sie sagt korrekt voraus, wie schnell Blasen aufsteigen und wie sie wackeln, und stimmt mit der realen Physik überein.

Kurz gesagt: Sie bauten einen neuen digitalen Motor, der endlich simulieren kann, wie Seifenblasen im 3D-Raum verhalten, wobei er den schwierigen Tanz zwischen der Form der Blase und den Seifenmolekülen, die auf und von ihrer Haut springen, bewältigt – alles ohne Genauigkeitsverlust oder einen Absturz des Computers.

Technische Zusammenfassung

Technisches Zusammenfassung: Eine hybride Volume-of-Fluid- und Phase-Field-Methode für interfaciale Strömungen mit löslichen Tensiden

Problemstellung
Die genaue numerische Simulation von interfacialen Strömungen, die lösliche Tenside beinhalten, stellt aufgrund der multiskaligen und multiphysikalischen Natur des Problems eine enorme Herausforderung dar. Im Gegensatz zu unlöslichen Tensiden unterliegen lösliche Spezies einem gekoppelten System, das Bulk-Konvektion–Diffusion, Interfacial-Transport und kinetische Austauschprozesse (Adsorption/Desorption) umfasst. Diese Wechselwirkungen erzeugen räumlich und zeitlich variierende Oberflächenspannungsfelder, die zu Marangoni-Spannungen führen, welche die Strömungstopologie, -stabilität und -bewegung erheblich verändern. Bestehende numerische Methoden weisen deutliche Einschränkungen auf:

• Grenzintegral-/Grenzelementmethoden: Auf Stokes-Strömungsregime beschränkt und haben Schwierigkeiten mit Topologieänderungen (Aufbrechen/Verschmelzen).
• Front-Tracking: Hohe Genauigkeit, aber komplexe Implementierung in 3D, erfordert eine explizite topologische Behandlung und steht vor Lastverteilungsproblemen im parallelen Rechnen.
• Level-Set-Methoden: Handhaben Topologieänderungen auf natürliche Weise, leiden jedoch oft unter Massenerhaltungsproblemen, die eine Neuinitialisierung und nachträgliche Korrekturen erfordern.
• Phase-Field-Methoden: Handhaben Topologieänderungen auf natürliche Weise, verlassen sich jedoch oft auf thermodynamische Gleichgewichtsannahmen, die durch starke Konvektion verletzt werden können, was zu Nicht-Erhaltung der Masse führt.
• Volume-of-Fluid (VoF)-Methoden: Hervorragend für Massenerhaltung und scharfe Grenzflächen, haben jedoch traditionell Schwierigkeiten, die für den Tensidtransport auf einer Grenzfläche mit Null-Dicke erforderlichen Oberflächendifferentialoperatoren zu definieren.

Es besteht ein spezifischer Bedarf an einer geometrischen VoF-basierten Methode, die Adaptive Mesh Refinement (AMR) unterstützt, lösliche Tenside mit dynamischem Bulk-Grenzflächen-Austausch handhabt und innerhalb eines skalierbaren, quelloffenen Rahmens operiert.

Methodik
Die Autoren stellen einen hybriden numerischen Rahmen vor, der im quelloffenen Code Basilisk implementiert ist und eine geometrische Volume-of-Fluid (VoF)-Methode mit einem diffusen Phase-Field-Ansatz kombiniert.

1. Grenzflächendarstellung:

   • VoF ($c$): Ein skalares Feld, das den Volumenanteil einer Flüssigkeit darstellt und für die konservative Grenzflächenverfolgung und die Impulslösung verwendet wird. Es gewährleistet eine exakte Massenerhaltung und unterstützt die geometrische Rekonstruktion.
   • Phase-Field ($\phi$): Eine Hilfsvariable, die sich über eine auf der Allen-Cahn-Gleichung basierende Accurate Conservative Diffuse-Interface (ACDI)-Formulierung entwickelt. Sie dient als glatter Träger für den Tensidtransport und liefert regularisierte Normalvektoren, Krümmungen und Oberflächengradienten, ohne dass explizite Oberflächenableitungen auf einer scharfen Grenzfläche berechnet werden müssen.
   • Kopplung: Das Phase-Field wird periodisch aus dem VoF-Feld unter Verwendung einer Vorzeichen-Abstandsfunktion reinitialisiert, um Konsistenz zu gewährleisten.

2. Steuernde Gleichungen:

   • Flüssigkeitsbewegung: Inkompressible Navier-Stokes-Gleichungen mit räumlich variierender Dichte und Viskosität. Oberflächenspannungskräfte werden als volumetrische Körperkräfte integriert: eine isotrope Kapillarkraft ($\sigma \kappa \mathbf{n} \delta_s$) und eine tangentialer Marangoni-Kraft ($\nabla_s \sigma \delta_m$).
   • Tensidtransport: Zwei skalare Felder werden gelöst: Bulk-Konzentration ($F$) und Grenzflächenkonzentration ($f$, eine volumetrische Darstellung von $\Gamma$).
     • Die Grenzflächentransportgleichung wird regularisiert, um Oberflächenableitungen zu vermeiden, indem die Oberflächekonzentration $\Gamma$ in ein volumetrisches Feld $f = \Gamma \delta_\phi$ umgewandelt wird, das innerhalb der diffusen Grenzfläche lokalisiert ist.
     • Kinetik: Adsorption und Desorption werden über regularisierte Quellen/Senken-Terme ($J$) modelliert, die an der Grenzfläche lokalisiert sind und Bulk- und Grenzflächenkonzentrationen verknüpfen. Das Modell unterstützt allgemeine Zustandsgleichungen (z. B. Langmuir-Isotherme), um die Oberflächenspannung mit der Konzentration in Beziehung zu setzen.

3. Numerische Diskretisierung:

   • Räumlich: Finite-Volumen-Diskretisierung auf dynamisch verfeinerten Quadtree/Octree-Gittern (AMR). Die Auflösung wird in der Nähe der Grenzfläche konzentriert, um steile Gradienten in Konzentration und Strömung aufzulösen.
   • Zeitlich: Operator-Splitting wird verwendet. Advektion wird explizit behandelt (geometrisches ungeteiltes Bell–Collela–Glaz-Schema), während Diffusion, Anti-Diffusion und Quellterme implizit unter Verwendung eines modifizierten Multigrid-Poisson-Lösers behandelt werden.
   • Stabilität: Eine Courant–Friedrichs–Lewy (CFL)-Bedingung wird basierend auf einer effektiven Geschwindigkeit durchgesetzt, die sowohl konvektive als auch anti-diffusive Beiträge umfasst, um Stabilität bei Vorhandensein starker Grenzflächengradienten zu gewährleisten.

Hauptbeiträge

• Hybrider Rahmen: Entwicklung einer Methode, die die massenerhaltenden und scharf-grenzflächen-Vorteile von VoF beibehält, während eine Phase-Field-Variable ausschließlich als glatter Träger für den Tensidtransport genutzt wird, wodurch die geometrischen Komplexitäten von Oberflächenableitungen auf stückweise linearen Grenzflächen vermieden werden.
• Fähigkeit für lösliche Tenside: Erweiterung vorheriger VoF-Phase-Field-Arbeiten (die auf unlösliche Tenside beschränkt waren) auf den vollständig löslichen Fall unter Einbeziehung gekoppelter Bulk-Advektion-Diffusion und kinetischer Austauschterme.
• Skalierbarkeit und Adaptivität: Der Rahmen ist vollständig adaptiv (AMR) und unterstützt 2D-, achsensymmetrische und 3D-Konfigurationen mit hoher paralleler Skalierbarkeit, was die Auflösung von Multiskalen-Gradienten sowohl in Bulk- als auch in Grenzflächenfeldern ermöglicht.
• Quelloffene Implementierung: Die Methode ist in Basilisk implementiert und bietet der Gemeinschaft ein reproduzierbares und zugängliches Werkzeug.

Ergebnisse und Validierung
Die Arbeit validiert die Methode durch eine Reihe von Benchmark-Tests und Anwendungsfällen:

1. Validierung ohne Strömung:

   • Expandierende Kugel (3D): Verifizierte die Wiederherstellung analytischer Lösungen für die Verdünnung unlöslicher Tenside während isotroper Expansion.
   • Rotierende Grenzfläche (2D): Bestätigte eine Konvergenz zweiter Ordnung für die Advektion-Diffusion unlöslicher Tenside auf einer bewegten Grenzfläche.
   • Adsorption/Desorption an flacher Grenzfläche: Validierte den gekoppelten Bulk-Grenzflächen-Löser gegen analytische Lösungen für Adsorption mit konstantem Fluss, reine Desorption und vollständige Langmuir-Kinetik. Die Methode erfasst korrekt den Übergang zum Gleichgewicht und die Bildung von Verarmungs-/Anreicherungsschichten im Bulk.

2. Achsensymmetrisch aufsteigende Blase:

   • Unlöslicher Grenzfall: Zeigte die Verringerung der terminalen Steiggeschwindigkeit aufgrund von Marangoni-Spannungen, wobei die Ergebnisse mit abnehmender Tensidkonzentration gegen den Grenzfall der sauberen Grenzfläche konvergierten.
   • Lösliche Dynamik: Systematische Variation der Biot-Zahl (Desorption) und der Damköhler-Zahl (Adsorption).
     • Hohe Biot-Zahl (schnelle Desorption): Das System nähert sich dem Grenzfall der sauberen Grenzfläche an, da Tenside schnell von der Grenzfläche entfernt werden.
     • Hohe Damköhler-Zahl (schnelle Adsorption): Das System nähert sich dem Grenzfall der oberflächeninkompressiblen (unlöslichen) Grenzfläche an, da die Grenzfläche schnell gesättigt wird.
   • Strömungsstruktur: Visualisierungen zeigten die Bildung von „Stagnationskappen" aus Tensid am Blasenrücken, die innere Rezirkulation unterdrücken und die Wake-Struktur modifizieren.

3. Dreidimensionale aufsteigende Blase:

   • Simulationen einer in 3D aufsteigenden Blase mit löslichen Tensiden erfassten komplexe Phänomene, einschließlich Formverformung, seitlicher Bahnabweichungen (helikoidale Bewegung) und der Bildung von hinteren Schleppe aufgrund von Desorption.
   • Die Methode erholte erfolgreich sowohl die sauberen als auch die unlöslichen asymptotischen Grenzfälle in 3D und demonstrierte Robustheit im Umgang mit instationären Strömungsmustern und verformenden Geometrien.

Bedeutung und Behauptungen
Die Arbeit behauptet, ein robuster und vielseitiger numerischer Rahmen zu bieten, der die Lücke zwischen klassischen Modellen für unlösliche Tenside und realistischen Formulierungen für lösliche Tenside schließt. Ihre Bedeutung liegt in:

• Vereinheitlichte Behandlung: Sie handhabt nichtlineare Kopplung, sich entwickelnde Grenzflächen und vollständig dreidimensionale Konfigurationen innerhalb einer einzigen quelloffenen Plattform.
• Physikalische Treue: Sie erfasst das komplexe Zusammenspiel zwischen Hydrodynamik, Bulk-/Grenzflächen-Transport und Marangoni-Spannungen, ohne sich auf Turbulenzschließungsmodelle oder empirische Näherungen zu verlassen.
• Praktische Anwendbarkeit: Die Methode eignet sich für Freiflächenströmungen, die Grenzflächenwiederverbindung beinhalten (z. B. Strahlaufbrechen, brechende Wellen) in Gegenwart von industriell und biologisch relevanten löslichen Tensiden.

Die Autoren vermerken bescheiden, dass der aktuelle Rahmen aufgrund inhärenter numerischer Diffusion in der Diskretisierung auf moderate Péclet-Zahlen beschränkt ist. Sie identifizieren die genaue Auflösung scharfer Grenzflächengradienten im Bereich hoher Péclet-Zahlen (konvektionsdominiert) als Herausforderung für zukünftige Arbeiten mit dem Ziel, Formulierungen zu entwickeln, die numerische Diffusion minimieren.