Publications

Anhand direkter numerischer Simulationen (DNS) werden detaillierte Untersuchungen in turbulenter Rayleigh-Bénard Konvektion in einer zylindrischen Zelle mit isothermen oberen und unteren Platten und adiabatischen Seitenwänden durchgeführt. Die Schwerpunkte der Untersuchungen sind die Statistiken der Temperatur T, ihrer Fluktuationen Theta, und ihrer Gradienten; die Statistik der thermischen Dissipationsrate Epsilon T und ihre Skalierung mit der Rayleigh-Zahl Ra; die Abhängigkeit des Wärmetransports vom Seitenverhältnis der Zelle γ und damit zusammenhängende Änderungen in der großskaligen Zirkulation (GSZ) sowie die Lagrangesche Teilchendynamik in Konvektion. Die Simulationsparameter sind Ra = 10 7 - 10 9, γ = 0.5 - 12 und die Prandtl-Zahl Pr = 0.7. Die Wahrscheinlichkeitsdichtefunktionen (WDF) von T und Theta weichen in allen Regionen der Zelle vom Gaußschen Fall ab. Anhand der Schiefe von ð Z Theta konnte die Rückkehr der kleinskaligen Turbulenz zur lokalen Isotropie mit zunehmender Ra im Innern der Zelle festgestellt werden. Ähnlich wie beim skalaren Mischen, weicht die WDF der thermischen Dissipationsrate der Temperaturfluktuationen, Epsilon Theta von der Log-Normalform ab. Die WDFs sind stets gestreckte exponentielle Verteilungen, deren Schweife mit wachsender Ra auf Grund zunehmender kleinskaliger Intermittenz weiter werden. Der Grad der Intermittenz ist stärker sowohl in der thermischen Grenzschicht (TGS) als auch im Volumen. Er ist stets stärker im Seitenwandbereich als im Innern der Zelle. Darüber hinaus liefert Epsilon Theta nicht nur im Volumen den dominanten Beitrag zur Gesamtdissipation, sondern trägt auch in der TGS signifikant bei. Das Potenzgesetz <Epsilon T> ˜ Ra Zeta ergibt immer einen negativen Exponenten Zeta, sowohl im Volumen als auch in von Plumes beherrschten Gebieten und im turbulenten Hintergrund. Das steht im Gegensatz zur Skalentheorie des Wärmetransports, stimmt aber gut mit Experimenten überein. Die Nusselt-Zahl, Nu, folgt dem Gesetz Nu = A(γ)×Ra ß (γ) mit einer Potenzgesetz-Abhängigkeit der Parameter A und β. Das Minimum der Kurve Nu(γ) liegt genau dort wo die GSZ einen Übergang von einer großen Rolle zu zwei hat. Nu(γ) variiert zwischen 3%-11% und wird geometrieunabhängig für γ≥ 8. Die Muster im vollen turbulenten Regime haben Ähnlichkeit mit Strukturen im schwach nichtlinearen Regime. Fünfeckige bzw. sechseckige im Rollen werden beobachtet, wenn Gamma ≥ 8. Die Lagrangesche Teilchendispersion in Konvektion zeigt einen Übergang vom ballistischen Regime zum Richardson-Regime, jedoch kein Taylor-Regime auf Grund die Endlichkeit der Konvektionszelle. Die Existenz des Richardsonregimes hängt sensitiv vom Anfangsabstand der Teilchen im Paar ab, ähnlich wie in homogen isotroper Turbulenz. Unser Interpolationsschema gibt die Nusseltzahlen im Lagrangeschen Bezugssystem richtig wieder. Die Statisitik der Komponenten des Beschleunigungsvektors ist sehr intermittent ähnlich zu isotroper Turbulenz. Alle drei Verteilungen fallen im Gegensatz zur Konvektion in einer unendlich ausgedehnten Strömgungsschicht zusammen.