Rayleigh-Bénard Convection in Open and Closed Rotating Cavities

2005 ◽  
Vol 129 (2) ◽  
pp. 305-311 ◽  
Author(s):  
Martin P. King ◽  
Michael Wilson ◽  
J. Michael Owen

Buoyancy effects can be significant in the rotating annular cavities found between compressor discs in gas-turbine engines, where Rayleigh numbers above 1012 are common. In some engines, the cavity is “closed” so that the air is confined between four rotating surfaces: two discs and inner and outer cylinders. In most engines, however, the cavity is “open” and there is an axial throughflow of cooling air at the center. For open rotating cavities, a review of the published evidence suggests a Rayleigh–Bénard type of flow structure, in which, at the larger radii, there are pairs of cyclonic and anti-cyclonic vortices. The toroidal circulation created by the axial throughflow is usually restricted to the smaller radii in the cavity. For a closed rotating annulus, solution of the unsteady Navier–Stokes equations, for Rayleigh numbers up to 109, show Rayleigh–Bénard convection similar to that found in stationary enclosures. The computed streamlines in the r-θ plane show pairs of cyclonic and anti-cyclonic vortices; but, at the larger Rayleigh numbers, the computed isotherms suggest that the flow in the annulus is thermally mixed. At the higher Rayleigh numbers, the computed instantaneous Nusselt numbers are unsteady and tend to oscillate with time. The computed time-averaged Nusselt numbers are in good agreement with the correlations for Rayleigh–Bénard convection in a stationary enclosure, but they are significantly higher than the published empirical correlations for a closed rotating annulus.

Author(s):  
Martin P. King ◽  
Michael Wilson ◽  
J. Michael Owen

Buoyancy effects can be significant in the rotating annular cavities found between compressor discs in gas-turbine engines, where Rayleigh numbers above 1012 are common. In some engines, the cavity is ‘closed’, so that the air is confined between four rotating surfaces: two discs and inner and outer cylinders. In most engines, however, the cavity is ‘open’, and there is an axial throughflow of cooling air at the centre. For open rotating cavities, a review of the published evidence suggests a Rayleigh-Be´nard type of flow structure, in which, at the larger radii, there are pairs of cyclonic and anti-cyclonic vortices. The toroidal circulation created by the axial throughflow is usually restricted to the smaller radii in the cavity. For a closed rotating annulus, solution of the unsteady Navier-Stokes equations, for Rayleigh numbers up to 109, show Rayleigh-Be´nard convection similar to that found in stationary enclosures. The computed streamlines in the r-θ plane show pairs of cyclonic and anti-cyclonic vortices; but, at the larger Rayleigh numbers, the computed isotherms suggest that the flow in the annulus is thermally mixed. At the higher Rayleigh numbers, the computed instantaneous Nusselt numbers are unsteady and tend to oscillate with time. The computed time-average Nusselt numbers are in good agreement with the correlations for Rayleigh-Be´nard convection in a stationary enclosure, but they are significantly higher than the published empirical correlations for a closed rotating annulus.


1980 ◽  
Vol 98 (2) ◽  
pp. 329-348 ◽  
Author(s):  
Bharat Jhaveri ◽  
G. M. Homsy

We consider the onset of Rayleigh–Bénard convection from random fluctuations arising within a fluid. In the specific case in which the fluctuations are thermodynamically determined, we reduce the problem to a random initial value problem for the Fourier modes. For the case of weak nonlinear convection, it is possible to truncate the number of modes and this truncated set is solved both by a Monte Carlo technique and by moment methods for various Rayleigh numbers. We find three stages in the evolution of ordered convection from random fluctuations which correspond to time intervals in which the fluctuations and the nonlinearity have different degrees of importance. It is shown that no simple moment truncation method will succeed and that the time for onset of convection is a mean over a distribution of times for which members of an ensemble exhibit appreciable convective transport.


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):  
Hannes J. Brauckmann ◽  
Bruno Eckhardt ◽  
Jörg Schumacher

Rayleigh–Bénard convection and Taylor–Couette flow are two canonical flows that have many properties in common. We here compare the two flows in detail for parameter values where the Nusselt numbers, i.e. the thermal transport and the angular momentum transport normalized by the corresponding laminar values, coincide. We study turbulent Rayleigh–Bénard convection in air at Rayleigh number Ra =10 7 and Taylor–Couette flow at shear Reynolds number Re S =2×10 4 for two different mean rotation rates but the same Nusselt numbers. For individual pairwise related fields and convective currents, we compare the probability density functions normalized by the corresponding root mean square values and taken at different distances from the wall. We find one rotation number for which there is very good agreement between the mean profiles of the two corresponding quantities temperature and angular momentum. Similarly, there is good agreement between the fluctuations in temperature and velocity components. For the heat and angular momentum currents, there are differences in the fluctuations outside the boundary layers that increase with overall rotation and can be related to differences in the flow structures in the boundary layer and in the bulk. The study extends the similarities between the two flows from global quantities to local quantities and reveals the effects of rotation on the transport. This article is part of the themed issue ‘Toward the development of high-fidelity models of wall turbulence at large Reynolds number’.


2018 ◽  
Vol 858 ◽  
pp. 437-473 ◽  
Author(s):  
B. Favier ◽  
J. Purseed ◽  
L. Duchemin

We study the evolution of a melting front between the solid and liquid phases of a pure incompressible material where fluid motions are driven by unstable temperature gradients. In a plane-layer geometry, this can be seen as classical Rayleigh–Bénard convection where the upper solid boundary is allowed to melt due to the heat flux brought by the fluid underneath. This free-boundary problem is studied numerically in two dimensions using a phase-field approach, classically used to study the melting and solidification of alloys, which we dynamically couple with the Navier–Stokes equations in the Boussinesq approximation. The advantage of this approach is that it requires only moderate modifications of classical numerical methods. We focus on the case where the solid is initially nearly isothermal, so that the evolution of the topography is related to the inhomogeneous heat flux from thermal convection, and does not depend on the conduction problem in the solid. From a very thin stable layer of fluid, convection cells appear as the depth – and therefore the effective Rayleigh number – of the layer increases. The continuous melting of the solid leads to dynamical transitions between different convection cell sizes and topography amplitudes. The Nusselt number can be larger than its value for a planar upper boundary, due to the feedback of the topography on the flow, which can stabilize large-scale laminar convection cells.


Sign in / Sign up

Export Citation Format

Share Document