Eddy thermal diffusivity model and mean temperature profiles in turbulent vertical convection

In dit artikel wordt een ruimteafhankelijk model voor wervelende thermische diffusiviteit voorgesteld om analytische resultaten voor het gemiddelde temperatuurprofiel in turbulente verticale natuurlijke convectie af te leiden, waarbij wordt aangetoond dat deze profielen door twee universele schalingsfuncties worden beschreven die in goede overeenstemming zijn met directe numerieke simulaties.

De kern

De Warmte-Verkeersregelaar: Hoe Turbulente Stroming Temperatuur Bepaalt

Stel je voor dat je twee enorme, oneindig hoge muren hebt. De ene muur is gloeiend heet (zoals een oven), en de andere is ijskoud (zoals een vriezer). Tussen deze muren zit lucht of water. Wat gebeurt er? De warme lucht wil omhoog, de koude lucht wil omlaag. Er ontstaat een dans van stroming, een wirwar van wervels en stromingen. Dit noemen we turbulente verticale convectie.

In dit wetenschappelijke artikel kijken Ho Yin Ng en Emily Ching naar precies dit fenomeen. Hun doel? Een simpele, maar krachtige formule vinden die vertelt hoe de temperatuur precies verloopt tussen die twee muren, of je nu met lucht (zoals in een kantoor) of met water (zoals in de oceaan) te maken hebt.

Hier is de kern van hun ontdekking, vertaald naar alledaagse taal:

1. Het Probleem: Een Chaos van Warmte

Wetenschappers weten al lang hoe dit werkt als de stroming rustig en voorspelbaar is (zoals een kalme rivier). Maar als de stroming turbulent wordt (zoals een wild stromende bergbeek), wordt het een chaos. De warmte wordt niet meer alleen door geleiding overgedragen, maar door de chaos van de wervels zelf.

Vroeger probeerden wetenschappers dit te modelleren met formules die leken op "inverse kubieke wortels" of "logaritmische lijnen". Het probleem? Die oude formules werkten goed voor lucht, maar faalden volledig als je het over water of andere vloeistoffen had. Het was alsof je probeerde een auto te besturen met een fietskaart: het werkt soms, maar niet overal.

2. De Oplossing: De "Verkeersregelaar" (Eddy Diffusiviteit)

De auteurs van dit artikel hebben een nieuw idee bedacht. Ze kijken naar iets dat ze "wervel-thermische diffusiviteit" noemen.

• De Analogie: Stel je voor dat de warmte een pakketje is dat van de warme muur naar de koude muur moet.
  • In de rustige zone vlak tegen de muur aan, is het pakketje vastgebonden aan de muur. Het moet voorzichtig en langzaam bewegen. Hier spelen de moleculen van het materiaal (de "moleculaire diffusiviteit") een grote rol.
  • In de chaotische zone in het midden, wordt het pakketje opgepakt door een enorme wervel en geslingerd. Hier is de "wervel-diffusiviteit" de baas.

De auteurs hebben een drie-laags model bedacht om dit te beschrijven:

1. De Muur-zone (Binnen): Hier is de stroming nog rustig. De warmte beweegt als een slak.
2. De Overgangszone (Midden): Hier begint het gevecht tussen rust en chaos.
3. De Centrum-zone (Buiten): Hier woedt de storm. De warmte wordt razendsnel rondgegooid.

Ze hebben een wiskundige "schakelaar" bedacht die precies zegt: "Op dit punt is de chaos net zo sterk als de rust, en op dat punt is de chaos de absolute baas."

3. De Grote Ontdekking: Twee Universele Regels

Het meest fascinerende aan hun werk is dat ze ontdekten dat je de hele temperatuurverdeling kunt beschrijven met slechts twee universele regels (functies):

• Regel 1 (Bij de muur): Dit is de "slak-regel". Het beschrijft hoe de temperatuur begint te dalen direct tegen de warme muur. Dit gedrag is voor alle vloeistoffen hetzelfde, als je het op de juiste manier meet.
• Regel 2 (In het midden): Dit is de "storm-regel". Dit beschrijft hoe de temperatuur zich gedraagt in het midden van de ruimte, waar de wervels het hardst waaien. Ook dit is universeel.

