Glassy phase transition in immiscible steady-state two-phase flow in porous media

Deze studie toont aan dat het gedrag van tweefasenstroom in poreuze media op macroscopische schaal succesvol kan worden voorspeld door de verdeling van druppels te mappen op een spin-glasmodel, waarbij de overgang naar een 'glazige' fase overeenkomt met een niet-lineair stromingsregime gekenmerkt door hysterese en sterke fluctuaties.

De kern

De Glazen Transformatie: Hoe Vloeistoffen in Poreuze Steen Net als een Verwarde Menigte Gedragen

Stel je voor dat je door een enorm, ondoordringbaar labyrint loopt. Dit labyrint is niet gemaakt van muren, maar van miljoenen kleine gaatjes in een steen (zoals zandsteen of grond). Door deze gaatjes stromen twee soorten vloeistoffen: water en olie (of een andere vloeistof die niet met water mengt).

Dit fenomeen heet twee-fasenstroom in poreuze media. Het klinkt als saaie natuurkunde, maar het is cruciaal voor dingen als het winnen van aardolie, het zuiveren van grondwater of het begrijpen van hoe regen door de bodem trekt.

Sinds meer dan een eeuw proberen wetenschappers een simpele regel te vinden die voorspelt hoe deze stromen zich gedragen op grote schaal (zoals in een heel olie-reservoir). Maar tot nu toe faalden ze. De regels op micro-niveau (in de kleine gaatjes) leken niets te zeggen over het gedrag op macro-niveau (in het hele reservoir).

Het Grote Geheim: De Spin-Glas Theorie

In dit paper ontdekken de auteurs (Santanu Sinha en zijn team) een verrassende manier om dit probleem op te lossen. Ze gebruiken een metafoor uit een heel ander vakgebied: de fysica van magneten.

Stel je een kamer vol met mensen voor. Iedereen heeft een klein bordje met een pijl: ofwel wijst hij naar links (-1) of naar rechts (+1).

• In een normale magnetische toestand (paramagnetisch) wijzen de mensen willekeurig, maar als je een sterke magneet (druk) op ze richt, wijzen ze allemaal in dezelfde richting.
• In een spin-glas toestand (de "glazen" fase) is het een chaos. Iedereen wil naar links, maar zijn buurman wil naar rechts. Ze zijn vastgekleefd in een verwarde positie. Ze kunnen niet bewegen, maar ze zijn ook niet netjes geordend. Het is een soort "bevroren chaos".

De auteurs hebben ontdekt dat de vloeistoffen in de steen zich precies zo gedragen als deze mensen met bordjes:

1. De Vloeistof als Spins: Ze hebben elke kleine druppel olie of water in de steen omgezet in een "spin" (een bordje).
2. De Machine Learning: Ze gebruikten een slim computerprogramma (een "Boltzmann Machine") om te leren hoe deze druppels zich gedragen. Het programma leerde de regels van de chaos in de steen en vertaalde ze naar de regels van de magneten.

De Drie Fasen van de Stroom

Het onderzoek toont aan dat er drie verschillende manieren zijn waarop deze vloeistoffen door de steen kunnen stromen, afhankelijk van hoe hard je ze duwt (de druk):

1. Fase I (De Bevroren Chaos / Het Glas):
   Als je de druk laag houdt, zijn de vloeistoffen vastgekleefd. Ze bewegen niet echt, of ze bewegen heel traag en onvoorspelbaar. Het is alsof de mensen in de kamer vastzitten in een menigte; ze kunnen niet vooruit, maar ze zijn ook niet stil. Dit is de glazen fase. Hier zie je grote schommelingen en "hysterie" (als je de druk verlaagt, komt het systeem niet terug naar de oude staat).

   • De ontdekking: De auteurs bewijzen dat deze "glazen" fase precies samenvalt met het punt waar de stroom van lineair (rechtlijnig) naar niet-lineair overgaat.

2. Fase II (De Power-Law Chaos):
   Als je de druk verhoogt, beginnen de vloeistoffen zich los te maken. Ze beginnen te bewegen, maar nog steeds in een chaotisch patroon. Ze vormen druppels die samensmelten en weer uit elkaar vallen. Dit is het gebied waar de stroom niet meer lineair is met de druk, maar een ingewikkelde kromme volgt.

3. Fase III (De Vlotte Stroom):
   Bij heel hoge druk bewegen alle vloeistoffen soepel. De chaos is weg, en de stroom volgt weer een simpele, rechte lijn.

