Boltzmann-Curtiss Description for Flows under Translational Nonequilibrium

Dit onderzoek toont aan dat de Boltzmann-Curtiss-benadering, die zowel rotatie- als translatievrijheidsgraden in acht neemt, in vergelijking met de Navier-Stokes-vergelijkingen aanzienlijk betere resultaten oplevert voor de beschrijving van schokgolven onder translatie- en rotatie-ongelijkwicht.

De kern

Hier is een uitleg van het onderzoek in eenvoudig Nederlands, met behulp van alledaagse vergelijkingen om de complexe fysica begrijpelijk te maken.

De Basis: Wat is er aan de hand?

Stel je voor dat je een auto rijdt op een drukke snelweg. Als het rustig is (lage snelheid), rijden alle auto's netjes in een rij en houden ze afstand. Dit is wat wetenschappers "evenwicht" noemen. De wiskunde die we normaal gebruiken om dit te beschrijven, heet de Navier-Stokes-vergelijkingen. Dit is als een simpele verkeersregelset die werkt als iedereen rustig rijdt.

Maar wat gebeurt er als er een enorme file ontstaat of een crash? Dan wordt het chaos. Auto's botsen, draaien, remmen en versnellen allemaal tegelijk. De snelheid is zo hoog dat de simpele regels niet meer werken. In de luchtvaart noemen we dit niet-evenwicht (nonequilibrium), vooral bij supersonische vliegtuigen (sneller dan het geluid) die door de lucht schieten en schokgolven creëren.

Het Probleem: De Simpele Regels Falen

De huidige wiskundige modellen (Navier-Stokes) gaan er vanuit dat moleculen (de deeltjes waar lucht uit bestaat) alleen maar als kleine balletjes bewegen. Ze vergeten echter dat moleculen ook kunnen draaien (zoals een tol) en trillen.

• De Analogie: Stel je voor dat je een groep mensen laat rennen.
  • De oude theorie (Navier-Stokes) kijkt alleen naar hoe snel ze rennen.
  • De nieuwe theorie (Boltzmann-Curtiss) kijkt ook naar hoe ze hun armen zwaaien en of ze om hun as draaien terwijl ze rennen.

Bij hoge snelheden (zoals in een schokgolf) is er niet genoeg tijd voor de moleculen om te "kalmeren" en te draaien. Ze blijven draaien terwijl ze alweer botsen. De oude theorie ziet dit niet en voorspelt daarom dat de schokgolf (de "muur" van druk) veel dunner is dan hij in werkelijkheid is. Het is alsof je een file voorspelt die slechts 1 meter lang is, terwijl hij in werkelijkheid 10 kilometer lang is.

De Oplossing: De "Vormende" Theorie (MCT)

De auteurs van dit paper, Mohamed Ahmed en zijn collega's, hebben een nieuwe manier bedacht om deze chaos te beschrijven. Ze gebruiken een theorie genaamd Morphing Continuum Theory (MCT), gebaseerd op de Boltzmann-Curtiss-vergelijking.

In plaats van moleculen als statische balletjes te zien, behandelen ze ze als kleine, levende objecten die zowel vooruit bewegen als rond hun eigen as draaien.

• De Creatieve Vergelijking:
  • Oude methode (Navier-Stokes): Je beschrijft een dansvloer alleen door te kijken naar hoe snel de mensen van links naar rechts lopen.
  • Nieuwe methode (Boltzmann-Curtiss): Je kijkt ook naar hoe ze draaien, springen en met elkaar interacteren. Je ziet dat ze soms vastlopen in een draaiende beweging voordat ze weer vooruit kunnen.

Door deze extra beweging (de rotatie) mee te nemen, krijgen ze een nieuw getal in hun wiskunde: een viskeuze koppeling. Dit is als een "demper" of "schokbreker" in de wiskunde die zorgt dat de schokgolf niet te dun wordt.

Wat hebben ze ontdekt?

Ze hebben dit getest met computersimulaties voor twee soorten gassen:

1. Argon: Een gas met enkelvoudige atomen (zoals een balletje dat alleen rolt).
2. Stikstof (N2): Een gas met twee atomen (zoals een dumbbell die rolt én draait).

Ze hebben gekeken naar schokgolven bij verschillende snelheden (van Mach 1.2 tot Mach 9).

