A HHO formulation for variable density incompressible flows where the density is purely advected

Dit artikel presenteert een Hybrid High-Order (HHO)-formulering voor incompressibele stromingen met variabele dichtheid die een exacte volumebewaring garandeert, waardoor de dichtheid puur wordt geadverteerd, en die geschikt is voor het simuleren van mengsels van onmengbare vloeistoffen met eigenschappen zoals druk-robustheid en hoge orde nauwkeurigheid.

De kern

Stel je voor dat je een gigantische, digitale soep aan het koken bent. In deze soep zweven twee soorten vloeistoffen die niet met elkaar willen mengen, zoals olie en water. De ene is zwaar (zoals een steen), de andere is licht (zoals een kurk). Als je deze soep laat staan, zakt de zware vloeistof naar beneden en stijgt de lichte naar boven. Dit proces heet de Rayleigh-Taylor-instabiliteit. Het is een heel chaotisch dansje van vloeistoffen die in elkaar draaien, strengen vormen en breken.

Deze wetenschappers (Botti en Massa) hebben een nieuwe, superkrachtige manier bedacht om dit soort vloeistofbewegingen op de computer te simuleren. Ze noemen hun methode HHO (Hybrid High-Order). Hier is hoe het werkt, vertaald naar alledaags taal:

1. De "Puzzel" van de Vloeistof

Stel je voor dat je de soep in duizenden kleine blokjes (een rooster) verdeelt om te meten wat er in elk blokje gebeurt.

• Het oude probleem: Veel oude methoden waren als een slechte fotograaf: ze maakten de foto wel scherp, maar vergeten dat de olie en water niet echt met elkaar vermengd moesten worden, of ze vergeten dat de druk in de soep de beweging beïnvloedt. Soms leek het alsof de zware olie plotseling lichter werd dan water, wat fysiek onmogelijk is.
• De nieuwe oplossing (HHO): Deze methode is als een meester-puzzelaar. Ze kijken niet alleen naar het binnenste van elk blokje, maar ook naar de randen waar de blokjes elkaar raken. Door deze extra informatie te gebruiken, kunnen ze de vloeistofbeweging zo nauwkeurig berekenen dat de "olie" echt als olie blijft en de "water" als water.

2. De "Onzichtbare Hand" van de Druk

Een van de coolste dingen aan hun methode is hoe ze omgaan met druk.

• De analogie: Stel je voor dat je een ballon opblaast. De lucht (de vloeistof) wil zich verplaatsen, maar de druk duwt er tegen. In veel computermodellen is de druk een "verkeerde vriend": als je de druk niet perfect berekent, gaat de hele berekening van de stroming fout.
• De HHO-magie: Hun methode is "druk-onafhankelijk" (pressure-robust). Het is alsof je een auto bouwt die perfect rijdt, ongeacht of je de bandenspanning net iets te hoog of te laag hebt. De druk wordt behandeld als een losse, onzichtbare hand die alleen de vorm van de vloeistof bepaalt, maar niet de snelheid verstoort. Dit maakt de simulatie veel stabieler en sneller.

3. De "Dichte Soep" die niet verdwijnt

In hun simulatie wordt de dichtheid (hoe zwaar de vloeistof is) puur "meegedragen" door de stroming.

• De analogie: Stel je voor dat je een zakje met rode verf in een rivier gooit. De rivier stroomt en het zakje wordt meegenomen. Het zakje mag niet uit elkaar vallen, het mag niet verdampen en het mag niet van kleur veranderen.
• Het resultaat: Hun computermodel zorgt ervoor dat dit zakje perfect intact blijft, zelfs als de rivier heel snel stroomt of heel turbulent wordt. Ze garanderen dat het volume exact behouden blijft, tot op het kleinste detail van de computer.

4. De "Tijdmachine" (ESDIRK)

Om te zien hoe de soep zich in de tijd ontwikkelt, gebruiken ze een speciale tijdmethode genaamd ESDIRK.

• De analogie: Normaal gesproken kijken we naar een film en nemen we een foto elke seconde. Als de film te snel gaat, mis je details. Deze nieuwe methode is als een filmcamera die niet alleen naar de huidige seconde kijkt, maar ook een beetje naar de toekomst en het verleden om de perfecte tussenfoto te maken. Hierdoor kunnen ze heel grote stappen in de tijd maken zonder dat de simulatie "uit elkaar valt" of onnauwkeurig wordt.

5. Waarom is dit zo slim? (De "Slimme Verpakking")

De grootste uitdaging bij zulke complexe berekeningen is dat de computer vaak vastloopt omdat er te veel gegevens zijn.

• De oplossing: De auteurs gebruiken een truc genaamd "static condensation".
• De analogie: Stel je voor dat je een heel groot kantoorgebouw moet inspecteren. In plaats van elke kamer (elk blokje in de rooster) te betreden en alles te meten, kijken ze alleen naar de deuropeningen (de randen). Ze berekenen wat er binnenin gebeurt op basis van wat er op de deuropening gebeurt.
• Het voordeel: Dit maakt de berekening enorm efficiënt. Je kunt heel fijne details zien (hoge precisie) zonder dat je computer het geheugen verliest. Het is alsof je een HD-film kunt streamen op een oude laptop.

