Prandtl and Rayleigh number dependence of the Reynolds number in turbulent thermal convection

Non-oberbeck-boussinesq effects in gaseous Rayleigh-Bénard convection

Non-Oberbeck-Boussinesq (NOB) effects are measured experimentally and calculated theoretically for strongly turbulent Rayleigh-Bénard convection of ethane gas under pressure where the material properties strongly depend on the temperature. Relative to the Oberbeck-Boussinesq case we find a decrease of the central temperature as compared to the arithmetic mean of the top- and bottom-plate temperature and an increase of the Nusselt number. Both effects are of opposite sign and greater magnitude than those for NOB convection in liquids like water.

Scaling in magnetohydrodynamic convection at high Rayleigh number

The theory of Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)] is extended to include the effect of a magnetic field on convection of an electrically conducting fluid. Different scaling laws are obtained depending on whether the bulk or the boundary layers make the major contribution to the dissipation. Scalings are obtained for both weak and strong magnetic fields. The predictions are shown to be in better agreement with experimental data than earlier theoretical models.

Azimuthal motion of the mean wind in turbulent thermal convection

We present an experimental study of the azimuthal motion of the mean wind in turbulent thermal convection. The experiments were conducted with cylindrical convection cells of unity aspect ratio and over the range of the Rayleigh number from 1 x 10(9) to 1 x 10(10). The azimuthal angle of the circulation plane of the mean wind was measured using both the particle image velocimetry and flow-visualization techniques. It is found that the azimuthal motion consists of erratic fluctuations and a time-periodic oscillation. The orientation of the wind is found to be "locked," i.e., it fluctuates about a preferred direction most of the time with all other orientations appearing as "transient states," and large excursions of the azimuthal angle often result in a net rotation which takes the wind back to the preferred orientation. The rate of erratic rotation of the circulation plane is found to have a strong dependence on Ra. Our result suggests that the oscillatory motion of the wind in its vertically oriented circulation plane and the orientational oscillation of the circulation plane itself have the same dynamic origin.

Bifurcation analysis of multiple steady flow patterns for Rayleigh-Bénard convection in a cubical cavity at Pr = 130

The bifurcation diagram of steady convective flow patterns inside a cubical cavity with adiabatic lateral walls heated from below and filled with silicone oil (Pr = 130) was determined for values of the Rayleigh number (Ra) up to 1.5 x 10(5). A continuation procedure based on the Galerkin spectral method was used to determine the steady convective solutions as a function of Ra. Bifurcations leading to either new steady or time-dependent solutions were identified and new steady solution branches were also continued. A total of fifteen steady solutions were tracked and the stability analysis predicted that six flow patterns were stable and that two, three, or even four of these patterns coexisted over certain ranges of Ra in the studied domain. Predicted flow patterns and transitions are in agreement with flow visualizations previously reported in the literature. The variation of the Nusselt number (Nu) as a function of Pr was investigated for three of the stable flow patterns identified: a x or y roll, a diagonal oriented roll and a pattern formed by four connected half rolls. It was found that whereas the Nusselt changes within the region 0.71 < or = Pr < or = 10 it tends to an asymptotic value with increasing Pr.

Scaling of the Reynolds number in turbulent thermal convection

A riddle in turbulent thermal convection is the apparent dispersion from 0.42 to 0.5 in the value of the scaling exponent of experimentally measured Reynolds number Re approximately Ragamma, where Ra is the Rayleigh number. The measured Re may be divided into two groups: one based on the circulation frequency of the mean wind and the other based on a directly measured velocity. With new experimental results we show that in frequency measurements the dispersion in gamma is a result of the evolution in the circulation path of the wind, and that in the velocity measurements it is caused by the inclusion of a counterflow in the mean velocity. When these factors are properly accounted for both groups give gamma=0.5, which may imply that a single mechanism is driving the flow for both low and high values of Ra.

Wind reversals in turbulent Rayleigh-Bénard convection

