The (Limited) Effect of Viscosity in Multiphase Turbulent Mixing

Probleemstelling
Meerfasig gas is alomtegenwoordig in astrofysische omgevingen zoals het interstellair medium (ISM), het circumgalactisch medium (CGM) en het intraclustermedium (ICM). De dynamica van deze systemen wordt bepaald door het samenspel tussen turbulentie, radiatief koelen en transportprocessen. Hoewel de rol van magnetische velden in het onderdrukken van de Kelvin-Helmholtz-instabiliteit (KHI) en het veranderen van menging in turbulente radiatieve menglagen (TRML's) is onderzocht, blijft de specifieke impact van viscositeit omstreden. Volledig geïoniseerde plasma's worden geacht zeer viskeus te zijn (Spitzer 1962), maar de netto effectieve viscositeit in astrofysische systemen is onzeker. Een kritieke vraag is of viscositeit de turbulentie en menging onderdrukt die noodzakelijk zijn voor de "universele" schaalrelaties tussen turbulente snelheid en koelingsrate (Fielding et al. 2020; Tan et al. 2021) om te gelden, met name in regimes waar het Karlovitz-getal ($Ka < 1$) laminair stroming suggereert.

Methodologie
De auteurs gebruikten geïdealiseerde numerieke simulaties met de Athena-code om de effecten van viscositeit te isoleren. De studie verliep in twee hoofdfasen:

1. Adiabatische 2D-simulaties: Om de kritieke viscositeit ($\eta_{Crit}$) te bepalen die nodig is om de KHI te onderdrukken, simuleerden de auteurs schuifstromen met variërende overdichtheden ($\chi = \rho_{cold}/\rho_{hot}$) en Mach-getallen ($M$). Zij gebruikten een opzet met een planair blok om de afhankelijkheid van overdichtheid te testen en een opzet met een planair vlak om de afhankelijkheid van het Mach-getal te testen, waarbij randartefacten werden vermeden. Zij leidden een analytische schatting af voor $\eta_{Crit}$ op basis van de vergelijking van de viskeuze tijdschaal en de KHI-groeitijdschaal, waarbij zij een harmonisch gemiddelde gebruikten voor kinematische viscositeit om impulsverschillen in meer-dichtheidsvloeistoffen te compenseren.
2. 3D Turbulente Radiatieve Menglagen (TRML's): Met behulp van de kritieke viscositeit die is afgeleid voor een overdichtheid van $\chi=100$, voerden de auteurs 3D-simulaties uit van TRML's met expliciete viscositeit. Zij varieerden het viscositeitsaandeel ten opzichte van $\eta_{Crit}$ (van 0% tot 50%) en het Damköhler-getal ($Da = t_{turb}/t_{cool}$) door de verstooringsgolflengte aan te passen. Dit stelde hen in staat om zowel regimes met zwak koelen ($Da < 1$) als met sterk koelen ($Da > 1$) te verkennen. De simulaties omvatten radiatief koelen (Townsend-algoritme) maar geen expliciete warmtegeleiding, zodat $Ka < 1$ bleef.

Belangrijkste Bijdragen en Resultaten

• Schaalwet voor Kritieke Viscositeit: De studie bevestigt dat de kritieke viscositeit die nodig is om de KHI te onderdrukken lineair schaalt met overdichtheid ($\eta_{Crit} \propto \chi$) in de limiet van hoge overdichtheid. Dit komt omdat het koude, dichte medium het grootste deel van de impuls draagt, waardoor sterkere wrijving nodig is om de instabiliteit te onderdrukken. De auteurs verifieerden ook dat bij hoge Mach-getallen compressibiliteitseffecten de KHI stabiliseren, zelfs zonder viscositeit, waardoor een kritiek Mach-getal ($M_{Crit}$) wordt gedefinieerd waarboven de instabiliteit wordt onderdrukt.
• Regime met Zwak Koelen ($Da < 1$): In dit regime zijn koeltijdschalen langer dan turbulente mengtijden. De auteurs vonden dat viscositeit turbulentie sterk onderdrukt, waardoor de stroming overgaat van turbulent naar laminair. Bijgevolg valt de universele relatie tussen turbulente snelheid ($u'$) en oppervlaktehelderheid ($Q$), die is afgeleid voor niet-viskeuze vloeistoffen, uiteen. Echter, tegenintuïtief blijft de totale bolometrische luminositeit grotendeels onaangetast. De auteurs schrijven dit toe aan een "flessenhalseffect" waarbij de lange koeltijd de totale energieverlies beperkt, ongeacht het mengmechanisme (diffusie versus convectie).
• Regime met Sterk Koelen ($Da > 1$): In dit regime zijn koeltijdschalen korter dan viskeuze tijdschalen. De auteurs ontdekten dat viscositeit een beperkt effect heeft op de dynamica van het systeem. Ondanks de aanwezigheid van viscositeit scherpt het koelproces de temperatuur- en snelheidsgradiënten aan het interface, wat het gladmakende effect van viscositeit effectief tegenwerkt. Dit stelt de KHI in staat zich te ontwikkelen en turbulentieniveaus te handhaven die vergelijkbaar zijn met niet-viskeuze gevallen. Als gevolg daarvan blijven de schaalrelaties tussen $u'$ en $Q$ intact en gedraagt het systeem zich effectief als niet-viskeus.
• Competitie tussen Viscositeit en Koelen: De auteurs leidden een analytische voorwaarde af ($Da > \zeta$, waarbij $\zeta$ het fractie van kritieke viscositeit is) voor wanneer koelen domineert boven viscositeit. Hun numerieke resultaten sluiten aan bij deze voorspelling. In het regime met sterk koelen leidt de toegenomen hoeveelheid gas met intermediaire temperatuur in viskeuze runs (door een iets bredere front) tot een bescheiden toename van de oppervlaktehelderheid (tot een factor van $\sim 2$) in vergelijking met niet-viskeuze runs, maar de fundamentele schaalwetten blijven gelden.

Betekenis en Implicaties
Het artikel stelt dat de rol van viscositeit in meerfasige turbulente menging beperkter is dan naïef wordt verwacht, met name in het regime met sterk koelen dat veel astrofysische omgevingen kenmerkt (bijvoorbeeld galactische uitstromingen, CGM, ICM).

• Observatieve Diagnostiek: De bevindingen suggereren dat gevestigde theoretische voorspellingen die lijnbreedten (turbulentie) relateren aan emissiviteit (koelen) nog steeds kunnen worden toegepast op observationele gegevens, zelfs in systemen waar viscositeit significant zou kunnen zijn, mits het systeem zich in het regime met sterk koelen bevindt.
• Subgrid-modellen: De resultaten bieden een fysische rechtvaardiging voor huidige meerfasige subgrid-modellen die worden gebruikt in grootschalige kosmologische simulaties. Deze modellen vertrouwen vaak op niet-viskeuze, kleinschalige simulaties om massatransferpercentages en koelvoorschriften te informeren. De auteurs tonen aan dat deze voorschriften geldig blijven, zelfs wanneer viscositeit aanwezig is, aangezien het koel-gedomineerde regime de viskeuze onderdrukking van turbulentie en menging effectief tenietdoet.
• Astrofysische Context: De studie benadrukt dat hoewel viscositeit de grootschalige morfologie kan beïnvloeden (bijvoorbeeld in door ram-druk gestroopte staarten), het de energie-dissipatie of massatransferpercentages in de snelle-koel-limiet fundamenteel niet verandert. Dit biedt vertrouwen in het interpreteren van de evolutie van meerfasig gas in galaxiehoven en clusters.

De auteurs wijzen op beperkingen met betrekking tot hun geïdealiseerde opzet, specifiek het gebruik van constante viscositeit in plaats van temperatuur-afhankelijke Spitzer-viscositeit, de afwezigheid van magnetische velden (die anisotrope viscositeit induceren) en het ontbreken van expliciete warmtegeleiding. Zij suggereren deze als avenues voor toekomstig werk om het begrip van transportprocessen in astrofysische plasma's verder te verfijnen.