Consistent multiple-relaxation-time lattice Boltzmann method for the volume averaged Navier-Stokes equations

In dit artikel wordt een consistent multiple-relaxation-time rooster-Boltzmann-methode ontwikkeld die de volume-gemiddelde Navier-Stokes-vergelijkingen met tweede-orde nauwkeurigheid herstelt door de holtefractie te ontkoppelen van de dichtheid en een straffingsbronterm in de momentruimte in te voeren om numerieke fouten te elimineren.

De kern

De "Slimme Simulatie": Hoe een nieuwe methode vloeistoffen door zand (en andere obstakels) beter laat stromen

Stel je voor dat je een enorme bak met water en duizenden kleine steentjes hebt. Je wilt precies weten hoe het water door de steentjes stroomt. Dit is lastig, want elke steen blokkeert het water op een andere manier. In de wereld van de computerwetenschappen noemen we dit het simuleren van vloeistof-vaste stof mengsels.

De auteurs van dit paper (Liu, Zhang, Min en Wu) hebben een nieuwe, slimme manier bedacht om dit te berekenen. Ze gebruiken een techniek die Lattice Boltzmann heet. Laten we dit uitleggen met een paar alledaagse metaforen.

1. Het Probleem: De "Grijze" Verwarring

Stel je voor dat je een digitale kaart maakt van een stad. Op sommige plekken zijn er huizen (de steentjes), en op andere plekken is er een park (het water).

• De oude methode: De computers die dit vroeger deden, gebruikten een simpele regel: "Als er een huis is, stroomt het water niet." Maar als je naar de rand van het huis kijkt, waar het water en de steen samenkomen, raakten ze in de war. De computer dacht: "Is dit nu water of steen?" Hierdoor ontstonden er spookstromen (in het Engels: spurious velocities). Het water leek vanzelf te gaan bewegen, zelfs als er niemand een kraan openhad. Alsof je een auto op een parkeerplaats ziet wegrijden zonder dat de motor aan staat.
• De druk-problematiek: In de oude methoden waren de druk en de dichtheid (hoe "vol" het water is) aan elkaar gekoppeld. Als je de dichtheid veranderde (bijvoorbeeld door meer steentjes toe te voegen), veranderde de druk onnatuurlijk. Dit zorgde voor rekenfouten, vooral als de overgang tussen water en steen heel scherp was.

2. De Oplossing: De "Meester-Regisseur" (MRT)

De auteurs hebben een nieuwe regisseur voor hun computerprogramma bedacht. Ze noemen dit een MRT-methode (Multiple-Relaxation-Time).

• De oude regisseur (SRT): Dit was een regisseur die alles met één handgreep deed. Als er iets misging, probeerde hij het met dezelfde simpele truc te fixeren. Dat werkte niet goed bij complexe situaties.
• De nieuwe regisseur (MRT): Deze regisseur heeft een team van specialisten. Hij kijkt naar de stroming en zegt: "Jij, specialist voor druk, doe dit. Jij, specialist voor snelheid, doe dat." Door de taken te verdelen, kan hij veel preciezer sturen.

3. De Twee Grote Verbeteringen

Verbetering A: De "Scheidingsmuur" (Decoupling)
In de oude methode waren de hoeveelheid steentjes (de "holte" of void fraction) en de waterdichtheid aan elkaar geplakt. Als je meer steentjes had, veranderde het water van gewicht.
De nieuwe methode bouwt een muur tussen deze twee. Ze zeggen: "De hoeveelheid steentjes bepaalt hoe veel ruimte er is, maar het gewicht van het water blijft gewoon water."

• Analogie: Stel je een bus voor. In de oude methode veranderde het gewicht van de bus als er meer mensen instapten, zelfs als de bus leeg was. In de nieuwe methode is de bus altijd even zwaar, ongeacht of er 1 of 100 mensen in zitten. Dit maakt de berekening veel rustiger en voorkomt die "spookstromen".

Verbetering B: De "Boete" (Penalty Source Term)
Soms maakt de computer kleine foutjes in de berekening van wrijving (viscositeit). De nieuwe methode heeft een slimme truc: een boete.
Als de computer een fout maakt in de wrijving, voegt de regisseur direct een kleine correctie toe (een "boete") om de fout te herstellen voordat het probleem groot wordt. Dit zorgt ervoor dat de stroming eruitziet zoals het in het echt zou moeten doen, zelfs als de overgang tussen water en steen heel scherp is (zoals een scherpe muur in een rivier).

4. Wat hebben ze bewezen?

De auteurs hebben hun nieuwe methode getest in drie situaties:

