Intermittency and non-universality of pair dispersion in isothermal compressible turbulence

Dit artikel toont via directe numerieke simulaties aan dat de statistieken van de paarverspreiding in isotherme, samendrukbare turbulentie niet-universeel zijn en sterk afhangen van de aard van de aandrijvingskracht en het Mach-getal, met name doordat verdubbelingstijden niet voldoen aan multifractale voorspellingen terwijl halveringstijden dat wel doen.

De kern

Titel: Waarom deeltjes in een storm niet altijd hetzelfde gedrag vertonen

Stel je voor dat je twee kleine, onzichtbare balletjes in een enorme, chaotische storm gooit. In de wereld van de natuurkunde noemen we dit "turbulentie". De vraag die deze wetenschappers zich stellen, is heel simpel: Hoe snel raken deze twee balletjes uit elkaar?

In een rustige, onbeweegbare lucht (of een heel soepel vloeibaar medium) weten we al lang hoe dit werkt. Het is een bekende wet van Richard Richardson: als je twee deeltjes loslaat, zullen ze exponentieel snel uit elkaar drijven, alsof ze door een onzichtbare kracht worden weggeblazen. Dit werkt perfect voor water of lucht die niet samendrukbaar is.

Maar wat gebeurt er in de ruimte? In sterrenwolken of bij het ontstaan van sterren is de lucht (gas) niet soepel, maar samendrukbaar. Het kan krimpen en uitdijen, en er ontstaan schokgolven (zoals een knal van een donder of een supersonische jet). Hier werkt de oude wet niet meer.

De auteurs van dit paper hebben gekeken naar wat er gebeurt in zo'n chaotisch, samendrukbaar gas. Ze hebben een digitale simulatie gemaakt (een soort super-computer storm) om te kijken hoe twee deeltjes zich gedragen.

Hier zijn de belangrijkste ontdekkingen, vertaald in alledaagse taal:

1. Het verschil tussen "uit elkaar drijven" en "naar elkaar toe vallen"

Stel je voor dat je twee deeltjes hebt. Er zijn twee scenario's:

• Scenario A (Verdubbelingstijd): Hoe lang duurt het voordat ze twee keer zo ver uit elkaar staan als nu?
• Scenario B (Halveringstijd): Hoe lang duurt het voordat ze de helft van hun huidige afstand tot elkaar hebben?

In een normale, soepele storm zou je denken dat deze twee tijden symmetrisch zijn. Als het snel is om uit elkaar te drijven, is het ook snel om weer dichterbij te komen.

De verrassing: In een samendrukbaar gas (zoals in de ruimte) is dit niet zo!

• Het gedrag van de deeltjes die uit elkaar drijven (verdubbeling) is heel anders dan dat van de deeltjes die naar elkaar toe vallen (halvering).
• Het is alsof je een ballon opblaast en weer leeglaat. Het opblazen (uitdijen) en het leeglopen (krimpen) volgen totaal andere regels.

2. De "Storm" maakt het verschil

De wetenschappers hebben twee soorten "stormen" nagebootst:

1. De Draaiende Storm (Solenoidal): Denk aan een reusachtige draaikolk of een tornado. Hier draait het gas om.
2. De Drukkende Storm (Irrotational): Denk aan een reusachtige zuiger die het gas in en uit duwt, zonder dat het draait.

De ontdekking:

• Als de storm draait (zoals een tornado), dan wordt het "uit elkaar drijven" bepaald door die draaiing. Maar het gedrag hangt ook sterk af van hoe hard de storm waait (de "Mach-getal", oftewel de snelheid ten opzichte van het geluid).
• Als de storm duwt en trekt (zonder te draaien), dan is het "uit elkaar drijven" heel anders. Het volgt geen bekende regels en hangt niet af van hoe hard de storm waait.

Dit betekent dat er geen universele wet is voor samendrukbaar gas. Je kunt niet zeggen "dit werkt altijd zo", want het hangt af van hoe de storm wordt veroorzaakt.

3. De "Schokgolven" als de echte boosdoeners

In deze stormen ontstaan er scherpe randen, schokgolven.

• Als de deeltjes dichterbij komen (halveringstijd), worden ze bijna altijd naar die schokgolven getrokken. Dit gedrag is voorspelbaar en volgt een vast patroon, ongeacht wat voor storm het is. Het is alsof alle deeltjes naar de randen van een krater worden getrokken.
• Maar als de deeltjes uit elkaar drijven (verdubbelingstijd), is het verhaal complexer. Ze worden beïnvloed door de draaiende wervels en de schokgolven. Dit maakt het gedrag onvoorspelbaar en afhankelijk van de specifieke aard van de storm.