Het mooie is: of je nu lucht, water of olie gebruikt, en of de temperatuurverschillen groot of klein zijn, deze twee regels werken altijd. Het is alsof ze hebben gevonden dat alle auto's, of het nu een Ferrari of een oude bestelbus is, op hetzelfde moment in de versnelling schakelen als ze een bepaalde snelheid bereiken.

4. Waarom is dit belangrijk?

Vroeger moesten wetenschappers voor elke nieuwe vloeistof of elke nieuwe temperatuur een nieuwe, ingewikkelde berekening maken. Met dit nieuwe model kunnen ze:

• Voorspellen hoe warmte zich verspreidt in gebouwen (voor betere ventilatie).
• Begrijpen hoe ijskappen in de poolzeeën smelten door de interactie met het water.
• Ontwerpen van efficiëntere industriële processen.

Ze hebben hun theorie getest met supercomputersimulaties (DNS-data) en bleek dat hun nieuwe "twee-regels-model" veel nauwkeuriger is dan de oude formules. Het past de data perfect, zelfs voor vloeistoffen die heel anders zijn dan lucht.

Samenvattend

Ng en Ching hebben een nieuwe "landkaart" getekend voor warmtestroming. Ze tonen aan dat hoewel turbulentie eruitziet als een chaotische storm, er onder die chaos een strakke, elegante orde schuilt. Door de ruimte op te delen in een "rustige muur-zone" en een "wilde centrum-zone", kunnen we de temperatuur overal en altijd voorspellen met slechts twee simpele formules.

Het is een beetje alsof ze de geheimtaal van de wind hebben vertaald naar een simpele instructie: "Bij de muur loop je langzaam, in het midden ren je als gek, en de overgang gebeurt op een heel specifiek punt."

Technische samenvatting

Probleemstelling

Het artikel richt zich op het begrijpen van turbulente verticale natuurlijke convectie in een vloeistof die is opgesloten tussen twee oneindige verticale wanden met verschillende temperaturen ($T_h$ en $T_c$). Hoewel laminaire convectie goed begrepen is, ontbreekt er nog steeds een volledig theoretisch inzicht in de turbulente regime, vooral wat betreft de gemiddelde temperatuurprofielen ($\bar{T}$) en snelheidsprofielen.

De huidige uitdagingen in de literatuur zijn:

• Bestaande modellen (zoals die van George & Capp, 1979, en Hölling & Herwig, 2005) voorspellen temperatuurprofielen die vaak afwijken van Direct Numerical Simulation (DNS) data, vooral voor verschillende Prandtl-getallen ($Pr$).
• Sommige modellen schenden de vereiste randvoorwaarden bij de wanden (waarbij de turbulentieflux en zijn afgeleiden nul moeten zijn).
• Er is behoefte aan een model dat de gemiddelde temperatuurprofielen voor een breed scala aan Rayleigh-getallen ($Ra$) en Prandtl-getallen universeel kan beschrijven.

Methodologie

De auteurs stellen een nieuw ruimte-afhankelijk model voor de wervelende thermische diffusiviteit ($K(x)$) voor. De aanpak omvat de volgende stappen:

1. Governing Equations: De gemiddelde warmtetransportvergelijking wordt geïntegreerd, waarbij de turbulentieflux ($\overline{u'T'}$) wordt gemodelleerd via een eddy-diffusiviteit benadering: $\overline{u'T'} = -K(x) \frac{d\bar{T}}{dx}$.
2. Drie-laags Model voor $K(x)$: Gezien de verschillende fysica in de buurt van de wanden (innerlijke laag) en in het midden van de kanaal (buitenste laag), wordt $K(x)$ gemodelleerd als een functie van drie gebieden:
   • Innerlijke laag (dicht bij de wand): Een kubische functie ($Ay^3$) die voldoet aan de randvoorwaarde dat de flux en zijn afgeleiden nul zijn bij de wand ($x=0$).
   • Buitenste laag (dicht bij het centrum): Een symmetrische kwadratische functie met een maximum in het midden ($x=H/2$).
   • Intermediaire laag: Een lineaire functie die de twee andere regio's continu verbindt.
