Compressible and immiscible fluids with arbitrary density ratio

Diese Arbeit präsentiert eine neue thermodynamikbasierte Theorie zur Modellierung kompressibler und unmischbarer Fluide mit beliebig hohen Dichteverhältnissen, welche die Einschränkungen der traditionellen Navier-Stokes- und Euler-Gleichungen adressiert, indem sie Dichteentwicklungs-Gleichungen herleitet, die die wahre Impulsentwicklung angemessen berücksichtigen.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, einen Tanz zwischen zwei sehr unterschiedlichen Partnern zu beschreiben: einem schweren, langsam beweglichen Elefanten (Wasser) und einer leichten, schnell beweglichen Feder (Luft). In der Welt der Fluidphysik sind diese beiden „unmischbar“, das heißt, sie vermischen sich nicht wie Öl und Wasser, aber sie haben eine unscharfe Grenze, an der sie aufeinandertreffen.

Lange Zeit haben Wissenschaftler darum gerungen, die „Tanzregeln“ (die Gleichungen) für Systeme zu formulieren, bei denen der Dichteunterschied riesig ist, wie etwa wenn der Elefant 1.000 Mal schwerer ist als die Feder. Die alten Regeln hatten einen entscheidenden Fehler: Sie behandelten das Gewicht des Elefanten so, als wäre es eine konstante, unveränderliche Zahl, selbst direkt an der unscharfen Grenze, an der der Elefant zur Feder wird. Das ist so, als würde man versuchen, eine Person, die sich in ein Gespenst verwandelt, zu beschreiben, indem man sagt: „Sie bleibt die ganze Zeit über ein fester Mensch“, was keinen Sinn ergibt.

Hier ist eine einfache Aufschlüsselung dessen, was dieses Paper vorschlägt, um dieses Problem zu lösen.

1. Das Problem mit den alten Regeln

Die traditionelle Art, die Berechnung der Fluidbewegung (unter Verwendung der Navier-Stokes- und Euler-Gleichungen), stützt sich auf eine Abkürzung, die als Boussinesq-Approximation bezeichnet wird.

• Die Analogie: Stellen Sie sich vor, Sie schieben einen Einkaufswagen. Wenn der Wagen voll mit Ziegeln ist, ist er schwer. Wenn er leer ist, ist er leicht. Die alten Regeln gingen davon aus, dass das Gewicht des Wagens niemals seine Form verändert, während Sie ihn schieben, selbst wenn der Wagen teilweise mit Ziegeln gefüllt und teilweise mit Luft leer ist. Es wird einfach angenommen, dass das Gewicht ein festes Durchschnittsgewicht ist.
• Der Fehler: In der Realität ändert sich das Gewicht, während sich der Wagen bewegt und die Ziegel sich verschieben (Diffusion). Die alten Regeln ignorierten die Tatsache, dass sich der „Impuls“ (Masse $\times$ Geschwindigkeit) ändert, weil sich die Masse selbst an der Grenze ändert. Sie nahmen auch an, dass Luft und Wasser einen perfekt vorhersehbaren Raum einnehmen, und ignorierten dabei, dass der Raum, den sie einnehmen, wenn sie sich an der Grenze vermischen, aufgrund des Drucks tatsächlich wackeln und sich verändern kann.

2. Der neue Ansatz: Energieminimierung

Anstatt zu raten, wie sich die Dichte ändert, geht der Autor von einem grundlegenden Prinzip aus: Die Natur versucht immer, so wenig Energie wie möglich zu verbrauchen.

• Die Analogie: Denken Sie an einen Ball, der einen Hügel hinunterrollt. Es ist ihm egal, welche spezifischen Schritte er macht; er will nur zum tiefsten Punkt (niedrigste Energie) gelangen. Der Autor nutzt dieses Konzept des „Energiehügels“, um neue Regeln dafür abzuleiten, wie Wasser und Luft interagieren.
• Die zentrale Neuerung: Der Autor führt das Konzept des „Überschussvolumens“ (Excess Volume) ein.
  • Stellen Sie sich vor, Sie haben einen Eimer Wasser und einen Eimer Luft. Wenn Sie sie zusammengießen, würden Sie vielleicht erwarten, dass das Gesamtvolumen genau die Summe der beiden Eimer ist. Aber auf mikroskopischer Ebene, wenn sie aufeinandertreffen, könnten sich die Moleküle etwas enger zusammenpacken oder lockerer verteilen, wodurch ein kleines bisschen „zusätzlicher“ oder „fehlender“ Raum entsteht.
  • Die alten Regeln gingen davon aus, dass dieser zusätzliche Raum überall Null sei. Dieses Paper sagt: „Nein, dieser zusätzliche Raum existiert, er ändert sich von Ort zu Ort und beeinflusst die Dichte.“

