2 results
Global and local statistics in turbulent convection at low Prandtl numbers
- Janet D. Scheel, Jörg Schumacher
-
- Journal:
- Journal of Fluid Mechanics / Volume 802 / 10 September 2016
- Published online by Cambridge University Press:
- 01 August 2016, pp. 147-173
-
- Article
- Export citation
-
Statistical properties of turbulent Rayleigh–Bénard convection at low Prandtl numbers $Pr$, which are typical for liquid metals such as mercury or gallium ($Pr\simeq 0.021$) or liquid sodium ($Pr\simeq 0.005$), are investigated in high-resolution three-dimensional spectral element simulations in a closed cylindrical cell with an aspect ratio of one and are compared to previous turbulent convection simulations in air for $Pr=0.7$. We compare the scaling of global momentum and heat transfer. The scaling exponent $\unicode[STIX]{x1D6FD}$ of the power law $Nu=\unicode[STIX]{x1D6FC}Ra^{\unicode[STIX]{x1D6FD}}$ is $\unicode[STIX]{x1D6FD}=0.265\pm 0.01$ for $Pr=0.005$ and $\unicode[STIX]{x1D6FD}=0.26\pm 0.01$ for $Pr=0.021$, which are smaller than that for convection in air ($Pr=0.7$, $\unicode[STIX]{x1D6FD}=0.29\pm 0.01$). These exponents are in agreement with experiments. Mean profiles of the root-mean-square velocity as well as the thermal and kinetic energy dissipation rates have growing amplitudes with decreasing Prandtl number, which underlies a more vigorous bulk turbulence in the low-$Pr$ regime. The skin-friction coefficient displays a Reynolds number dependence that is close to that of an isothermal, intermittently turbulent velocity boundary layer. The thermal boundary layer thicknesses are larger as $Pr$ decreases and conversely the velocity boundary layer thicknesses become smaller. We investigate the scaling exponents and find a slight decrease in exponent magnitude for the thermal boundary layer thickness as $Pr$ decreases, but find the opposite case for the velocity boundary layer thickness scaling. A growing area fraction of turbulent patches close to the heating and cooling plates can be detected by exceeding a locally defined shear Reynolds number threshold. This area fraction is larger for lower $Pr$ at the same $Ra$, but the scaling exponent of its growth with Rayleigh number is reduced. Our analysis of the kurtosis of the locally defined shear Reynolds number demonstrates that the intermittency in the boundary layer is significantly increased for the lower Prandtl number and for sufficiently high Rayleigh number compared to convection in air. This complements our previous findings of enhanced bulk intermittency in low-Prandtl-number convection.
Local boundary layer scales in turbulent Rayleigh–Bénard convection
- Janet D. Scheel, Jörg Schumacher
-
- Journal:
- Journal of Fluid Mechanics / Volume 758 / 10 November 2014
- Published online by Cambridge University Press:
- 08 October 2014, pp. 344-373
-
- Article
- Export citation
-
We compute fully local boundary layer scales in three-dimensional turbulent Rayleigh–Bénard convection. These scales are directly connected to the highly intermittent fluctuations of the fluxes of momentum and heat at the isothermal top and bottom walls and are statistically distributed around the corresponding mean thickness scales. The local boundary layer scales also reflect the strong spatial inhomogeneities of both boundary layers due to the large-scale, but complex and intermittent, circulation that builds up in closed convection cells. Similar to turbulent boundary layers, we define inner scales based on local shear stress that can be consistently extended to the classical viscous scales in bulk turbulence, e.g. the Kolmogorov scale, and outer scales based on slopes at the wall. We discuss the consequences of our generalization, in particular the scaling of our inner and outer boundary layer thicknesses and the resulting shear Reynolds number with respect to the Rayleigh number. The mean outer thickness scale for the temperature field is close to the standard definition of a thermal boundary layer thickness. In the case of the velocity field, under certain conditions the outer scale follows a scaling similar to that of the Prandtl–Blasius type definition with respect to the Rayleigh number, but differs quantitatively. The friction coefficient $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}c_{\epsilon }$ scaling is found to fall right between the laminar and turbulent limits, which indicates that the boundary layer exhibits transitional behaviour. Additionally, we conduct an analysis of the recently suggested dissipation layer thickness scales versus the Rayleigh number and find a transition in the scaling. All our investigations are based on highly accurate spectral element simulations that reproduce gradients and their fluctuations reliably. The study is done for a Prandtl number of $\mathit{Pr}=0.7$ and for Rayleigh numbers that extend over almost five orders of magnitude, $3\times 10^5\le \mathit{Ra} \le 10^{10}$, in cells with an aspect ratio of one. We also performed one study with an aspect ratio equal to three in the case of $\mathit{Ra}=10^8$. For both aspect ratios, we find that the scale distributions depend on the position at the plates where the analysis is conducted.