Turbulence Modelling of Mixing Layers under Anisotropic Strain

✨

Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Each language version is independently generated for its own context, not a direct translation.

Stel je een kom voor met twee vloeistoffen van verschillende kleuren, zoals olie en water, die op een tafel staan. Als je de tafel schudt, wordt de grens tussen hen rommelig en beginnen ze te mengen. Dit is vergelijkbaar met wat er gebeurt in een turbulente menglaag in de fysica: twee vloeistoffen met verschillende dichtheden worden tegen elkaar geduwd, waardoor een chaotische, wervelende warboel ontstaat.

Dit artikel gaat over het begrijpen van wat er gebeurt wanneer je niet alleen de tafel schudt, maar ook de hele kamer waarin het mengen plaatsvindt, uitrekt of samendrukt.

Hieronder volgt een uiteenzetting van het verhaal van het artikel, met behulp van eenvoudige analogieën:

1. De Setting: De Kamer Uitrekken

In veel realistische scenario's – zoals een ster die explodeert (supernova) of een kernfusiebom die wordt gecomprimeerd – is de ruimte waarin de vloeistoffen mengen niet gewoon stilstaand. De ruimte zelf rekt uit of krimpt.

De Analogie: Stel je de menglaag voor als een stuk deeg dat wordt gekneed. Meestal bestuderen wetenschappers hoe het deeg mengt als je het gewoon rondduwt. Maar in dit artikel vragen de auteurs zich af: "Wat gebeurt er als iemand, terwijl je het deeg kneedt, ook de tafel waarop het deeg ligt, uitrekt in de lengte of samendrukt in de breedte?"
Het Probleem: De "rek" (spanning) is niet overal even groot. Als je een elastiekje trekt, wordt het langer in de ene richting maar dunner in de andere. Dit heet anisotrope spanning. De meeste computermodellen die worden gebruikt om deze menging te voorspellen, gaan ervan uit dat de rek in elke richting gelijk is (zoals het opblazen van een perfecte ballon), wat niet overeenkomt met de werkelijkheid.

2. Het Hulpmiddel: Het "K-L"-Model

Om te voorspellen hoe de vloeistoffen mengen, gebruiken de auteurs een computerprogramma genaamd het K-L-turbulentiemodel.

De Analogie: Denk aan dit model als een receptenboek voor het voorspellen van chaos. Het bevat twee hoofdingrediënten die het bijhoudt:
1. Hoeveel energie er in de wervels zit (Turbulente Kinetische Energie).
2. Hoe groot de wervels zijn (Turbulente Lengteschaal).
Het model probeert te raden hoe groot de wervels zullen worden naarmate de vloeistoffen mengen. Het lastige deel is een regel in het recept die de "bulk compressie"-term wordt genoemd. Deze regel vertelt het model hoe de grootte van de wervels verandert wanneer de hele kamer wordt samengedrukt of uitgerekt.

3. Het Experiment: Drie Verschillende Regels Testen

De auteurs draaiden computersimulaties om te zien welke "regel" voor de bulk compressie het beste werkte wanneer de kamer in specifieke richtingen werd uitgerekt. Ze testten drie versies van het recept:

De "Gemiddelde" Regel: Deze gaat ervan uit dat de rek in alle richtingen gelijk is (de standaardinstelling).
De "Lengterichting"-Regel: Deze gaat ervan uit dat de grootte van de wervels verandert op basis van alleen hoe sterk de kamer wordt uitgerekt langs de mengrichting.
De "Zijwaartse" Regel: Deze gaat ervan uit dat de grootte van de wervels verandert op basis van hoe sterk de kamer wordt uitgerekt overdwars de mengrichting (loodrecht op de stroming).

4. De Resultaten: De "Zijwaartse" Regel Wint

De auteurs vergeleken hun computervoorspellingen met zeer gedetailleerde, hoogoplossende simulaties (die fungeren als een "perfecte" referentie).