3. Die neuen „Tanzregeln“ (Die Ergebnisse)

Indem der Autor den sich ändernden „zusätzlichen Raum“ berücksichtigt und die Energieminimierung nutzt, leitet er einen neuen Satz von Gleichungen ab, die im Wesentlichen drei Dinge bewirken:

A. Eine neue Art, Schall zu hören (Schallgeschwindigkeit)
Das Paper zeigt, dass die Schallgeschwindigkeit nicht einfach eine Zufallszahl ist, sondern direkt daraus resultiert, wie sich die Energie des Fluids ändert, wenn es zusammengedrückt wird.

• Die Metapher: Denken Sie an den Schall als eine Welle in einer Menschenmenge. Die Geschwindigkeit dieser Welle hängt davon ab, wie dicht die Menschen (Moleküle) gepackt sind und wie viel Energie sie haben. Die neue Formel berechnet diese Geschwindigkeit ganz natürlich, ohne dass man ihr vorher sagen muss, wie hoch sie ist. Sie legt sogar nahe, dass in einem Gas die Schallgeschwindigkeit in etwa der durchschnittlichen Geschwindigkeit entspricht, mit der die Gasmoleküle umherhüpfen.

B. Eine neue Regel für Druck und Geschwindigkeit (Bernoulli-Gesetz)
Sie haben wahrscheinlich schon vom Bernoulli-Prinzip gehört: „Wenn ein Fluid schneller fließt, sinkt sein Druck.“

• Die Wendung: Die alte Regel funktioniert großartig für Wasser, das in einem Rohr fließt, aber sie bricht zusammen, wenn man einen riesigen Sprung in der Dichte hat (wie Wasser, das auf Luft trifft). Der Autor erstellt ein verallgemeinertes Bernoulli-Gesetz.
• Die Metapher: Stellen Sie sich einen Fluss vor, der in einen Wasserfall fließt. Die alte Regel besagt, dass die Energie gleich bleibt. Die neue Regel sagt: „Warte, während das Wasser zu Nebel wird, geht etwas Energie verloren oder transformiert sich, weil die ‚Dichte/Beschaffenheit‘ des Wassers sich ändert.“ Die neue Gleichung berücksichtigt diesen Energieübergang und ist daher auch dann genau, wenn sich die Natur des Fluids von schwer zu leicht ändert.

C. Der „Buckel“ in der Dichte
Dies ist vielleicht das visuellste Ergebnis.

• Die alte Sichtweise: Wenn man sich die Grenze zwischen Wasser und Luft ansieht, sagten die alten Modelle, dass die Dichte einfach glatt von „schwerem Wasser“ zu „leichter Luft“ gleiten würde, wie eine abfallende Rampe.
• Die neue Sichtweise: Die Mathematik des Autors sagt einen Buckel voraus. Wenn man die Grenze überquert, steigt die Dichte tatsächlich kurzzeitig leicht an, bevor sie auf das Luftniveau abfällt.
• Die Metapsaher: Stellen Sie sich eine Menschenmenge (Wasser) vor, die versucht, durch eine Tür in einen leeren Flur (Luft) zu schlüpfen. Während sie durch die Tür zwängen, könnten sie für einen kurzen Moment dichter zusammenrücken, bevor sie sich verteilen. Die neue Theorie sagt diesen „Packungs-Buckel“ voraus, was mit fortgeschrittenen Computersimulationen (der sogenannten Dichtefunktionaltheorie) bereits bestätigt wurde, die die alten einfachen Modelle jedoch übersehen hatten.

Zusammenfassung

Dieses Paper schlägt einen neuen Weg vor, die physikalischen Gesetze für Fluide zu formulieren, die sich in ihrem Gewicht stark unterscheiden (wie Wasser und Luft).

1. Es hört auf vorzutäuschen, dass das Gewicht an der Grenze konstant ist.
2. Es gibt zu, dass sich der Raum, den Moleküle einnehmen, ändert (Überschussvolumen).
3. Es nutzt das Prinzip der „geringsten Energie“, um neue Regeln abzuleiten, die erklären, wie Schall reist, wie sich der Druck verändert und warum die Dichte an der Wasser-Luft-Grenze tatsächlich einen kleinen „Höcker“ aufweist.

