Nonlinear Bénard convection with rotation

1973 ◽  
Vol 57 (3) ◽  
pp. 433-458 ◽  
Author(s):  
J. C. Morgan

The equations for nonlinear Bénard convection with rotation for a layer of fluid, thickness d, are derived using the Glansdorff & Prigogine (1964) evolutionary criterion as used by Roberts (1966) in his paper on non-rotational Bénard convection. The parameters of the problem in this case are the Rayleigh number R = αgΔθd/vK, the Taylor number T = 4d4Ω23/v2 and the Prandtl number Pr = v/K, where α is the coefficient of volume expansion, g the acceleration due to gravity, Δθ the temperature difference between the horizontal surfaces, v the kinematic viscosity, K the thermal diffusivity and Ω3 the rotation rate about the vertical direction. The asymptotic solution for two-dimensional cells (rolls) is investigated for large Rayleigh numbers and large Taylor numbers. For rolls the convection equations are found t o be independent of the Prandtl number. However, the solutions depend upon the Prandtl number for another reason. The rotational problem differs from the non-rotational one in that the Rayleigh number and the horizontal wavenumber a of the convection are now functions of the Taylor number. These are taken to be R ∼ ρTα′ and α ∼ ATβ, where α′ and β are positive numbers. Thermal layers develop as R becomes large with ρ or T becoming large. The order in which ρ and T are allowed to increase is important since the horizontal wavenumber a also increases with T and the convection equations can be reduced in this case. A liquid of large Prandtl number such as water has v [Gt ] K. Since R ∼ O (1/vK) and T ∼ O(1/v2), ρ will be greater than T for a given (large) Δθ and Ω3. Similarly, for a liquid of small Prandtl number such as mercury v [Lt ] K, and T is greater than ρ for a given Δθ and Ω3. For rigid-rigid horizontal boundaries with ρ large and then T large the ρ thermal layer has the same structure as for the non-rotating problem. As T → ∞ three types of thermal layers are possible: a linear Ekman layer, a nonlinear Ekman layer and a Blasius-type thermal layer. When the horizontal boundaries are both free the ρ thermal layer is again of the same structure as for non-rotating BBnard convection. As T → ∞ a nonlinear Ekman layer and a Blasius-type thermal layer are possible.When T is large and then ρ made large the differential equations governing the convection are reduced from eighth order to sixth order owing to a becoming large as T → ∞. There are Ekman layers as T → ∞, when the horizontal boundaries are both rigid. The ρ thermal layers now have a different structure from the non-rotating problem for both rigid-rigid and free-free horizontal boundaries. The equation for small amplitude convection near to the marginal case is derived and the solution for free-free horizontal boundaries is obtained.

A recent study by Cross et al . (1980) has described a class of finite-amplitude phase-winding solutions of the problem of two-dimensional Rayleigh-Bénard convection in a shallow fluid layer of aspect ratio 2 L (≫ 1) confined laterally by rigid side-walls. These solutions arise at Rayleigh numbers R = R 0 + O ( L -1 ) where R 0 is the critical Rayleigh number for the corresponding infinite layer. Nonlinear solutions of constant phase exist for Rayleigh numbers R = R 0 + O ( L -2 ) but of these only the two that bifurcate at the lowest value of R are stable to two-dimensional linearized disturbances in this range (Daniels 1978). In the present paper one set of the class of phase-winding solutions is found to be stable to two-dimensional disturbances. For certain values of the Prandtl number of the fluid and for stress-free horizontal boundaries the results predict that to preserve stability there must be a continual readjustment of the roll pattern as the Rayleigh number is raised, with a corresponding increase in wavelength proportional to R - R 0 . These solutions also exhibit hysteresis as the Rayleigh number is raised and lowered. For other values of the Prandtl number the number of rolls remains unchanged as the Rayleigh number is raised, and the wavelength remains close to its critical value. It is proposed that the complete evolution of the flow pattern from a static state must take place on a number of different time scales of which t = O(( R - R 0 ) -1 ) and t = O(( R - R 0 ) -2 ) are the most significant. When t = O(( R - R 0 ) -1 ) the amplitude of convection rises from zero to its steady-state value, but the final lateral positioning of the rolls is only completed on the much longer time scale t = O(( R - R 0 ) -2 ).


Author(s):  
Ojas Satbhai ◽  
Subhransu Roy ◽  
Sudipto Ghosh

