A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

Die Autoren stellen eine scharfe und konservative 3D-VOF-Methode vor, die auf einem eingebetteten Grenzflächen-Framework basiert, um Kontaktliniendynamiken mit Hysterese auf komplexen Geometrien durch eine modifizierte Advektionsschemata für Mischzellen, eine Umverteilungsstrategie zur Beseitigung der CFL-Einschränkung und ein neuartiges Kontaktwinkel-Imponierungsverfahren präzise zu simulieren.

Das Wesentliche

🌊 Der digitale Tropfen: Wie man Wassertropfen auf komplexen Oberflächen simuliert

Stellen Sie sich vor, Sie lassen einen Wassertropfen auf eine glatte Tischplatte fallen. Das ist einfach: Er breitet sich rund aus. Aber was passiert, wenn Sie den Tropfen auf eine unregelmäßige, gewölbte oder sogar durchlöcherte Oberfläche fallen lassen? Vielleicht auf einen Kegel, eine gewellte Metallplatte oder ein feines Netz?

Genau das ist die Herausforderung, der sich die Autoren dieses Papers stellen. Sie haben einen neuen, hochpräzisen digitalen Werkzeugkasten entwickelt, um zu berechnen, wie sich Flüssigkeiten auf solchen komplexen 3D-Oberflächen verhalten.

Hier ist die Geschichte ihrer Lösung, aufgeteilt in drei einfache Teile:

1. Das Problem: Der "verlorene" Tropfen in kleinen Ecken

In Computersimulationen wird die Welt normalerweise in kleine Würfel (Gitterzellen) zerlegt. Wenn eine flache Wand durch diese Würfel schneidet, entstehen seltsame, winzige Ecken – sogenannte "geschnittene Zellen" (cut cells).

• Die alte Methode: Stellen Sie sich vor, Sie versuchen, Wasser durch ein Sieb zu gießen, aber das Sieb hat Löcher, die so klein sind, dass der Computer den Durchfluss nicht genau messen kann. Die alten Methoden haben hier oft "vergessen", wie viel Wasser eigentlich wo ist. Der Tropfen schien plötzlich zu verschwinden oder sich an falschen Stellen zu sammeln.
• Die neue Lösung (Der "Umlade-Truck"): Die Autoren haben eine neue Regel erfunden. Wenn eine Zelle zu klein ist, um das Wasser genau zu messen, nehmen sie das Wasser nicht weg. Stattdessen sagen sie: "Okay, dieser kleine Bereich ist voll, aber wir haben zu viel Wasser reingepackt. Wir geben das Übermaß sofort an den nächsten Nachbarn weiter."
  • Die Analogie: Stellen Sie sich einen riesigen Umzug vor. Wenn ein kleines Zimmer (die geschnittene Zelle) nicht genug Platz für alle Möbel hat, schieben die Umzugshelfer die überschüssigen Kisten sofort in das nächste Zimmer, damit nichts verloren geht. So bleibt die Gesamtmenge an Wasser (die Masse) immer perfekt erhalten, egal wie winzig die Ecken sind.

2. Der Kontaktwinkel: Der "mürrische" Tropfen

Wenn ein Tropfen auf eine Oberfläche trifft, bildet er einen bestimmten Winkel. Das nennt man den Kontaktwinkel.

• Auf einer hydrophoben (wasserabweisenden) Oberfläche steht der Tropfen wie eine Perle (hoher Winkel).
• Auf einer hydrophilen (wasserliebenden) Oberfläche flacht er ab (niedriger Winkel).

Das Problem bei komplexen Oberflächen (z. B. einer schiefen oder gekrümmten Wand) ist, dass der Computer oft den Winkel falsch berechnet. Er versucht, den Tropfen wie auf einer flachen Ebene zu berechnen, was zu Fehlern führt.

