Development and application of a multiphase Lagrangian structure function model in anisotropic turbulence

Technische Samenvatting: Ontwikkeling en Toepassing van een Multifasisch Lagrangiaans Structuurfunctiemodel in Anisotrope Turbulentie

Probleemstelling
De energetische respons van traagheidspartikels op turbulente stroming is cruciaal voor het begrijpen van meerfasendynamica en het ontwikkelen van modellen met een gereduceerde orde voor de dispersie van deeltjes. Hoewel talrijke studies de energie van deeltjes in homogene isotrope turbulentie (HIT) hebben onderzocht, vertonen praktische stromingen vaak ruimtelijke inhomogeniteiten en anisotropieën, met name nabij vaste wanden. In wandgebonden turbulente stromingen, zoals turbulente wandlagen (TBL), wordt grootschalige anisotropie gegenereerd door gemiddelde afschuiving (mean shear) en randvoorwaarden. Een belangrijke onopgeloste vraag is hoe deze grootschalige anisotropie wordt overgedragen naar de deeltalfase en of dit voortduurt in de formeel isotrope kleinschalige bereiken. Bovendien vereist de complexe koppeling tussen deeltelijk traagheid, gravitationele bezinking en stromingsanisotropie in deze omgevingen verder onderzoek om de modellering van bezinkende deeltjes met traagheid nauwkeurig te kunnen uitvoeren.

Methodologie
De auteurs hanteren een gecombineerde theoretische en numerieke aanpak om de schaalafhankelijke deeltjesenergetica in een turbulente wandlaag onder condities van netto depositie te onderzoeken.

1. Theoretisch Kader:

   • De studie maakt gebruik van een kinetische theorie-benadering om exacte, zij het niet-gesloten, behoudswetten af te leiden voor de tweede-orde structuurfuncties (SoSF) van de deeltalfase bij willekeurige schalen.
   • De deeltjesdynamica wordt gemodelleerd met behulp van de Maxey-Riley vergelijkingen voor kleine, sferische, niet-interagerende deeltjes in het limiet van een grote deeltel-tot-vloeistof dichtheidsverhouding.
   • De auteurs definiëren een metriek voor de deeltalfase-anisotropie, $\eta^{(p)}$, gebaseerd op het verschil tussen longitudinale en verticale tweede-orde structuurfuncties.
   • Er wordt een asymptotische analyse uitgevoerd op de afgeleide continuümvergelijkingen om de balans van termen die de energiefluxen en anisotropie beheersen te verhelderen, waarbij wordt uitgegaan van machtswet-gedragingen voor relatieve snelheidsmomenten en de radiale distributiefunctie.

2. Numerieke Simulaties:

   • Direct Numerical Simulations (DNS) van een eenrichtingsgekoppelde deeltjesbeladen wandgebonden turbulente stroming worden uitgevoerd bij een matig wrijvingsReynoldsgetal ($Re_\tau = 315$).
   • De simulaties bestrijken de logaritmische laag en de bufferlaag van de TBL.
   • De studie varieert het Stokesgetal ($St^+$) om deeltelijk traagheid te beoordelen en de bezinkingsparameter ($Sv^+$) om de effecten van gravitationele drift te beoordelen.
   • Zowel Eulerische (vloeistof) als Lagrangiaanse (deeltje) perspectieven worden gebruikt, met een specifieke focus op "slab-wise" statistieken waarbij de verticale separatie nul is om de lokale verticale dynamica te isoleren.

Belangrijkste Bijdragen

• Afleiding van Behoudswetten: Het artikel leidt een exacte behoudsvergelijking af voor deeltjesenergetica op willekeurige schalen, waarbij termen voor transport, koppeling aan eerste-orde momenten en energiefluxen expliciet worden geïdentificeerd. Dit biedt een continuümkader voor het analyseren van schaalafhankelijke energieoverdracht in anisotrope meerfasische stromingen.
• Kwantificering van Anisotropieoverdracht: Het werk kwantificeert systematisch hoe grootschalige anisotropie uit de dragende vloeistof wordt "gedoneerd" aan de deeltalfase en hoe deeltelijk traagheid en bezinking deze overdracht modificeren.
• Asymptotische Analyse: De auteurs leveren een asymptotisch argument om de DNS-data te interpreteren, waarbij ze de waargenomen niet-monotone gedraging van anisotropie-metrieken koppelen aan de competitie tussen grootschalige stromingsanisotropie en kleinschalige inter-schaal energieoverdracht.

Resultaten

• Traagheid en Anisotropie: Bij kleine Stokesgetallen vertoont de deeltalfase een bijna-isotropie op kleine schalen, die de vloeistof weerspiegelt. Echter, naarmate de traagheid toeneemt ($St^+ = 10$), behoudt de deeltalfase significante anisotropie, zelfs op kleine schalen, wat aangeeft dat deeltelijk traagheid het geheugen van de grootschalige stromingsanisotropie behoudt.
• Bezinkingseffecten: De invloed van de bezinkingsparameter ($Sv^+$) is sterk afhankelijk van de stroomregio. In de bufferlaag ($z^+ = 30$), waar de onderliggende stroming inherent anisotroop is, leidt een toenemende $Sv^+$ tot een complexe modificatie van structuurfuncties waarbij de verticale component toeneemt terwijl de longitudinale component afneemt. In de logaritmische laag ($z^+ = 90$), waar de kleinschalige stroming isotroop is, drijven hoge bezinkingssnelheden ($Sv^+$) de deeltalfase naar een isotrope staat op kleine schalen.
• Schaal-Niet-Lokaliteit: De deeltelstructuurfuncties vertonen een vlakkere schaling (benaderend $r^{2/3}$) in het dissipatieve bereik vergeleken met de verwachte $r^2$ schaling van de vloeistof. Dit reflecteert de schaal-niet-lokaliteit van traagheidsdeeltjes, primair gedreven door traagheid en gemoduleerd door bezinking.
• Gedrag van de Anisotropie-metriek: De deeltel-anisotropie-metriek $\eta^{(p)}$ vertoont niet-monotoon gedrag naarmate de scheidingsschaal $r$ afneemt. Het bereikt een lokaal minimum voordat het weer toeneemt naarmate $r/z \to 0$. Dit suggereert een competitie tussen de door de grootschalige stroming opgelegde anisotropie en de anisotropie die wordt gegenereerd door de inter-schaal energieoverdracht binnen de deeltalfase zelf.

Betekenis en Claims
Het artikel claimt cruciaal inzicht te bieden in de schaalafhankelijke energiefluxen van traagheidsdeeltjes in anisotrope turbulentie. Door aan te tonen dat deeltelijk traagheid grootschalige anisotropie kan behouden in de formeel isotrope kleine schalen, benadrukken de auteurs een beperking in modellen die lokale isotropie voor de deeltalfase aannemen. De afgeleide behoudswetten en de geïdentificeerde mechanismen van anisotropieoverdracht hebben directe implicaties voor de ontwikkeling van continuümmodellen voor bezinkende deeltjes met traagheid, met name met betrekking tot hoe twee-weg gekoppelde dynamica beïnvloed kan worden door deze niet-lokale anisotrope effecten. Het werk benadrukt dat voor nauwkeurige modellering de gelijktijdige overweging van deeltelijk traagheid en gravitationele drift noodzakelijk is, aangezien hun effecten niet simpelweg additief zijn maar impliciet gekoppeld zijn binnen de deeltjesdynamica.