Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBoundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution
382
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Results on the Prandtl-Blasius type kinetic and thermal boundary layer thicknesses in turbulent Rayleigh-Benard convection in a broad range of Prandtl numbers are presented. By solving the laminar Prandtl-Blasius boundary layer equations, we calculate the ratio of the thermal and kinetic boundary layer thicknesses, which depends on the Prandtl number Pr only. It is approximated as $0.588Pr^{-1/2}$ for $Prll Pr^*$ and as $0.982 Pr^{-1/3}$ for $Pr^*llPr$, with $Pr^*= 0.046$. Comparison of the Prandtl--Blasius velocity boundary layer thickness with that evaluated in the direct numerical simulations by Stevens, Verzicco, and Lohse (J. Fluid Mech. 643, 495 (2010)) gives very good agreement. Based on the Prandtl--Blasius type considerations, we derive a lower-bound estimate for the minimum number of the computational mesh nodes, required to conduct accurate numerical simulations of moderately high (boundary layer dominated) turbulent Rayleigh-Benard convection, in the thermal and kinetic boundary layers close to bottom and top plates. It is shown that the number of required nodes within each boundary layer depends on Nu and Pr and grows with the Rayleigh number Ra not slower than $simRa^{0.15}$. This estimate agrees excellently with empirical results, which were based on the convergence of the Nusselt number in numerical simulations.
rate research
Read More
We report an experimental study aiming to clarify the role of boundary conditions (BC) in high Rayleigh number $10^8 < {rm{Ra}} < 3 times 10^{12}$ turbulent thermal convection of cryogenic helium gas. We switch between BC closer to constant heat flux (CF) and constant temperature (CT) applied to the highly conducting bottom plate of the aspect ratio one cylindrical cell 30 cm in size, leading to dramatic changes in the temperature probability density function and in power spectral density of the temperature fluctuations measured at the bottom plate, while the dynamic thermal behaviour of the top plate and bulk convective flow remain unaffected. Within our experimental accuracy, we find no appreciable changes in Reynolds number Re(Ra) scaling, in the dimensionless heat transfer efficiency expressed via Nusselt number Nu(Ra) scaling, nor in the rate of direction reversals of large scale circulation.
Recently, in Zhang et al. (2020), it was found that in rapidly rotating turbulent Rayleigh-Benard convection (RBC) in slender cylindrical containers (with diameter-to-height aspect ratio $Gamma=1/2$) filled with a small-Prandtl-number fluid ($Pr approx0.8$), the Large Scale Circulation (LSC) is suppressed and a Boundary Zonal Flow (BZF) develops near the sidewall, characterized by a bimodal PDF of the temperature, cyclonic fluid motion, and anticyclonic drift of the flow pattern (with respect to the rotating frame). This BZF carries a disproportionate amount ($>60%$) of the total heat transport for $Pr < 1$ but decreases rather abruptly for larger $Pr$ to about $35%$. In this work, we show that the BZF is robust and appears in rapidly rotating turbulent RBC in containers of different $Gamma$ and in a broad range of $Pr$ and $Ra$. Direct numerical simulations for $0.1 leq Pr leq 12.3$, $10^7 leq Ra leq 5times10^{9}$, $10^{5} leq 1/Ek leq 10^{7}$ and $Gamma$ = 1/3, 1/2, 3/4, 1 and 2 show that the BZF width $delta_0$ scales with the Rayleigh number $Ra$ and Ekman number $Ek$ as $delta_0/H sim Gamma^{0} Pr^{{-1/4, 0}} Ra^{1/4} Ek^{2/3}$ (${Pr<1, Pr>1}$) and the drift frequency as $omega/Omega sim Gamma^{0} Pr^{-4/3} Ra Ek^{5/3}$, where $H$ is the cell height and $Omega$ the angular rotation rate. The mode number of the BZF is 1 for $Gamma lesssim 1$ and $2 Gamma$ for $Gamma$ = {1,2} independent of $Ra$ and $Pr$. The BZF is quite reminiscent of wall mode states in rotating convection.
We perform direct numerical simulations of rotating Rayleigh--Benard convection of fluids with low ($Pr=0.1$) and high ($Pr=5$) Prandtl numbers in a horizontally periodic layer with no-slip top and bottom boundaries. At both Prandtl numbers, we demonstrate the presence of an upscale transfer of kinetic energy that leads to the development of domain-filling vortical structures. Sufficiently strong buoyant forcing and rotation foster the quasi-two-dimensional turbulent state of the flow, despite the formation of plume-like vertical disturbances promoted by so-called Ekman pumping from the viscous boundary layer.
The shape of velocity and temperature profiles near the horizontal conducting plates in turbulent Rayleigh-B{e}nard convection are studied numerically and experimentally over the Rayleigh number range $10^8lesssim Ralesssim3times10^{11}$ and the Prandtl number range $0.7lesssim Prlesssim5.4$. The results show that both the temperature and velocity profiles well agree with the classical Prandtl-Blasius laminar boundary-layer profiles, if they are re-sampled in the respective dynamical reference frames that fluctuate with the instantaneous thermal and velocity boundary-layer thicknesses.
On the basis of (i) Particle Image Velocimetry data of a Turbulent Boundary Layer with large field of view and good spatial resolution and (ii) a mathematical relation between the energy spectrum and specifically modeled flow structures, we show that the scalings of the streamwise energy spectrum $E_{11}(k_{x})$ in a wavenumber range directly affected by the wall are determined by wall-attached eddies but are not given by the Townsend-Perry attached eddy models prediction of these spectra, at least at the Reynolds numbers $Re_{tau}$ considered here which are between $10^{3}$ and $10^{4}$. Instead, we find $E_{11}(k_{x}) sim k_{x}^{-1-p}$ where $p$ varies smoothly with distance to the wall from negative values in the buffer layer to positive values in the inertial layer. The exponent $p$ characterises the turbulence levels inside wall-attached streaky structures conditional on the length of these structures.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
