Prandtl number dependence of rotating internally heated convection

Dit onderzoek toont aan dat het Prandtl-getal de warmtetransportmechanismen in roterende, intern verwarmde convectie sterk beïnvloedt door de dynamiek van de stabiele onderste laag te reguleren, terwijl de globale gemiddelde temperatuur voornamelijk wordt bepaald door de onstabiele bovenste grenslaag.

De kern

Stel je voor dat je een grote, diepe pan soep op het vuur zet. Maar in plaats van dat de hitte alleen van onderaf komt (zoals bij een gewone pan), wordt de soep overal in de pot tegelijk verwarmd. Dit noemen wetenschappers "intern verwarmde convectie".

In deze studie kijken twee onderzoekers naar wat er gebeurt met die soep als je de dikte van de vloeistof verandert. In de natuurkunde noemen we deze "dikte" de Prandtl-getal.

• Lage Prandtl-getal: Denk aan water of gesmolten ijzer. Dit is dun en vloeibaar; warmte verspreidt zich heel snel, maar de vloeistof zelf is niet erg "stroperig".
• Hoge Prandtl-getal: Denk aan honing of lava. Dit is dik en stroperig; warmte verspreidt zich traag, en de vloeistof wil niet snel bewegen.

De onderzoekers hebben gekeken naar twee scenario's:

1. De pan staat stil (geen draaiing).
2. De pan draait (zoals de aarde die om haar as draait).

Hier is wat ze ontdekten, vertaald naar alledaagse beelden:

1. Als de pan stil staat (Geen rotatie)

Stel je voor dat de hitte de soep laat koken. De bovenkant van de soep wordt onstabiel en begint te borrelen (dit is de "onstabiele laag"). De onderkant blijft echter koel en rustig, alsof er een deksel op zit (dit is de "stabiele laag").

• Bij dunne vloeistoffen (Lage Prandtl-getal, zoals water):
  De turbulentie van het koken bovenin is zo krachtig dat het de onderkant ook gaat "schudden". Het is alsof je een grote lepel door de hele pan roert. De onderkant wordt actief, de soep wordt goed gemengd. De onderzoekers noemen dit een "symmetrie-herstel": de rustige onderkant wordt wakker geschud door de chaos bovenin.

• Bij dikke vloeistoffen (Hoge Prandtl-getal, zoals honing):
  Hier is de onderkant heel anders. Omdat de vloeistof zo stroperig is, kan de turbulentie van bovenin de onderkant niet bereiken. De onderkant wordt een "dode zone". Het is alsof er een dikke laag honing op de bodem ligt die volledig stil blijft, terwijl er bovenop een stormachtige kookpartij gaande is. De warmte wordt dan vooral via de bovenkant afgevoerd.

Het verrassende resultaat:
Ondanks dat de onderkant bij de dikke vloeistof volledig stil ligt en bij de dunne vloeistof wild beweegt, is de gemiddelde temperatuur van de hele soep bijna hetzelfde! De temperatuur wordt vooral bepaald door wat er bovenin gebeurt. De "dikte" van de vloeistof verandert dus hoe de soep binnenin beweegt, maar niet hoe heet de pot gemiddeld is.

2. Als de pan draait (Rotatie)

Nu draai je de pan (of de aarde). Dit voegt een nieuwe kracht toe: de Corioliskracht. Dit zorgt ervoor dat stromingen gaan draaien en kolommen vormen, net zoals een tornado.

• Wat gebeurt er met de warmte?
  Rotatie helpt de warmte om sneller naar boven te komen, maar alleen als de vloeistof niet te dun is.

  • Bij dunne vloeistoffen (water): De rotatie helpt niet echt om de warmte beter af te voeren. De warmte ontsnapt te snel via de zijkanten voordat de draaiing het kan regelen.
  • Bij dikke vloeistoffen (honing): Hier werkt de rotatie als een pomp. Het duwt de warmte efficiënter naar boven. Dit heet "Ekman-pomping". De soep koelt daardoor sneller af en wordt efficiënter gemengd.

• De "Dode Zone" onder rotatie:
  Bij de dikke vloeistoffen blijft de onderkant vaak nog steeds rustig, maar de rotatie zorgt ervoor dat de kolommen van warmte die van boven komen, sterker worden. Het is alsof de draaiing de "torens" van warmte versterkt, terwijl de bodem eronder nog steeds in slaap blijft liggen.