Der Autor behauptet, dass dieser neue Rahmen für jedes Gemisch von Fluiden funktioniert, unabhängig davon, wie unterschiedlich ihre Gewichte sind, und öffnet damit die Tür für genauere Computersimulationen von Phänomenen wie fallendem Regen, brechenden Wellen oder aufsteigenden Blasen.

Technische Zusammenfassung

Technische Zusammenfassung: Kompressible und unmischbare Fluide mit beliebigen Dichteverhältnissen: Diffuse Grenzflächen-Theorie

Problemstellung
Aktuelle numerische Strömungsmechanik-Modelle (CFD) haben Schwierigkeiten, unmischbare Fluidsysteme mit hohen Dichteverhältnissen, wie etwa Wasser-Luft (Verhältnis $\approx 1000$), genau zu simulieren. Bestehende Ansätze stützen sich primär auf zwei Methoden, die beide signifikante Einschränkungen aufweisen:

1. Die Boussinesq-Approximation: Diese nimmt eine konstante Dichte ($\rho$) innerhalb der Grenzfläche an und nutzt die Gleichung $\rho(du/dt) = -\nabla p + \dots$. Diese Formulierung versäumt es, die „wahre“ Impulsentwicklung $d(\rho u)/dt$ zu berücksichtigen, wodurch Dichteänderungen und der daraus resultierende Impulstransfer vernachlässigt werden. Zudem können dadurch Gravitations- und Auftriebskräfte nicht auf natürliche Weise ohne ad-hoc hinzugefügte Terme erzeugt werden, da diese aus Dichteunterschieden resultieren.
2. Die Interpolationsmethode: Diese definiert die Dichte als eine lineare Funktion der Volumenkonzentration ($\rho = \rho_w \phi + \rho_a(1-\phi)$). Dieser Ansatz erzeugt einen fundamentalen Konflikt zwischen Inkompressibilität ($d\rho/dt = 0$) und Diffusion ($d\phi/dt \neq 0$). Zudem wird ein einheitliches Überschussvolumen (oft Null) über das gesamte Gebiet vorausgesetzt, was der physikalischen Realität widerspricht, in der das lokale Überschussvolumen vom lokalen Druck abhängt. Folglich können die Standard-Kontinuitätsgleichungen und die lineare Interpolation die nicht-monotonen Dichteprofile, die an der Wasser-Luft-Grenzfläche beobachtet werden, nicht erfassen.

Methodik
Die Autoren schlagen ein neuartiges theoretisches Rahmenwerk vor, das auf dem thermodynamischen Prinzip der Energieminimierung basiert und sich von der traditionellen Ableitung der Navier-Stokes- und Euler-Gleichungen unterscheidet. Die Methodik umfasst:

• Reformulierung von Volumen und Dichte: Anstatt ein Überschussvolumen von Null anzunehmen, führt das Modell ein räumlich variierendes Überschussvolumen ($v_e$) innerhalb der Grenzfläche ein. Dieses Volumen wird zwischen den Fluidbestandteilen umverteilt, was zu partiellen Dichten ($\rho'_w, \rho'_a$) führt, die nicht mehr konstant sind, sondern von lokalem Druck und Mischung abhängen. Dies erfordert zwei gekoppelte Evolutionsgleichungen: eine für die Volumenkonzentration ($\phi$) und eine für die Dichte ($\rho$).
• Dissipations-Erhaltungssatz: Die Autoren etablieren ein Theorem, das den ersten Hauptsatz der Thermodynamik (Erhaltung) mit dem zweiten Hauptsatz (Dissipation) verknüpft. Durch die Definition der Gesamtenergie des Systems ($E$) als Summe aus freier Energie ($F$), Druckenergie ($P$), kinetischer Energie ($K$) und Wandfreier Energie ($W$) leiten sie kinetische Gleichungen durch Minimierung der Energiedissipation ab.
• Modifizierte Gibbs-Duhem-Beziehung: Um mathematische Inkonsistenzen zwischen räumlichen Gradienten und funktionalen Ableitungen in bestehenden Modellen zu lösen, wird eine modifizierte isotherme Gibbs-Duhem-Beziehung abgeleitet. Diese berücksichtigt explizit geschwindigkeitsabhängige Entropie und Trägheitsarbeit, was zu einer rigorosen Definition des Drucks führt, der Beiträge aus freier Energie und Geschwindigkeitsentropie enthält.
• Definition der kinetischen Energie: Die kinetische Energie wird als $K = \int \rho u \cdot u \, d\Omega$ (linearer Impulsgehalt $\zeta = \rho u$) formuliert und nicht wie üblich als $\frac{1}{2}\rho u^2$. Diese Unterscheidung ermöglicht die Trennung von Impulsbilanz (Änderung von sowohl $u$ als auch $\rho$) und Geschwindigkeitsdissipation (Änderung nur von $u$).