De Vondst: De standaard "Gemiddelde" regel was acceptabel, maar niet geweldig. De "Lengterichting"-regel maakte de voorspellingen daadwerkelijk slechter.
De Winnaar: De "Zijwaartse" Regel (met gebruikmaking van de dwarsspanning) was het meest accuraat.
Waarom? De auteurs leggen uit dat wanneer je een menglaag uitrekt, de grote "eddies" (wervels) zich verschillend gedragen afhankelijk van de richting. Het blijkt dat de grootte van deze wervels gevoeliger is voor hoe de ruimte verandert zijwaarts (transversaal) dan voor hoe deze verandert in lengterichting. Door de zijwaartse rek te gebruiken om de grootte van de wervels in het recept aan te passen, voorspelde het model de mengbreedte en energie veel nauwkeuriger.

5. Het Grotere Geheel: Een Nieuw "Driedelig" Recept

Het artikel keek ook naar hoe deze complexe vergelijkingen kunnen worden vereenvoudigd tot een "Drijfvermogen-Drag"-model (een eenvoudigere manier om over de menging na te denken).

Ze realiseerden zich dat de "breedte van het mengsel" en de "grootte van de wervels" eigenlijk reageren op verschillende krachten. De breedte rekt uit met de trek in lengterichting, maar de grootte van de wervels reageert op de zijwaartse druk.
De Conclusie: Om de beste voorspelling te krijgen, heb je een model nodig dat deze twee dingen apart behandelt. In plaats van één regel voor alles, heb je een driedelig model nodig dat de breedte en de grootte van de wervels onafhankelijk van elkaar ontwikkelt.

Samenvatting

Kortom, dit artikel gaat over het repareren van een computermodel dat wordt gebruikt om te voorspellen hoe vloeistoffen mengen wanneer de ruimte eromheen wordt vervormd. De auteurs ontdekten dat de standaardmanier waarop wordt berekend hoe de "wervels" krimpen of groeien, verkeerd was voor deze specifieke omstandigheden. Door de regel te wijzigen zodat wordt gekeken naar hoe de ruimte zijwaarts wordt uitgerekt in plaats van dit gewoon te middelen, maakten ze het model veel nauwkeuriger. Dit helpt wetenschappers complexe gebeurtenissen zoals sterexplosies of fusie-energie-experimenten beter te begrijpen, waarbij vloeistoffen voortdurend op ongelijke manieren worden samengedrukt en uitgerekt.

Each language version is independently generated for its own context, not a direct translation.

1. Probleemstelling

Turbulente menglagen, zoals die gegenereerd worden door Rayleigh-Taylor (RTI) en Richtmyer-Meshkov (RMI) instabiliteiten, zijn fundamenteel voor fysica met hoge energiedichtheid, waaronder traagheidsopsluitingfusie (ICF) en astrofysische fenomenen zoals supernova's. In deze scenario's ondervindt de stroming vaak anisotrope vervormingssnelheden als gevolg van radiale beweging in convergente geometrieën (implosies) of divergente geometrieën (explosies).

Huidige Reynolds-gegemiddelde Navier-Stokes (RANS) turbulentiemodellen, met name het veelgebruikte K-L-model (een tweevergelijkingmodel dat oplost voor turbulente kinetische energie $K$ en turbulente lengteschaal $L$ ), vertrouwen doorgaans op sluitingsaannames die zijn afgeleid voor isotrope vervorming of onsamendrukbare stromingen. De standaard sluiting voor de term voor bulkcompressie in de transportvergelijking voor de lengteschaal neemt aan dat de lengteschaal evolueert op basis van de gemiddelde isotrope vervormingssnelheid (divergentie van snelheid gedeeld door 3). Echter, in realistische anisotrope scenario's (bijvoorbeeld sferische implosies) verschillen de vervormingssnelheden in de radiale (axiale) en omtrekrichting (transversale) aanzienlijk. Het artikel adresseert het gebrek aan theoretische definitie en nauwkeurige modellering van hoe bulkcompressie de turbulente lengteschaal beïnvloedt onder deze anisotrope omstandigheden.