3. Analytische Afleiding: Door dit model voor $K(x)$ in de geïntegreerde warmtetransportvergelijking te substitueren, leiden de auteurs analytische uitdrukkingen af voor het gemiddelde temperatuurprofiel $\bar{T}(x)$.
4. Validatie: Het model wordt gevalideerd met behulp van uitgebreide DNS-data (Howland et al., 2022) voor $1 \le Pr \le 100$ en $10^6 \le Ra \le 10^9$.

Belangrijkste Bijdragen

• Universele Schalingsfuncties: De analyse toont aan dat de gemiddelde temperatuurprofielen voor verschillende $Ra$ en $Pr$ worden beschreven door twee universele schalingsfuncties:
  • $F_i$ voor de innerlijke regio (naast de wanden).
  • $F_o$ voor de buitenste regio (dicht bij het centrum).
• Karakteristieke Scales: De auteurs introduceren nieuwe karakteristieke temperatuur- en lengteschalen ($T_i, l_i$ voor de binnenlaag en $T_o, H$ voor de buitenlaag) die afhankelijk zijn van twee modelparameters ($A$ en $C_m$). Deze parameters bepalen de karakteristieke snelheden voor warmtetransport in de respectievelijke regio's.
• Correctie van Bestaande Modellen: Het artikel demonstreert dat eerdere modellen (zoals die van George & Capp) falen omdat ze de afhankelijkheid van het Prandtl-getal in de schalingsfuncties niet correct meenemen of verkeerde randvoorwaarden gebruiken. Het nieuwe model lost dit op door expliciet rekening te houden met moleculaire diffusiviteit in de binnenlaag.

Resultaten

1. Overeenkomst met DNS: De analytische resultaten tonen uitstekende overeenkomst met de DNS-data. De gemiddelde temperatuurprofielen voor alle geteste $Pr$-waarden (1 tot 100) en $Ra$-waarden collapse naar de voorspelde universele curves $F_i$ en $F_o$.
2. Nusselt-getal Scaling: Het model voorspelt correct dat het Nusselt-getal ($Nu$) schaalt als $Nu \sim Ra^{1/3}$ in de hoge-$Ra$ limiet, wat consistent is met eerdere theoretische en numerieke bevindingen.
3. Foutanalyse:
   • De voorspellingen van het nieuwe model voor $Nu$ hebben relatief kleine fouten (meestal < 10%) ten opzichte van DNS-data.
   • In tegenstelling hiermee leiden modellen die $K(x)$ als puur kubisch over het hele kanaal beschouwen, tot aanzienlijk grotere fouten (tot 47% voor lage $Pr$).
   • De voorspellingen van George & Capp (inverse kubische wortel) en Li et al. (lineaire relatie in een specifieke variabele) sluiten niet goed aan bij de data, vooral niet voor hoge $Pr$-waarden.
4. Parameter Afhankelijkheid: De parameters $A$ en $C_m$ blijken te schalen met $Ra$ en $Pr$ op een manier die de karakteristieke snelheden in de binnen- en buitenlaag bepaalt.

Significantie

Deze studie biedt een robuust en analytisch raamwerk voor het begrijpen van turbulente verticale convectie. De belangrijkste implicaties zijn:

• Verbeterde Voorspellingskracht: Het model levert een nauwkeurigere beschrijving van temperatuurprofielen dan bestaande literatuur, wat essentieel is voor engineering-toepassingen zoals ventilatie in gebouwen en het modelleren van ijs-oceaan interacties in poolgebieden.
• Fundamenteel Inzicht: Het benadrukt het belang van het onderscheiden van de fysica in de wandnabijheid (waar moleculaire diffusiviteit belangrijk is) en de kernregio (waar turbulentie dominant is).
• Unificatie: Door het gebruik van twee universele schalingsfuncties, wordt een brug geslagen tussen verschillende stromingsregimes en vloeistofeigenschappen, wat een stap voorwaarts is in de fundamentele theorie van turbulente natuurlijke convectie.

Kortom, de auteurs hebben een succesvol gesloten model ontwikkeld dat de complexiteit van turbulente warmtetransport in verticale kanalen nauwkeurig beschrijft zonder de noodzaak van empirische aanpassingen voor elke specifieke $Ra$ of $Pr$.