Linear Stability Analysis of Two-phase, Two-Component Flow in Porous Media

Technische Samenvatting: Lineaire Stabiliteitsanalyse van Tweefasige, Tweekomponentige Stroming in Poreuze Media

Probleemstelling
Viskeuze vingervorming-instabiliteiten tijdens de verplaatsing van vloeistoffen in poreuze media compromitteren aanzienlijk de efficiëntie van toepassingen zoals enhanced oil recovery, CO2-sequestratie en grondwaterremediatie. Hoewel er uitgebreide literatuur bestaat over lineaire stabiliteitsanalyse voor volledig immiscibele en volledig miscibele verplaatsingen, blijft het tussenliggende geval van gedeeltelijk miscibele stroming — waarbij beperkte massaoverdracht tussen fasen plaatsvindt — grotendeels onverkend. Specifiek ontbreekt een analytisch kader dat rekening houdt met de complexe interactie tussen zwaartekracht, capillaire krachten, mechanische dispersie en interfasiale massaoverdracht in een tweefasig, tweekomponentig systeem. Deze studie richt zich op de voorspelling van onstabiele perturbatiemodi in dergelijke gedeeltelijk miscibele systemen, met de focus op het scenario waarbij een gasvormige vloeistof een vloeibare vloeistof verplaatst.

Methodologie
De auteurs ontwikkelen een wiskundig model voor tweedimensionale stroming in een isotroop poreus medium, waarbij zij uitgaan van thermodynamisch evenwicht voor een binair mengsel (componenten "a" en "b"). De methodologie verloopt via de volgende stappen:

1. Goverende Vergelijkingen: Het model leidt massa-balansvergelijkingen af voor een tweefasig, tweekomponentig systeem, waarbij de wet van Darcy, fractionele stroming, capillaire druk en mechanische dispersie worden geïntegreerd. Een belangrijke vereenvoudigende aanname is dat de dichtheid van de vloeibare fase constant is in de eenfasige regio, terwijl de viscositeit varieert met de samenstelling.
2. Basisstaat Formulering: De basisstaat wordt gedefinieerd als een eendimensionale, bewegende golfoplossing (schokfront) afgeleid van de Buckley-Leverett-benadering, aangepast voor massaoverdracht. De schokconfiguratie wordt berekend met de methode van karakteristieken, waarbij rekening wordt gehouden met de discontinuïteit in de stromingsfuncties bij de overgang van de tweefasige regio naar de zuivere vloeistofregio.
3. Lineaire Stabiliteitsanalyse: Kleine-amplitude, golfachtige perturbaties worden op de basisstaat gesuperponeerd. De goverende vergelijkingen worden gelineariseerd om een eigenwaardeprobleem te vormen waarbij de eigenwaarde de groeirate van de perturbatie ($\sigma$) vertegenwoordigt en de eigenfuncties de verstoringsprofielen weergeven.
4. Numerieke Oplossing: Een kritische uitdaging die wordt geïdentificeerd, is de discontinuïteit in de afgeleiden van de eigenfuncties bij het overgangspunt tussen tweefasige en eenfasige stroming. Om dit op te lossen, leiden de auteurs specifieke sprongcondities (jump conditions) af voor de afgeleiden van de eigenfuncties bij dit interface. Het eigenwaardeprobleem wordt numeriek opgelost met behulp van de Matched Initial Value Problem (MIVP) methode. Dit houdt in dat de lineaire gewone differentiaalvergelijkingen vanuit de ver-veld grenzen naar het overgangspunt worden geïntegreerd, waarbij de oplossingen worden samengevoegd met behulp van de afgeleide sprongcondities.

Belangrijkste Bijdragen