• Die alte Methode: Sie versuchte, den Tropfenrand mit einer einfachen, geraden Linie zu verlängern. Das ist wie wenn man versucht, die Konturen eines Berges mit einem Lineal zu zeichnen – es passt einfach nicht.
• Die neue Lösung (Der "Paraboloid-Verfeinerer"): Die Autoren haben eine Methode entwickelt, die den Tropfenrand nicht mit einer geraden Linie, sondern mit einer gekrümmten Parabel (wie eine sanfte Hängebrücke) modelliert.
  • Die Analogie: Stellen Sie sich vor, Sie wollen einen Zaun um einen krummen Garten bauen. Die alte Methode würde gerade Bretter verwenden, die Lücken lassen. Die neue Methode passt die Bretter so an, dass sie der Kurve des Gartens perfekt folgen. Zudem berücksichtigen sie, dass ein Tropfen manchmal "klebt" (Hysterese) – er bewegt sich nicht sofort, wenn die Oberfläche rau ist, sondern wartet, bis der Druck groß genug ist. Ihr Algorithmus simuliert dieses "Kleben" und "Loslassen" sehr realistisch.

3. Der Test: Von der Kugel bis zum Netz

Um zu beweisen, dass ihr Werkzeug funktioniert, haben sie es an extremen Tests geprüft:

• Der rotierende Tropfen: Ein Tropfen, der sich um eine Kugel dreht. Die neue Methode hat gezeigt, dass kein einziger Wassertropfen verloren ging, während die alten Methoden hier Fehler machten.
• Der Tropfen auf dem Kegel: Ein Tropfen, der sich allein durch die Krümmung des Kegels nach oben bewegt (wie eine Schnecke). Die Simulation zeigte genau das gleiche Verhalten wie in der echten Physik.
• Der Tropfen auf dem Netz: Ein Tropfen, der auf ein feines Drahtgeflecht fällt. Hier muss der Tropfen entscheiden, ob er durch die Löcher fällt oder oben bleibt. Die Simulation hat das Verhalten perfekt nachgebildet.

🏆 Das Fazit in einem Satz

Die Autoren haben einen neuen, extrem genauen "digitalen Wasser-Simulator" gebaut, der nicht nur auf flachen Flächen funktioniert, sondern auch auf den wildesten, krummsten und kleinsten 3D-Oberflächen, ohne dabei auch nur einen einzigen Wassertropfen zu verlieren oder die Physik zu verfälschen.

Warum ist das wichtig?
Dies hilft Ingenieuren, bessere Tintenstrahldrucker zu entwickeln, effizientere Brennstoffzellen zu bauen oder zu verstehen, wie Wasser auf mikroskopisch kleinen Strukturen in der Natur (wie auf Lotusblättern) läuft. Es ist ein großer Schritt vorwärts für die digitale Welt der Flüssigkeiten.

Technische Zusammenfassung

Hier ist eine detaillierte technische Zusammenfassung des Artikels „A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries" von Huang et al. (2026) auf Deutsch.

1. Problemstellung

Die Dynamik von Kontaktlinien (der Schnittstelle zwischen Flüssigkeit, Gas und Festkörper) ist für zahlreiche ingenieurwissenschaftliche Prozesse entscheidend, wie z. B. beim Tintenstrahldruck oder bei der Mikrofluidik. Während die Simulation auf einfachen, ebenen Substraten gut etabliert ist, stellt die Modellierung auf komplexen, gekrümmten oder beliebig geformten Festkörpergrenzflächen erhebliche numerische Herausforderungen dar.

Die Hauptprobleme, die in dieser Arbeit adressiert werden, sind:

• Massenerhaltung in Mischzellen: Bei der Verwendung von Cartesianischen Gittern mit eingebetteten Grenzen (Embedded Boundary Method, EBM) entstehen „Mischzellen" (Mixed Cells), die von der Festkörpergrenze durchschnitten werden. Herkömmliche VOF-Verfahren (Volume-of-Fluid) führen hier oft zu lokalen Massenerhaltungsfehlern und physikalisch nicht sinnvollen Ansammlungen von Fluid im Festkörper.
• Zeitschrittbeschränkungen: Kleine abgeschnittene Zellen (Small Cut Cells) erzwingen extrem kleine Zeitschritte aufgrund der CFL-Bedingung (Courant-Friedrichs-Lewy), was die Recheneffizienz drastisch reduziert.
• Kontaktwinkel-Bedingungen: Die genaue Auferlegung von Kontaktwinkelbedingungen auf gekrümmten Oberflächen ist schwierig. Herkömmliche Methoden (z. B. lineare Extrapolation oder einfache Ghost-Cell-Methoden) führen bei extremen Kontaktwinkeln oder komplexen Geometrien zu Fehlern in der Krümmungsschätzung und damit zu ungenauen Oberflächenspannungskräften.
• Hysterese: Die Modellierung des Kontaktwinkel-Hysterese-Effekts (Unterschied zwischen voranschreitendem und zurückweichendem Winkel) in 3D auf komplexen Oberflächen ist bisher kaum robust gelöst.