Direct numerical simulations for low Prandtl number fluid (Pr = 0.0216) are used to study the steady-state Rayleigh–Bénard convection (RB) in a two-dimensional unit aspect ratio box. The steady-state RB convection is characterized by analyzing the time-averaged temperature-field, and flow field for a wide range of Rayleigh number (2.1 × 105 ⩽ Ra ⩽ 2.1 × 108). It is seen that the time-averaged and space-averaged Nusselt number (Nuh¯) at the hot-wall monotonically increases with the increase in Rayleigh number (Ra) and the results show a power law scaling Nuh¯∝Ra0.2593. The current Nusselt number results are compared with the results available in the literature. The complex flow is analyzed by studying the frequency power spectra of the steady-state signal of the vertical velocity at the midpoint of the box for different Ra and probability density function of dimensionless temperature at various locations along the midline of the box.


2018 ◽  
Vol 837 ◽  
Author(s):  
Giovanni Fantuzzi

We prove the first rigorous bound on the heat transfer for three-dimensional Rayleigh–Bénard convection of finite-Prandtl-number fluids between free-slip boundaries with an imposed heat flux. Using the auxiliary functional method with a quadratic functional, which is equivalent to the background method, we prove that the Nusselt number $\mathit{Nu}$ is bounded by $\mathit{Nu}\leqslant 0.5999\mathit{R}^{1/3}$ uniformly in the Prandtl number, where $\mathit{R}$ is the Rayleigh number based on the imposed heat flux. In terms of the Rayleigh number based on the mean vertical temperature drop, $\mathit{Ra}$, we obtain $\mathit{Nu}\leqslant 0.4646\mathit{Ra}^{1/2}$. The scaling with Rayleigh number is the same as that of bounds obtained with no-slip isothermal, free-slip isothermal and no-slip fixed-flux boundaries, and numerical optimisation of the bound suggests that it cannot be improved within our bounding framework. Contrary to the two-dimensional case, therefore, the $\mathit{Ra}$-dependence of rigorous upper bounds on the heat transfer obtained with the background method for three-dimensional Rayleigh–Bénard convection is insensitive to both the thermal and the velocity boundary conditions.


2017 ◽  
Vol 830 ◽  
Author(s):  
Pranav Joshi ◽  
Hadi Rajaei ◽  
Rudie P. J. Kunnen ◽  
Herman J. H. Clercx

This experimental study focuses on the effect of horizontal boundaries with pyramid-shaped roughness elements on the heat transfer in rotating Rayleigh–Bénard convection. It is shown that the Ekman pumping mechanism, which is responsible for the heat transfer enhancement under rotation in the case of smooth top and bottom surfaces, is unaffected by the roughness as long as the Ekman layer thickness $\unicode[STIX]{x1D6FF}_{E}$ is significantly larger than the roughness height $k$. As the rotation rate increases, and thus $\unicode[STIX]{x1D6FF}_{E}$ decreases, the roughness elements penetrate the radially inward flow in the interior of the Ekman boundary layer that feeds the columnar Ekman vortices. This perturbation generates additional thermal disturbances which are found to increase the heat transfer efficiency even further. However, when $\unicode[STIX]{x1D6FF}_{E}\approx k$, the Ekman boundary layer is strongly perturbed by the roughness elements and the Ekman pumping mechanism is suppressed. The results suggest that the Ekman pumping is re-established for $\unicode[STIX]{x1D6FF}_{E}\ll k$ as the faces of the pyramidal roughness elements then act locally as a sloping boundary on which an Ekman layer can be formed.


1988 ◽  
Vol 190 ◽  
pp. 451-469 ◽  
Author(s):  
D. R. Jenkins

The relationship between observations of cellular Rayleigh-Bénard convection using shadowgraphs and theoretical expressions for convection planforms is considered. We determine the shadowgraphs that ought to be observed if the convection is as given by theoretical expressions for roll, square or hexagonal planforms and compare them with actual experiments. Expressions for the planforms derived from linear theory, valid for low supercritical Rayleigh number, produce unambiguous shadowgraphs consisting of cells bounded by bright lines, which correspond to surfaces through which no fluid flows and on which the vertical component of velocity is directed downwards. Dark spots at the centre of cells, indicating regions of hot, rising fluid, are not accounted for by linear theory, but can be produced by adding higher-order terms, predominantly due to the temperature dependence of a material property of the fluid, such as its viscosity.


Sign in / Sign up

Export Citation Format

Share Document