Zentrale Beiträge
Die Ableitung liefert eine generalisierte Dichteevolutionsgleichung und drei wesentliche theoretische Fortschritte:

1. Generalisiertes Bernoulli-Gesetz: Das Standard-Bernoulli-Prinzip ($p + \frac{1}{2}\rho u^2 = \text{konst}$) wird als Spezialfall gezeigt. Die generalisierte Form berücksichtigt Änderungen der Geschwindigkeitsentropie und der Dichteentwicklung an der Grenzfläche: $p + \rho u^2 = \text{konstant}$ (unter spezifischen stationären, unverdrängten Bedingungen), wobei anerkannt wird, dass an Grenzflächen aufgrund von Dichteunterschieden Energiesprünge auftreten.
2. Universelle Schallgeschwindigkeitsformulierung: Eine allgemeine Expression für die Schallgeschwindigkeit wird bascheierend auf der Ableitung der freien Energie nach der Dichte hergeleitet: $c = \sqrt{d\Psi/d\rho}$. Für ideale Gase entspricht dies näherungsweise der mittleren thermischen Geschwindigkeit der Moleküle.
3. Generalisierte Zustandsgleichung (EOS): Die Theorie zeigt, dass die Dichteevolutionsgleichung natürlich eine Zustandsgleichung für Fluide generiert (z. B. $p \sim \rho \ln \rho$ für ideale Gase), ohne dass eine vorgegebene Zustandsgleichung erforderlich ist, was im Gegensatz zu Standardmethoden der diffusen Grenzfläche steht.

Ergebnisse und Validierung
Das Paper validiert die Theorie durch drei Fallstudien:

• Schallausbreitung: Die abgeleitete Schallgeschwindigkeit ist konsistent mit der Störungstheorie und sagt korrekt die Dispersionsrelation für Schallwellen in Luft voraus.
• Ideale Gas-Zustandsgleichung: Im stationären Zustand reproduziert das Modell das ideale Gasgesetz, was bestätigt, dass die Dichteevolutionsgleichung thermodynamische Zustandsbeziehungen inhärent erfasst.
• Dynamik der Wasser-Luft-Grenzfläche: Eine analytische Lösung für eine bewegte, flache Wasser-Luft-Grenzfläche in einem Rohr zeigt, dass das Modell ein nicht-monotones Dichteprofil vorhersagt. Konkret steigt die Dichte vom Bulk-Wasserwert zu einem maximalen Peak im Zentrum der Grenzfläche an, bevor sie zum Bulk-Gaswert abfällt. Dieses Ergebnis stimmt mit Erkenntnissen der Dichtefunktionaltheorie (DFT) [11, 12] überein und widerspricht den monoton fallenden Profilen, die durch lineare Interpolationsmethoden vorhergesagt werden. Das Modell reproduziert zudem korrekt die Young-Laplace-Gleichung für gekrümmte Grenzflächen.

Bedeutung
Das Paper beansprucht, einen allgemeinen Rahmen für unmischbare Fluide mit beliebig hohen Dichteverhältnissen zu bieten, der sowohl für kompressible als auch inkompressible Regime anwendbar ist. Indem die Dichteevolution als dynamische Variable behandelt wird, die mit dem Überschussvolumen und der viskosen Dissipation gekoppelt ist, löst das Modell die Widersprüche auf, die in der Boussinesq-Approximation und den linearen Interpolationsmethoden bestehen. Die Arbeit legt nahe, dass viskose Dissipation die Dichte inhärent verändert, was die strikte Annahme der Inkompressibilität in dissipativen Strömungen infrage stellt. Die Autoren postulieren, dass diese Theorie ein neues Fenster für die CFD öffnet und eine physikalisch konsistente Beschreibung des Impuls- und Energietransfers an Fluid-Fluid-Grenzflächen bietet, an denen die Dichteverhältnisse extrem sind.