Compressible and immiscible fluids with arbitrary density ratio

Dit artikel presenteert een nieuwe thermodynamische theorie voor het modelleren van samendrukbare en immiscibele vloeistoffen met arbitrair hoge dichtheidsverhoudingen, waarbij de beperkingen van traditionele Navier-Stokes- en Euler-vergelijkingen worden aangepakt door dichtheidsevolutievergelijkingen af te leiden die de werkelijke momentumevolutie correct verantwoorden.

De kern

Stel je voor dat je een dans probeert te beschrijven tussen twee zeer verschillende partners: een zware, traag bewegende olifant (water) en een lichte, snel bewegende veer (lucht). In de wereld van de vloeistofmechanica zijn deze twee "immiscibel", wat betekent dat ze niet mengen zoals olie en water, maar ze hebben wel een vage grens waar ze elkaar ontmoeten.

Al een lange tijd worstelen wetenschappers met het schrijven van de "dansregels" (de vergelijkingen) voor systemen waar het dichtheidsverschil enorm is, zoals wanneer de olifant 1.000 keer zwaarder is dan de veer. De oude regels hadden een groot gebrek: ze behandelden het gewicht van de olifieant alsof het een constante, onveranderlijke waarde was, zelfs precies op de vage grens waar de olifant een veer wordt. Dit is alsof je probeert te beschrijven dat iemand in een geest verandert door te zeggen: "Ze blijven de hele tijd een solide mens," wat geen zin heeft.

Hier is een eenvoudige uiteenzetting van wat dit artikel voorstelt om dat probleem op te lossen.

1. Het probleem met de oude regels

De traditionele manier om te berekenen hoe vloeistoffen bewegen (met behulp van de Navier-Stokes- en Euler-vergelijkingen) vertrouwt op een afkorting die de Boussinesq-benadering wordt genoemd.

• De analogie: Stel je voor dat je een winkelwagentje duwt. Als het wagentje vol bakstenen zit, is het zwaar. Als het leeg is, is het licht. De oude regels gingen ervan uit dat als je een wagentje duwt dat deels gevuld is met bakstenen en deels met lucht, het gewicht van het wagentje nooit verandert terwijl je duwt. Het gaat er simpelweg vanuit dat het gewicht een vast gemiddelde is.
• Het gebrek: In werkelijkheid, naarmate het wagentje beweegt en de bakstenen verschuiven (diffusie), verandert het gewicht. De oude regels negeerden het feit dat het "impulsmoment" (massa $\times$ snelheid) verandert omdat de massa zelf verandert bij de grens. Ze gingen er ook vanuit dat de lucht en het water een perfect voorspelbare hoeveelheid ruimte innemen, waarbij ze negeerden dat wanneer ze op de grens mengen, de ruimte die ze innemen daadwerkelijk kan trillen en veranderen door druk.

2. De nieuwe aanpak: Energieminimalisatie

In plaats van te gokken hoe de dichtheid verandert, begint de auteur vanuit een fundamenteel principe: De natuur probeert altijd de minste hoeveelheid energie te gebruiken.

• De analogie: Denk aan een bal die een heuvel afrolt. De bal geeft niet om de specifieke stappen die hij neemt; hij wil gewoon naar de bodem (de laagste energie). De auteur gebruikt dit "energieheuvel"-concept om nieuwe regels af te leiden voor hoe het water en de lucht met elkaar interageren.
• De belangrijkste innovatie: De auteur introduceert een concept genaamd "Excess Volume" (overschotvolume).
  • Stel je voor dat je een emmer water en een emmer lucht hebt. Als je ze bij elkaar giet, zou je verwachten dat het totale volume exact de som is van de twee emmers. Maar op microscopisch niveau, wanneer ze elkaar ontmoeten, kunnen de moleculen dichter op elkaar gepakt zitten of juist losser, wat een beetje "extra" of "ontbrekende" ruimte creëert.
  • De oude regels gingen ervan uit dat deze extra ruimte overal nul was. Dit artikel zegt: "Nee, die extra ruimte bestaat, deze verandert van plaats naar plaats, en het beïnvloedt de dichtheid."