Waarom is dit belangrijk?

Deze ontdekking is cruciaal voor de sterrenkunde.
Sterren worden geboren in enorme wolken van gas en stof. Deze wolken zijn vol turbulentie en schokgolven.

• Als we willen begrijpen hoe gassen zich mengen in deze wolken (bijvoorbeeld voor de vorming van nieuwe sterren of het verspreiden van chemische stoffen), kunnen we de oude, simpele wetten van Richardson niet meer gebruiken.
• De auteurs zeggen: "We moeten onze modellen opnieuw bekijken." De manier waarop dingen mengen in de ruimte is ingewikkelder en minder voorspelbaar dan we dachten.

Samenvattend in één zin:

In een gewone storm gedragen deeltjes zich voorspelbaar, maar in een storm die kan krimpen en uitdijen (zoals in de ruimte), gedragen deeltjes die uit elkaar drijven zich totaal anders dan deeltjes die naar elkaar toe vallen, en dit gedrag hangt af van de precieze aard van de storm, waardoor er geen enkele "one-size-fits-all" wet voor bestaat.

Technische samenvatting

Titel: Intermitentie en non-universaliteit van koppelverspreiding in isotherme comprimeerbare turbulentie

Auteurs: Sadhitro De, Dhrubaditya Mitra en Rahul Pandit.

---

1. Het Probleem

De basis van ons begrip van turbulentie rust op de wet van Richardson voor de verspreiding van deeltjesparen in incompressibele stromingen, waarbij de gemiddelde kwadratische scheiding $\langle R^2(t) \rangle$ schaalt met $t^3$. Deze wet geldt binnen het inertiebereik van incompressibele stromingen. Echter, veel astrofysische systemen (zoals moleculaire wolken) worden gekenmerkt door comprimeerbare turbulentie (van transsonisch tot supersonisch).

Bestaande theorieën, zoals het multifractale model, voorspellen dat de dynamische schalingsexponenten voor het verdubbelen van de afstand tussen deeltjes (doubling time) en het halveren daarvan (halving time) identiek zouden moeten zijn en voldoen aan een specifieke "brugrelatie" met de structuurfuncties van de snelheid. De auteurs stellen de vraag of deze wetten en voorspellingen gelden voor comprimeerbare, isotherme ideale gassen, en of er sprake is van intermitentie (niet-uniformiteit in de statistische verdeling) en universaliteit (onafhankelijkheid van de manier waarop turbulentie wordt aangedreven).

2. Methodologie

De auteurs hebben Directe Numerieke Simulaties (DNS) uitgevoerd van een tweedimensionale, isotherme, comprimeerbare ideale gasstroom.

• Governing Equations: De Navier-Stokes vergelijkingen voor een isotherm gas ($P = c_s^2 \rho$), inclusief een bulkviscositeit om schokgolven op te lossen.
• Aandrijving: De turbulentie wordt aangedreven door grote-schaal krachten op een golfgetal $k_f = 3$. Er zijn twee soorten aandrijvingen onderzocht:
  1. Solenoidaal (divergentievrij): $\nabla \cdot \mathbf{f} = 0$.
  2. Irrotationaal (rotatievrij): $\nabla \times \mathbf{f} = 0$.
• Regimes: Simulaties zijn uitgevoerd in zowel transsonische ($Ma \approx 1$) als supersonische ($Ma > 1$) regimes.
• Numerieke Oplossing: Een pseudo-spectrale code op een $4096^2$ rooster met een exponentiële filter om aliasing te voorkomen en schokgolven stabiel te houden.
• Lagrange Deeltjes: $N^2$ passieve tracers werden uniform verdeeld en gevolgd. Hun beweging wordt bepaald door de lokale snelheid.
• Analyse Methode: In plaats van alleen te kijken naar de tijdsafhankelijkheid van $\langle R^2(t) \rangle$ (wat een te klein schaalbereik bleek te hebben voor robuuste metingen), gebruikten de auteurs de exit-time methode. Hierbij wordt de tijd $\tau$ gemeten die nodig is voor een paar om een drempel te passeren:
  • Verdubbelingstijd ($\tau_D$): Tijd om van $R$ naar $1.5R$ te gaan ($\sigma > 1$).
  • Halveringstijd ($\tau_H$): Tijd om van $R$ naar $0.75R$ te gaan ($\sigma < 1$).
    De statistieken van deze tijden worden geanalyseerd via de negatieve momenten: $\langle \tau^{-p} \rangle \sim R^{-\chi_p}$.

3. Belangrijkste Bijdragen en Resultaten

A. Breuk van de Multifractale Voorspelling