• Het Resultaat: De oude methode (Navier-Stokes) gaf een schokgolf die veel te smal was. De nieuwe methode (Boltzmann-Curtiss) gaf een schokgolf die precies leek op de echte metingen en op de zeer nauwkeurige (maar rekenkracht-geheime) simulaties die bekend staan als DSMC.
• De Analogie: Als je een foto maakt van een schokgolf, ziet de oude theorie eruit als een dunne, scherpe lijn. De nieuwe theorie ziet eruit als een zachte, gebogen boog die precies overeenkomt met wat je in de echte wereld ziet.

Waarom is dit belangrijk?

1. Betere Voorspellingen: Voor vliegtuigen die sneller dan het geluid vliegen, is het cruciaal om te weten hoe warm het wordt en hoeveel druk er op het vliegtuig staat. Als je de schokgolf te smal voorspelt, kun je het vliegtuig verkeerd ontwerpen. Het kan oververhitten of instabiel worden.
2. Rekenkracht: De meest nauwkeurige methode (DSMC) is als het proberen te berekenen van het gedrag van elke individuele muis in een bos. Dat kost een enorme hoeveelheid computerkracht en tijd. De oude methode (Navier-Stokes) is als het tellen van de bomen; dat is snel, maar onnauwkeurig.
   • De nieuwe methode (MCT) is als het tellen van de bomen, maar dan met een slimme schatting van waar de muisjes zitten. Het is sneller dan de dure methode, maar nauwkeuriger dan de simpele methode.

Conclusie

Kortom: De auteurs hebben een nieuwe wiskundige "bril" ontworpen om naar snelle luchtstromen te kijken. Door te erkennen dat luchtmoleculen niet alleen bewegen, maar ook draaien, kunnen ze schokgolven veel beter voorspellen. Dit helpt ingenieurs om veiligere en efficiëntere supersonische vliegtuigen te bouwen, zonder dat ze dagenlang op een supercomputer hoeven te wachten voor het antwoord.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Boltzmann-Curtiss Description for Flows under Translational Nonequilibrium" in het Nederlands.

Probleemstelling

Traditionele continuümtheorieën, zoals de Navier-Stokes (NS) vergelijkingen, blijken ontoereikend voor stromingen onder extreme niet-evenwichtscondities, zoals die voorkomen in schokgolven bij supersonische en hypersonische snelheden.

• Beperkingen van Navier-Stokes: De NS-vergelijkingen zijn afgeleid uit de Chapman-Enskog-expansie van de klassieke Boltzmann-vergelijking, die uitgaat van een lichte afwijking van het translatie-evenwicht. Ze verwaarlozen vaak de rotatie-vrijheidsgraden van moleculen en nemen aan dat de bulkviscositeit nul is (Stokes-hypothese).
• Fysische realiteit: In schokgolven is de karakteristieke lengteschaal vergelijkbaar met de gemiddelde vrije weglengte van gasdeeltjes. Hierdoor zijn intermoleculaire botsingen onvoldoende om alle energiemodi (translatie, rotatie, vibratie) direct in evenwicht te brengen. Dit leidt tot thermische en chemische niet-evenwichtstoestanden.
• Gevolg: NS-simulaties voorspellen vaak schokgolven die veel te dun zijn en onnauwkeurige dichtheids- en temperatuurprofielen, in tegenstelling tot experimentele data en Direct Simulation Monte Carlo (DSMC) resultaten.

Methodologie

De auteurs presenteren een nieuw theoretisch raamwerk gebaseerd op de Boltzmann-Curtiss-verdeling, dat uitbreidt op de klassieke kinetische theorie door de rotatie-vrijheidsgraden expliciet te modelleren.

1. Boltzmann-Curtiss Vergelijking:

   • In plaats van gasmoleculen als puntdeeltjes te behandelen, worden ze beschouwd als structuren met een eindige grootte, inclusief variabelen voor interne vibratie en spin (gyratie).
   • De vergelijking omvat zowel lineaire impuls ($p_i$) als hoekimpuls ($M_i$).
   • De evenwichtsverdelingsfunctie ($f_0$) omvat zowel de verstoorde snelheid ($v'$) als de verstoorde gyratie ($\omega'$).

2. Eerste-orde Benadering:

   • De auteurs gebruiken een eerste-orde benadering van de Boltzmann-Curtiss-vergelijking om de constitutieve relaties af te leiden.
   • Dit leidt tot de Morphing Continuum Theory (MCT), die een koppeling introduceert tussen translatie en rotatie in de spannings-tensor.

