Model of incompressible turbulent flows via a kinetic theory

Dit artikel presenteert een uitgebreid kinetisch theoriemodel voor incompressibele turbulente stromingen dat, door aanpassingen voor wandgebonden stromingen en een consistente relaxatietijd, zowel lineaire als niet-lineaire wervelviscositeitsmodellen reproduceert en experimentele data nauwkeurig voorspelt met een fysiek onderbouwde, minder empirische benadering.

De kern

De "Atomaire" Blik op Turbulentie: Een Nieuwe Manier om Stroming te Begrijpen

Stel je voor dat je naar een drukke menigte mensen kijkt die door een station rennen.

• De oude manier (conventionele modellen): Je kijkt naar de menigte als één grote, vloeibare massa. Je zegt: "Gemiddeld bewegen ze snel naar rechts, maar er is wat chaos." Dit werkt goed voor grote stromingen, maar als je de details wilt zien (bijvoorbeeld hoe iemand struikelt of hoe de chaos precies ontstaat), moet je veel gokwerk doen en "empirische" (op ervaring gebaseerde) regels bedenken. Het is alsof je probeert een auto te repareren zonder te weten hoe de motor er van binnen uitziet; je draait maar aan bouten tot het loopt.
• De nieuwe manier (dit artikel): In plaats van naar de menigte als één blok te kijken, kijken we naar elk individu in de menigte. We kijken naar hun snelheid, hun richting en hoe ze met elkaar botsen. Dit is de kinetische theorie. Het is alsof we een super-snelheidscamera gebruiken die elke beweging van elke persoon vastlegt.

De auteurs van dit artikel (Xin, Guo en Chen) hebben een nieuwe wiskundige formule bedacht die deze "individuele" blik op turbulente lucht- en waterstromingen mogelijk maakt. Hier is hoe ze het hebben aangepakt:

1. Het Probleem: De "Gok" in de Oude Modellen

Turbulente stroming (zoals wind in een storm of water in een rivier) is chaotisch. De oude modellen gebruiken een concept genaamd "wervelviscositeit".

• De Metafoor: Stel je voor dat je boter op brood smeert. De oude modellen zeggen: "De boter is dikker als je harder smeert." Maar ze weten niet waarom de boter dikker wordt. Ze gebruiken een getal dat ze handmatig moeten instellen (een "ad-hoc" parameter). Als je de situatie verandert, moet je dat getal vaak opnieuw uitvinden. Het is als een sleutel die maar op één slot past.

2. De Oplossing: De "Kinetische" Benadering

De auteurs gebruiken een theorie die oorspronkelijk voor gassen is bedacht (waar moleculen botsen), maar ze passen dit toe op grote wervels in een vloeistof.

• De Metafoor: In plaats van te raden hoe de boter zich gedraagt, kijken ze naar de moleculen zelf. Ze zeggen: "Als we weten hoe deze deeltjes botsen en hoe ze energie uitwisselen, kunnen we de stroming voorspellen zonder te hoeven gokken."
• Ze hebben een nieuwe formule bedacht die de "ontspanningstijd" (hoe snel een wervel weer normaal wordt na een botsing) beter berekent. Hierdoor komen ze uit op getallen die veel logischer zijn dan de oude, willekeurige instellingen.

3. De Uitbreiding: Van Open Veld naar de Muur

Het originele model werkte alleen voor stromingen in het "open veld" (zoals wind ver weg van de grond). Maar in de echte wereld stromen we vaak langs muren (zoals een vliegtuigvleugel of een pijp).

• Het Probleem: Dicht bij een muur gedraagt de stroming zich heel anders. De wrijving is groot en de wervels worden klein en gedempt.
• De Oplossing: De auteurs hebben een "laag-Reynolds" versie van hun model gemaakt.
  • De Metafoor: Stel je voor dat je door een drukke gang loopt. In het midden van de gang (het "open veld") ren je vrij. Maar als je dicht bij de muur loopt, moet je voorzichtig zijn, je moet je schrap zetten en je bewegingen zijn anders. Het nieuwe model heeft een speciale "muur-regel" bedacht die precies weet hoe je je moet gedragen als je die muur raakt, zonder dat je de hele gang hoeft te herbouwen.

4. De Test: De "Couette-stroming"

Om te bewijzen dat hun nieuwe formule werkt, hebben ze een simpele maar lastige test gedaan: twee platen die tegen elkaar in bewegen met water ertussen (een "Couette-stroming").