Waarom is dit belangrijk?

Dit onderzoek is niet alleen over soep; het gaat over planeten en sterren.

• De kern van de aarde is gesmolten ijzer (dun, lage Prandtl-getal).
• De mantel van de aarde is rotsachtig lava (dik, hoge Prandtl-getal).
• De atmosfeer van Jupiter is gas (zeer lage Prandtl-getal).

De onderzoekers laten zien dat we niet kunnen zeggen "rotatie maakt alles beter gemengd". Het hangt af van hoe "stroperig" het materiaal is.

• In een dunne kern (ijzer) zorgt rotatie voor andere effecten dan in een dikke mantel (lava).
• Als je probeert te begrijpen hoe de aarde warmte afvoert of hoe de atmosfeer van een planeet werkt, moet je rekening houden met deze "dikte". Een model dat werkt voor water, faalt voor lava.

Kortom:
De "dikte" van een vloeistof bepaalt of de onderkant van een verwarmde laag actief wordt of in een droomtoestand blijft. Rotatie kan de warmte-afvoer verbeteren, maar alleen als de vloeistof dik genoeg is om die rotatiekracht te "grijpen". Het is een delicate dans tussen warmte, beweging en de dikte van het materiaal.

Technische samenvatting

Probleemstelling

Het artikel onderzoekt de invloed van het Prandtl-getal ($Pr$) op penetratief intern verwarmde convectie (IHC), zowel in niet-roterende als roterende regimes. Intern verwarmde convectie is een fundamenteel proces in de atmosfeer en binnenste lagen van planetaire en stellaire objecten (zoals de aardkern en mantel), waar warmte volumetrisch wordt gegenereerd in plaats van via een grenslaag.

Een centrale uitdaging in het modelleren van deze systemen is de enorme variatie in de vloeistofeigenschappen, specifiek het Prandtl-getal (de verhouding tussen viskeuze en thermische diffusie), dat varieert van $O(10^{-2})$ in vloeibare ijzerkernen tot $O(10^{23})$ in planetaire mantels. Hoewel het effect van rotatie en $Pr$ op klassieke Rayleigh-Bénard-convectie (RBC) goed bestudeerd is, is de interactie tussen interne verwarming, rotatie en een breed scala aan $Pr$-waarden nog onvoldoende onderzocht. In IHC ontstaat een unieke dynamiek met een onstabiel gestratificeerde bovenlaag en een stabiel gestratificeerde onderlaag, wat de stromingsdynamiek ingewikkelder maakt dan bij RBC.

Methodologie

De auteurs hebben driedimensionale Directe Numerieke Simulaties (DNS) uitgevoerd met de open-source code AFiD.

• Systeem: Een vloeistoflaag met diepte $d$, begrensd door twee horizontale platen met gelijke temperatuur ($T=0$), waarbij de vloeistof uniform intern wordt verwarmd. Het systeem roteert om de verticale as met een hoeksnelheid $\Omega$.
• Vergelijkingen: De Navier-Stokes-vergelijkingen onder de Boussinesq-benadering, genormaliseerd met de schaal van Roberts (1967).
• Parameterbereik:
  • Prandtl-getal: $Pr \in [0.1, 100]$ (met specifieke focus op $0.1, 0.3, 1, 3, 10, 30, 100$).
  • Rayleigh-getal: $R \in [3.16 \times 10^5, 10^{10}]$.
  • Ekman-getal: $E \in [10^{-6}, \infty]$ (waarbij $E=\infty$ staat voor het niet-roterende geval).
• Oplossing: De simulaties zijn opgezet met periodieke randvoorwaarden in horizontale richting en niet-glijdende/isotherme randvoorwaarden aan de platen. De resolutie is strikt getest om ervoor te zorgen dat dissipatievergelijkingen binnen 1% nauwkeurig zijn.

Belangrijkste Bijdragen

1. Uitgebreid Parameterbereik: Dit is de eerste studie die DNS toepast op roterende IHC over een breed bereik van $Pr$, inclusief het laag-$Pr$-regime dat computationeel zeer veeleisend is.
2. Scheiding van Dynamische Regimes: Het artikel onderscheidt duidelijk het gedrag van de onstabiele bovenlaag (die de globale temperatuur bepaalt) van de stabiele onderlaag (die sterk afhankelijk is van $Pr$).
3. Rotatie-effecten bij Laag $Pr$: Het onthult dat de bekende "Ekman-pomping" die warmtetransport in RBC verbetert, niet optreedt bij lage $Pr$ in IHC, terwijl het wel aanwezig is bij hoge $Pr$.