3. Nieuw Bulkviscositeitsmodel:

   • Een cruciale afleiding is dat de relaxatietijd ($\tau$) in dit model de tijd is die nodig is om zowel translatie- als rotatie-evenwicht te bereiken.
   • Hieruit volgt een nieuwe uitdrukking voor de bulkviscositeit ($\mu_B$). In tegenstelling tot de NS-theorie (waar $\mu_B = 0$), leidt dit model tot een relatie waarbij de verhouding tussen bulkviscositeit en schuifviscositeit ongeveer $1/3$is ($\lambda' = -1/3 \mu$).
   • Deze term is inherent aan de spannings-tensor en houdt automatisch rekening met afwijkingen van zowel translatie- als rotatie-evenwicht.

4. Numerieke Implementatie:

   • Er zijn simulaties uitgevoerd voor monatomische gassen (Argon) en diatomische gassen (Stikstof/N2) in een Mach-getalbereik van 1,2 tot 9.
   • Een finite volume solver met een tweede-orde flux-splitting schema (KNP) werd gebruikt.
   • De resultaten werden vergeleken met experimentele metingen, DSMC-simulaties en klassieke NS-oplossingen.

Belangrijkste Bijdragen

• Theoretische Uitbreiding: De paper introduceert een rigorieuze afleiding van de bulkviscositeit vanuit de Boltzmann-Curtiss-verdeling, in plaats van deze als een aanpasbare parameter te behandelen. Dit koppelt de bulkviscositeit direct aan de relaxatietijd van rotatie- en translatiebewegingen.
• MCT Framework: Het toepassen van de Morphing Continuum Theory op schokgolven, waarbij de spannings-tensor een extra koppelterm bevat die de interactie tussen lineaire en hoekimpuls beschrijft.
• Verbeterde Voorspellingskracht: Het aantonen dat het nieuwe model, zonder de complexiteit van hogere-orde vergelijkingen (zoals Burnett-vergelijkingen die numeriek instabiel kunnen zijn), betere resultaten oplevert dan de standaard NS-vergelijkingen.

Resultaten

De numerieke simulaties tonen significante verbeteringen ten opzichte van de Navier-Stokes-vergelijkingen:

• Argon (Monatomisch):

  • De NS-vergelijkingen voorspellen schokgolven die tot 60% te dun zijn bij hoge Mach-getallen.
  • Het Boltzmann-Curtiss-model voorspelt een schokdikte die binnen 10% van de experimentele data en DSMC-resultaten ligt.
  • De voorspelde dichtheidsprofielen, normale spanningen en warmtestromen komen veel nauwkeuriger overeen met DSMC-data dan de NS-oplossingen.
  • De resultaten zijn bijna identiek aan die van de complexere Burnett-vergelijkingen.

• Stikstof (Diatomisch):

  • Bij Mach-getallen waar zowel translatie- als rotatie-niet-evenwicht optreedt (bijv. Mach 5 en 10), voorspelt het nieuwe model nauwkeurigere dichtheidsprofielen dan NS.
  • De temperatuurprofielen tonen een afwijking ten opzichte van DSMC, maar dit wordt toegeschreven aan het ontbreken van experimentele temperatuurdata in schokgolven voor validatie. Het model slaagt er echter wel in om de fysica van de rotatie-relaxatie correct te incorporeren.

Betekenis en Conclusie

De studie toont aan dat de Navier-Stokes-vergelijkingen falen bij grote afwijkingen van thermisch evenwicht omdat ze de rotatie-vrijheidsgraden en de bijbehorende bulkviscositeit niet correct modelleren.

• Validiteit: De vergelijkingen afgeleid van de Boltzmann-Curtiss-verdeling zijn geldig voor een breder scala aan niet-evenwichtscondities dan die van de Maxwell-Boltzmann-verdeling.
• Efficiëntie: Het MCT-framework biedt een compromis tussen nauwkeurigheid en rekentijd. Het is nauwkeuriger dan NS en veel efficiënter dan DSMC of Burnett-vergelijkingen. In eerder onderzoek bleek MCT bijvoorbeeld een orde van grootte minder rekentijd te vereisen dan NS voor DNS van turbulente stromingen, terwijl het toch de fysica van schokgolven correct vastlegde.
• Toekomst: De resultaten suggereren dat het gebruik van een correct afgeleide bulkviscositeit, gebaseerd op de Boltzmann-Curtiss-theorie, leidt tot een aanzienlijke verbetering in het modelleren van hypersonische niet-evenwichtsstromingen met continuümoplossers.