Immiscible two-phase flow in porous media: a statistical mechanics approach

Dit artikel presenteert een overzicht van een statistisch-mechanistische aanpak voor immiscibele tweefasenstroming in poreuze media, die een hanteerbare continuüm-beschrijving biedt met nieuwe emergente variabelen zoals de meebewegende snelheid, en zo de beperkingen van de traditionele relatieve doorlaatbaarheidstheorie overbrugt.

De kern

De Stroom van Twee Vloeistoffen door een Zwam: Een Nieuwe Wiskundige Bril

Stel je voor dat je een enorme, dichte spons hebt. Nu gooi je niet alleen water door die spons, maar ook olie. De olie en het water willen niet mengen (ze zijn "onmengbaar"), en ze vechten om de kleine gaatjes in de spons. Soms stopt de olie, soms stroomt het water, en soms haken ze vast aan elkaar.

Voor wetenschappers is dit een enorme hoofdpijn. De oude manier om dit te beschrijven (de "Darcy-theorie") werkt goed voor simpele situaties, maar faalt als het echt ingewikkeld wordt. Het is alsof je probeert het verkeer in een drukke stad te beschrijven met alleen maar de snelheid van één auto, terwijl je eigenlijk duizenden auto's hebt die botsen, inhalen en vastlopen.

In dit artikel presenteren de auteurs, Alex Hansen en Santanu Sinha, een nieuwe manier om naar dit probleem te kijken. Ze gebruiken een idee uit de statistische mechanica (een tak van de natuurkunde die normaal gesproken gaat over atomen en hitte), maar dan toegepast op vloeistoffen die niet warmte uitwisselen.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het Probleem: De "Schaal" is de Sleutel

Stel je voor dat je door een bos loopt.

• Op micro-niveau (de poriën): Je ziet elke boom, elke tak en elke pad. Je ziet hoe een vogel (olie) probeert langs een tak (water) te vliegen. Dit is heel complex en chaotisch.
• Op macro-niveau (de Darcy-schaal): Je kijkt vanuit een helikopter. Je ziet alleen een groen veld. Je ziet niet de individuele bomen, maar je ziet wel dat het veld als geheel water vasthoudt of doorlaat.

Het oude probleem was: hoe vertaal je het gedrag van de individuele bomen (de micro-wereld) naar het gedrag van het hele veld (de macro-wereld)? De oude theorieën deden dit met "gokjes" (fenomenologische vergelijkingen) die vaak niet klopten.

2. De Oplossing: De "Informatie-thermodynamica"

De auteurs gebruiken een slimme truc bedacht door Edwin Jaynes in de jaren '50. In plaats van te kijken naar hitte en energie (zoals in normale thermodynamica), kijken ze naar informatie en onzekerheid.

• De Analogie: Stel je voor dat je een doos vol met rode en blauwe knikkers hebt die door een trillende machine worden geschud. Je weet niet precies waar elke knikker zit, maar je kent de kansen.
• In hun theorie kijken ze niet naar de krachten tussen de vloeistoffen, maar naar de kansverdeling van hoe de vloeistoffen zich in de poriën verdelen. Ze vragen: "Wat is de meest waarschijnlijke manier waarop deze vloeistoffen zich kunnen verdelen, gegeven wat we weten?"

Door deze "onzekerheid" (entropie) te maximaliseren, krijgen ze een nieuwe wiskundige structuur die er precies uitziet als thermodynamica, maar dan voor stroming.

3. Nieuwe Begrippen: De "Agitatie" en de "Meewegende Snelheid"

In deze nieuwe theorie ontstaan er nieuwe variabelen die in de oude theorie ontbraken.

A. Agitatie (Agiture) in plaats van Temperatuur
In normale thermodynamica is temperatuur een maat voor hoe snel de deeltjes trillen. In deze stroming is er geen hitte, maar er is wel druk.

• De Vergelijking: Stel je voor dat je de spons hard schudt. Hoe harder je schudt (hoe hoger de druk), hoe meer de vloeistoffen "onrustig" worden en hoe makkelijker ze door de obstakels komen.
• Ze noemen dit Agitatie. Het is een maat voor hoe "geagiteerd" of onrustig het systeem is door de stroming. Hoe meer agitatie, hoe makkelijker de vloeistoffen zich herschikken.

B. De Co-bewegende Snelheid (Co-moving Velocity)
Dit is misschien wel het belangrijkste nieuwe idee.