2. Methodik

Die Autoren entwickeln ein neues numerisches Framework innerhalb des Open-Source-Lösers Basilisk, das auf einem gekoppelten geometrischen VOF-Ansatz und der Embedded Boundary Method (EBM) basiert.

A. Geometrische VOF-Advektion in Mischzellen

• Rekonstruktion: Ein modifizierter Flood-Fill-Algorithmus (3D-Erweiterung des 2D-Ansatzes) wird verwendet, um die Flüssigkeits-Gas-Grenzfläche (Γl) innerhalb der unregelmäßigen Polyeder in Mischzellen zu rekonstruieren. Dies ermöglicht eine präzise Berechnung des Flüssigkeitsvolumens trotz der Schnittmenge mit dem Festkörper (Γs).
• Konservative Flux-Korrektur: Um lokale Massenerhaltung sicherzustellen, wird eine Korrektur der Advektionsgeschwindigkeit eingeführt. Der effektive Advektionsbreiten-Parameter ($u_n^*$) wird so angepasst, dass das durch die Festkörpergeometrie blockierte Volumen kompensiert wird. Dies verhindert systematische Über- oder Unterfüllung von Zellen.
• Redistributionsstrategie (Flux Redistribution): Um die strenge CFL-Bedingung durch kleine Cut-Cells zu umgehen, wird eine Redistributionsstrategie entwickelt.
  1. Der Zeitschritt wird unabhängig von der Fluidfraktion ($c_s$) gewählt.
  2. Falls eine Zelle nach der Advektion überfüllt oder unterfüllt ist, wird das überschüssige bzw. fehlende Volumen konservativ auf benachbarte Zellen entlang der Strömungsrichtung umverteilt.
  3. Dies eliminiert die CFL-Einschränkung für Cut-Cells vollständig, während die globale und lokale Volumenerhaltung gewahrt bleibt.

B. Kontaktwinkel-Bedingung mittels Höhenfunktion (Height Function - HF)

• Paraboloid-Anpassung (Paraboloid Fitting): Anstatt einer linearen Extrapolation wird eine parabolische Anpassung der Grenzfläche verwendet.
  1. Pre-Fitting: Zuerst wird eine vorläufige Parabeloberfläche basierend auf den Höhenfunktionen (HF) in den Fluidzellen angepasst, um eine stabile geometrische Approximation zu erhalten.
  2. Constraint-Fitting: Anschließend wird eine finale Parabeloberfläche angepasst, die explizit die vorgeschriebenen Kontaktwinkelbedingungen am Kontaktpunkt erfüllt. Dies geschieht durch Lösung eines nichtlinearen Gleichungssystems (Newton-Verfahren).
• Vertikale Höhenfunktion (VHF): Für extremste Kontaktwinkel ($\theta < 45^\circ$ oder $\theta > 135^\circ$) oder fast horizontale Oberflächen, wo die Standard-HF-Methode versagt, wird eine vertikale Höhenfunktion eingeführt. Diese erlaubt die Definition von HF-Werten auch in degenerierten Geometrien, indem die Schnittlinie der Grenzfläche mit dem Festkörper direkt genutzt wird.
• Hysterese-Modellierung: Ein iterativer Bisektionsalgorithmus wird innerhalb eines Hysterese-Fensters $[\theta_{rec}, \theta_{adv}]$ angewendet. Der Kontaktwinkel wird so angepasst, dass die Kontaktlinie an ihrer vorherigen Position „festgeklebt" (pinned) bleibt, solange der lokale Winkel nicht die Schwellenwerte erreicht.