De "Glassy" Transformatie

Het meest spannende deel van dit paper is de conclusie over Fase I.
De auteurs zeggen: "Dit is geen gewone stroming. Dit is een glazen toestand."

Net zoals glas (zoals vensterglas) eigenlijk een vloeistof is die zo traag is dat het lijkt op een vast stof, is deze stroming in de steen een vloeistof die vastzit in een verwarde, bevroren staat. De overgang van deze "glazen" staat naar de chaotische stroom (Fase II) is een glazen overgang.

Waarom is dit belangrijk?

Voor honderden jaren dachten ingenieurs dat ze deze stromingen met simpele formules konden beschrijven. Dit paper zegt: "Nee, het is veel complexer."
Het laat zien dat als je probeert olie of water uit de grond te halen, je soms in een "glazen" valstrik terechtkomt waar de vloeistoffen vastzitten in een verwarde menigte.

Samenvattend in één zin:
De auteurs hebben ontdekt dat vloeistoffen die door poreuze steen stromen, zich op een heel specifieke manier gedragen die precies hetzelfde is als de verwarde magneten in een "spin-glas", en dat dit verklaart waarom stroming soms vastloopt in een chaotische, bevroren staat voordat het weer gaat vloeien.

Het is een briljante brug geslagen tussen de wereld van vloeistoffen in de grond en de wereld van magneten in de kwantumwereld, met behulp van slimme computers die de taal van de natuur hebben vertaald.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Glassy phase transition in immiscible steady-state two-phase flow in porous media" in het Nederlands.

Probleemstelling

Twee-fasenstroom in poreuze media is een alomtegenwoordig verschijnsel dat al meer dan een eeuw wordt bestudeerd. Desondanks ontbreekt er een succesvolle theorie die stroming op macroscopische schaal (de zogenaamde Darcy-schaal) kan voorspellen op basis van een "microscopisch" model. Hoewel fractale patronen op porieniveau in de jaren 80 veel aandacht kregen, bieden deze weinig inzicht in het gedrag op grotere schaal, wat cruciaal is voor toepassingen zoals hydrologie en reservoirtechniek.

Een specifiek aandachtspunt is de overgang tussen verschillende stromingsregimes. In de literatuur worden drie regimes onderscheiden op basis van het verband tussen stroomsnelheid en drukgradiënt:

1. Regime I: Lineair gedrag (onderverdeeld in Ia, waar interfaces bevroren zijn, en Ib, waar interfaces bewegen maar geen transport veroorzaken).
2. Regime II: Een niet-lineaire machtswet (power-law).
3. Regime III: Terugkeer naar lineair gedrag bij zeer hoge capillaire getallen.

De aard van de overgang tussen regime Ib en regime II is onduidelijk. Het doel van dit onderzoek is om te bepalen of deze overgang een echte faseovergang is, en zo ja, van welk type.

Methodologie

De auteurs hanteren een innovatieve benadering door het stromingsprobleem te mappen op een statistisch mechanisch model:

1. Dynamisch Poreus Netwerk (DPN) Model:

   • Ze gebruiken een 2D diamantrooster van 40x40 links (poriën) met periodieke randvoorwaarden.
   • Het model simuleert onmengbare twee-fasenstroom (bevochtigend en niet-bevochtigend) waarbij zowel viskeuze als capillaire krachten worden meegenomen.
   • De stroming wordt gestuurd door een globaal drukverschil, en de configuraties worden vastgelegd in een stationaire toestand.

2. Mapping naar een Spin-Glas Model:

   • De verzadigingsconfiguratie van de vloeistoffen in het netwerk wordt gemappt op een spin-systeem ($\sigma_i = \pm 1$). Een spin is +1 als de niet-bevochtigende vloeistof de link domineert, en -1 anderszins.
   • Met behulp van het Jaynes Maximum Entropy Principe wordt een configuratiekansverdeling $P(\{\sigma\})$ afgeleid die overeenkomt met een spin-glas Hamiltoniaan:
     $$H(\{\sigma\}) = -\sum_i h_i \sigma_i - \sum_{i<j} J_{ij} \sigma_i \sigma_j$$
     Hierbij zijn $h_i$ lokale veldkoppelingen en $J_{ij}$ spin-spin koppelingen.