3. De nieuwe "Dansregels" (De resultaten)

Door rekening te houden met deze veranderende "extra ruimte" en door energieminimalisatie te gebruiken, leidt de auteur een nieuwe set vergelijkingen af die drie hoofdzaken doen:

A. Een nieuwe manier om geluid te horen (Geluidsnelheid)
Het artikel laat zien dat de snelheid van het geluid niet zomaar een willekeurig getal is; het komt rechtstreeks voort uit hoe de energie van de vloeistof verandert wanneer deze wordt samengedrukt.

• De metafoor: Denk aan geluid als een rimpeling in een menigte. De snelheid van die rimpeling hangt af van hoe dicht de mensen (moleculen) op elkaar gepakt zitten en hoeveel energie ze hebben. De nieuwe formule berekent deze snelheid op een natuurlijke manier, zonder dat vooraf moet worden verteld wat het is. Het suggereert zelfs dat in een gas de geluidssnelheid ongeveer gelijk is aan de gemiddelde snelheid waarmee de gasmoleculen rondstuiteren.

B. Een nieuwe regel voor druk en snelheid (Wet van Bernoulli)
Je hebt waarschijnlijk wel eens gehoord van het principe van Bernoulli: "Wanneer een vloeistof sneller beweegt, daalt de druk."

• De draai: De oude regel werkt uitstekend voor water dat door een pijp stroomt, maar het loopt vast wanneer er een enorme sprong in dichtheid is (zoals water dat de lucht raakt). De auteur creëert een Gegeneraliseerde Wet van Bernoulli.
• De metafoor: Stel je een rivier voor die een waterval in stroomt. De oude regel zegt dat de energie gelijk blijft. De nieuwe regel zegt: "Wacht even, naarmate het water in nevel verandert (lucht), gaat er wat energie verloren of wordt deze getransformeerd omdat de 'volheid' van het water verandert." De nieuwe vergelijking houdt rekening met deze energieverschuiving, waardoor het nauwkeurig blijft, zelfs wanneer de vloeistof van aard verandert van zwaar naar licht.

C. De "Bump" in de dichtheid
Dit is misschien wel het meest visuele resultaat.

• Het oude beeld: Als je naar de grens tussen water en lucht kijkt, zeiden de oude modellen dat de dichtheid gewoon vloeiend zou afnemen van "zwaar water" naar "lichte lucht", als een helling die naar beneden loopt.
• Het nieuwe beeld: De wiskunde van de auteur voorspelt een bump (bult). Terwijl je de grens oversteekt, gaat de dichtheid eigenlijk eerst iets omhoog voordat deze daalt naar het niveau van de lucht.
• De metafoor: Stel je een menigte mensen (water) voor die probeert door een deur naar een lege gang (lucht) te glippen. Terwijl ze door de deuropening persen, kunnen ze voor een fractie van een seconde dichter op elkaar gepakt raken voordat ze zich verspreiden. De nieuwe theorie voorspelt deze "packing bump", wat overeenkomt met wat geavanceerde computersimulaties (Density Functional Theory) al hebben gezien, maar wat de oude eenvoudige modellen misten.

Samenvatting

Dit artikel stelt een nieuwe manier voor om de natuurwetten te schrijven voor vloeistoffen die erg verschillen in gewicht (zoals water en lucht).

1. Het stopt met de suggestie dat het gewicht constant is bij de grens.
2. Het geeft toe dat de ruimte die moleculen innemen verandert (excess volume).
3. Het gebruikt het principe van "laagste energie" om nieuwe regels af te leiden die verklaren hoe geluid reist, hoe druk verandert en waarom de dichtheid bij de water-luchtgrens daadwerkelijk een kleine "heuvel" vertoont.

