Heteroclinic behavior in rotating Rayleigh-Bénard convection

We investigate numerically the appearance of heteroclinic behavior in a three-dimensional, buoyancy-driven fluid layer with stress-free top and bottom boundaries, a square horizontal periodicity with a small aspect ratio, and rotation at low to moderate rates about a vertical axis. The Prandtl number is 6.8. If the rotation is not too slow, the skewed-varicose instability leads from stationary rolls to a stationary mixed-mode solution, which in turn loses stability to a heteroclinic cycle formed by unstable roll states and connections between them. The unstable eigenvectors of these roll states are also of the skewed-varicose or mixed-mode type and in some parameter regions skewed-varicose like shearing oscillations as well as square patterns are involved in the cycle. Always present weak noise leads to irregular horizontal translations of the convection pattern and makes the dynamics chaotic, which is verified by calculating Lyapunov exponents. In the nonrotating case, the primary rolls lose, depending on the aspect ratio, stability to traveling waves or a stationary square pattern. We also study the symmetries of the solutions at the intermittent fixed points in the heteroclinic cycle. Received 10 June 1999     

Thermal-noise effect on the transition to Rayleigh-Bénard convection

We report measurements of fluctuation and roll patterns near the transition to Rayleigh-Bénard convection which are consistent with a fluctuation-induced first-order transition, as predicted by Swift and Hohenberg. Above onset, we find convection rolls with noise-induced fluctuations, time-dependent amplitude modulation and roll undulation, and homogeneous dislocation nucleation.     

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection

In turbulent thermal convection in cylindrical samples with an aspect ratio Γ≡D/L (D is the diameter and L the height), the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number 1/Ro(c), with 1/Ro(c)∝1/Γ. We present a Ginzburg-Landau-like model that explains the existence of a bifurcation at finite 1/Ro(c) as a finite-size effect. The model yields the proportionality between 1/Ro(c) and 1/Γ and is consistent with several other measured or computed system properties.     

Rayleigh-Bénard convection in large-aspect-ratio domains

The coarsening and wave number selection of striped states growing from random initial conditions are studied in a nonrelaxational, spatially extended, and far-from-equilibrium system by performing large-scale numerical simulations of Rayleigh-Bénard convection in a large-aspect-ratio cylindrical domain with experimentally realistic boundaries. We find evidence that various measures of the coarsening dynamics scale in time with different power-law exponents, indicating that multiple length scales are required in describing the time dependent pattern evolution. The translational correlation length scales with time as t0.12, the orientational correlation length scales as t0.54, and the density of defects scale as t(-0.45). The final pattern evolves toward the wave number where isolated dislocations become motionless, suggesting a possible wave number selection mechanism for large-aspect-ratio convection.     

Dislocation dynamics in Rayleigh-Bénard convection

Theoretical results on the dynamics of dislocations in Rayleigh-Bénard convection are reported both for a Swift-Hohenberg model and the Oberbeck-Boussinesq equations. For intermediate Prandtl numbers the motion of dislocations is found to be driven by the superposition of two independent contributions: (i) the Peach-Koehler force and (ii) an advection force on the dislocation core by its self-generated mean flow. Their competition allows to explain the experimentally observed bound dislocation pairs.     

Large scale structures in Rayleigh-Bénard convection at high Rayleigh numbers

Direct numerical simulations of Rayleigh-Bénard convection in a plane layer with periodic boundary conditions at Rayleigh numbers up to 10(7) show that flow structures can be objectively classified as large or small scale structures because of a gap in spatial spectra. The typical size of the large scale structures does not always vary monotonically as a function of the Rayleigh number but broadly increases with increasing Rayleigh number. A mean flow (whose average over horizontal planes differs from zero) is also excited but is weak in comparison with the large scale structures. The large scale circulation observed in experiments should therefore be a manifestation of the large scale structures identified here.     

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.     

New perspectives in turbulent Rayleigh-Bénard convection

Recent experimental, numerical and theoretical advances in turbulent Rayleigh-Bénard convection are presented. Particular emphasis is given to the physics and structure of the thermal and velocity boundary layers which play a key role for the better understanding of the turbulent transport of heat and momentum in convection at high and very high Rayleigh numbers. We also discuss important extensions of Rayleigh-Bénard convection such as non-Oberbeck-Boussinesq effects and convection with phase changes.     

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.     

Chaotic travelling rolls in Rayleigh-Bénard convection

In this paper we investigate two-dimensional (2D) Rayleigh-Bénard convection using direct numerical simulation in Boussinesq fluids with Prandtl number P = 6.8 confined between thermally conducting plates. We show through the simulation that in a small range of reduced Rayleigh number r (770 < r < 890) the 2D rolls move chaotically in a direction normal to the roll axis. The lateral shift of the rolls may lead to a global flow reversal of the convective motion. The chaotic travelling rolls are observed in simulations with free-slip as well as no-slip boundary conditions on the velocity field. We show that the travelling rolls and the flow reversal are due to an interplay between the real and imaginary parts of the critical modes.     

New perspectives in turbulent Rayleigh-Bénard convection

Recent experimental, numerical and theoretical advances in turbulent Rayleigh-Bénard convection are presented. Particular emphasis is given to the physics and structure of the thermal and velocity boundary layers which play a key role for the better understanding of the turbulent transport of heat and momentum in convection at high and very high Rayleigh numbers. We also discuss important extensions of Rayleigh-Bénard convection such as non-Oberbeck-Boussinesq effects and convection with phase changes.     

Transition to the ultimate state of turbulent Rayleigh-Bénard convection

Measurements of the Nusselt number Nu and of a Reynolds number Re(eff) for Rayleigh-Bénard convection (RBC) over the Rayleigh-number range 10(12)?Ra?10(15) and for Prandtl numbers Pr near 0.8 are presented. The aspect ratio Γ≡D/L of a cylindrical sample was 0.50. For Ra?10(13) the data yielded Nu∝Ra(γ(eff)) with γ(eff)?0.31 and Re(eff)∝Ra(ζ(eff)) with ζ(eff)?0.43, consistent with classical turbulent RBC. After a transition region for 10(13)?Ra?5×10(14), where multistability occurred, we found γ(eff)?0.38 and ζ(eff)=ζ?0.50, in agreement with the results of Grossmann and Lohse for the large-Ra asymptotic state with turbulent boundary layers which was first predicted by Kraichnan.     

Reorientation of the large-scale circulation in turbulent Rayleigh-Bénard convection

We present measurements of the orientation theta0(t) of the large-scale circulation (LSC) of turbulent Rayleigh-Bénard convection in cylindrical cells of aspect ratio 1. Theta0(t) undergoes irregular reorientations. It contains reorientation events by rotation through angles delta theta with a monotonically decreasing probability distribution p(delta theta), and by cessations (where the LSC stops temporarily) with a uniform p(delta theta). Reorientations have Poissonian statistics in time. The amplitude of the LSC and the magnitude of the azimuthal rotation rate have a negative correlation.