2. Methodologie

De auteurs voerden een systematisch onderzoek uit met behulp van één-dimensionale RANS-simulaties om de effecten van vervormingsrichting op de evolutie van menglagen te isoleren.

Numeriek Kader: De simulaties maakten gebruik van het K-L-turbulentiemodel geïmplementeerd in de FLAMENCO-code. Het model lost transportvergelijkingen op voor Favre-gegemiddelde turbulente kinetische energie ( $\tilde{K}$ ) en turbulente lengteschaal ( $L$ ).
Testgevallen: Het onderzoek baseerde zich op het "quarter-scale $\theta$ -group" Richtmyer-Meshkov-instabiliteitsgeval (zwaar-naar-licht configuratie, Atwood-getal $A_t=0.5$ , schok-Machgetal $M_a \approx 1.84$ ).
Toepassing van Vervorming: Om de effecten van anisotropie te isoleren, werden uniforme normale vervormingssnelheden toegepast op het simulatiedomein in twee verschillende richtingen:
- Axiale Vervorming ( $S_A$ ): Vervorming toegepast in de richting van het mengen (loodrecht op het interface).
- Transversale Vervorming ( $S_T$ ): Vervorming toegepast in de richtingen evenwijdig aan het interface (homogeen vlak).
- Zowel profielen met constante snelheid als met constante vervormingssnelheid werden getest.
Variaties in Sluiting: Drie verschillende sluitingen voor de term voor bulkcompressie ( $\bar{\rho} L S_\phi$ $\overset{ρ}{ˉ} L S_{ϕ}$ ) in de lengteschaalvergelijking werden vergeleken:
1. Isotrope Sluiting (Standaard): Schaal $L$ met de gemiddelde vervormingssnelheid ( $S_I = \nabla \cdot \mathbf{u} / 3$ ).
2. Axiale Sluiting: Schaal $L$ met de axiale vervormingssnelheid ( $S_A$ ).
3. Transversale Sluiting: Schaal $L$ met de transversale vervormingssnelheid ( $S_T$ ).
Validatie: De RANS-resultaten werden vergeleken met data van hoge fideliteit Implicit Large Eddy Simulation (ILES) uit eerdere studies (Pascoe et al.), die dienden als grondwaarheid.

3. Belangrijkste Bijdragen

Identificatie van Tekortkomingen in Sluiting: Het artikel toont aan dat de standaard isotrope sluitingsaanname faalt in het nauwkeurig voorspellen van de evolutie van turbulente menglagen onder anisotrope vervorming, met name wat betreft de turbulente kinetische energie (TKE) en de integraalbreedte.
Voorstel van Transversale Sluiting: De auteurs stellen een nieuwe sluitingsstrategie voor en valideren deze, waarbij de evolutie van de turbulente lengteschaal onder bulkcompressie wordt gedreven door de transversale vervormingssnelheid ( $S_T$ ) in plaats van de gemiddelde isotrope vervorming.
Afleiding van een Drie-vergelijking Drijfkracht-Weerstandsmodel: Via zelfgelijkvormigheidsanalyse hebben de auteurs een vereenvoudigd drijfkracht-weerstandsmodel afgeleid. Zij vonden dat een enkele lengteschaal niet gelijktijdig de schaling van de integraalbreedte (die volgt op axiale vervorming) en de turbulente lengteschaal (die volgt op transversale vervorming) kan vastleggen. Derhalve ontwikkelden zij een drie-vergelijkingmodel dat de integraalbreedte en de turbulente lengteschaal afzonderlijk doet evolueren.
Generalisatie naar Andere Modellen: Het onderzoek breidt de bevindingen uit naar andere populaire twee-vergelijkingmodellen (K- $\epsilon$ en K- $\omega$ ), waarbij equivalente modificaties voor hun termen voor bulkcompressie worden afgeleid op basis van de transversale vervormings-sluiting.