• Het Oude Denkbeeld: Je dacht dat als water en olie door de spons stromen, ze gewoon hun eigen snelheid hebben.
• Het Nieuwe Inzicht: De auteurs ontdekken dat de twee vloeistoffen niet onafhankelijk van elkaar bewegen. Ze "trekken" elkaar mee.
• De Vergelijking: Denk aan een drukke voetgangersbrug. Als je loopt, loop je niet alleen op je eigen tempo. Als er een grote groep mensen voor je loopt, moet je misschien sneller lopen om ze in te halen, of juist trager als je vastloopt. Maar er is ook een "meewegende stroom": als de hele menigte beweegt, beweeg jij mee, zelfs als je zelf niet extra kracht zet.
• In de theorie is er een extra snelheidcomponent, de co-moving velocity. Dit is de snelheid waarmee de ene vloeistof de andere "meesleept" of "beïnvloedt". Het is alsof de olie en het water een danspartner zijn die elkaar voortdurend duwen en trekken.

4. Waarom is dit belangrijk?

Deze theorie laat zien dat het gedrag van vloeistoffen in porieuze media (zoals grond, olie-reservoirs of filters) niet zomaar "rommelig" is. Het volgt strikte wiskundige regels, net als gas in een fles of water in een rivier.

• Glasachtige Fasen: Ze ontdekten dat er een punt is waarop de stroming "vastvriest" (zoals glas), zelfs als er nog steeds stroming is. Dit verklaart waarom sommige processen in de natuur zo onvoorspelbaar en "hysteretisch" zijn (waarbij de uitkomst afhangt van de geschiedenis).
• Toekomst: Met deze nieuwe "thermodynamische" bril kunnen ingenieurs betere modellen maken voor het winnen van olie, het zuiveren van water of het begrijpen van grondwaterstromen. Ze hoeven niet meer te gokken, maar kunnen berekenen op basis van de fundamentele wetten van informatie en onzekerheid.

Kortom:
De auteurs hebben een brug gebouwd tussen de chaotische wereld van individuele vloeistofdruppels in een spons en de grote, overzichtelijke wereld van stroming. Ze gebruiken de wiskunde van "kans en onzekerheid" om te bewijzen dat er een diepe, elegante orde schuilgaat in wat eruitzag als puur chaos. En ze ontdekten dat vloeistoffen in een spons elkaar meer "meevoeren" dan we ooit dachten.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Immiscible two-phase flow in porous media: a statistical mechanics approach" in het Nederlands.

Titel: Immiscibele tweefasestroming in poreuze media: een statistisch-mechanistische benadering

Auteurs: Alex Hansen en Santanu Sinha (PoreLab, NTNU, Noorwegen)

---

1. Het Probleem: Het Schaalvergrotingsprobleem (Scale-up)

De kernuitdaging in de fysica van immiscibele (niet-mengbare) tweefasestroming in poreuze media is het vinden van een juiste beschrijving van de stroming op schalen die groot genoeg zijn om het medium als een continuüm te beschouwen (de Darcy-schaal), terwijl deze beschrijving direct gekoppeld blijft aan de fysica op de porie-schaal.

• Huidige status: De bestaande aanpak is gebaseerd op fenomenologische vergelijkingen, specifiek de relatieve permeabiliteitstheorie (Muskat, Leverett, jaren '30-'40). Deze theorie splitst de stroming in twee fasen met eigen drukken ($p_w, p_n$) en snelheden, gekoppeld via een capillaire drukfunctie ($p_c$) en relatieve permeabiliteiten ($k_{rw}, k_{rn}$).
• Beperkingen: De centrale aanname is dat deze parameters uitsluitend afhangen van de verzadiging ($S_w$). Dit is een sterke, vaak onjuiste aanname die de complexiteit van de interface-dynamica en de geschiedenis van het systeem negeert. Pogingen om hierboven uit te komen hebben vaak geleid tot onbeheersbare complexiteit zonder praktische toepasbaarheid.

2. Methodologie: Generaliseerde Statistische Mechanica

De auteurs passen de generalisatie van statistische mechanica voor niet-thermische systemen toe, zoals oorspronkelijk voorgesteld door Edwin T. Jaynes in de jaren '50. In plaats van thermodynamische entropie gebaseerd op moleculaire beweging, gebruiken ze de Shannon-informatie-entropie als maatstaf voor onwetendheid over het systeem.

De aanpak:

1. Representatief Elementair Gebied (REA): Het systeem wordt gemodelleerd als een cilindrisch monster met een REA loodrecht op de stromingsrichting.
2. Configuraties: De toestand van het systeem wordt beschreven door een configuratie $X$ (verdeling van vloeistoffen en snelheidsveld).
3. Maximale Entropie: Er wordt een configuratie-entropie gedefinieerd ($S = -\int p(X) \ln p(X) dX$). Deze wordt gemaximaliseerd onder de beperkingen van bekende gemiddelden:
   • Totale volumestroom ($Q_p$)
   • Oppervlakte van de natte vloeistof ($A_w$)
   • Totale porie-oppervlakte ($A_p$)
4. Afleiding van een Thermodynamisch Formalisme: Door de entropie te maximaliseren, ontstaat er een verdeling van waarschijnlijkheid die leidt tot een formalisme dat wiskundig identiek is aan de thermodynamica, maar voor stroming.