Wat hebben ze bewezen?

Ze hebben hun methode getest op twee scenario's:

1. Een rustige stroming: Hier bleek hun methode extreem nauwkeurig te zijn, zelfs met grove roosters.
2. De Rayleigh-Taylor instabiliteit (de zware en lichte vloeistof):
   • Bij een gematigd verschil in gewicht (lage Atwood-getal) konden ze met een grof rooster en hoge precisie dezelfde resultaten halen als andere methoden met een heel fijn rooster. Dit bespaart enorm veel rekentijd.
   • Bij een groot verschil in gewicht (hoge Atwood-getal, heel zware olie boven heel licht water) was hun methode de enige die stabiel bleef zonder dat de dichtheid negatief werd (wat fysiek onmogelijk is). Ze hielden de "olie" en "water" gescheiden, zelfs in het meest chaotische deel van de simulatie.

Kortom: Deze wetenschappers hebben een nieuwe, super-efficiënte en zeer nauwkeurige manier bedacht om te simuleren hoe vloeistoffen met verschillende gewichten zich mengen en bewegen. Het is alsof ze een nieuwe, onbreekbare wet hebben bedacht voor hoe computers vloeistoffen moeten "voelen", waardoor we in de toekomst veel realistischere voorspellingen kunnen doen over weer, branden of zelfs het ontwerp van nieuwe materialen.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het paper "A HHO formulation for variable density incompressible flows where the density is purely advected" van Botti en Massa, in het Nederlands.

Probleemstelling

Het paper richt zich op de numerieke simulatie van incompressibele stromingen met variabele dichtheid, specifiek voor mengsels van niet-mengbare vloeistoffen. De uitdaging bij dergelijke problemen ligt in het behouden van:

1. Exacte volumebewaring: De dichtheid moet puur worden vervoerd (geconveerd) zonder numerieke diffusie of bronnen/senken, wat essentieel is voor de fysieke correctheid van de interface tussen vloeistoffen.
2. Stabiliteit: Vooral in regimes met dominante convectie en bij hoge dichtheidsverhoudingen (zoals bij Rayleigh-Taylor-instabiliteit).
3. Druk-robustheid: De fouten in snelheid en dichtheid mogen niet afhankelijk zijn van de nauwkeurigheid van de drukoplossing.
4. Hoog-orde nauwkeurigheid: De behoefte aan methoden die zowel in de ruimte als in de tijd hoge-orde nauwkeurigheid bieden voor complexe stromingspatronen.

Bestaande methoden (zoals DG of klassieke Galerkin-methoden) hebben vaak moeite met het combineren van exacte incompressibiliteit, druk-robustheid en efficiënte implementatie op ongestructureerde roosters.

Methodologie

De auteurs stellen een nieuwe Hybrid High-Order (HHO) formulering voor, gekoppeld aan Explicit Singly Diagonal Implicit Runge-Kutta (ESDIRK) tijdsintegratie.

1. Ruimtelijke discretisatie (HHO):

• Hybride variabelen: De methode gebruikt polynomen op zowel de mesh-cellen als op de randen (skelet) van de cellen voor snelheid, druk en dichtheid.
• H(div)-conformiteit: Door een hybride drukruimte te gebruiken en lokale reconstructie-operatoren toe te passen, wordt een snelheidsveld gerealiseerd dat puntsgewijs divergentievrij is en H(div)-conform. Dit is cruciaal voor de exacte behoudswet van volume.
• Puur advectie van dichtheid: Omdat het snelheidsveld divergentievrij is, kan de massabehoudsvergelijking worden herschreven als de materiële tijdsafgeleide van de dichtheid ($\partial_t \rho + \mathbf{u} \cdot \nabla \rho = 0$). Dit garandeert dat de dichtheid puur wordt vervoerd.
• Stabilisatie:
  • Voor convectie worden schuif-symmetrische formuleringen gebruikt, aangevuld met upwind-stabilisatie voor massatransport en momentum.
  • Er is een specifieke stabilisatie voor de tijdsafgeleide van de dichtheid toegevoegd om robuustheid te garanderen wanneer de normale snelheidscomponent op een rand nul is.
• Randvoorwaarden: Zwakke oplegging van Dirichlet- en Neumann-randvoorwaarden.

2. Tijdsintegratie (ESDIRK):