The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in turbulent Rayleigh-Bénard convection is theoretically analyzed. The force and thermal balance on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers in the range 10(-2) < or = Pr < or = 10(3) and 10(7) < or = Ra < or = 10(12), the model has the following features: (i) chaotic reversals may be exhibited at Ra > or = 10(7); (ii) the Reynolds number based on the root mean square velocity scales as Re(rms) approximately Ra([0.41...0.47]) (depending on Pr), and as Re(rms) approximately Pr(-[0.66...0.76]) (depending on Ra); and (iii) the mean reversal frequency follows an effective scaling law omega/(nu L(-2)) approximately Pr(-(0.64 +/- 0.01))Ra(0.44 +/- 0.01). The phase diagram of the model is sketched, and the observed transitions are discussed.

Three-dimensional flow structures and dynamics of turbulent thermal convection in a cylindrical cell

The technique of particle image velocimetry is used to study the velocity field of turbulent Rayleigh-Bénard convection in an aspect-ratio-1 cylindrical cell filled with water. By measuring the two-dimensional (2D) velocity vector map in different vertical cross sections of the cell, we investigate the 3D structures and dynamics of turbulent thermal convection. The experiment reveals how thermal plumes synchronize their emissions and organize their motions spatially between the top and bottom plates, leading to an oscillatory motion in the bulk region of the fluid with a period equal to twice the plume's cell-crossing time. From the measured instantaneous velocity vector map, we find the phase relationship between the velocity components along different directions and at different positions in a 2D plane. These phase relations illustrate how the convecting fluid in different regions of the cell interact with each other and generate a synchronized and coherent motion in a closed system.

Fluctuations of temperature gradients in turbulent thermal convection

Broad theoretical arguments are proposed to show, formally, that the magnitude G of the temperature gradients in turbulent thermal convection at high Rayleigh numbers obeys the same advection-diffusion equation that governs the temperature fluctuation T , except that the velocity field in the new equation is substantially smoothed. This smoothed field leads to a -1 scaling of the spectrum of G in the same range of scales for which the spectral exponent of T lies between -7/5 and -5/3 . This result is confirmed by measurements in a confined container with cryogenic helium gas as the working fluid for Rayleigh number Ra=1.5x 10(11) . Also confirmed is the logarithmic form of the autocorrelation function of G . The anomalous scaling of dissipation-like quantities of T and G are identical in the inertial range, showing that the analogy between the two fields is quite deep.

Plume motion and large-scale circulation in a cylindrical Rayleigh-Bénard cell

We used the time correlation of shadowgraph images to determine the angle Theta of the horizontal component of the plume velocity above (below) the center of the bottom (top) plate of a cylindrical Rayleigh-Bénard cell of aspect ratio Gamma identical with D/L=1 (D is the diameter and L approximately 87 mm is the height) in the Rayleigh-number range 7 x 10(7)</=R</=3 x 10(9) for a Prandtl number sigma=6. We expect that Theta gives the direction of the large-scale circulation. It oscillates time periodically. Near the top and bottom plates Theta(t) has the same frequency but is anticorrelated.

Particle image velocimetry measurement of the velocity field in turbulent thermal convection

The spatial structure of the velocity field in turbulent Rayleigh-Bénard convection in water has been measured using the particle image velocimetry technique, with the Rayleigh number Ra varying from 9 x 10(8) to 9 x 10(11) and the Prandtl number remaining approximately constant (Pr approximately 4). The study provides a direct confirmation that a rotatory mean wind indeed persists for the highest value of Ra reached in the experiment. The measurement reveals that the mean flow in the central region of the convection cell is of the shape of a coherent elliptical rotating core for Ra below 1 x 10(10). Above this Ra, the orientation of the elliptical core changes by a 90 degrees angle and an inner core rotating at a lower rate inside the original bulk core emerges. It is further found that the rotation frequencies of the inner core and the outer shell have distinct scalings with Ra; the scaling exponent for the outer-shell is 0.5 and it is 0.4 for the inner core. From the measured rms and skewness distributions of the velocity field, we find that velocity fluctuations at the cell center are neither homogeneous nor isotropic. The turbulent energy production fields further reveal that the mean wind is not driven by turbulent fluctuations associated with Reynolds stress.