Resultaten

1. Niet-roterende IHC

• Globale Temperatuur ($\langle T \rangle$): De gemiddelde temperatuur is opmerkelijk ongevoelig voor $Pr$ en wordt voornamelijk bepaald door de dynamiek van de onstabiele bovenste grenslaag. De schaling volgt $\langle T \rangle \sim R^{-0.2}$.
• Stromingsmorfologie:
  • Laag $Pr$($0.1 - 1$): De hoge thermische diffusie zorgt ervoor dat turbulentie vanuit de bulk de stabiele onderlaag "agiteert". Dit leidt tot een "symmetrieherstel" waarbij de onderlaag actief wordt en turbulentie vertoont.
  • Hoge $Pr$($10 - 100$): De hoge viscositeit en lage thermische diffusie onderdrukken fluctuaties in de onderlaag. De stroming wordt gedomineerd door coherente, verticale plumes die voornamelijk in de bovenste helft actief zijn. Bij $Pr=100$ ontstaat een volledig "dode zone" (quasi-still) aan de onderkant.
• Warmtetransport: De verticale convectieve warmtestroom $\langle wT \rangle$ is afhankelijk van $Pr$. Bij lage $Pr$ wordt meer warmte via de onderkant afgevoerd door de turbulentie, wat resulteert in een lagere efficiëntie in het transport van warmte naar de bovenkant vergeleken met hoge $Pr$-gevallen.

2. Roterende IHC

• Versterking van Convectie: Rotatie verhoogt de verticale convectieve warmtestroom $\langle wT \rangle$ over het hele bereik van $Pr$.
• Efficiëntie en Ekman-pomping:
  • De globale koelingsefficiëntie (gemeten als reductie in $\langle T \rangle$) verbetert alleen voor $Pr \geq 1$. Dit komt door het optreden van Ekman-pomping, waarbij de Corioliskracht kolommen vormt die warmte efficiënter transporteren.
  • Bij **lage $Pr$** ($<1$) is er geen verbetering in de globale efficiëntie. De hoge thermische diffusie laat warmte ontsnappen voordat deze verticaal kan worden getransporteerd door de kolomstructuren; de Ekman-pomping is hier dus ineffectief.
• Grenslaagdynamiek: Rotatie introduceert een temperatuurgradiënt in de bulk en verandert de dikte van de grenslagen. De bovenste viskeuze grenslaag volgt de klassieke schaling $\delta_\nu \sim E^{1/2}$, maar de interactie met de stabiele onderlaag maakt de overgang naar het geostrofische regime complexer dan in RBC.
• Dissipatie: Bij hoge rotatie verschuift de dominante bron van thermische dissipatie van de bovenste grenslaag naar de bulk, wat wijst op een fundamentele verandering in het transportmechanisme.

Betekenis en Conclusie

De studie toont aan dat hoewel intern verwarmde convectie (IHC) op de bovenste grenslaag schalingen deelt met Rayleigh-Bénard-convectie, de interne stratificatie zorgt voor een unieke en kritische gevoeligheid voor het Prandtl-getal.

• Geofysische Relevantie: De bevindingen zijn cruciaal voor het begrijpen van warmtetransport in planetaire binnenkernen en mantels. Bijvoorbeeld, in systemen met hoge $Pr$ (zoals de aardmantel) kan de onderste laag van het convectieve systeem effectief "verdoezeld" worden door stabiele stratificatie, waardoor de effectieve diepte van convectie dynamisch wordt bepaald door $Pr$.
• Rotatie-effect: Het bevestigt dat rotatie de warmtetransportmechanismen fundamenteel verandert, maar dat dit alleen leidt tot een efficiënter globaal transport bij voldoende hoge $Pr$-waarden. Bij lage $Pr$ (zoals in metalen kernen) blijft het effect van rotatie beperkt tot lokale herordening zonder globale efficiëntiewinst.

Kortom, dit werk biedt een nieuw inzicht in hoe vloeistofeigenschappen en rotatie samenwerken om de thermische evolutie van hemellichamen te sturen, en benadrukt dat modellen die $Pr$ negeren of aannemen dat IHC identiek is aan RBC, tekortschieten in het voorspellen van transport in deze complexe systemen.