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

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh-Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh-Bénard convection. Our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

Elliptical Instability and Multiple-Roll Flow Modes of the Large-Scale Circulation in Confined Turbulent Rayleigh-Bénard Convection

The large-scale circulation (LSC) of fluid is one of the main concepts in turbulent thermal convection as it is known to be important in global heat and mass transport in the system. In turbulent Rayleigh-Bénard convection (RBC) in slender containers, the LSC is formed of several dynamically changing convective rolls that are stacked on top of each other. The present study reveals the following two important facts: (i) the mechanism which causes the twisting and breaking of a single-roll LSC into multiple rolls is the elliptical instability and (ii) the heat and momentum transport in RBC, represented by the Nusselt (Nu) and Reynolds (Re) numbers, is always stronger (weaker) for smaller (larger) number n of the rolls in the LSC structure. Direct numerical simulations support the findings for n=1,…,4 and the diameter-to-height aspect ratio of the cylindrical container Γ=1/5, the Prandtl number Pr=0.1 and Rayleigh number Ra=5×10^{5}. Thus, Nu and Re are, respectively, 2.5 and 1.5 times larger for a single-roll LSC (n=1) than for a LSC with n=4 rolls.

Robustness of heat transfer in confined inclined convection at high Prandtl number

We investigate the dependency of the magnitude of heat transfer in a convection cell as a function of its inclination by means of experiments and simulations. The study is performed with a working fluid of large Prandtl number, Pr≃480, and at Rayleigh numbers Ra≃10^{8} and Ra≃5×10^{8} in a quasi-two-dimensional rectangular cell with unit aspect ratio. By changing the inclination angle (β) of the convection cell, the character of the flow can be changed from moderately turbulent, for β=0^{∘}, to laminar and steady at β=90^{∘}. The global heat transfer is found to be insensitive to the drastic reduction of turbulent intensity, with maximal relative variations of the order of 20% at Ra≃10^{8} and 10% at Ra≃5×10^{8}, while the Reynolds number, based on the global root-mean-square velocity, is strongly affected with a decay of more than 85% occurring in the laminar regime. We show that the intensity of the heat flux in the turbulent regime can be only weakly enhanced by establishing a large-scale circulation flow by means of small inclinations. However, in the laminar regime the heat is transported solely by a slow large-scale circulation flow which exhibits large correlations between the velocity and temperature fields. For inclination angles close to the transition regime in-between the turbulentlike and laminar state, a quasiperiodic heat-flow bursting phenomenon is observed.

Temperature Fluctuation Profiles in Turbulent Thermal Convection: A Logarithmic Dependence versus a Power-Law Dependence

We report an experimental measurement of the rms temperature (σ_{T}) profiles in two regions inside a large aspect ratio (Γ=4.2) rectangular convection cell. It is found that, in the region where the boundary layer is sheared by a large-scale wind, σ_{T} has a power-law dependence on the vertical distance (z) from the plate, whereas in the region where plumes are abundant, σ_{T} has a logarithmic dependence on z. The power-law profile may be understood by balancing the inertia force and the viscous force in the equations of motion, and the logarithmic profile may be understood in terms of the balance between the buoyancy and the inertia forces. When normalized by a convective temperature scale, θ_{*}, the profiles of σ_{T} collapse onto a single curve for different values of the Rayleigh number. This shows that the convective temperature first proposed by Deardorff is the suitable temperature scale outside the thermal boundary layer for both logarithmic and power-law profiles. Our finding suggests a strong connection between plumes and the logarithmic rms temperature profile. The present Letter reveals that multiple force balance mechanisms can coexist in the bulk of highly turbulent flows.

Jump rope vortex in liquid metal convection

Understanding large-scale circulations (LSCs) in turbulent convective systems is important for the study of stars, planets, and in many industrial applications. The canonical model of the LSC is quasi-planar with additional horizontal sloshing and torsional modes [Brown E, Ahlers G (2009) J Fluid Mech 638:383-400; Funfschilling D, Ahlers G (2004) Phys Rev Lett 92:194502; Xi HD et al. (2009) Phys Rev Lett 102:044503; Zhou Q et al. (2009) J Fluid Mech 630:367-390]. Using liquid gallium as the working fluid, we show, via coupled laboratory-numerical experiments in tanks with aspect ratios greater than unity ([Formula: see text]), that the LSC takes instead the form of a "jump rope vortex," a strongly 3D mode that periodically orbits around the tank following a motion much like a jump rope on a playground. Further experiments show that this jump rope flow also exists in more viscous fluids such as water, albeit with a far smaller signal. Thus, this jump rope mode is an essential component of the turbulent convection that underlies our observations of natural systems.

