Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

Hoogresolutie grote-wervel-simulaties tonen aan dat, hoewel passieve scalar-menging in gestratificeerde homogene turbulentie zich in de transversale richting vergelijkbaar gedraagt met niet-gestratificeerde gevallen, stratificatie verticale menging ernstig remt door het onderdrukken van grootschalig roeren, wat leidt tot onderscheidende groeigrenzen en fluctuatie-intensiteiten die specifieke modelleringbenaderingen voor scalair flux onderbouwen.

De kern

Stel je voor dat je een druppel inkt ziet uitlopen in een glas water. Als je het water zachtjes roert, verspreidt de inkt zich gelijkmatig. Zo denken wetenschappers doorgaans over menging in vloeistoffen: turbulentie werkt als een reusachtige lepel, die alles door elkaar roert totdat het uniform is.

Maar wat gebeurt er als dat water gelaagd is? Stel je voor dat het water onderaan zwaar en zout is, terwijl het water bovenaan licht en zoet is. Dit heet stratificatie. In de echte wereld komt dit voor in de oceaan (waar diep water dichter is) en in de atmosfeer (waar de lucht dunner en lichter wordt naarmate je hoger komt).

Dit artikel beschrijft een high-tech computerexperiment dat de volgende vraag stelt: Hoe verandert deze gelaagdheid de manier waarop een "vlek" (zoals vervuiling of rook) zich verspreidt door een turbulente vloeistof?

Hier is het verhaal van wat ze ontdekten, opgesplitst in eenvoudige concepten.

De Opzet: Twee Soorten "Vlekken"

De onderzoekers creëerden een virtuele vloeistof die aanvankelijk chaotisch draaide (turbulent). Vervolgens introduceerden ze twee verschillende "vlekken" (passieve scalairen) om te zien hoe ze zich gedroegen:

1. De Horizontale Vlek: Een laag inkt die zijwaarts uitgespreid was, zoals een plat vel papier dat op het water drijft.
2. De Verticale Vlek: Een laag inkt die verticaal uitgespreid was, zoals een verticale muur van kleur.

Ze voerden twee simulaties uit: één in normaal water (niet-gelaagd) en één in "gelaagd" water (gestratificeerd).

De Grote Ontdekking: Het Verschil tussen "Boven-Onder" en "Zij-aan-Zij"

1. De Zij-aan-Zij Vlek (Transversale Laag)

Wat er gebeurde: Toen de vlek zijwaarts uitgespreid was, verhinderde de gelaagdheid de verspreiding niet. Sterker nog, het verspreidde zich iets sneller in het gelaagde water dan in het normale water.
De Analogie: Stel je een menigte mensen voor die in een gang rennen. Als de vloer perfect vlak is (niet-gelaagd), rennen ze in alle richtingen. Als de vloer een lichte, onzichtbare helling heeft (gestratificeerd), rennen ze zijwaarts nog steeds prima, misschien zelfs iets energiever. De "inkt" verspreidt zich breed in beide gevallen.
De Vangst: Hoewel de algehele verspreiding vergelijkbaar was, was de "inkt" in het gelaagde water meer "spits". In plaats van een gladde gradiënt had het scherpere, meer gezaagde randen. Het was meer "intermitterend", wat betekent dat er zakken met pure inkt en zakken met puur water waren, met minder een gladde overgang in het midden.

2. De Boven-Onder Vlek (Verticale Laag)

Wat er gebeurde: Hier gebeurde de magie (en de beperking). In het normale water verspreidde de verticale vlek zich gemakkelijk boven en onder, net zoals de zijwaartse. Maar in het gelaagde water stopte de verspreiding bijna volledig.
De Analogie: Stel je voor dat je een dikke milkshake probeert te roeren met een lepel. Als je de lepel probeert op en neer te bewegen, verzetten de lagen van de shake zich tegen je. Het zware materiaal wil onderaan blijven, en het lichte materiaal wil bovenaan blijven. Het "roerende" beweging wordt platgedrukt.
Het Resultaat: De verticale vlek groeide een klein beetje in het allerbegin, maar daarna stuitte het op een "plafond". Het kon niet breder worden omdat de stabiele lagen van de vloeistof werkten als een deksel, waardoor de turbulentie niet kon mengen in verticale richting. De vloeistof kon nog steeds zijwaarts draaien, maar kon niet boven en onder mengen.

Waarom Is Dit Belangrijk? (Het "Waarom" Achter het "Wat")

De onderzoekers ontdekten dat de vloeistof zich in verticale richting gedraagt als een veer. Zodra de turbulentie probeert een zware laag omhoog te duwen of een lichte laag omlaag, trekt de zwaartekracht het terug. Dit stopt de "roerende" beweging.

De vloeistof kan echter nog steeds zijwaarts draaien. De "verticale lengte" van de turbulentie wordt dus vergrendeld op een specifieke grootte (bepaald door hoe sterk de zwaartekracht en de vloeistoflagen zijn), en de vlek kan niet groter worden dan die maat.

Het "Recept" voor Voorspelling

