Energy-conserving model of hexagonal pattern in Rayleigh-Bénard convection

We have constructed an energy-conserving sixteen mode dynamical system to model hexagonal pattern in Rayleigh-Bénard convection of Boussinesq fluids with symmetric stress-free thermally conducting boundaries. The model shows stable roll pattern at the onset of convection. Hexagon is found to appear in the system via sausage and (or) stationary rhombus patterns. Both up and down hexagons arise periodically or chaotically with roll, sausage and rhombus patterns. Hexagonal patterns exist for all values of the Prandtl number, 1 ≤ Pr ≤ 5 explored here. However the pattern is more prominent for small Pr and k < kc , where k denotes the wave number. The plot of Nusselt number matches with previous theoretical result. In dissipationless limit, the total energy and the unavailable energy are constants though the kinetic energy, the potential energy and the available energy vary with time. The derived model does not diverge for large values of Rayleigh number Ra.

Steady Rayleigh–Bénard convection between no-slip boundaries

The central open question about Rayleigh–Bénard convection – buoyancy-driven flow in a fluid layer heated from below and cooled from above – is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $Ra\to \infty$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed ‘classical’ $Nu \sim Ra^{1/3}$ or ‘ultimate’ $Nu \sim Ra^{1/2}$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $10^{14}$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$ , we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $Ra\to \infty$ , it cannot achieve the ultimate scaling.

Topographically Controlled Marangoni-Rayleigh-Bénard Convection in Liquid Metals

Convection induced in a layer of liquid with a top free surface by a distribution of heating elements at the bottom can be seen as a variant of standard Marangoni–Rayleigh–Bénard Convection where in place of a flat boundary at constant temperature delimiting the system from below, the underlying thermal inhomogeneity reflects the existence of a topography. In the present work, this problem is investigated numerically through solution of the governing equations for mass, momentum and energy in their complete, three-dimensional time-dependent and non-linear form. Emphasis is given to a class of liquids for which thermal diffusion is expected to dominate over viscous effects (liquid metals). Fixing the Rayleigh and Marangoni number to 104 and 5 × 103, respectively, the sensitivity of the problem to the geometrical, kinematic and thermal boundary conditions is investigated parametrically by changing: the number and spacing of heating elements, their vertical extension, the nature of the lateral boundary (solid walls or periodic boundary) and the thermal behavior of the portions of bottom wall between adjoining elements (assumed to be either adiabatic or at the same temperature of the hot blocks). It is shown that, like the parent phenomena, this type of thermal flow is extremely sensitive to the specific conditions considered. The topography can be used to exert a control on the emerging flow in terms of temporal response and patterning behavior.

Jump rope vortex flow in liquid metal Rayleigh–Bénard convection in a cuboid container of aspect ratio

We study the topology and the temporal dynamics of turbulent Rayleigh–Bénard convection in a liquid metal with a Prandtl number of 0.03 located inside a box with a square base area and an aspect ratio of $\varGamma = 5$ . Experiments and numerical simulations are focused on the Rayleigh number range $6.7 \times 10^{4} \leqslant Ra \leqslant 3.5 \times 10^{5}$ , where a new cellular flow regime has been reported previously (Akashi et al., Phys. Rev. Fluids, vol. 4, 2019, 033501). This flow structure shows symmetries with respect to the vertical planes crossing at the centre of the container. The dynamic behaviour is dominated by strong three-dimensional oscillations with a period length that corresponds to the turnover time. Our analysis reveals that the flow structure in the $\varGamma = 5$ box corresponds in key features to the jump rope vortex structure, which has recently been discovered in a $\varGamma = 2$ cylinder (Vogt et al., Proc. Natl Acad. Sci. USA, vol. 115, 2018, pp. 12674–12679). While in the $\varGamma = 2$ cylinder a single jump rope vortex occurs, the coexistence of four recirculating swirls is detected in this study. Their approach to the lid or the bottom of the convection box causes a temporal deceleration of both the horizontal velocity at the respective boundary and the vertical velocity in the bulk, which in turn is reflected in Nusselt number oscillations. The cellular flow regime shows remarkable similarities to properties commonly attributed to turbulent superstructures.