scholarly journals Axially homogeneous Rayleigh–Bénard convection in a cylindrical cell

2011 ◽  
Vol 691 ◽  
pp. 52-68 ◽  
Author(s):  
Laura E. Schmidt ◽  
Enrico Calzavarini ◽  
Detlef Lohse ◽  
Federico Toschi ◽  
Roberto Verzicco

AbstractPrevious numerical studies have shown that the ‘ultimate regime of thermal convection’ can be attained in a Rayleigh–Bénard cell when the kinetic and thermal boundary layers are eliminated by replacing both lateral and horizontal walls with periodic boundary conditions (homogeneous Rayleigh–Bénard convection). Then, the heat transfer scales like $\mathit{Nu}\ensuremath{\sim} {\mathit{Ra}}^{1/ 2} $ and turbulence intensity as $\mathit{Re}\ensuremath{\sim} {\mathit{Ra}}^{1/ 2} $, where the Rayleigh number $\mathit{Ra}$ indicates the strength of the driving force (for fixed values of $\mathit{Pr}$, which is the ratio between kinematic viscosity and thermal diffusivity). However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh–Bénard convection in a laterally confined geometry – a small-aspect-ratio vertical cylindrical cell – and show evidence of the ultimate regime as $\mathit{Ra}$ is increased: in spite of the lateral confinement and the resulting kinetic boundary layers, we still find $\mathit{Nu}\ensuremath{\sim} \mathit{Re}\ensuremath{\sim} {\mathit{Ra}}^{1/ 2} $ at $\mathit{Pr}= 1$. Further, it is shown that the system supports solutions composed of modes of exponentially growing vertical velocity and temperature fields, with $\mathit{Ra}$ as the critical parameter determining the properties of these modes. Counter-intuitively, in the low-$\mathit{Ra}$ regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions, which can dominate the dynamics and lead to very high and unsteady heat transfer. As $\mathit{Ra}$ is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the root-mean-square velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased.

2014 ◽  
Vol 16 (5) ◽  
pp. 053042 ◽  
Author(s):  
Pavel Urban ◽  
Pavel Hanzelka ◽  
Věra Musilová ◽  
Tomáš Králík ◽  
Marco La Mantia ◽  
...  

2014 ◽  
Vol 758 ◽  
pp. 344-373 ◽  
Author(s):  
Janet D. Scheel ◽  
Jörg Schumacher

AbstractWe 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.


2018 ◽  
Vol 857 ◽  
pp. 374-397 ◽  
Author(s):  
Hadi Rajaei ◽  
Kim M. J. Alards ◽  
Rudie P. J. Kunnen ◽  
Herman J. H. Clercx

Background rotation causes different flow structures and heat transfer efficiencies in Rayleigh–Bénard convection. Three main regimes are known: rotation unaffected, rotation affected and rotation dominated. It has been shown that the transition between rotation-unaffected and rotation-affected regimes is driven by the boundary layers. However, the physics behind the transition between rotation-affected and rotation-dominated regimes are still unresolved. In this study, we employ the experimentally obtained Lagrangian velocity and acceleration statistics of neutrally buoyant immersed particles to study the rotation-affected and rotation-dominated regimes and the transition between them. We have found that the transition to the rotation-dominated regime coincides with three phenomena; suppressed vertical motions, strong penetration of vortical plumes deep into the bulk and reduced interaction of vortical plumes with their surroundings. The first two phenomena are used as confirmations for the available hypotheses on the transition to the rotation-dominated regime while the last phenomenon is a new argument to describe the regime transition. These findings allow us to better understand the rotation-dominated regime and the transition to this regime.


Sign in / Sign up

Export Citation Format

Share Document