3. Boltzmann Machine Learning (BLM):

   • Om de parameters ($h_i$ en $J_{ij}$) te bepalen die de DPN-simulaties het beste reproduceren, wordt een inverse Ising-modellering (Boltzmann Machine Learning) gebruikt.
   • De auteurs minimaliseren de Kullback-Leibler-divergentie tussen de data uit de DPN-simulaties en Monte Carlo (MC) simulaties van het spin-model.
   • Ze verifiëren het model door te controleren of het niet alleen tweepunts-correlaties, maar ook driedpunts-correlaties correct voorspelt.

4. Fasediagram Constructie:

   • Voor een breed scala aan capillaire getallen ($Ca$) en niet-bevochtigende verzadigingen ($S_n$) worden thermodynamische grootheden gemeten: magnetisatie ($m$), Edwards-Anderson ordeparameter ($q$), uniforme susceptibiliteit ($\chi_m$) en spin-glas susceptibiliteit ($\chi_{sg}$).

Belangrijkste Resultaten

• Validatie van het Model: Het trainde spin-model reproduceert de stromingspatronen uit de DPN-simulaties zeer nauwkeurig, inclusief hogere-orde correlaties. Dit bevestigt dat de dissipatieve dynamiek van de stroming effectief kan worden beschreven door een evenwichtsstatistische mechanica van een spin-glas.
• Identificatie van de Glasovergang:
  • Bij lage capillaire getallen ($Ca$) vertoont het systeem een spin-glas fase (SG), gekenmerkt door een nul totale magnetisatie ($m \approx 0$) maar een hoge lokale magnetisatie ($q > 0$).
  • Bij hoge $Ca$ bevindt het systeem zich in een paramagnetische fase (P) met zowel $m \approx 0$ als $q \approx 0$.
  • De overgang tussen deze fasen wordt gemarkeerd door een piek in de spin-glas susceptibiliteit ($\chi_{sg}$).
• Koppeling aan Darcy-schaal:
  • De kritische lijn die de paramagnetische fase scheidt van de spin-glas fase in het fasediagram ($Ca$ vs. $S_n$) valt exact samen met de overgang van regime Ib naar regime II op de Darcy-schaal.
  • Regime Ib (de lineaire regime met bewegende interfaces) correspondeert met de glasachtige stromingsfase.
  • Regime II (de niet-lineaire machtswet) correspondeert met de paramagnetische fase.
• Fysische Interpretatie:
  • De glasachtige fase (Ib) wordt gekenmerkt door sterke fluctuaties, hysterese en een breed scala aan tijdschalen. Dit komt overeen met het "kinetisch gevangen" gedrag van vloeistofclusters (ganglions) die moeilijk mobiliseerbaar zijn door capillaire krachten.
  • De overgang naar regime II is een echte kritische faseovergang (glasovergang), waarbij de stroming overgaat van een bevroren, fluctuerende toestand naar een toestand waar interfaces geleidelijk mobiliseren en de stroming niet-lineair wordt.

Bijdrage en Significantie

1. Theoretische Doorbraak: Het artikel biedt een succesvolle brug tussen microscopische poreuze stroming en macroscopische Darcy-beschrijvingen door gebruik te maken van concepten uit de gecondenseerde materie-fysica (spin-glas). Het lost het probleem op van het voorspellen van Darcy-gedrag op basis van poreuze schaal-fysica.
2. Nieuwe Classificatie: Het identificeert regime Ib als een echte "glasachtige stromingsfase" en de overgang Ib-II als een glasovergang. Dit verklaart waarom dit regime hysterese en grote fluctuaties vertoont, kenmerken die typisch zijn voor glasachtige systemen.
3. Methodologische Innovatie: Het toont aan dat machine learning (Boltzmann Machine Learning) en het Maximum Entropy Principe effectief kunnen worden gebruikt om effectieve Hamiltonianen af te leiden voor complexe, niet-evenwichtssystemen zoals twee-fasenstroom.
4. Praktische Implicaties: Voor toepassingen in reservoirengineering en bodemkunde biedt dit inzicht in de voorwaarden waaronder stroming onvoorspelbaar wordt (hysterese) en wanneer de overgang naar een meer gestructureerde, niet-lineaire stroming plaatsvindt. Dit is cruciaal voor het modelleren van vloeistoftransport in complexe geologische formaties.

Kortom, het paper bewijst dat de complexe dynamiek van onmengbare twee-fasenstroom in poreuze media fundamenteel kan worden begrepen als een glasovergang, waarbij de "glasachtige" toestand overeenkomt met een regime van bevroren maar fluctuerende vloeistofpatronen.