Smooth and rough boundaries in turbulent Taylor-Couette flow

We examine the torque required to drive the smooth or rough cylinders in turbulent Taylor-Couette flow. With rough inner and outer walls the scaling of the dimensionless torque G is found to be consistent with pure Kolmogorov scaling G approximately Re2. The results are interpreted within the Grossmann-Lohse theory for the relative role of the energy dissipation rates in the boundary layers and in the bulk; as the boundary layers are destroyed through the wall roughness, the torque scaling is due only to the bulk contribution. For the case of one rough and one smooth wall, we find that the smooth cylinder dominates the dissipation rate scaling, i.e., there are corrections to Kolmogorov scaling. A simple model based on an analogy to electrical circuits is advanced as a phenomenological organization of the observed relative drag functional forms. This model leads to a qualitative prediction for the mean velocity profile within the bulk of the flow.

Phenomenology of Rayleigh-Taylor turbulence

I analyze the advanced mixing regime of the Rayleigh-Taylor incompressible turbulence in the small Atwood number Boussinesq approximation. The prime focus of my phenomenological approach is to resolve the temporal behavior and the small-scale spatial correlations of velocity and temperature fields inside the mixing zone, which grows as proportional, variant t(2). I show that the "5/3"-Kolmogorov scenario for velocity and temperature spectra is realized in three spatial dimensions with the viscous and dissipative scales decreasing in time, proportional, variant t(-1/4). The Bolgiano-Obukhov scenario is shown to be valid in two dimensions with the viscous and dissipative scales growing, proportional, variant t(1/8).

Nusselt number measurements for turbulent Rayleigh-Bénard convection

We present high-precision measurements of the Nusselt number N as a function of the Rayleigh number R for a cylindrical sample of water (Prandtl number sigma=4.4) of height L approximately equal to 50 cm and aspect ratio Gamma identical with D/L approximately equal to 1 (D is the diameter) for 3 x 10(9)< or =R< or =6 x 10(10). For R approximately 3 x 10(9) the data are consistent with existing results for acetone (sigma=4.0, R< or =3 x 10(9)). There the measurements are also consistent with a model by Grossmann and Lohse (GL). As R increases, the measurements fall below the GL prediction. Near R=6 x 10(10) the prediction is 8% above the data.

Ultimate state of thermal convection

The ultimate regime of thermal convection, the so-called Kraichnan regime [R. H. Kraichnan, Phys. Fluids 5, 1374 (1962)]], hitherto has been elusive. Here numerical evidence for that regime is presented by performing simulations of the bulk of turbulence only, eliminating the thermal and kinetic boundary layers and replacing them with periodic boundary conditions.

Temperature oscillations in turbulent Rayleigh-Bénard convection

A systematic study of temperature oscillations in turbulent thermal convection was carried out in two aspect-ratio-one convection cells filled with water. Temperature correlation functions and local velocity fluctuations were measured over varying Rayleigh numbers and spatial positions across the entire cell. These measurements fully characterize the spatial structure of the temperature oscillation and reveal the mixing and emission dynamics of the thermal plumes near the conducting surface. A sharp transition from a random chaotic state to a correlated turbulent state of finite coherence time is found when the Rayleigh number becomes larger than a critical value Ra(c) approximately equal 5 x 10(7). Above Ra(c) the measured temperature correlation functions show a well-defined oscillation with a finite coherence time. The oscillation period is found to be twice as large as the cell crossing time. The experiment demonstrates how the thermal plumes in a closed cell organize themselves both in space and time and generate coherent oscillations in a turbulent environment.