1. Gelijkmatige stroming: Water door een heel gelijkmatig zandbed. Dit werkt al goed bij oude methoden, maar hun methode deed het ook perfect.
2. Onregelmatige stroming: Water door een zandbed waar de steentjes willekeurig verspreid zijn. Hier faalden de oude methoden vaak door de "spookstromen", maar de nieuwe methode hield het water rustig en realistisch.
3. De "Schok" test: Een situatie waar er plotseling een muur van steentjes staat (van 0% water naar 100% steen). De oude methoden maakten hier enorme rekenfouten. De nieuwe methode gleed er soepel overheen, alsof het water een slimme zwemmer is die weet hoe hij om een obstakel moet zwemmen zonder te struikelen.

Conclusie

Kortom: Deze paper introduceert een slimmere, nauwkeurigere manier om te simuleren hoe vloeistoffen door materialen met gaten (zoals zand, schuim of bloedcellen) stromen.

• Vroeger: De computer maakte veel ruis en "spookbewegingen", vooral bij scherpe overgangen.
• Nu: Met deze nieuwe "MRT-regisseur" en de "scheidingsmuur" tussen druk en dichtheid, krijgen we een heel rustig, stabiel en accuraat beeld van de stroming.

Dit is belangrijk voor ingenieurs die bijvoorbeeld willen weten hoe medicijnen door het lichaam stromen, hoe olie uit de grond komt, of hoe batterijen het beste gemaakt kunnen worden. Het maakt de simulaties betrouwbaarder en sneller.

Technische samenvatting

Titel: Consistente multiple-relaxation-time lattice Boltzmann-methode voor de volume-gemiddelde Navier-Stokes-vergelijkingen

Auteurs: Yang Liu, Xuan Zhang, Jingchun Min, Xiaomin Wu
Affiliatie: Tsinghua Universiteit en Beijing Institute of Technology, China.

---

1. Het Probleem

De volume-gemiddelde Navier-Stokes-vergelijkingen (VANSE) vormen de hoeksteen voor het simuleren van vloeistof-vastestof meerfasensystemen (zoals gefluïdiseerde bedden, bloedstroom en poreuze media). Hoewel de Lattice Boltzmann-methode (LBM) populair is voor dergelijke simulaties vanwege zijn efficiëntie en geschiktheid voor complexe grensvlakken, kampen bestaande dichtheidsgebaseerde LBM-schema's voor VANSE met ernstige tekortkomingen:

• Spurious Velocities (Schijnbare snelheden): Bestaande methoden genereren ongewenste numerieke snelheden in gebieden met grote gradiënten in het holtefractie (porositeitsveld, $\phi$). Dit komt door numerieke fouten in de discretisatie van de correctiekracht voor de drukgradiënt.
• Inconsistentie: De herstelde macroscopische vergelijkingen zijn niet volledig consistent met de theoretische VANSE, vooral bij niet-uniforme holtefractievelden.
• Beperkingen van SRT: De recente drukgebaseerde LBM-methode (Fu et al.) loste sommige problemen op, maar de gebruikte Single-Relaxation-Time (SRT) operator beperkt de toepasbaarheid bij hoge viscositeiten en kan instabiliteit veroorzaken.
• Noodzaak: Er is een dringende behoefte aan een robuust, dichtheidsgebaseerd schema dat consistent is met VANSE, geschikt is voor grote spatiotemporele gradiënten in de holtefractie en vrij is van significante schijnbare snelheden.

2. Methodologie

De auteurs ontwikkelen een nieuw Multiple-Relaxation-Time (MRT) Lattice Boltzmann-schema (MRTLB-VANSE) binnen een dichtheidsgebaseerd kader. De kern van de methode bestaat uit drie innovatieve stappen:

1. Ontkoppeling van Holtefractie en Dichtheid:

   • In plaats van de gebruikelijke toestandsvergelijking $p = c_s^2 \phi \rho$ te gebruiken, wordt een voorlopige toestandsvergelijking geïntroduceerd: $p = c_s^2 \kappa \rho$, waarbij $\kappa$ een constante is (dicht bij de minimale holtefractie).
   • De evenwichtsverdeling (equilibrium distribution) wordt aangepast om deze voorlopige vergelijking te incorporeren. Hierdoor wordt de holtefractie $\phi$ ontkoppeld van de dichtheid $\rho$ in de drukgradiëntterm.
   • Dit elimineert de noodzaak om de gradiënt van de holtefractie ($\nabla \phi$) te discretiseren in de krachtterm, wat een belangrijke bron van numerieke ruis en onfysische oscillaties bij discontinuïteiten wegneemt.