Locally induced laminar convection in liquid nitrogen and silicone oils

We present an experimental study of a laminar convective phenomenon induced by a centimetric heater totally immersed in a liquid pool (Rayleigh number ranging from 10(4) to 10(7)). This local heating is observed to induce a laminar convection that differs from the classical Rayleigh-Bénard cells created by heating the whole bottom of the fluid: the convection pattern is no more periodic. In order to obtain a complete map of the velocity field, we use Particle Image Velocimetry technique. The vertical velocity between the counter-rotating convective cells is used as the relevant physical parameter to describe the phenomenon. The potential cooling applications of this problem lead us to choose liquid nitrogen as an experimental fluid. We thus compare the results obtained for various temperature gradients in liquid nitrogen with experiments performed at room temperature with silicone oils of various viscosities. The theoretical law for the maximal vertical velocity from classical Rayleigh-Bénard experiments is adapted to the specific geometry investigated by using a new definition for the characteristic wavelength. This length is studied and appears to be dependent on the liquid properties. We finally obtain a remarkable agreement between theory and experimental data.

Ability of a low-dimensional model to predict geometry-dependent dynamics of large-scale coherent structures in turbulence

We test the ability of a general low-dimensional model for turbulence to predict geometry-dependent dynamics of large-scale coherent structures, such as convection rolls. The model consists of stochastic ordinary differential equations, which are derived as a function of boundary geometry from the Navier-Stokes equations [Brown and Ahlers, Phys. Fluids 20, 075101 (2008); Phys. Fluids 20, 105105 (2008)]. We test the model using Rayleigh-Bénard convection experiments in a cubic container. The model predicts a mode in which the alignment of a convection roll stochastically crosses a potential barrier to switch between diagonals. We observe this mode with a measured switching rate within 30% of the prediction.

Boundary layer structure in turbulent Rayleigh–Bénard convection

The structure of the boundary layers in turbulent Rayleigh–Bénard convection is studied by means of three-dimensional direct numerical simulations. We consider convection in a cylindrical cell at aspect ratio one for Rayleigh numbers of $\mathit{Ra}= 3\ensuremath{\times} 1{0}^{9} $ and $3\ensuremath{\times} 1{0}^{10} $ at fixed Prandtl number $\mathit{Pr}= 0. 7$. Similar to the experimental results in the same setup and for the same Prandtl number, the structure of the laminar boundary layers of the velocity and temperature fields is found to deviate from the prediction of Prandtl–Blasius–Pohlhausen theory. Deviations decrease when a dynamical rescaling of the data with an instantaneously defined boundary layer thickness is performed and the analysis plane is aligned with the instantaneous direction of the large-scale circulation in the closed cell. Our numerical results demonstrate that important assumptions of existing classical laminar boundary layer theories for forced and natural convection are violated, such as the strict two-dimensionality of the dynamics or the steadiness of the fluid motion. The boundary layer dynamics consists of two essential local dynamical building blocks, a plume detachment and a post-plume phase. The former is associated with larger variations of the instantaneous thickness of velocity and temperature boundary layer and a fully three-dimensional local flow. The post-plume dynamics is connected with the large-scale circulation in the cell that penetrates the boundary region from above. The mean turbulence profiles taken in localized sections of the boundary layer for each dynamical phase are also compared with solutions of perturbation expansions of the boundary layer equations of forced or natural convection towards mixed convection. Our analysis of both boundary layers shows that the near-wall dynamics combines elements of forced Blasius-type and natural convection.

Chaotic gas turbine subject to augmented Lorenz equations

Inspired by the chaotic waterwheel invented by Malkus and Howard about 40 years ago, we have developed a gas turbine that randomly switches the sense of rotation between clockwise and counterclockwise. The nondimensionalized expressions for the equations of motion of our turbine are represented as a starlike network of many Lorenz subsystems sharing the angular velocity of the turbine rotor as the central node, referred to as augmented Lorenz equations. We show qualitative similarities between the statistical properties of the angular velocity of the turbine rotor and the velocity field of large-scale wind in turbulent Rayleigh-Bénard convection reported by Sreenivasan et al. [Phys. Rev. E 65, 056306 (2002)]. Our equations of motion achieve the random reversal of the turbine rotor through the stochastic resonance of the angular velocity in a double-well potential and the force applied by rapidly oscillating fields. These results suggest that the augmented Lorenz model is applicable as a dynamical model for the random reversal of turbulent large-scale wind through cessation.

Effect of weak rotation on large-scale circulation cessations in turbulent convection