De auteur beweert dat dit nieuwe kader werkt voor elke mengeling van vloeistoffen, ongeacht hoe verschillend hun gewichten zijn, wat de deur opent naar nauwkeurigere computersimulaties van zaken als vallende regen, brekende golven of opstijgende bellen.

Technische samenvatting

Technische Samenvatting: Samendrukbare en Onmengbare Vloeistoffen met Willekeurige Dichtheidsverhoudingen: Diffuse Interface Theorie

Probleemstelling
Huidige computationele vloeistofdynamica (CFD) modellen worstelen met het nauwkeurig simuleren van onmengbare vloeistofsystemen met hoge dichtheidsverhoudingen, zoals water-lucht (verhouding $\approx 1000$). Bestaande benaderingen vertrouwen op twee primaire methoden, die beide aanzienlijke beperkingen kennen:

1. De Boussinesq-benadering: Veronderstelt een constante dichtheid ($\rho$) binnen de interface, waarbij gebruik wordt gemaakt van de vergelijking $\rho(du/dt) = -\nabla p + \dots$. Deze formulering houdt geen rekening met de "werkelijke" impulsontwikkeling $d(\rho u)/dt$, waardoor dichtheidsveranderingen en de resulterende momentumoverdracht worden genegeerd. Bovendien kan deze methode niet op natuurlijke wijze zwaartekracht- en drijfvermogenskrachten genereren zonder ad-hoc toevoegingen, aangezien deze voortkomen uit dichtheidsverschillen.
2. De Interpolatiemethode: Definieert dichtheid als een lineaire functie van volumconcentratie ($\rho = \rho_w \phi + \rho_a(1-\phi)$). Deze benadering creëert een fundamenteel conflict tussen onsamendrukbaarheid ($d\rho/dt = 0$) en diffusie ($d\phi/dt \neq 0$). Bovendien gaat het uit van een uniforme overtollige volume (vaak nul) over het gehele domein, wat in strijd is met de fysieke realiteit waar lokaal overtollig volume afhangt van de lokale druk. Bijgevolg slagen standaard continuïteitsvergelijkingen en lineaire interpolatie er niet in om de niet-monotone dichtheidsprofielen te vangen die worden waargenomen bij het water-lucht interface.

Methodologie
De auteurs stellen een nieuw theoretisch kader voor dat gebaseerd is op het thermodynamische principe van energiedisperatie (minimalisatie), wat verschilt van de traditionele afleiding van de Navier-Stokes- en Euler-vergelijkingen. De methodologie omvat:

• Herformulering van Volume en Dichtheid: In plaats van een nul overtollig volume aan te nemen, introduceert het model een ruimtelijk variërend overtollig volume ($v_e$) binnen de interface. Dit volume wordt herverdeeld tussen de vloeistofcomponenten, wat leidt tot partiële dichtheden ($\rho'_w, \rho'_a$) die niet langer constant zijn maar afhankelijk zijn van lokale druk en menging. Dit vereist twee gekoppelde evolutie-vergelijkingen: één voor de volumconcentratie ($\phi$) en één voor de dichtheid ($\rho$).
• Dissipatie-Conserveringsprincipe: De auteurs stellen een stelling vast die de eerste wet van de thermodynamica (conservering) verbindt met de tweede wet (dissipatie). Door de totale systeemenergie ($E$) te definiëren als de som van vrije energie ($F$), drukenergie ($P$), kinetische energie ($K$) en wandvrije energie ($W$), leiden zij kinetische vergelijkingen af door de energie-dissipatie te minimaliseren.
• Gemodificeerde Gibbs-Duhem Relatie: Om wiskundige incompatibiliteiten tussen ruimtelijke gradiënten en functionele afgeleiden in bestaande modellen op te lossen, is een gemodificeerde isotherme Gibbs-Duhem relatie afgeleid. Deze houdt expliciet rekening met snelheid-afhankelijke entropie en inertieel werk, wat leidt tot een rigoureuze definitie van druk die bijdragen bevat van vrije energie en snelheid-entropie.
• Definitie van Kinetische Energie: De kinetische energie wordt geformuleerd als $K = \int \rho u \cdot u \, d\Omega$ (lineaire momentumdichtheid $\zeta = \rho u$) in plaats van de conventionele $\frac{1}{2}\rho u^2$. Dit onderscheid maakt het mogelijk om momentumbalans (verandering van zowel $u$ als $\rho$) te scheiden van snelheid-dissipatie (verandering van enkel $u$).

