Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media

Deze studie presenteert een rigoureus afgeleid macroscopisch model voor de stroming van twee niet-mengbare vloeistoffen in inhomogeen poreus medium, waarbij de volume-averagetechniek wordt gebruikt om de Navier-Stokes- en Cahn-Hilliard-vergelijkingen op poreus schaal te upscalen naar Darcy-schaal, inclusief een formele integratie van natgedrag en numerieke validatie.

De kern

De Grote Reis van Vloeistoffen: Hoe een Nieuw Model de Stroom door Poreuze Steen Begrijpt

Stel je voor dat je door een enorm, ingewikkeld labyrint loopt. Dit labyrint is niet gemaakt van muren, maar van miljoenen kleine steentjes met gaatjes ertussen. Dit noemen we porieus materiaal (zoals zand, steen of zelfs een spons). Nu laat je twee soorten vloeistof door dit labyrint stromen: bijvoorbeeld water en olie. Ze mogen niet mengen, maar ze duwen elkaar wel weg.

De vraag die wetenschappers al jaren bezighoudt, is: Hoe gedragen deze vloeistoffen zich als we niet naar elk individueel gaatje kijken, maar naar het hele labyrint als één groot geheel?

Dit artikel van Chunhua Zhang en zijn team geeft een antwoord op die vraag met een slimme nieuwe manier van kijken. Hier is de uitleg in simpele taal:

1. Het Probleem: Te Klein of Te Groot?

Er zijn twee manieren om naar dit probleem te kijken:

• De Microscopische Manier (Pore-schaal): Je kijkt door een supersterke microscoop naar één klein gaatje. Je ziet precies hoe de olie en het water langs de wanden glijden, hoe ze plakkend zijn (natmakend) en hoe ze botsen. Dit is heel accuraat, maar als je een heel olie-reservoir wilt simuleren, heb je miljarden van deze gaatjes. Dat is voor computers te zwaar om te berekenen.
• De Macroscopische Manier (Darcy-schaal): Je kijkt naar het hele labyrint als één grote, vage vloeistof. Je gebruikt simpele regels (zoals "water stroomt sneller door droge grond"). Dit is makkelijk te rekenen, maar het mist de details. Het weet bijvoorbeeld niet precies hoe de "plakkerigheid" (natmakendheid) van de steen de stroom beïnvloedt.

De uitdaging: Hoe maak je een brug tussen deze twee werelden? Hoe vertaal je de complexe details van de kleine gaatjes naar een simpele regel voor het grote geheel, zonder de waarheid te verliezen?

2. De Oplossing: De "Gemiddelde Chef"

De auteurs gebruiken een techniek die ze volume-averaging noemen. Stel je voor dat je een grote kom met een mengelmoes van water, olie en steentjes hebt. In plaats van elke druppel te tellen, nemen ze een "gemiddelde chef" die kijkt naar wat er gemiddeld gebeurt in een stukje van de kom.

Ze hebben een wiskundig recept (het Navier-Stokes-Cahn-Hilliard model) gebruikt om de beweging van de vloeistoffen te beschrijven. Vervolgens hebben ze dit recept "opgeblazen" (upscaled) naar het grote niveau.

3. Het Nieuwe Trucje: De "Plakkracht" in de Formule

Het echte genie van dit onderzoek zit in hoe ze omgaan met natmakendheid (wetting).

• Oude manier: In de oude simpele modellen werd de "plakkracht" van de vloeistof aan de steen vaak als een losse, geschatte waarde toegevoegd. Het was alsof je een recept zegt: "Voeg een beetje plakkracht toe," zonder te weten waarom.
• Nieuwe manier: De auteurs hebben bewezen hoe je die plakkracht rechtstreeks uit de wiskunde van de kleine gaatjes haalt. Ze hebben een nieuwe term in hun formule bedacht die we chemisch potentieel noemen.
  • Vergelijking: Stel je voor dat de chemische potentieel een kompas is. In de oude modellen keek het kompas alleen naar de druk. In dit nieuwe kompas wijst de naald ook naar de "plakkerigheid" van de wanden. Als de wanden water liefhebben (hydrofiel), wijst het kompas in die richting. Als ze water haten (hydrofoob), wijst het juist de andere kant op.

Dit betekent dat hun nieuwe formule automatisch begrijpt hoe de steen de vloeistof beïnvloedt, zonder dat je handmatig moet gokken.

4. Wat hebben ze bewezen?

Om te laten zien dat hun nieuwe "gemiddelde chef" het goed doet, hebben ze drie proeven gedaan:

1. De Stroomtest: Ze lieten vloeistof door een rechte buis met steentjes stromen. Hun nieuwe formule gaf precies hetzelfde antwoord als de bekende wiskundige theorieën. (De chef kan goed tellen).
2. De Schokgolf (Buckley-Leverett): Ze lieten water een olieveld binnenstomen. Ze zagen hoe de water-front zich verplaatste. Hun model voorspelde precies hoe snel dit ging en hoe scherp de overgang was, net zoals in de echte wereld.
3. De Vingers (Viscous Fingering): Dit is het coolste deel. Als je dunne olie door dikke olie duwt, ontstaan er vaak vreemde, vingerachtige patronen (alsof de olie "pootjes" heeft).
   • Ze lieten zien dat als de steen hydrofiel is (water-liefhebbend), deze "vingers" worden onderdrukt en de vloeistof stroomt rustiger.
   • Als de steen hydrofoob is (water-hatend), worden de vingers wilder en onvoorspelbaarder.
   • Hun model kon dit gedrag perfect nabootsen, wat laat zien dat ze de "natmakendheid" echt goed hebben gevangen.

5. Waarom is dit belangrijk voor de wereld?

Dit onderzoek is niet alleen leuk wiskunde. Het heeft grote gevolgen voor:

• Olie- en gaswinning: Om meer olie uit de grond te halen, moet je weten hoe water de olie door de steen duwt.
• CO2-opslag: Om CO2 veilig in de grond op te slaan, moet je begrijpen hoe het zich verplaatst door poreus gesteente.
• Grondwaterzuivering: Om vervuiling te verwijderen, moet je weten hoe chemicaliën door de bodem reizen.

Conclusie:
De auteurs hebben een brug gebouwd tussen de microscopische chaos van individuele gaatjes en de makkelijke regels voor grote projecten. Ze hebben een nieuwe "taal" ontwikkeld die de plakkracht van de steen automatisch in de berekening verwerkt. Hierdoor kunnen ingenieurs in de toekomst veel nauwkeuriger voorspellen hoe vloeistoffen zich gedragen in de ondergrond, zonder dat ze een supercomputer nodig hebben om elk klein gaatje te simuleren.

Het is alsof ze een vertaler hebben gevonden die de complexe gedachten van een individuele druppel water perfect kan vertalen naar een helder commando voor het hele leger.

Technische samenvatting

Titel: Opschaling van het Navier-Stokes-Cahn-Hilliard-model voor incompressibele meerfasestroming in inhomogeen poreus medium

Auteurs: Chunhua Zhang, Peiyao Liu, Cheng Peng, Lian-Ping Wang, Zhaoli Guo.

---

1. Probleemstelling

Tweefasestromingen door poreuze media komen veel voor in toepassingen zoals grondwaterzuivering, verbeterde oliewinning en CO2-opslag. Het modelleren van deze stromingen is uitdagend vanwege de complexe geometrie van het poreuze medium en de enorme schaalverschillen:

• Poreuze schaal: De stroming wordt beheerst door de volledige Navier-Stokes-vergelijkingen en de fase-interface-evolutie door de Cahn-Hilliard-vergelijking. Directe numerieke simulatie (DNS) is hier mogelijk, maar te rekenintensief voor grote industriële schalen.
• Darcy-schaal (Macroscopisch): Voor engineeringtoepassingen zijn gemiddelde modellen nodig (Darcy-schaal). Bestaande empirische modellen (zoals de Darcy-wet met relatieve permeabiliteit en capillaire druk) hebben echter beperkingen:
  • Ze zijn vaak gebaseerd op experimentele fits in plaats van strikte afleidingen uit fysische wetten.
  • Ze behandelen capillaire druk en relatieve permeabiliteit vaak als functies van verzadiging alleen, wat de hysterese en het vochtgedrag (wetting) onvoldoende vastlegt.
  • Er is een gebrek aan een theoretisch raamwerk dat poreuze dynamiek consistent koppelt aan macroscopische beschrijvingen, met name voor wat betreft het vochtgedrag aan de wanden.

2. Methodologie

De auteurs gebruiken de volume-averagetechniek (volume averaging method) om de poreuze schaalvergelijkingen op te schalen naar de Darcy-schaal.

• Poreuze Schaal Model: Het systeem wordt beschreven door het gekoppelde Navier-Stokes-Cahn-Hilliard (NSCH) systeem voor twee onmengbare, incompressibele vloeistoffen.
  • De fase-interface wordt gemodelleerd met een ordeparameter $\phi$ (diffuse interface theorie).
  • De vrije-energiefunctional omvat volumetrische energie, interfaciale energie en wandvrije energie (voor vochtgedrag).
  • Randvoorwaarden omvatten een niet-slip voorwaarde voor snelheid en een Neumann-voorwaarde voor het chemisch potentieel, waarbij de vochtgrensvoorwaarde (wetting boundary condition) expliciet wordt meegenomen.