We investigate the effect of weak rotation on the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection, using the theory for cessations in a low-dimensional stochastic model of the flow previously studied. We determine the cessation frequency of the LSC as a function of rotation, and calculate the statistics of the amplitude and azimuthal velocity fluctuations of the LSC as a function of the rotation rate for different Rayleigh numbers. Furthermore, we show that the tails of the reorientation PDF remain unchanged for rotating systems, while the distribution of the LSC amplitude and correspondingly the cessation frequency are strongly affected by rotation. Our results are in close agreement with experimental observations.

Rare fluctuations and large-scale circulation cessations in turbulent convection

In turbulent Rayleigh-Bénard convection, a large-scale circulation (LSC) develops in a nearly vertical plane and is maintained by rising and falling plumes detaching from the unstable thermal boundary layers. Rare but large fluctuations in the LSC amplitude can lead to extinction of the LSC (a cessation event), followed by the reemergence of another LSC with a different (random) azimuthal orientation. We extend previous models of the LSC dynamics to include momentum and thermal diffusion in the azimuthal plane and calculate the tails of the probability distributions of both the amplitude and azimuthal angle. Our analytical results are in very good agreement with the experimental data.

Coherent oscillations of turbulent Rayleigh-Bénard convection in a thin vertical disk

A well-defined oscillation is observed in the power spectrum of several fluctuating signals in turbulent Rayleigh-Bénard convection occurring in a thin vertical disk filled with water. The experiment reveals that the coherent oscillations are produced by periodic emission of thermal plumes, which gives rise to periodic pulses of forcing, resulting in a pulsed large-scale circulation in the thin cell. The experimental results agree well with the theoretical predictions made from two coupled nonlinear delayed equations.

Transition phenomena in unstably stratified turbulent flows

We study experimentally and theoretically the transition phenomena caused by external forcing from Rayleigh-Bénard convection with large-scale circulation (LSC) to the limiting regime of unstably stratified turbulent flow without LSC, where the temperature field behaves like a passive scalar. In the experiments we use the Rayleigh-Bénard apparatus with an additional source of turbulence produced by two oscillating grids located near the sidewalls of the chamber. When the frequency of the grid oscillations is larger than 2 Hz, the LSC in turbulent convection is destroyed, and the destruction of the LSC is accompanied by a strong change of the mean temperature distribution. However, in all regimes of the unstably stratified turbulent flow the ratio [(ℓ{x}∇{x}T)²+(ℓ{y}∇{y}T)² + (ℓ{z}∇{z}T)²]/<θ²> varies slightly (even in the range of parameters where the behavior of the temperature field is different from that of the passive scalar). Here ℓ{i} are the integral scales of turbulence along the x,y,z directions, and T and θ are the mean and fluctuating parts of the fluid temperature. At all frequencies of the grid oscillations we have detected long-term nonlinear oscillations of the mean temperature. The theoretical predictions based on the budget equations for turbulent kinetic energy, turbulent temperature fluctuations, and turbulent heat flux, are in agreement with the experimental results.

Structure of viscous boundary layers in turbulent Rayleigh-Bénard convection

Highly resolved profiles of the mean velocity in turbulent Rayleigh-Bénard convection in air are presented and discussed. The present work extends our recently performed experiments at constant aspect ratio [Phys. Rev. Lett. 99, 234504 (2007)] to variable aspect ratios. The experiments cover a range of Rayleigh numbers 10(9)<Ra<10(12) and aspect ratios 1.13<Gamma<11.3 whereas the Prandtl number is fixed at Pr=0.7. The major finding of the present work is that the profiles of the mean horizontal velocity and its fluctuations are virtually invariant against the variation in Ra or Gamma if the wall distance is scaled by the displacement thickness of the boundary layer. Furthermore we have studied typical length scales of the boundary layer and their scaling with Ra and Gamma. Regarding a potential transition of the heat transport toward the ultimate regime we found that the boundary layer Reynolds number remains below Redelta=250 which is significantly lower than the predicted limit.

Enhanced heat transport by turbulent two-phase Rayleigh-Bénard convection

We report measurements of turbulent heat transport in samples of ethane (C2H6) heated from below while the applied temperature difference DeltaT straddled the liquid-vapor coexistence curve T(phi)(P). When the sample top temperature T(t) decreased below T(phi), droplet condensation occurred and the latent heat of vaporization H provided an additional heat-transport mechanism. The effective conductivity lambda(eff) increased linearly with decreasing T(t), and reached a maximum value lambda(eff)(*) that was an order of magnitude larger than the single-phase lambda(eff). As P approached the critical pressure, lambda(eff)(*) increased dramatically even though H vanished. We attribute this phenomenon to an enhanced droplet-nucleation rate as the critical point is approached.