3. Belangrijkste Bijdragen en Concepten

A. Niet-thermische Thermodynamica op Darcy-schaal

De methologie levert een nieuw formalisme op met geëmergeerde variabelen die analoog zijn aan thermodynamische grootheden:

• Agiture ($\theta$): Een variabele analoog aan temperatuur, gedefinieerd als $\theta \propto |\nabla p|$ (de grootte van het drukgradiënt). Het is een maat voor de "agitatietemperatuur" of de intensiteit van de stromingsdynamica.
• Stroom-afgeleide ($\mu$): Analoge aan chemisch potentiaal, gekoppeld aan het natte oppervlak (en dus de verzadiging).
• Stroom-druk ($\pi$): Een nieuwe intensieve variabele gekoppeld aan de porie-oppervlakte, zonder equivalent in gewone thermodynamica.

B. De Glas-overgang (Glass Transition)

Het artikel identificeert een fase-overgang in de stroming die overeenkomt met een glas-overgang:

• Bij lage capillaire getallen ($Ca$) is de stroming "bevroren" (fasen Ia/Ib), waarbij interfaces vastzitten of slechts intermitterend bewegen.
• Bij hogere $Ca$ wordt de stroming vloeibaar.
• De overgang wordt gekenmerkt door een piek in de spin-glas susceptibiliteit en een verandering in de relatie tussen snelheid en drukgradiënt van lineair naar een machtswet ($v \sim |\nabla p|^\beta$). Dit verklaart de waargenomen hysterese en grote fluctuaties in de stroming.

C. De Co-bewegende Snelheid (Co-moving Velocity, $v_m$)

Dit is de meest cruciale nieuwe variabele. De auteurs tonen aan dat de gemiddelde snelheid ($v_p$) niet simpelweg de som is van de individuele fasensnelheden ($v_w, v_n$) vermenigvuldigd met hun verzadiging. Er is een extra term nodig:
$$v_w = \hat{v}_w - S_n v_m$$
$$v_n = \hat{v}_n + S_w v_m$$
Waarbij $\hat{v}$ de "thermodynamische snelheden" zijn en $v_m$ de co-bewegende snelheid.

• $v_m$ beschrijft hoe de fasen ten opzichte van elkaar bewegen als de verzadiging verandert.
• De auteurs tonen aan dat $v_m$ een lineaire relatie heeft met de stroom-afgeleide $\mu$ (en dus met de verzadigingsafhankelijkheid van de stroomsnelheid).

4. Resultaten

1. Fasediagram: Er wordt een Darcy-schaal fasediagram gepresenteerd dat een glasachtige fase en een niet-glasachtige fase scheidt door een kritieke lijn. Deze lijn komt exact overeen met het begin van het machtswet-gedrag in de stroming.
2. Verband met gewone thermodynamica: De auteurs trekken een verrassende parallel tussen de co-bewegende snelheid in poreuze media en het co-molaire volume in mengsels van vloeistoffen (zoals water en ethanol). Net als het co-molaire volume corrigeert de co-bewegende snelheid voor de interactie tussen componenten die niet lineair optellen.
3. Oplossing voor de Relatieve Permeabiliteit: Het nieuwe formalisme biedt een manier om de stroming te beschrijven zonder de problematische aanname dat relatieve permeabiliteit en capillaire druk alleen van de verzadiging afhangen. De "driving forces" worden de gradiënten van agiture en stroom-afgeleide, wat de noodzaak van een aparte capillaire drukfunctie ($p_c$) in de vergelijkingen elimineert.

5. Betekenis en Conclusie

• Fundamenteel: Het artikel toont aan dat immiscibele tweefasestroming in poreuze media een uitstekend voorbeeld is van een niet-thermisch systeem dat toch de wiskundige structuur van de thermodynamica kan aannemen. Dit biedt een robuust theoretisch raamwerk dat de pore-schaal fysica koppelt aan de Darcy-schaal.
• Praktisch: De introductie van de co-bewegende snelheid ($v_m$) is essentieel voor het nauwkeurig modelleren van stroming, vooral in de glasachtige fase die relevant is voor veel praktische toepassingen (zoals oliewinning of CO2-opslag). Het verklaart waarom traditionele modellen vaak falen bij het voorspellen van hysterese en niet-lineariteiten.
• Toekomst: De auteurs suggereren dat de volgende stap het toepassen van niet-evenwichtsthermodynamica is om tijdsafhankelijke (transiënte) stromingen te beschrijven, waarbij de drukgradiënten worden vervangen door gradiënten in agiture en stroom-afgeleide.

Kortom, dit werk biedt een paradigmaverschuiving: van fenomenologische, vaak geforceerde vergelijkingen naar een fundamenteel, statistisch-mechanistisch onderbouwd formalisme dat de complexiteit van poreuze media op een beheersbare en fysiek consistente manier beschrijft.