A Shakhov-based Bhatnagar-Gross-Krook model for polyatomic molecules and for atomic as well as polyatomic mixtures

Dit artikel beschrijft de uitbreiding van het Shakhov-BGK-model in de open-source PICLas-code voor polyatomaire moleculen en mengsels, waarbij de nauwkeurigheid van de nieuwe aanpak voor niet-evenwichtstoestanden en schokgolven wordt geverifieerd door vergelijking met DSMC- en ESBGK-metodes.

De kern

De "Verkeersregelaar" voor Moleculen: Een Nieuwe Manier om Ruimtevaart te Simuleren

Stel je voor dat je een enorme stad moet simuleren waar miljoenen auto's rondrijden. In de stad (zoals in de atmosfeer dicht bij de aarde) rijden ze in dichte files. Maar naarmate je hoger komt (zoals in de ruimte), wordt het steeds stiller en rijden de auto's ver uit elkaar.

Wetenschappers gebruiken computers om te voorspellen hoe deze auto's (moleculen) zich gedragen. Dit is cruciaal voor het ontwerpen van raketten en satellieten. Maar hier zit een probleem:

• De oude methode (DSMC): Dit is alsof je elke auto individueel volgt. In de dichte stad werkt dit perfect, maar het kost enorm veel rekenkracht. In de ruimte is het ook lastig, omdat je dan duizenden lege plekken moet simuleren.
• De snellere methode (CFD): Dit is alsof je alleen naar het totale verkeer kijkt (gemiddelde snelheid, druk). Dit is supersnel, maar faalt als de auto's te ver uit elkaar staan of als ze chaotisch gedrag vertonen (zoals in een schokgolf).

De auteurs van dit artikel hebben een nieuwe, slimme "hybride" methode ontwikkeld die de beste van beide werelden combineert. Ze noemen dit het Shakhov-BGK-model.

1. Het Probleem: De "Moeilijke" Moleculen

In de ruimte zijn er twee soorten "auto's":

1. Atomen: Simpele bolletjes (zoals Helium of Argon). Ze bewegen alleen maar heen en weer.
2. Moleculen: Complexere constructies (zoals Stikstof $N_2$ of Koolstofdioxide $CO_2$). Deze kunnen niet alleen bewegen, maar ook draaien en trillen (zoals een springveer).

De oude simpele modellen wisten niet goed hoe ze met die draai- en trillende energie moesten omgaan. Ze gaven vaak een verkeerd antwoord over hoe snel de moleculen opwarmen of afkoelen.

2. De Oplossing: Een Slimme "Verkeersregelaar"

De auteurs hebben hun model (dat in de open-source software PICLas zit) verbeterd. Ze hebben een nieuwe "verkeersregelaar" bedacht die:

• Zowel simpele atomen als complexe moleculen begrijpt.
• Zelfs mengsels kan simuleren (bijvoorbeeld een mix van Stikstof en Koolstofdioxide, zoals op Mars).
• De "interne energie" (het draaien en trillen) correct meerekent.

De Metafoor van de Dansvloer:
Stel je een dansvloer voor.

• Bij de oude methode (Standard BGK) moesten de dansers (moleculen) altijd precies in een rechte lijn dansen. Als ze moesten draaien of trillen, werd de dans onnatuurlijk.
• Bij de nieuwe Shakhov-methode krijgen de dansers een speciale instructie: "Als de temperatuur stijgt, draai je net iets anders dan normaal." Hierdoor blijft de dans (de fysica) realistisch, zelfs als de dansers heel snel bewegen of als er een plotselinge drukgolf komt.

3. Waarom is dit beter dan de vorige versie?

Er was al een vergelijkbare methode genaamd ESBGK. Die was al best goed, maar had een zwak punt: het zag schokgolven (zoals de luchtstroom voor een raket die met supersnelheid vliegt) vaak te "wazig". Het was alsof je door een wazige bril keek; je zag dat er iets aan de hand was, maar de randen waren niet scherp.

De nieuwe Shakhov-methode werkt als een scherp objectief.

• In hun tests lieten ze een raketkop (een afgeronde kegel) door de lucht vliegen.
• De nieuwe methode zag de schokgolf (de muur van samengeperste lucht voor de raket) veel scherper en nauwkeuriger dan de oude methode.
• Het resultaat kwam bijna exact overeen met de "gouden standaard" (DSMC), maar dan veel sneller berekend.

4. De Praktijk: Wat hebben ze getest?