Origin of the temperature oscillation in turbulent thermal convection

We report an experimental study of the three-dimensional spatial structure of the low-frequency temperature oscillations in a cylindrical Rayleigh-Bénard convection cell. Through simultaneous multipoint temperature measurements it is found that, contrary to the popular scenario, thermal plumes are emitted neither periodically nor alternately, but randomly and continuously, from the top and bottom plates. We further identify a new flow mode-the sloshing mode of the large-scale circulation (LSC). This sloshing mode, together with the torsional mode of the LSC, are found to be the origin of the oscillation of the temperature field.

Azimuthal motion, reorientation, cessation, and reversal of the large-scale circulation in turbulent thermal convection: a comparative study in aspect ratio one and one-half geometries

We report a systematic experimental study of the orientation and the flow strength of the large-scale circulation (LSC) in water-filled cylindrical Rayleigh-Bénard convection cells with aspect ratios 2.3, 1, and 0.5 by both direct velocity measurement and the indirect multithermal-probe measurement. Unlike its weak effect in the system's global heat transport, the aspect ratio Gamma is found to play an important role in the dynamics of the azimuthal motion of the LSC. It is found that in larger Gamma geometries the azimuthal motion of the LSC's vertical plane is confined in smaller azimuthal region than that in smaller Gamma geometries. The twisting motion between top and bottom parts of the LSC observed in the Gamma=1 geometry is found to be absent in the Gamma=1/2 case. It is found that in the Gamma=1/2 geometry the orientational change mid R:Deltavarphimid R: through a reorientation has an exponential distribution, in contrast to the power-law distribution for the Gamma=1 case. Despite the difference in orientational change, the occurrence of the reorientations is a Poisson process in both geometries. Using the conditional average of the time interval between adjacent cessations or reversals on the rebound flow strength, we demonstrate the possibility to empirically predict when the next cessation or reversal will most likely occur if the rebound flow strength of the preceding cessation or reversal is given.

Non-Oberbeck-Boussinesq effects in turbulent thermal convection in ethane close to the critical point

As shown in earlier work [Ahlers, J. Fluid Mech. 569, 409 (2006)], non-Oberbeck-Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Bénard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point, the origin of non-Oberbeck-Boussinesq corrections is very different, as will be shown in the present paper. Namely, the main origin of NOB corrections then lies in the strong temperature dependence of the isobaric thermal expansion coefficient beta(T). More precisely, it is the nonlinear T dependence of the density rho(T) in the buoyancy force that causes another type of NOB effect. We demonstrate this through a combination of experimental, numerical, and theoretical work, the last in the framework of the extended Prandtl-Blasius boundary-layer theory developed by Ahlers as cited above. The theory comes to its limits if the temperature dependence of the thermal expension coefficient beta(T) is significant. The measurements reported here cover the ranges 2.1<or similar to Pr<or similar to 3.9 and 5x10(9)<or similar to Ra<or similar to 2x10(12) and are for cylindrical samples of aspect ratios 1.0 and 0.5.

Numerical and experimental investigation of structure-function scaling in turbulent Rayleigh-Bénard convection

Direct numerical simulation and stereoscopic particle image velocimetry of turbulent convection are used to gather spatial data for the calculation of structure functions. We wish to add to the ongoing discussion in the literature whether temperature acts as an active or passive scalar in turbulent convection, with consequences for structure-function scaling. The simulation results show direct confirmation of the scalings derived by Bolgiano and Obukhov for turbulence with an active scalar for both velocity and temperature statistics. The active-scalar range shifts to larger scales when the forcing parameter (Rayleigh number) is increased. Furthermore, a close inspection of local turbulent length scales (Kolmogorov and Bolgiano lengths) confirms conjectures from earlier studies that the oft-used global averages are not suited for the interpretation of structure functions. In the experiment, a characterization of the domain-filling large-scale circulation of confined convection is carried out for comparison with other studies. The measured velocity fields are also used to calculate velocity structure functions, further confirming the Bolgiano-Obukhov scalings when interpreted with the local turbulent length scales found in the simulations. An extended self-similarity analysis shows that the relative scalings are different for the Kolmogorov and Bolgiano-Obukhov regimes.

Mean velocity profile in confined turbulent convection

In this Letter, highly resolved measurements of the horizontal velocity inside the boundary layer of turbulent Rayleigh-Bénard convection are reported. They were performed in a cylindrical box with an aspect ratio Gamma=1.13 which was filled with air with a Prandtl number Pr=0.7. The horizontal velocity was measured along the central axis close to the cooling plate in a range of Rayleigh numbers between Ra=10;{11} and Ra=10;{12} using a two-dimensional laser Doppler velocimeter. We demonstrate that the profile of the mean velocity strongly differs from that of classical shear flows like the Blasius shape of a laminar flat plate boundary layer or a turbulent logarithmic velocity profile with standard coefficients.