• Nieuw Kader: Dit werk presenteert, naar weten van de auteurs, de eerste lineaire stabiliteitsanalyse van een echt tweefasig, tweekomponentig systeem dat expliciet rekening houdt met interfasiale massaoverdracht.
• Sprongcondities: De afleiding van sprongcondities voor de afgeleiden van eigenfuncties bij de faseovergangsdiscontinuïteit maakt een nauwkeurige oplossing van het eigenwaardeprobleem mogelijk in systemen waar stromingsfuncties discontinu zijn.
• Uitgebreide Parameterstudie: De studie onderzoekt systematisch de effecten van initiële samenstelling, dispersiviteiten (longitudinaal en transversaal), miscibiliteitsniveaus (via binaire interactiecoëfficiënten), viscosieratio's en zwaartekracht op de stabiliteit van het verplaatsingsfront.

Resultaten

• Stabiliserend Effect van Massaoverdracht: De belangrijkste bevinding is dat massaoverdracht een overwegend stabiliserend effect heeft op de verplaatsing. Het vermindert de maximale groeirate ($\sigma_{max}$) door de evenwichtsviscositeitsratio ($M_e$) te verlagen en de schokeigenschappen te wijzigen (het verlagen van de schoksnelheid $v_s$ en het verhogen van de gasverzadiging $S_{gs}$). Deze stabilisatie is bijzonder uitgesproken bij hoge viscositeitscontrasten.
• Interacties met Zwaartekracht: Bij opwaartse verplaatsingen dempt massaoverdracht de door zwaartekracht geïnduceerde instabiliteit. Het dichtheidscontrast dat de gravitationele instabiliteit aandrijft, versterkt ook de massaoverdracht, wat het destabiliserende effect van de zwaartekracht tegengaat. Bij neerwaartse verplaatsingen werken zowel zwaartekracht als massaoverdracht als stabiliserende mechanismen.
• Niet-lineariteit van Dispersie: De studie identificeert een niet-lineaire relatie tussen de dimensieloze longitudinale dispersiecoëfficiënt ($D^*_{\ell,xx}$) en instabiliteit. Er bestaat een kritische "meest gevaarlijke" waarde voor $D^*_{\ell,xx}$ waarbij zowel de groeirate als de afkapgetal (cutoff wavenumber, $n_{cut}$) worden gemaximaliseerd. Dit suggereert een complexe interactie waarbij intermediaire niveaus van mechanische dispersie, relatief ten opzichte van capillaire krachten, de instabiliteit kunnen vergroten vergeleken met zowel zeer lage als zeer hoge dispersieregimes.
• Afkapgetal (Cutoff Wavenumber): Hoewel massaoverdracht de groeirate aanzienlijk vermindert, is het effect op het afkapgetal minder uitgesproken en complexer, waarbij in specifieke gevallen soms hogere $n_{cut}$ waarden worden bereikt voor gedeeltelijk miscibele gevallen vergeleken met immiscibele gevallen.

Betekenis en Claims
Het artikel claimt de kloof te overbruggen tussen de goed bestudeerde limietgevallen van immiscibele en miscibele verplaatsing. Door aan te tonen dat massaoverdracht een stabiliserend mechanisme is, suggereren de resultaten dat gedeeltelijk miscibele condities robuustere verplaatsingsfronten kunnen bieden dan puur immiscibele scenario's, met name in omgevingen met een hoog viscositeitscontrast. De identificatie van de kritische dispersiecoëfficiënt waarde benadrukt een voorheen ongerapporteerde complexe interactie tussen capillaire krachten en mechanische dispersie in de context van viskeuze vingervorming. De auteurs merken op dat hoewel hun model de fysica van multicomponent-systemen vereenvoudigt, het een noodzakelijk kader biedt voor het begrijpen van de fundamentele stabiliteitsmechanismen in gedeeltelijk miscibele stromingen. De studie erkent beperkingen, zoals de aanname van een constante vloeistofdichtheid in de eenfasige regio en vaste relatieve permeabiliteitscurven, en stelt deze voor als gebieden voor toekomstige verfijning.