Ze hebben drie soorten tests gedaan om hun nieuwe regelaar te bewijzen:

1. De Supersonische Couette-stroom: Twee wanden die langs elkaar schuiven met gas ertussen. Ze keken of de temperatuur correct werd voorspeld bij verschillende mengsels (alleen stikstof, een mix van helium en argon, en een mix van stikstof-atomen en -moleculen). Resultaat: Het werkt perfect, zelfs bij zware mengsels.
2. De Hypersonische Raketkop: Een raketkop die met extreme snelheid door de atmosfeer vliegt.
   • Test 1: Alleen Stikstof ($N_2$).
   • Test 2: Een mix van $CO_2$ en $N_2$ (zoals de atmosfeer van Mars).
   • Test 3: Een mix van drie soorten ($N_2$, $O_2$ en $N$).
   • Resultaat: In alle gevallen zag de nieuwe methode de schokgolf en de hitteverdeling veel nauwkeuriger dan de oude methode. Vooral bij de complexe $CO_2$-mix (waar de moleculen veel kunnen trillen) was het verschil groot.

Conclusie in één zin

De auteurs hebben een slimmer, sneller en nauwkeuriger rekenmodel gemaakt dat de complexe dans van moleculen in de ruimte beter nabootst dan ooit tevoren. Hierdoor kunnen ingenieurs in de toekomst sneller en goedkoper ontwerpen voor satellieten en raketten, zonder dat ze duizenden jaren aan rekenkracht nodig hebben.

Het is alsof ze van een handgetekende schets zijn gegaan naar een 3D-simulatie met een hoge resolutie, die toch net zo snel draait als de oude schets.

Technische samenvatting

Titel: Shakhov-BGK voor polyatomische mengsels

Auteurs: M. Pfeiffer en F. Tuttas (Instituut voor RuimteSystemen, Universiteit van Stuttgart)
Context: Uitbreiding van de open-source deeltjescode PICLas voor het modelleren van niet-evenwichtstoestanden in de interne vrijheidsgraden van polyatomische moleculen en mengsels.

---

1. Het Probleem

Numerieke simulaties van gasdynamica in ruimtevaarttoepassingen, micro- en nanoschaalstromingen en vacuümtechnologie stellen grote eisen. Deze toepassingen omvatten vaak grote dichtheidsgradiënten die variëren van het continue regime tot het vrije moleculaire stromingsregime.

• DSMC (Direct Simulation Monte Carlo): Hoewel DSMC zeer nauwkeurig is voor niet-evenwichtsstromingen, wordt de rekenkosten aanzienlijk hoger in de overgangs- en continue regimes.
• CFD (Computational Fluid Dynamics): Traditionele CFD-methoden zijn efficiënt in het continue regime, maar het koppelen ervan aan DSMC is problematisch vanwege fundamenteel verschillende benaderingen en statistische ruis aan de grenzen.
• Bestaande BGK-modellen: De Bhatnagar-Gross-Krook (BGK) methode benadert de botsingsintegraal via een relaxatieproces, wat minder rekenkracht vereist dan DSMC. Echter, standaard BGK-modellen hebben vaak een onjuiste Prandtl-getal en kunnen de warmteflux niet correct modelleren. Bestaande uitbreidingen voor mengsels (zoals ESBGK) zijn vaak complex (meerdere relaxatietermen) of minder nauwkeurig in het vastleggen van schokgolven.

2. Methodologie

De auteurs hebben de Shakhov-BGK (SBGK) methode geïmplementeerd en uitgebreid in de open-source deeltjescode PICLas. De kern van de methode is als volgt:

• Enkele Relaxatieterm: In plaats van complexe modellen met meerdere relaxatietermen voor elke soortcombinatie, gebruikt dit model een enkele relaxatieterm per soort. Dit vereenvoudigt de berekening aanzienlijk en verlaagt de rekenkosten, terwijl het toch de juiste macroscopische eigenschappen behoudt.
• Doelverdeling (Target Distribution):
  • Voor atomen wordt de doelfunctie puur translationeel.
  • Voor moleculen (diatomisch en polyatomisch) wordt de doelfunctie gefactoreerd in translationele, rotatie- en vibratiecomponenten. Vibratie wordt gemodelleerd als een gekwantiseerd harmonisch oscillator (SHO) met vibratiekwantumgetallen.