Measured oscillations of the velocity and temperature fields in turbulent Rayleigh-Bénard convection in a rectangular cell

Temperature and velocity oscillations have been found in a rectangular Rayleigh-Bénard convection cell, in which one large-scale convection roll exists. At Rayleigh number Ra=8.9x10(11) and Prandtl number Pr=4, temperature oscillations can be observed in most parts of the system and the oscillation period remains almost constant, tT=74+/-2 s. Velocity oscillation can only be found in its horizontal component vy (perpendicular to the large-scale circulation plane) near the cell sidewall, its oscillation period is also constant, tv=65+/-2 s, at these positions. Temperature and velocity oscillations have different Ra dependences, which are, respectively, indicated by the Péclect number PeT=0.55Ra0.47 and Pev=0.28Ra0.50. In comparison to the case of a cylindrical cell, we find that velocity oscillations are affected by the system geometry.

Large-scale circulation model for turbulent Rayleigh-Bénard convection

A model for the large-scale-circulation (LSC) of turbulent Rayleigh-Bénard convection in cylindrical samples is presented. It consists of two physically motivated stochastic ordinary differential equations, one each for the strength and the azimuthal orientation of the LSC. Stochastic forces represent phenomenologically the influence of turbulent fluctuations. Consistent with measurements, the model yields an azimuthally meandering LSC with occasional rotations, and with more rare cessations. As in experiment, cessations have a uniform distribution of LSC orientation changes.

Morphological evolution of thermal plumes in turbulent Rayleigh-Bénard convection

An experimental study of the morphological evolution of thermal plumes in turbulent thermal convection is presented. Individual sheetlike plumes are extracted and their area, circumference, and "heat content" are found to all exhibit log-normal distributions. As the sheetlike plumes move across the plate they collide and convolute into spiraling swirls. These swirls then spiral away from the plates to become mushroomlike plumes which are accompanied by strong vertical vorticity. The measured profiles of plume numbers and of vertical vorticity quantify the morphological transition of sheetlike plumes to mushroomlike ones and the mixing and merging or clustering of mushroomlike plumes. The fluctuating vorticity is found to have the same exponential distribution and scaling behavior as the fluctuating temperature.

Breakdown of wind in turbulent thermal convection

We report experiments on turbulent Rayleigh-Bénard convection of air in a cylindrical large-scale facility with a diameter of 7 meters and Rayleigh numbers up to Ra approximately 10(12). The facility is used to explore the structure of the large-scale circulation for continuously varying aspect ratios between Gamma approximately 1 and Gamma approximately 10. Using autocorrelation functions derived from high-resolution time series of temperature and velocity near the cooling plate we demonstrate that the well-known single-roll structure (often called "wind") breaks down when the aspect ratio increases beyond the critical value Gamma(1) = 1.68+/-0.22. We further show that at Gamma(2) = 3.66+/-0.46 a second transition from an oscillatory two-roll structure to an unstable multi-roll structure takes place. The value of Gamma(2) represents a lower bound for the aspect ratio that is necessary to reach homogeneous convection--a turbulent state that is free from the influence of lateral walls.

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.

Oscillating large-scale circulation in turbulent Rayleigh-Bénard convection

We report on the dynamics and structure of the turbulent velocity field in a high-Rayleigh-number (Ra = 5.9 x 10(8))thermal convection cell with an aspect ratio of 4. Spectral density functions (measured with laser Doppler velocimetry) indicated the existence of a large-scale periodic component. The long-time mean flow field (measured with particle image velocimetry) revealed that the large-scale circulation in the aspect-ratio-4 cell consists of two corotating rolls. The periodicity in the flow could be traced back to the alternating growth and decay of these rolls.

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.

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.

Measurements of the local convective heat flux in turbulent Rayleigh-Bénard convection

A systematic study of the local convective heat transport in turbulent thermal convection is carried out in small-aspect-ratio cells filled with water. The local convective heat flux is obtained from the simultaneous velocity and temperature measurements over varying Rayleigh numbers and spatial positions across the entire convection cell. Large fluctuations of the local convective heat flux are found mainly in the vertical direction and they are determined primarily by the thermal plumes in the system. The experiment reveals the spatial distribution of the local convective heat flux in a closed cell and thus settles a long-debated issue on how heat is transported in small-aspect-ratio cells.