Belangrijkste Bijdragen
De afleiding levert een gegeneraliseerde dichtheidsvergelijking op en drie belangrijke theoretische vooruitgang:

1. Gegeneraliseerde Wet van Bernoulli: Het standaard Bernoulli-principe ($p + \frac{1}{2}\rho u^2 = \text{const}$) wordt getoond als een speciaal geval. De gegeneraliseerde vorm houdt rekening met snelheid-entropie veranderingen en dichtheidsontwikkeling bij de interface: $p + \rho u^2 = \text{constant}$ (onder specifieke stationaire onge unimpactieve condities), waarbij erkend wordt dat er energie-sprongen optreden over interfaces door dichtheidsverschillen.
2. Universele Geluidssnelheid Formulering: Een algemene expressie voor de geluidssnelheid is afgeleid op basis van de afgeleide van de vrije energie ten opzichte van de dichtheid: $c = \sqrt{d\Psi/d\rho}$. Voor ideale gassen benadert dit de gemiddelde thermische snelheid van moleculen.
3. Gegeneraliseerde Toestandsvergelijking (EOS): De theorie demonstreert dat de dichtheidsontwikkelingsvergelijking op natuurlijke wijze een EOS voor vloeistoffen genereert (bijv. $p \sim \rho \ln \rho$ voor ideale gassen) zonder dat een voorgeschreven EOS vereist is, wat contrasteert met standaard diffuse-interface methoden.

Resultaten en Validatie
Het artikel valideert de theorie aan de hand van drie casestudies:

• Geluidspropagatie: De afgeleide geluidssnelheid is consistent met perturbatieanalyse en voorspelt correct de dispersierelatie van geluidsgolven in lucht.
• Ideale Gas EOS: In de stationaire toestand reproduceert het model de ideale gaswet, wat bevestigt dat de dichtheidsontwikkelingsvergelijking inherent thermodynamische toestandsrelaties vastlegt.
• Dynamica van het Water-Lucht Interface: Een analytische oplossing voor een bewegend plat water-lucht interface in een buis laat zien dat het model een niet-monotoon dichtheidsprofiel voorspelt. Specifiek neemt de dichtheid toe van de bulkwaterwaarde naar een maximale piek in het centrum van de interface voordat deze afneemt naar de bulkgaswaarde. Dit resultaat komt overeen met bevindingen uit de Density Functional Theory (DFT) [11, 12] en spreekt de monotone afname tegen die door lineaire interpolatiemethoden wordt voorspeld. Het model reproduceert ook correct de Young-Laplace vergelijking voor gebogen interfaces.

Betekenis
Het artikel beweert een algemeen kader te bieden voor onmengbare vloeistoffen met willekeurig hoge dichtheidsverhoudingen, toepasbaar op zowel samendrukbare als onsamendrukbare regimes. Door de ontwikkeling van dichtheid te behandelen als een dynamische variabele gekoppeld aan overtollig volume en viskeuze dissipatie, lost het model de tegenstrijdigheden op die inherent zijn aan de Boussinesq-benadering en lineaire interpolatiemethoden. Het werk suggereert dat viskeuze dissipatie de dichtheid inherent verandert, wat de strikte onsamendrukbaarheidsveronderstelling uitdaagt in dissipative stromingen. De auteurs stellen dat deze theorie een nieuw venster opent voor CFD, door een fysiek consistente beschrijving te bieden van momentum- en energieoverdracht bij vloeistof-vloeistof interfaces waar de dichtheidsverhoudingen extreem zijn.