• Correctie voor het Prandtl-getal: Het Shakhov-model past de doelfunctie aan om het correcte Prandtl-getal ($Pr$) te verkrijgen. Dit wordt gedaan door een correctieterm voor de warmteflux in te bouwen. De parameter $\alpha$ in de formule wordt dynamisch berekend op basis van de massafractie, dichtheidsfractie en interne vrijheidsgraden van het mengsel.
• Transporteigenschappen: Om de viscositeit en warmtegeleidingsvermogen van mengsels nauwkeurig te bepalen, wordt gebruikgemaakt van de eerste benadering gebaseerd op botsingsintegralen (collision integrals) volgens de Chapman-Enskog-theorie, in plaats van de simpelere Wilke-mengregel. Dit is essentieel voor mengsels met grote massaverschillen.
• Relaxatiefrequentie: De relaxatiefrequentie $\nu$ wordt bepaald door de mengselviscositeit en de translationele temperatuur. De interne energie-relaxatie volgt het Landau-Teller-formulier om de overgang tussen translationele, rotatie- en vibratie-temperaturen correct te modelleren.

3. Belangrijkste Bijdragen

1. Uitbreiding naar Polyatomische Mengsels: Voor het eerst is een SBGK-model met één relaxatieterm succesvol geïmplementeerd voor mengsels die atomen en polyatomische moleculen (met rotatie en vibratie) bevatten.
2. Behoud van Eigenschappen: Het model garandeert behoud van massa, impuls en energie. Het reproduceert de Maxwelliaanse verdeling in thermodynamisch evenwicht.
3. Efficiëntie: Door gebruik te maken van één relaxatieterm en alleen de momenten van het totale mengsel te hoeven evalueren (in plaats van per soort), is het model computatieel veel efficiënter dan multi-relaxatiemodels.
4. Nauwkeurige Schokoplossing: Het model is specifiek ontworpen om schokgolven nauwkeuriger te vangen dan het ESBGK-model, wat cruciaal is voor hypersonische stromingen.

4. Resultaten en Validatie

Het model is gevalideerd via twee hoofdtestcases en vergeleken met DSMC (de "gouden standaard") en het ESBGK-model:

• Case 1: Supersonische Couette-stroming:

  • Getest met zuivere $N_2$, een Ar-He mengsel, en een $N_2$-$N$ mengsel.
  • De resultaten voor temperatuurprofielen toonden uitstekende overeenkomst met DSMC.
  • Bij het Ar-He mengsel (groot massaverschil) was er een lichte afwijking, wat consistent is met eerdere bevindingen bij vergelijkbare modellen. De $N_2$-$N$ case toonde bijna perfecte overeenkomst.

• Case 2: Hypersonische stroming rond een afgeronde kegel (70°):

  • Getest met zuivere $N_2$, een $CO_2$-$N_2$ mengsel (veel interne vrijheidsgraden), en een $N_2$-$O_2$-$N$ mengsel.
  • Schokgolfstructuur: Het SBGK-model ving de start van de temperatuurstijging in de schokgolf aanzienlijk nauwkeuriger op dan het ESBGK-model. Het ESBGK-model neigde tot een te vroeg begin van de stijging, wat resulteerde in een bredere (en minder accurate) schokgolf.
  • Interne Vrijheidsgraden: De modellering van rotatie- en vibratietemperaturen (niet-evenwicht) kwam zeer goed overeen met DSMC, zelfs in de complexe schokregio.
  • Warmteflux: De berekende warmteflux op het oppervlak van de kegel stemde uitstekend overeen met zowel DSMC als experimentele metingen.

5. Betekenis en Conclusie

Dit werk biedt een krachtig en efficiënt alternatief voor DSMC in regimes waar de Knudsen-getallen laag tot gemiddeld zijn (overgangs- en continue regimes), maar waar niet-evenwichtseffecten (zoals in schokgolven) nog steeds belangrijk zijn.

• Kernvoordeel: Het SBGK-model combineert de rekenefficiëntie van BGK-methoden met de nauwkeurigheid van DSMK voor schokgolven en interne energie-relaxatie.
• Toekomstperspectief: De auteurs zien een koppeling met DSMC voor multi-schaalproblemen als een logische volgende stap. Ook is de implementatie van chemische reacties in dit model een geplande uitbreiding.

Samenvattend bewijst dit artikel dat het Shakhov-BGK-model met één relaxatieterm een robuuste, nauwkeurige en computatie-efficiënte methode is voor het simuleren van complexe polyatomische gasmengsels in ruimtevaarttoepassingen.