4. Resultaten

Integraalbreedte ( $W$ ):
- Onder axiale vervorming wordt de breedte van de menglaag direct uitgerekt of gecomprimeerd door het snelheidsgradiënt van de gemiddelde stroming. Alle sluitingen vingen deze trend op, maar de transversale sluiting leverde de meest nauwkeurige voorspelling van de turbulente groeisnelheid op, waarbij de Gemiddelde Absolute Percentage Fout (MAPE) met een factor vier werd verkleind ten opzichte van de isotrope sluiting.
- Onder transversale vervorming toonde de isotrope sluiting bijna geen variatie in breedte (onjuist het opheffen van schuifproductie en compressie-effecten), terwijl de transversale sluiting correct de lichte toename in groeisnelheid voorspelde die in ILES-data werd waargenomen.
Turbulente Kinetische Energie (TKE):
- De transversale sluiting presteerde aanzienlijk beter dan de andere methoden in het voorspellen van TKE.
- De axiale sluiting presteerde het slechtst, waarbij vaak het verkeerde teken van verandering voor TKE onder expansie werd voorspeld.
- De isotrope sluiting overschatte de invloed van vervorming op TKE.
Zelfgelijkvormigheidsanalyse:
- De analyse onthulde dat onder anisotrope vervorming de verhouding tussen de piek van de turbulente lengteschaal en de integraalbreedte ( $\beta_W$ ) niet constant is, waardoor de standaard zelfgelijkvormigheidsaanname die in veel K-L-calibraties wordt gebruikt, wordt geschonden.
- Dit maakte de ontwikkeling noodzakelijk van het drie-vergelijking drijfkracht-weerstandsmodel, waarbij de integraalbreedte schaalt met axiale vervorming ( $S_A$ ) en de turbulente lengteschaal schaalt met transversale vervorming ( $S_T$ ).
Model-equivalentie: De afgeleide modificaties voor het K-L-model werden succesvol gemapt naar de K- $\epsilon$ en K- $\omega$ modellen, wat suggereert dat de term voor bulkcompressie in deze modellen ook moet worden aangepast om rekening te houden met transversale vervorming in plaats van alleen gemiddelde divergentie.

5. Betekenis

Dit werk biedt een kritieke vooruitgang in de modellering van samendrukbare turbulente menging in toepassingen met hoge energiedichtheid.

Verbeterde Voorspellende Capaciteit: Door de sluiting voor bulkcompressie te corrigeren om anisotropie in aanmerking te nemen, kan het K-L-model (en zijn equivalenten) nu de groei van menglagen en TKE in implosie/explosie-scenario's nauwkeuriger voorspellen, wat centraal staat in Traagheidsopsluitingfusie (ICF) en astrofysische simulaties (bijvoorbeeld supernova's).
Fysisch Inzicht: Het onderzoek verduidelijkt de onderscheiden fysische rollen van axiale en transversale vervorming: axiale vervorming drijft voornamelijk de geometrische uitrekking van de menglaag aan, terwijl transversale vervorming de evolutie van de turbulente wervels (lengteschaal) en dissipatie regelt.
Praktische Toepassing: De voorgestelde modificaties zijn relatief eenvoudig te implementeren in bestaande RANS-codes, maar bieden aanzienlijke verbeteringen in nauwkeurigheid zonder de rekenkosten van LES of DNS. Het afgeleide drie-vergelijking drijfkracht-weerstandsmodel biedt een rekenkundig efficiënt hulpmiddel voor technische toepassingen waar volledige RANS-simulaties te duur zijn.

Concluderend stelt het artikel dat voor anisotrope vervormingsomgevingen de turbulente lengteschaal moet worden gemodelleerd als evoluerend met de transversale vervormingssnelheid, een bevinding die de langdurige aanname van isotrope schaling in standaard turbulentiesluitingen uitdaagt.