Het artikel probeerde ook een eenvoudig wiskundig "recept" te creëren om te voorspellen hoe deze vlekken zich zouden verspreiden, zonder dat een supercomputer nodig was.

• Als je de vorm van de vlek kent: Je kunt een simpele formule met één getal gebruiken om te voorspellen hoe snel het zich zijwaarts verspreidt. Dit werkt zeer goed.
• Als je de vorm niet kent: Dan moet je de vorm raden (aannemende dat het eruitziet als een gladde kromme). Als je dit doet, heb je een formule met twee getallen nodig. Dit werkt uitstekend nadat de vlek de tijd heeft gehad om een ritme te vinden met de draaiende vloeistof.

De Conclusie

• Zijwaartse menging: Stabiele lagen (zoals in de diepe oceaan of de bovenste atmosfeer) stoppen zijwaartse menging niet; ze maken het misschien zelfs iets intenser en gezaagder.
• Verticale menging: Stabiele lagen werken als een rem. Ze voorkomen dat de vloeistof bijna volledig boven en onder mengt.
• Het onderscheid tussen "Roeren" en "Mengen": De vloeistof kan nog steeds "roeren" (bewegen) zijwaarts, maar kan niet "mengen" (verblenden) verticaal, omdat de lagen zich verzetten tegen het verwisselen.

De auteurs merken op dat hun experiment een specifiek type vloeistofeigenschap gebruikte (een Prandtl-getal van 0,7). Ze waarschuwen dat als de vloeistof "dikker" was of andere eigenschappen had (Prandtl-getal > 1), de resultaten zouden kunnen veranderen vanwege een "omgekeerd" effect waarbij de menging zijn eigen opwaartse kracht creëert. Maar voor de omstandigheden die ze testten, geldt de regel "zijwaarts is vrij, boven-onder is geblokkeerd".

Technische samenvatting

Technische Samenvatting: Evolutie van Passieve Scalar Menglagen in Gestratificeerde en Ongestratificeerde Homogene Turbulentie

Probleemstelling
De dispersie van deeltjes en scalairen in stabiel gestratificeerde omgevingen, zoals de bovenste troposfeer, stratosfeer en oceanische gemengde lagen, vormt een kritieke uitdaging in de atmosferische en oceanische wetenschappen. Hoewel de effecten van stabiele stratificatie op het dichtheidsveld en het verval van turbulentie goed gedocumenteerd zijn, blijft de dispersie van passieve scalairen (bijvoorbeeld aerosolen, verontreinigingen, microplastics) minder goed begrepen. Specifiek is het onduidelijk hoe een passieve scalar menglaag (PSML) evolueert wanneer deze wordt geïntroduceerd in homogene isotrope turbulentie (HIT) die vervolgens verval ondergaat onder invloed van een stabiliserend dichtheidsgradiënt. De centrale vraag is hoe de introductie van opwaartse krachtschalen en anisotropie de zelfgelijkvormige evolutie van scalarstatistieken, zoals waargenomen in ongestratificeerde HIT, verandert.

Methodologie
De auteurs hebben gebruikgemaakt van hoogopgeloste large-eddy simulaties (LES) om het verval van zowel ongestratificeerde als stabiel gestratificeerde homogene turbulentie te onderzoeken. De simulaties maakten gebruik van de niet-divergente Navier-Stokes-vergelijkingen met de niet-hydrostatische Boussinesq-benadering.

• Stromingsconfiguratie: Er werden twee simulaties uitgevoerd. De eerste was ongestratificeerde HIT. De tweede was identiek wat betreft de beginvoorwaarden, maar bevatte een uniforme stabiliserende dichtheidsgradiënt ($g \neq 0$), waardoor de stroming na ongeveer één opwaartse krachttijd overging naar sterk gestratificeerde turbulentie.
• Passieve Scalairen: Er werden twee passieve scalar menglagen geïnitieerd: één met een gemiddelde gradiënt in de verticale richting ($\xi_z$) en één in de transversale richting ($\xi_y$). Deze werden gemodelleerd als mengfracties ($0 \le \xi \le 1$).
• Numerieke Details: Om diffusieve effecten die afwijkingen van zelfgelijkvormigheid veroorzaken te minimaliseren, gebruikten de simulaties hyperviscositeit en hyperdiffusiviteit met gelijke coëfficiënten. Het moleculaire Prandtl-getal werd ingesteld op $Pr = 0.7$. De Reynoldsgetallen waren hoog ($Re_h > 15.000$ en $Re_\lambda > 450$ voor het ongestratificeerde geval), wat de aanwezigheid van inertieel-convectieve en inertieel-diffusieve bereiken garandeerde. De resolutie was voldoende om de opwaartse krachtlengteschaal op te lossen ($L_b/\Delta \gg 1$).
• Analyse: De studie volgde de evolutie van scalarprofielbreedtes, variantie, hogere orde momenten (scheefheid, kurtosis), fluxen en mengefficiëntie. Daarnaast werden wervel-diffusiviteitsmodellen voor de transversale scalarflux geëvalueerd onder twee voorwaarden: wanneer het gemiddelde scalarprofiel bekend is en wanneer het moet worden aangenomen.