• Het probleem wordt opgelost als een stelsel van differentiaal-algebraïsche vergelijkingen (DAE's) met index 2.
• Er worden ESDIRK-schema's gebruikt, die speciaal zijn ontworpen voor DAE's en implicititeit bieden voor stabiliteit bij stijve problemen.
• De methode is volledig impliciet en behoudt de hoge-orde nauwkeurigheid in de tijd.

3. Algebraïsche oplossing:

• Static Condensation: Cellen-variabelen worden lokaal geëlimineerd, zodat alleen de randvariabelen (face unknowns) globaal gekoppeld zijn. Dit vermindert de grootte van het globale stelsel en het geheugengebruik aanzienlijk.
• Preconditionering: Voor de lineaire oplosseerders binnen Newton-methode wordt een p-multigrid strategie gebruikt, wat de convergentie versnelt, vooral bij hoge polynoomgraden ($k$).

Belangrijkste Bijdragen

• Exacte volumebewaring: De methode garandeert cell-by-cell behoud van volume tot op machineprecisie, wat leidt tot een puur advectief gedrag van de dichtheid.
• Druk-robustheid: De formulering is bewezen druk-robust; irrotationale krachten beïnvloeden alleen het drukveld en niet de snelheids- of dichtheidsfouten.
• Stabiliteit: De methode is stabiel in convectie-gedomineerde regimes dankzij upwind-stabilisatie en is in staat om dichtheidsgrenzen te behouden bij lage polynoomgraden ($k=0$).
• Efficiëntie: Door static condensation en p-multigrid preconditionering is de methode schaalbaar en memory-efficiënt, zelfs voor hoge-orde benaderingen.
• Validatie: Uitgebreide numerieke validatie voor zowel ruimtelijke als temporele convergentie, inclusief complexe fysische scenario's.

Resultaten

De auteurs hebben de methode getest op drie verschillende scenario's:

1. Kovasznay-stroming (Steady State):

   • Toont optimale convergentiepercentages voor snelheid, druk en dichtheid.
   • De variant HHO-$\pi^{k+1}$ (waarbij de geconveerde snelheid niet wordt geprojecteerd) levert bij $k=0$ tweede-orde nauwkeurigheid op voor snelheid, terwijl de standaardvariant ($k$) slechts eerste-orde levert. Bij $k>0$ bereiken beide varianten de theoretische convergentiepercentages.
   • De divergentie van het snelheidsveld is tot op machineprecisie nul.

2. Guermond en Quartapelle Manufactured Solution (Unsteady):

   • Validatie van zowel ruimtelijke als temporele convergentie met ESDIRK-schema's (3-staps en 6-staps).
   • De resultaten bevestigen de druk-robustheid: snelheidsfouten zijn onafhankelijk van de drukfouten.
   • Temporele convergentiepercentages komen overeen met de theoretische orde van de gebruikte ESDIRK-schema's (2e orde voor TR-BDF2, 4e orde voor het 6-staps schema).

3. Rayleigh-Taylor Instabiliteit (RTI):

   • Laag Atwood-getal ($At=0.5$): Simulaties bij $Re=1000$ en $Re=5000$.
     • Vergelijking tussen hoge-orde (k=6) op een grof rooster en lage-orde (k=1) op een fijn rooster toont aan dat hoge-orde methoden zeer nauwkeurige resultaten leveren met aanzienlijk minder vrijheidsgraden (14x kleiner Jacobiaan).
     • Bij hoge $Re$ en hoge $k$ treden er spurious oscillaties op bij scherpe interfaces, wat wijst op de noodzaak van oscillatiecontrolestrategieën voor extreme regimes.
   • Hoog Atwood-getal ($At=0.75$): Simulatie bij $Re=1000$.
     • De methode is bounds-preserving (behoudt grenzen) bij $k=0$, wat betekent dat de dichtheid nooit negatief wordt, zelfs niet bij hoge dichtheidsverhoudingen (7:1).
     • Hoewel de ruimtelijke resolutie bij $k=0$ beperkt is voor zeer fijne details, wordt de hoofdstructuur van de plume correct gereproduceerd.

Significantie

Dit werk biedt een krachtig en robuust numeriek raamwerk voor de simulatie van niet-mengbare vloeistoffen met variabele dichtheid. De combinatie van HHO-methoden (voor ruimtelijke nauwkeurigheid en behoudseigenschappen) met ESDIRK-tijdsintegratie (voor stabiliteit en hoge-orde tijdsnauwkeurigheid) is uniek.

De belangrijkste implicaties zijn:

• Het vermogen om hoge dichtheidsverhoudingen en hoge Reynolds-getallen te simuleren zonder instabiliteit, dankzij de druk-robustheid en de mogelijkheid om dichtheidsgrenzen te behouden.
• De efficiëntie door static condensation en p-multigrid maakt hoge-orde simulaties haalbaar op complexe roosters.
• De methode vormt een solide basis voor toekomstige uitbreidingen naar meer complexe multi-fase stromingen, zoals gekoppeld met het Cahn-Hilliard-Navier-Stokes systeem.

Kortom, de paper presenteert een state-of-the-art methode die de beperkingen van bestaande discretisaties voor variabele dichtheid overwint en een nieuwe standaard zet voor nauwkeurige en efficiënte CFD-simulaties van mengsels.