2. Penalty Source Term in Momentruimte:

   • De introductie van de voorlopige toestandsvergelijking veroorzaakt een bias in de derde-orde momenten, wat leidt tot fouten in de viskeuze spanningsensor en een schending van de Galileïse invariantie.
   • Om dit op te lossen, wordt een strafbronterm (penalty source term) toegevoegd in de momentruimte van de MRT-collisieoperator. Deze term corrigeert specifiek de fouten in de viskeuze spanningsensor en garandeert dat de herstelde macroscopische impulsvergelijking Galileïsch invariant is.

3. MRT-Collisieoperator:

   • Het gebruik van de MRT-operator (in plaats van SRT) biedt betere numerieke stabiliteit, vooral bij lage relaxatietijden (hoge viscositeit), en maakt de implementatie van de bovengenoemde bronterm mogelijk.

3. Belangrijkste Bijdragen

• Theoretische consistentie: Het artikel biedt een volledige Chapman-Enskog-analyse die bewijst dat het nieuwe schema de VANSE exact herleidt met tweede-orde nauwkeurigheid, inclusief de correcte drukgradiënt en viskeuze termen.
• Eliminatie van Schijnbare Snelheden: Door de ontkoppeling van $\phi$ en $\rho$ in de evenwichtsverdeling, blijven de deeltjesverdelingen isotroop tijdens het stromingsproces, zelfs in niet-uniforme velden. Dit onderdrukt schijnbare snelheden aanzienlijk, zelfs bij discontinuïteiten in de holtefractie.
• Verbeterde Stabiliteit: Het schema overwint de stabiliteitsproblemen van SRT-schema's bij hoge viscositeiten en werkt betrouwbaar voor holtefracties kleiner dan 0,5, waar eerdere methoden vaak divergeerden.

4. Resultaten en Validatie

De methode is uitgebreid gevalideerd via numerieke tests en vergelijkingen met analytische oplossingen:

• Uniforme Poreuze Stroming: De methode reproduceert nauwkeurig de analytische oplossingen voor Darcy-Brinkman stroming (Poiseuille en Couette stroming) in homogene media.
• Niet-uniforme Deeltjesstroming: In een 1D-test met ruimtelijke en tijdsafhankelijke fluctuaties in de holtefractie, toonde de methode perfecte overeenstemming met de analytische oplossing, wat de nauwkeurigheid van de drukcorrectie bevestigt.
• Discontinue Velden (No-Flow Test): In een scenario met een plotselinge discontinuïteit in de holtefractie (zonder externe stroom), genereerde het MRTLB-VANSE-schema extreem lage schijnbare snelheden. Dit was een groot verbetering ten opzichte van eerdere dichtheidsgebaseerde methoden (Zhang et al., Bukreev et al., Blais et al.), die aanzienlijke onfysische snelheden vertoonden.
• Convergentie (MMS): Met de "Method of Manufactured Solutions" (MMS) werd aangetoond dat de methode tweede-orde convergentie behaalt voor zowel divergentievrije als niet-divergentievrije stromingen, ongeacht de grootte van de holtefractiegradiënt of de viscositeit.
• Vergelijking met SRT: In vergelijking met de recente drukgebaseerde SRT-methode (Fu et al.), presteerde het MRTLB-VANSE-schema beter bij hoge viscositeiten, waar de SRT-methode zijn nauwkeurigheid verloor.

5. Significantie en Toekomstperspectief

Dit werk is significant omdat het een robuust, consistent en nauwkeurig numeriek gereedschap biedt voor het simuleren van complexe vloeistof-vastestof interacties.

• Toepassingsgebied: De methode is ideaal voor industriële toepassingen zoals mineralenoplossing, koksverbranding, matrixdeformatie en capillaire/viskeuze fractuur in poreuze media.
• Schaalbaarheid: Omdat de LBM van nature paralleliseerbaar is, is deze methode zeer geschikt voor grootschalige simulaties.
• Toekomst: De auteurs wijzen erop dat het uitbreiden naar 3D eenvoudig is en dat het koppelen aan concentratie- en temperatuurvelden, evenals onmengbare meerfasemodellen (kleurgradiënt, pseudo-potentiaal), veelbelovende paden opent voor het simuleren van complexe vloeistof-gas-vastestof processen.

Conclusie: De voorgestelde MRTLB-VANSE-methode lost de inherente inconsistenties en numerieke instabiliteiten van eerdere dichtheidsgebaseerde benaderingen op, en biedt een betrouwbare standaard voor toekomstige coarse-grained simulaties van meerfasensystemen.