Computational homogenization of unsteady flows in a periodic porous medium

1. Probleemstelling

Het artikel behandelt de uitdaging van het modelleren van onstabiliseerde filtratie van viskeuze, onsamendrukbare vloeistoffen in periodieke poreuze media.

• Microschaal: De stroming wordt beheerst door de onstabiliseerde Stokes-vergelijkingen. Directe numerieke simulatie van individuele poriën is computatief onhaalbaar voor macroscopische domeinen.
• Macroschaal: De standaard Darcy-wet is ontoereikend voor onstabiliseerde stromingen waarbij viskeuze diffusietijdschalen vergelijkbaar zijn met observatietijden. In deze gevallen vertoont de stroming geheugeneffecten, wat een integro-differentiële Darcy-wet vereist.
• De Kernuitdaging: De macroscopische filtratiesnelheid hangt af van de geschiedenis van de drukgradiënt via een tensor-geheugenkern $K(t)$. Het oplossen van de resulterende niet-lokale integro-differentiële vergelijking is computatief duur omdat het het opslaan van de oplossingengeschiedenis op alle vorige tijdstappen vereist, wat leidt tot hoge geheugen- en tijdscomplexiteit ($O(N^2)$).

2. Methodologie

De auteur stelt een computationele homogenisatietechnologie voor die het niet-lokale macroscopische probleem transformeert in een lokaal systeem van differentiaalvergelijkingen. De aanpak bestaat uit drie hoofdfasen:

A. Microschaal Karakterisering (Celproblemen)

In plaats van onstabiliseerde Stokes-problemen direct op de cel op te lossen voor elke macroscopische tijdstap, maakt de methode gebruik van een spectrale aanpak:

1. Stabiele Probleem: Oplossing van hulpelijke randwaardeproblemen om de standaard (statische) doorlatendheidstensor $\bar{K}$ te bepalen.
2. Spectraal Probleem: Oplossing van een partiële eigenwaardeprobleem (Stokes-spectraal probleem) op de periodiciteitscel $Y_f$:
   $$-\Delta_y \phi + \nabla_y \eta = \lambda \phi, \quad \nabla_y \cdot \phi = 0$$
   Dit levert eigenwaarden $\lambda_k$ en eigenfuncties $\phi_k$ op.
3. Kernrepresentatie: De componenten van de dynamische doorlatendheidstensor $K_{ij}(t)$ worden uitgedrukt als een oneindige som van exponentiële functies afgeleid van de spectrale data:
   $$K_{ij}(t) = \sum_{k=1}^{\infty} a_i^k a_j^k \exp(-\lambda_k t)$$
   waarbij $a_i^k$ coëfficiënten zijn die zijn afgeleid van de eigenfuncties.

B. Benadering en Lokalisatie

Om het probleem computatief hanteerbaar te maken:

• Truncatie: De oneindige som wordt afgekapt tot de eerste $m$ termen.
• Delta-Benadering: De staart van de reeks (termen $k > m$) wordt benaderd als een Dirac-deltafunctie $\delta(t)$, die bijdraagt aan de statische doorlatendheidstensor.
• Lokale Transformatie: Het niet-lokale convolutie-integraal wordt vervangen door een systeem van $m$ hulpelijke gewone differentiaalvergelijkingen (ODE's). Nieuwe variabelen $c_k(x,t)$ worden geïntroduceerd:
  $$\frac{\partial c_k}{\partial t} + \lambda_k c_k - p = 0$$
  De macroscopische vergelijking wordt een gekoppeld systeem van de drukvergelijking en deze $m$ ODE's, waardoor de geheugeneffecten effectief worden gelokaliseerd.

C. Numerieke Implementatie

• Ruimtelijke Discretisatie: Finite Element Method (FEM) met gemengde elementen (Hood-Taylor: kwadratische snelheid, lineaire druk) voor de celproblemen en lineaire elementen voor het macroscopische domein.
• Tijdsdiscretisatie: Twee-niveau impliciete differentieschema's met een gewichtsparameter $\sigma$.
  • $\sigma = 1/2$: Tweede-orde nauwkeurigheid.
  • $\sigma \geq 1/2$: Ondinget voorwaarde gestabiel.

3. Belangrijkste Bijdragen

1. Spectraal Gebaseerde Kernbenadering: Het artikel introduceert een methode om de componenten van de geheugenkern direct te berekenen uit de eigenwaarden en eigenfuncties van het Stokes-spectrale probleem op de periodiciteitscel, waardoor de noodzaak om onstabiliseerde celproblemen op te lossen bij elke macroscopische tijdstap wordt vermeden.
2. Niet-lokaal-naar-lokaal Transformatie: Het slaagt erin de integro-differentiële Darcy-wet te transformeren in een lokaal systeem van partiële differentiaalvergelijkingen (PDE's) gekoppeld aan ODE's. Dit reduceert de computatieve complexiteit van kwadratisch (in tijdstappen) naar lineair, aangezien de oplossing op stap $n+1$ alleen afhankelijk is van stap $n$.
3. Stabiliteitsanalyse: De auteur levert rigoureuze a priori-schattingen die de onvoorwaardelijke stabiliteit van de voorgestelde twee-niveau differentieschema's voor $\sigma \geq 1/2$ bewijzen, waardoor wordt gegarandeerd dat de numerieke oplossing begrensd blijft en consistent is met het continue probleem.
4. Efficiënte Filterstrategie: De studie toont aan dat de geheugenkern nauwkeurig kan worden benaderd met een zeer klein aantal eigenwaarden (bijvoorbeeld $m=3$ tot $m=10$), omdat hogere-orde eigenwaarden snel afnemen, wat aanzienlijke computatiewinst mogelijk maakt zonder in te leveren op nauwkeurigheid.

4. Numerieke Resultaten

De methodologie werd getest op een 2D-benchmarkprobleem met een periodiek poreus medium dat elliptische obstakels bevat.

• Celprobleem: De doorlatendheidstensor werd berekend voor verschillende aspectverhoudingen van de ellipsen ($\gamma$). Resultaten toonden aan dat naarmate de ellips verlengt, de niet-diagonale doorlatendheidscomponenten toenemen, wat anisotropie aangeeft.
• Eigenwaarde Convergentie: Eigenwaarden berekend op verfijnde netten convergeerden snel. De eerste 100 eigenwaarden toonden ongeveer lineaire groei.
• Invloed van Geheugeneffecten:
  • Simulaties toonden een snelle "heraanpassing" van het drukprofiel op het initiële tijdstip als gevolg van geheugeneffecten.
  • Het systeem convergeert snel naar de stabiele oplossing, die overeenkomt met de standaard Darcy-oplossing (zonder geheugen) zodra de overgangsfase voorbij is.
  • Efficiëntie: Het gebruik van slechts $m=3$ eigenwaarden was voldoende om de essentiële geheugeneffecten voor de testcase te vangen, met een verwaarloosbare fout vergeleken met het gebruik van $m=100$.
• Computatiekosten: De complexiteit van het voorgestelde algoritme is slechts iets hoger dan die van een standaard lokaal filtratieprobleem, aangezien het het oplossen van $m$ eenvoudige hulpelijke lokale evolutionaire problemen expliciet vereist.

5. Betekenis

Dit werk biedt een robuust en efficiënt raamwerk voor het simuleren van onstabiliseerde stromingen in poreuze media waarbij inertie- en geheugeneffecten significant zijn (bijvoorbeeld snelle drukveranderingen, hoge-frequentie oscillaties).

• Praktijktoepasbaarheid: Door een niet-lokaal probleem om te zetten in een lokaal probleem, maakt de methode grootschalige 3D-simulaties van onstabiliseerde filtratie haalbaar op standaard high-performance computersystemen.
• Theoretische Strenge: De combinatie van homogenisatietheorie, spectrale analyse en stabiliteitsbewijzen biedt een wiskundig onderbouwde alternatief voor empirische geheugenmodellen.
• Veelzijdigheid: De aanpak kan worden uitgebreid naar andere soorten poreuze media en stromingsregimes waar geheugenkernen exponentiële afname vertonen, mits de spectrale eigenschappen van de onderliggende microstructuur kunnen worden berekend.