Het multifractale model voorspelt dat de dynamische exponenten voor verdubbeling ($\chi^D_p$) en halvering ($\chi^H_p$) gelijk moeten zijn en voldoen aan de relatie $\chi_p = p - \zeta_p$, waarbij $\zeta_p$ de exponenten zijn van de Eulerische snelheidsstructuurfuncties.

• Resultaat: De auteurs vinden dat $\chi^D_p \neq \chi^H_p$. De multifractale voorspelling geldt alleen voor de halveringstijden, maar faalt voor de verdubbelingstijden. Dit is een fundamenteel verschil met incompressibele turbulentie.

B. Non-Universaliteit en Aandrijvingsafhankelijkheid

De exponenten voor verdubbelingstijden zijn niet universeel; ze hangen sterk af van hoe de turbulentie wordt aangedreven en van het Mach-getal ($Ma$).

• Solenoidale Aandrijving: Voor deze stromingen geldt een nieuwe relatie: $\chi^D_p = p - \zeta^s_p$, waarbij $\zeta^s_p$ de schalingsexponenten zijn van de solenoidale component van de snelheid. Deze exponenten hangen sterk af van het Mach-getal.
• Irrotationale Aandrijving: Voor deze stromingen voldoet $\chi^D_p$ aan geen enkele bekende multifractale relatie en is deze onafhankelijk van het Mach-getal. De groei van de afstand wordt hier niet alleen bepaald door de solenoidale component, maar ook door gebieden met $\nabla \cdot \mathbf{u} > 0$ (expansie).

C. De Rol van Schokgolven en Structuren

Door de Legendre-transformatie toe te passen op de exponenten, kunnen de fractale dimensies $D(h)$ van de structuren die de verspreiding domineren, worden bepaald.

• Halveringstijden ($\chi^H_p$): Worden gecontroleerd door schokgolven, die 1-dimensionale structuren zijn met een Hôlder-exponent $h \to 0$. Dit gedrag is universeel en onafhankelijk van de aandrijving.
• Verdubbelingstijden ($\chi^D_p$):
  • Bij solenoidale aandrijving in supersonische regimes ($S2$) worden de exponenten bepaald door dunne, langwerpige vlekken van hoge vorticiteit nabij schokgolven (niet de schokgolven zelf).
  • Bij transsonische regimes ($S1$) zijn de structuren minder ruimte-invullend dan schokgolven.
  • Dit verklaart waarom de verdubbelingstijden niet-universeel zijn: ze worden beïnvloed door de specifieke interactie tussen de aandrijving, het Mach-getal en de vorming van vorticiteitsstructuren.

D. Richardson's Wet

Hoewel er een machtswet-groei werd gevonden voor $\langle R^2(t) \rangle$, was het schaalbereik te klein om de exponent $z$ (waarbij $\langle R^2 \rangle \sim t^{2z}$) robuust te meten. De studie concludeert dat de klassieke Richardson-wet ($z=3/2$) niet direct en universeel generaliseerbaar is naar comprimeerbare turbulentie zonder rekening te houden met het onderscheid tussen verdubbeling en halvering.

4. Betekenis en Conclusie

Deze studie levert een cruciale bijdrage aan het begrip van turbulentie in astrofysische contexten (zoals sterrenvorming in moleculaire wolken).

1. Fundamenteel Inzicht: De wet van Richardson moet worden herzien voor comprimeerbare stromingen. Het onderscheid tussen het groeien van een paar (verdubbeling) en het krimpen (halvering) is essentieel; ze vertonen verschillende dynamische multiscalings.
2. Non-Universaliteit: In tegenstelling tot de incompressibele theorie, waar inertiebereik-fluctuaties universeel zijn, hangen de statistieken van koppelverspreiding in comprimeerbare turbulentie af van de aard van de externe kracht en het Mach-getal.
3. Beperking van Bestaande Modellen: Het huidige multifractale model is ontoereikend om de statistieken van koppelverspreiding in comprimeerbare turbulentie volledig te beschrijven. Een nieuwe theoretische raamwerk is nodig dat de asymmetrie tussen expansie en compressie in acht neemt.
4. Astrofysische Implicatie: Mengprocessen in grote astrofysische systemen kunnen niet worden beschreven als een schaalafhankelijke effectieve viscositeit (zoals in incompressibele turbulentie), maar vereisen een heroverweging van hoe schokgolven en vorticiteit de transport en menging van gassen beïnvloeden.

De auteurs speculeren dat deze kwalitatieve bevindingen ook geldig zijn voor driedimensionale turbulentie, wat de relevantie voor de bredere astrofysica verder onderstreept.