• Opschaling (Upscaling):
  • Er wordt een representatief elementair volume (REV) gedefinieerd dat voldoet aan schaalbeperkingen ($l \ll r_0 \ll L$, waarbij $l$ de poriegrootte is en $L$ de macroscopische schaal).
  • De vergelijkingen worden gemiddeld over het REV. De grootheden worden ontbonden in een intrinsiek gemiddeld deel en een ruimtelijke afwijking (Gray's decomposition).
  • Sluiting (Closure): De afgeleide vergelijkingen bevatten ongesloten termen (afwijkingen). Deze worden gesloten door:
    1. Lokale evenwichtshypothese (chemisch potentieel is constant binnen een REV).
    2. Het oplossen van lokale "closure-problemen" op het REV om de afwijkingen ($\hat{\phi}, \hat{u}$) te relateren aan de gemiddelde grootheden.
    3. Het modelleren van de wrijvingskracht (drag force) als een lineair (Darcy) en kwadratisch (Forchheimer) term.
• Numerieke Implementatie: De afgeleide macroscopische vergelijkingen worden opgelost met een Lattice Boltzmann Methode (LBM) die is aangepast voor de REV-schaal.

3. Belangrijkste Bijdragen

1. Formele Afleiding van Capillaire Effecten: In tegenstelling tot eerdere werken die capillaire druk empirisch invoeren, leiden de auteurs het gemiddelde chemisch potentieel strikt af. Hierin wordt de vochtgrensvoorwaarde (wetting boundary condition) direct geïntegreerd. Dit resulteert in een term die de capillaire krachten en het vochtgedrag intrinsiek koppelt aan de macroscopische verzadiging.
2. Unificatie van Schalen: Het model vormt een brug tussen de poreuze schaal (NSCH) en de Darcy-schaal. Het reduceert tot de klassieke Darcy- of Brinkman-vergelijkingen onder specifieke aannames, maar behoudt de thermodynamische consistentie van het faseveldmodel.
3. Enkelfluid-model: Het resultaat is een enkelfluid-model (single-fluid model) voor twee fasen, in plaats van een tweefluid-model. Dit elimineert de noodzaak voor empirische interface-randvoorwaarden en biedt een meer fundamentele beschrijving van de interfaciale dynamiek.
4. Inclusie van Inhomogeniteit: Het model houdt rekening met variaties in porositeit ($\nabla \epsilon$), wat essentieel is voor inhomogene media.

4. Resultaten en Validatie

De auteurs hebben het model gevalideerd via verschillende numerieke tests:

• Poiseuille-stroming in homogeen medium: De simulatie van een eenfasestroming in een kanaal met poreus medium toonde uitstekende overeenstemming met de analytische oplossing (Brinkman-vergelijking).
• Poiseuille-stroming in een vloeistof-poreus kanaal: Voor een kanaal dat deels gevuld is met poreus medium, toonde het model een goede overeenstemming met analytische oplossingen voor de snelheidsprofielen en de stress-sprong aan het interface.
• Buckley-Leverett probleem: Een 1D simulatie van verdringing van een niet-vochtende fase door een vochtende fase. De numerieke schokfronten kwamen goed overeen met de analytische oplossing (Welge-methode), hoewel er lichte numerieke oscillaties werden waargenomen nabij het front (een bekend kenmerk van centrale-differentie methoden zonder flux-limiters).
• Viskeuze vingers (Viscous Fingering): Simulaties van instabiliteit in poreuze media toonden aan dat het model het effect van vochtigheid (wettability) op de vingerpatronen kan vastleggen. Hydrofiele media (lage contacthoek) disperseren de vingers en remmen de vorming ervan af, wat consistent is met experimentele waarnemingen.

5. Betekenis en Conclusie

Dit werk biedt een fundamenteel theoretisch raamwerk voor het modelleren van meerfasestroming in poreuze media. De belangrijkste doorbraak is de rigoureuze afleiding van capillaire krachten en vochtgedrag vanuit de microscopische fysische wetten, in plaats van het gebruik van empirische correlaties.

• Het model lost het probleem op van het "koppelen" van poreuze details aan macroscopische vergelijkingen.
• Het biedt een nieuwe theoretische weg om vochtgedrag (wettability) en specifieke oppervlakte-oppervlakten in een verenigd formalisme op te nemen.
• Hoewel sommige sluitingsparameters (zoals de effectieve viscositeit en dispersiecoëfficiënten) in deze studie zijn vereenvoudigd, legt het de basis voor toekomstig onderzoek om deze parameters nauwkeuriger te bepalen voor complexe geometrieën en om de relatie met traditionele modellen (zoals de Leverett J-functie) verder te verduidelijken.

Kortom, deze studie levert een schaal-overbruggend model dat de thermodynamische consistentie van faseveldmodellen combineert met de praktische toepasbaarheid van Darcy-achtige wetten voor complexe meerfasestromingen.