Belangrijkste Resultaten

1. Transversale Menging ($\xi_y$): In de richting loodrecht op de zwaartekracht is de evolutie van de passieve scalar in het gestratificeerde geval globaal vergelijkbaar met het ongestratificeerde geval, hoewel deze iets sneller verspreidt.

   • De scalarfluctuaties zijn intenser in het gestratificeerde geval, met een piekvariantie die ongeveer 15% hoger ligt.
   • De interface tussen turbulente en niet-turbulente stroming is in het gestratificeerde geval intermitterender.
   • Zodra het scalarveld quasi-evenwicht bereikt met het snelheidsveld, vallen de genormaliseerde fluxprofielen voor zowel gestratificeerde als ongestratificeerde gevallen samen, wat wijst op vergelijkbare onderliggende transportmechanismen ondanks de stratificatie.

2. Verticale Menging ($\xi_z$): Stratificatie remt verticale menging drastisch.

   • De verticale laagbreedte ($\delta_z$) groeit aanvankelijk, maar stopt na ongeveer één opwaartse krachttijd ($t^* \approx 3$) snel met uitbreiden.
   • Dit staken van de groei wordt toegeschreven aan de onderdrukking van grootschalig verticaal roeren door opwaartse krachten.
   • De uiteindelijke breedte van de laag is evenredig met de verticale integraallengteschaal van de horizontale snelheid ($L_v$), die beperkt is om een verticaal Froude-getal van de orde één te behouden ($Fr_v \sim O(1)$).
   • Hoewel het gemiddelde profiel monotoon groeit door moleculaire diffusie, onthult de fluxanalyse een "omgekeerde flux"-fenomeen waarbij roeren soms scalarvezels terugveegt naar hun oorspronkelijke kant, wat consistent is met intern golfachtig gedrag.

3. Mengefficiëntie: De cumulatieve menging (dissipatie van scalarvariantie) in het gestratificeerde geval overtreft op latere tijdstippen uiteindelijk die van het ongestratificeerde geval. Dit wordt toegeschreven aan het feit dat menging volgt op roeren, en dat de horizontale roersnelheid in de gestratificeerde stroming evenredig is met de horizontale lengteschaal, die groter wordt dan in het ongestratificeerde geval.

4. Modellering van Flux:

   • Bekend Profiel: Als het gemiddelde scalarprofiel bekend is, voorspelt een wervel-diffusiviteitsmodel met één constante ($F_y = -c_1 \delta_y v' \partial \xi / \partial y$) de flux effectief voor zowel gestratificeerde als ongestratificeerde gevallen, met $c_1 \approx 0.4$.
   • Aangenomen Profiel: Als de profielform (verondersteld een foutfunctie te zijn) en de scalarlengteschaal moeten worden geschat uit snelheidsstatistieken, is een model met twee constanten vereist. Dit model presteert goed op latere tijdstippen wanneer de scalar in quasi-evenwicht is met het snelheidsveld, waardoor de scalarlengteschaal kan worden geschaald aan de kinetische energiedissipatiesnelheid.

Betekenis en Aanspraken
Het artikel claimt een gedetailleerde karakterisering te bieden van passieve scalar menging in vervallende gestratificeerde turbulentie, waardoor de kloof tussen theoretische voorspellingen van anisotrope diffusie en numerieke waarnemingen wordt overbrugd.

• Anisotropie: De studie bevestigt dat stratificatie een scherp onderscheid creëert in menging: horizontale menging blijft actief en wordt zelfs versterkt, terwijl verticale menging effectief wordt onderdrukt zodra de stroming sterk gestratificeerd wordt.
• Modelleerbaarheid: De auteurs tonen aan dat standaard wervel-diffusiviteitsbenaderingen levensvatbaar blijven voor gestratificeerde stromingen, mits de scalar in quasi-evenwicht is met het snelheidsveld. Zij benadrukken dat het aannemen van een bekend profiel de modelleereisen aanzienlijk vereenvoudigt ten opzichte van het aannemen van een profielform.
• Beperkingen en Toekomstig Werk: De auteurs merken bescheiden op dat hun resultaten specifiek zijn voor een Prandtl-getal van 0,7. Zij stellen expliciet dat recente literatuur suggereert dat hogere Prandtl-getallen ($Pr > 1$) op kleine schalen omgekeerde opwaartse krachtfuxen induceren, wat de vervalsnelheden van de stroming en bijgevolg de mengsnelheden van passieve scalairen zou kunnen veranderen. Derhalve zijn de huidige bevindingen niet direct van toepassing op scenario's met hoge Prandtl-getallen zonder verder onderzoek.

De studie concludeert dat stratificatie weliswaar de verticale dynamiek van een menglaag fundamenteel verandert, maar dat de transversale menging veel kenmerken van ongestratificeerde turbulentie behoudt, zij het met verhoogde intermitterendheid en fluctuatiedichtheid.