• Het Resultaat: Hun model voorspelde de snelheid van het water en de wrijving op de platen bijna perfect, in overeenstemming met super-computersimulaties (DNS) en echte experimenten.
• De Nuance: Het model is heel goed in het voorspellen van de gemiddelde stroming en de schuifkrachten. Het is echter nog niet 100% perfect in het voorspellen van de exacte richting van de kleine trillingen dicht bij de muur. Het is alsof je de snelheid van de auto perfect voorspelt, maar de trillingen van de motor nog niet helemaal precies.

5. Waarom is dit belangrijk?

Dit onderzoek is een grote stap voorwaarts omdat het:

1. Minder gokwerk vereist: De formules komen voort uit de natuurwetten zelf, niet uit handmatige instellingen.
2. Meer details laat zien: Het kan complexe effecten voorspellen die de oude modellen missen (zoals hoe stroming zich gedraagt bij plotselinge veranderingen).
3. Een brug slaat: Het verbindt de wereld van de "grote stromingen" (wat ingenieurs nodig hebben) met de wereld van de "microscopische deeltjes" (waar de natuurkunde echt gebeurt).

Samenvattend:
De auteurs hebben een nieuwe "bril" ontworpen om naar turbulente stroming te kijken. In plaats van naar de stroming te kijken als een vaag, wazig geheel, kijken ze door deze bril naar de individuele bewegingen die de chaos veroorzaken. Hierdoor kunnen ze de stroming nauwkeuriger voorspellen, met minder gissen en meer natuurkunde. Het is een stap richting een universele "taal" voor stroming die in elke situatie werkt, van een kraan tot een orkaan.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Model of incompressible turbulent flows via a kinetic theory" in het Nederlands.

Titel: Model van incompressibele turbulente stromingen via kinetische theorie

Auteurs: Ziyang Xin, Zhaoli Guo en Hudong Chen

---

1. Probleemstelling

Turbulente stromingen zijn overal in de natuur en techniek aanwezig, maar hun complexe aard maakt voorspellende modellering uitdagend. De huidige standaard, de Reynolds-gegemiddelde Navier-Stokes (RANS) vergelijkingen, vertrouwt vaak op ad hoc sluitingsmodellen (zoals lineaire wervelviscositeit) die empirische parameters vereisen. Deze modellen hebben fundamentele beperkingen:

• Ze gaan uit van een schaalenscheiding tussen de gemiddelde en fluctuerende stroming, wat niet altijd geldt.
• Ze kunnen niet-Newtonse effecten (zoals niet-lineaire spanningen) en complexe anisotropieën in de nabijheid van wanden vaak niet goed beschrijven.
• Bestaande kinetische benaderingen voor turbulentie (gebaseerd op de Boltzmann-vergelijking) hebben tot nu toe te kampen met onrealistische transportcoëfficiënten (bijv. een te hoge turbulent Prandtl-getal) en ontbrekende behandelingen voor wandgebonden stromingen (wall-bounded flows).

Het doel van dit onderzoek is om een theoretisch onderbouwd, fysisch consistent kinetisch model voor incompressibele turbulentie te ontwikkelen dat minder afhankelijk is van empirische parameters en zowel vrij als wandgebonden stromingen kan modelleren.

2. Methodologie

De auteurs bouwen voort op een eerder door Chen et al. (2023) ontwikkeld kinetisch model, gebaseerd op de Klimontovich-vergelijking voor vloeistofelementen. De kern van de methodologie omvat:

• Kinetische Representatie: In plaats van moleculaire thermische fluctuaties, worden turbulente fluctuaties gemodelleerd via een snelheidsverdelingsfunctie (VDF) voor vloeistofelementen. De gemiddelde vergelijkingen worden afgeleid uit een ensemble-gemiddelde Boltzmann-BGK (Bhatnagar-Gross-Krook) vergelijking.
• Chapman-Enskog (CE) Expansie: De auteurs voeren een CE-expansie uit om de relatie tussen het kinetische model en traditionele macroscopische sluitingsmodellen (zoals lineaire en niet-lineaire wervelviscositeit) te analyseren. Dit gebeurt tot de tweede orde.
• Herdefinitie van de Relaxatietijd ($\tau$): Een cruciale bijdrage is de herkeuze van de relaxatietijd. In plaats van de eerder gebruikte waarde, kiezen de auteurs voor $\tau = K / (7\epsilon)$ (waarbij $K$ de turbulente kinetische energie en $\epsilon$ de dissipatiesnelheid is). Deze keuze zorgt ervoor dat de afgeleide transportcoëfficiënten (zoals de turbulent Prandtl-getal) beter overeenkomen met gevestigde theorieën.
• Uitbreiding naar Wandgebonden Stromingen:
  • HR-BGK (High-Reynolds): Voor stromingen buiten de viskeuze sublaag wordt gebruik gemaakt van wandfuncties (logaritmische wet) en een niet-evenwicht extrapolatie randvoorwaarde.
  • LR-BGK (Low-Reynolds): Voor de nabijheid van de wand wordt een model ontwikkeld dat dempingsfuncties (van Driest) en een expliciete term voor moleculaire diffusie van $K$ bevat. Hierbij wordt een volledig diffuus Maxwell-randvoorwaarde toegepast om de niet-glijdende wandconditie en nul-turbulentie-energie bij de wand te garanderen.