3. Wichtige Beiträge

1. Erster vollständig geometrischer und konservativer 3D-VOF-Algorithmus für Kontaktlinien auf beliebig komplexen 3D-Oberflächen, der strikte lokale Massenerhaltung garantiert.
2. Entwicklung einer Redistributionsstrategie, die die CFL-Beschränkung durch kleine Cut-Cells in 3D-Simulationen vollständig aufhebt, ohne die Genauigkeit zu beeinträchtigen.
3. Neuartige 3D-HF-Methode mit Paraboloid-Anpassung, die robust Kontaktwinkel auf gekrümmten und eingebetteten Oberflächen auferlegt und signifikant genauer ist als lineare Extrapolationsmethoden.
4. Integration von Kontaktwinkel-Hysterese in das 3D-Framework, was realistischere Benetzungsdynamiken ermöglicht.

4. Ergebnisse und Validierung

Die Methode wurde durch eine Reihe anspruchsvoller Benchmark-Tests validiert:

• Volumenerhaltung: Tests mit rotierenden Flüssigkeitskugelschalen (Concentric Sphere & Conical Cutout) zeigten, dass die korrigierte Advektion lokale Fehler auf Maschinengenauigkeit ($\sim 10^{-15}$) reduziert, während herkömmliche Methoden signifikante Fehler aufweisen. Die Redistributionsmethode ermöglichte lange Simulationszeiten ohne die sonst erforderlichen extrem kleinen Zeitschritte.
• Tropfenausbreitung auf ebenen und gekrümmten Flächen: Die Simulation von Tropfen auf ebenen Platten (horizontal und geneigt) und auf sphärischen Oberflächen bestätigte die hohe Genauigkeit der Kontaktwinkelauferlegung. Im Gegensatz zu linearen Methoden zeigte die Paraboloid-Methode eine perfekte Rotationssymmetrie und war unabhängig von der Position der eingebetteten Grenze im Gitter.
• Schubgetriebene Tropfenbewegung mit Hysterese: Simulationen von Tropfen unter Scherströmung zeigten das korrekte Verhalten von Anheften (Pinning) und Lösen (Depinning) der Kontaktlinie, übereinstimmend mit experimentellen Daten und früheren Studien.
• Komplexe Geometrien: Die Methode bewährte sich bei der Simulation von Tropfen auf wellenförmigen Substraten, beim Aufprall auf Platten mit Öffnungen (Orifice) und beim spontanen Klettern auf konischen sowie mesh-förmigen Oberflächen. In allen Fällen wurden die Kontaktlinien-Dynamiken stabil und genau aufgelöst.

5. Bedeutung

Diese Arbeit stellt einen bedeutenden Fortschritt in der numerischen Strömungsmechanik dar. Sie schließt die Lücke zwischen der hohen Genauigkeit scharfer Grenzflächenmethoden (Sharp Interface) und der Fähigkeit, komplexe Geometrien effizient zu behandeln.

• Praktische Relevanz: Das Verfahren ist essenziell für die genaue Simulation von Benetzungsphänomenen in realen Anwendungen, wo Oberflächen selten perfekt eben sind (z. B. Mikrofluidik-Chips, Beschichtungsprozesse, biologische Systeme).
• Effizienz: Durch die Beseitigung der Cut-Cell-Zeitschrittbeschränkung werden 3D-Simulationen mit komplexen Festkörpern erstmals praktikabel und recheneffizient.
• Robustheit: Die Kombination aus konservativer Advektion und der robusten Kontaktwinkel-Modellierung macht das Werkzeug zu einer zuverlässigen Basis für zukünftige Forschungen im Bereich der Mehrphasenströmungen.

Zusammenfassend bietet das vorgestellte Framework ein leistungsfähiges, genaues und konservierendes Werkzeug zur Untersuchung von Benetzungsphänomenen auf beliebig komplexen 3D-Oberflächen.