• Numerieke Implementatie: Het model wordt getest op turbulente vlakke Couette-stroming. Er wordt gebruik gemaakt van gereduceerde verdelingsfuncties en een gediskretiseerde snelheidsmethode om de kosten te verlagen.

3. Belangrijkste Bijdragen

1. Theoretische Analyse en Coëfficiënten: Door de relaxatietijd aan te passen, worden transportcoëfficiënten afgeleid die kwantitatief overeenkomen met standaard $K-\epsilon$ modellen ($C_\mu \approx 0.0816$, $Pr_T = 0.7$). Dit is een verbetering ten opzichte van eerdere kinetische modellen die onrealistische waarden opleverden.
2. Niet-lineaire Sluiting en Hoge Orde Effecten: De tweede-orde CE-expansie levert niet-lineaire wervelviscositeitsmodellen en gesloten uitdrukkingen voor de turbulentie-energieflux ($Q$) op. Dit omvat termen met materiële afgeleiden en hogere orde gradiënten die zelden in traditionele modellen worden meegenomen.
3. Nieuwe Koppeling voor Wandgebonden Stroming: Het artikel introduceert een systematische manier om kinetische theorie toe te passen op wandgebonden turbulentie, inclusief een laag-Reynolds-getal model dat de viskeuze sublaag direct kan oplossen zonder uitsluitend op empirische wandfuncties te vertrouwen.
4. Fysisch Inzicht: Het model toont aan dat gemiddelde turbulente stroming zich gedraagt als een verdunde gasstroom bij een eindig Knudsen-getal, waardoor niet-Newtonse effecten en anisotropieën worden gevangen die door lineaire modellen worden gemist.

4. Resultaten

Het model is gevalideerd tegen experimentele data en Directe Numerieke Simulaties (DNS) voor turbulente Couette-stroming bij verschillende Reynolds-getallen ($Re = 1300$ en $Re = 10133$):

• Gemiddelde Snelheidsprofielen: Zowel het HR-BGK als het LR-BGK model voorspellen de gemiddelde snelheidsprofielen zeer nauwkeurig, inclusief het logaritmische gebied.
• Wrijvingscoëfficiënt: De voorspelde huidwrijvingscoëfficiënten ($C_f$) komen uitstekend overeen met experimentele en DNS-data over een breed bereik van Reynolds-getallen.
• Reynolds-schuifspanning: De verdeling van de Reynolds-schuifspanning wordt nauwkeurig weergegeven, inclusief de piekwaarde en locatie.
• Anisotropie en Normale Spanningen: Hoewel het model in staat is om de anisotropie van de normale spanningen (niet-Newtonse effecten) te vangen, wordt de anisotropie in de nabijheid van de wand (bufferlaag) iets onderschat ten opzichte van DNS. Dit wordt toegeschreven aan de beperkingen van het enkel-relaxatietijd BGK-model in deze complexe regio.
• Niet-evenwicht Eigenschappen: Analyse van de verdelingsfuncties toont duidelijke afwijkingen van de lokale evenwichtsverdeling (Gaussisch), vooral voor impuls- en energietransport, wat de fysische basis van het model bevestigt.

5. Significantie

Deze studie biedt een fundamenteel nieuwe benadering voor turbulente modellering:

• Minder Empirisme: Het model reduceert de afhankelijkheid van ad hoc parameters door deze af te leiden uit kinetische theorie-principes.
• Universeelheid: Het biedt een raamwerk dat zowel lineaire als niet-lineaire effecten, en zowel vrij als wandgebonden stromingen, in één theoretisch kader beschrijft.
• Toekomstperspectief: Het model legt de basis voor het modelleren van complexere fenomenen (zoals samendrukbare stromingen en 3D-effecten) en suggereert dat kinetische theorie een krachtig alternatief is voor traditionele RANS-modellen, vooral in situaties waar schaalenscheiding faalt of niet-Newtonse effecten dominant zijn.

Kortom, dit werk valideert dat een kinetische theorie-benadering voor turbulentie niet alleen theoretisch consistent is, maar ook numeriek accuraat presteert in het voorspellen van complexe wandgebonden stromingen.