Large wavenumber convection in the rotating annulus

Abstract

The annulus model considers convection between concentric cylinders with sloping endwalls. It is used as a simplified model of convection in a rapidly rotating sphere. Large azimuthal wavenumbers are preferred in this problem, and this has been exploited to develop an asymptotic approach to nonlinear convection in the annulus. The problem is further reduced because the Taylor-Proudman constraint simplifies the dependence in the direction of the rotation vector, so that a nonlinear system dependent only on the radial variable and time results. As Rayleigh number is increased a sequence of bifurcations is found, from steady solutions to periodic solutions and 2-tori, typically ending in chaotic behaviour. Both the magnetic (MHD convection) and non-magnetic problem has been considered, and in the non-magnetic case our bifurcation sequence can be compared with those found by previous two-dimensional numerical simulations.     

On simple models for simulation of nonlinear processes in convection and turbulence

The mechanism of nonlinear interaction in hydrodynamics is studied with dynamical systems having finite degrees of freedom. The equations are assumed to have the same integrals of motion and main features as those peculiar to hydrodynamical equations. The simplest system of this kind is a triplet (a system described by three parameters). Its equations of motion coincide with the Euler equations in the theory of the gyroscope. The forced motion of a triplet is treated theoretically. A real hydrodynamical system controlled by the equations of motion of a triplet was devised and verified in the laboratory.

The simplest theoretical model of baroclinic motion which provides a basis for studies of of forced heat convection in an ellipsoidal cavity was also constructed. Under certain conditions, the addition of rotation causes a regime of motion analogous to the Rossby regime in a rotating annulus.

More complicated models constructed from a large number of interacting triplets can simulate the cascade process of energy transformation in developed turbulence.     

Nonlinear Convection in a Rotating Annulus with a Finite Gap

Thermal convection in a fluid-filled gap between the two corotating, concentric cylindrical sidewalls with sloping curved ends driven by radial buoyancy was first studied by Busse (Busse, F.H., "Thermal instabilities in rapidly rotating systems", J. Fluid Mech . 44 , 441-460 (1970)). The annulus model captures the key features of rotating convection in full spherical geometry and has been widely employed to study convection, magnetoconvection and dynamos in planetary systems, usually in connection with the small-gap approximation neglecting the effect of azimuthal curvature of the annulus. This article investigates nonlinear thermal convection in a rotating annulus with a finite gap through numerical simulations of the full set of nonlinear convection equations. Three representative cases are investigated in detail: a large-gap annulus with the ratio of the radii ( s i and s o ) of the sidewalls ξ = s i / o s = 0.1, a medium-gap annulus with ξ = 0.35 and a small-gap annulus with ξ = 0.8. Near the onset of convection, the effect of rapid rotation through the sloping ends forces the first (Hopf) bifurcation in the form of small-scale, steadily drifting rolls (thermal Rossby waves). At moderately large Rayleigh numbers, a variety of different convection patterns are found, including mixed-mode steadily drifting, quasi-periodic (vacillating) and temporally chaotic convection in association with various temporal and spatial symmetry-breaking bifurcations. Our extensive simulations suggest that competition between nonlinear and rotational effects with increasing Rayleigh number leads to an unusual sequence of bifurcation characterized by enlarging the spatial scale of convection.     

Thermal convection in differentially rotating systems

Abstract

The onset of convection in a cylindrical fluid annulus is analyzed in the case when the cylindrical walls are rotating differentially, a temperature gradient in the radial direction is applied, and the centrifugal force dominates over gravity. The small gap approximation is used and no-slip conditions on the cylindrical walls are assumed. It is found that over a considerable range of the parameter space either convection rolls aligned with the axis of rotation or rolls in the perpendicular (azimuthal) direction are preferred. It is shown that by a suitable redefinition of parameters, results for finite amplitude Taylor vortices and for convection rolls in the presence of shear can be applied to the present problem. Weakly nonlinear results for transverse rolls in a Couette flow indicate the possibility of subcritical bifurcation for Prandtl numbers P less than 0.82. Heat and momentum transports are derived as functions of P and the problem of interaction between transverse and longitudinal rolls is considered. The relevance of the analysis for problems of convection in planetary and stellar atmospheres is briefly discussed.     

Two-and three-dimensional linear convection in a rotating annulus

Abstract

An exceptional case to the model-independent theory of Knobloch (1995) is presented, by investigating a rotating cylindrical annulus of height H and side wall radii r _o and r _i, with non-slip, perfectly thermally conducting side walls and thermally insulating stress-free ends. Radial heating permits the possibility of either two- or three-dimensional convective solutions being the preferred mode. An analytical solution is obtained for the two-dimensional case and a numerical solution for the three-dimensional solution, which is also applied to the two-dimensional solution. It is shown that both two- and three-dimensional solutions can be realized depending on the aspect ratio, γ = H/d, where d = r _o-r _i is the thickness of the annulus, the radii ratio λ = r _i/r _o and the rotation rate of the model. For γ = O(1) and λ = 0.4, the preferred convective solution is three-dimensional when the Taylor number, T < 10² and two-dimensional for T > 10². For small aspect ratios, γ ? 1, the preferred mode is two-dimensional for all rotation rates.     

A benchmark comparison of numerical methods for infinite Prandtl number thermal convection in two-dimensional Cartesian geometry

Abstract

A comparison is made between seven different numerical methods for calculating two-dimensional thermal convection in an infinite Prandtl number fluid. Among the seven methods are finite difference and finite element techniques that have been used to model thermal convection in the Earth's mantle. We evaluate the performance of each method using a suite of four benchmark problems, ranging from steady-state convection to intrinsically time-dependent convection with recurring thermal boundary layer instabilities. These results can be used to determine the accuracy of other computational methods, and to assist in the development of new ones.     

Laboratory experiments in a baroclinic annulus with heating and cooling on the horizontal boundaries

Abstract

Experiments have been performed in a cylindrical annulus with horizontal temperature gradients imposed upon the horizontal boundaries and in which the vertical depth was smaller than the width of the annulus. Qualitative observations were made by the use of small, suspended, reflective flakes in the liquid (water).

Four basic regimes of flow were observed: (1) axisymmetric flow, (2) deep cellular convection, (3) boundary layer convective rolls, and (4) baroclinic waves. In some cases there was a mix of baroclinic and convective instabilities present. As a “mean” interior Richardson number was decreased from a value greater than unity to one less than zero, axisymmetric baroclinic instability of the Solberg type was never observed. Rather, the transition was from non-axisymmetric baroclinic waves, to a mix of baroclinic and convective instability, to irregular cellular convection.     

Topology of Rayleigh–Bénard convection and magnetoconvection in plane layer

ABSTRACT

The present study aims to link the dynamics of geophysical fluid flows with their vortical structures in physical space and to study the transition of these structures due to the control parameters. The simulations are carried in a rectangular box filled with liquid gallium for three different cases, namely, Rayleigh–Bénard convection (RBC), magnetoconvection (MC) and rotating magnetoconvection (RMC). The physical setup and material properties are similar to those considered by Aurnou and Olson in their experimental work. The simulated results are validated with theoretical results of Chandrasekhar and experimental results of Aurnou and Olson. The results are also topologically verified with the help of Euler number given by Ma and Wang. For RBC, the onset is obtained at Ra greater than 1708 and at this Ra, the symmetric rolls are orientated in/along a horizontal axis. As the value of Ra increases further, the width of the horizontal rolls starts to amplify. It is observed that these two-dimensional rolls are nothing but the cross-sections of three-dimensional (3D) cylindrical rolls with wave structures. When the vertically imposed magnetic field is added to RBC, the onset of convection is delayed due to the effect of Lorentz force on the thermal buoyancy force. The presence of 3D rectangular structures is highlighted and analysed. When the magnetically influenced rectangular box rotates about vertical axis at low rotation rates in magnetoconvection model, the onset of convection gets further delayed by magnetic field, which is in general agreement with the theoretical predictions. The critical Ra increases linearly with magnetic field intensity. Coherent thermal oscillations are detected near the onset of convection, at moderate rotation rates.     

Nonlinear convection in an imposed horizontal magnetic field

Abstract

Nonlinear two-dimensional magnetoconvection, with a Boussinesq fluid driven across the field-lines, is taken as a model for giant-cell convection in the sun and late-type stars. A series of numerical experiments shows the sensitivity of the horizontal scale of convection to the applied field and to the Rayleigh number R. Overstable oscillations occur in cells as broad as they are deep, but increasing R leads to steady motions of much greater wavelength. Purely geometrical effects can cause oscillation: this work implies that strong horizontal field will in general lead to time-dependent convection.     

On the transition between axisymmetric and non-axisymmetric flow in a rotating liquid annulus subject to a horizontal temperature gradient: Hysteresis effects at moderate Taylor number and baroclinic waves beyond the eady cut-off at high Taylor number

Abstract

The transition between axisymmetric and non-axisymmetric régimes of flow in a rotating annulus of liquid subject to horizontal temperature gradient is known from previous experimental studies to depend largely on two dimensionless parameters. These are Θ, which is proportional to the impressed density contrast Δρ and inversely proportional to the square of the angular speed of rotation ω, and  (Taylor number), which is proportional to ω² /v² where v is the coefficient of kinematic viscosity. At moderate values of , around 10⁷, the critical value of Θ above which axisymmetric flow is found to OCCUT and below which non-axisymmetric fully-developed baroclinic waves (sloping convection) occur, is fairly insensitive to . Though sharp, the transition exhibits marked hysteresis when the upper surface of the liquid is free (but not when the upper surface is in contact with a rigid lid), and it is argued on the basis of the experimental evidence supported by various results of baroclinic instability theory that both the sharpness of the transition and the hysteresis phenomenon are consequences of the combined effects of potential vorticity gradients and viscosity on the process of sloping convection.

We also present some new experiments on fully-developed baroclinic waves, conducted in a large rotating annulus using liquids of very low viscosity (di-ethyl ether), thus attaining values of  as high as 10⁹ to 10¹⁰. The transition from axisymmetric to non-axisymmetric flow is found to lose its sharpness at such high values of , and it is argued that this occurs because viscosity is no longer able to inhibit instabilities at wavelengths less than the so-called ‘Eady short-wave cut-off’, which owe their existence to potential vorticity gradients in the main body of the fluid.     

Transitions to broad cells in a nonlinear thermal convection system

Abstract

This paper explores the properties of a two-dimensional, Boussinesq convection model with an ad hoc term in the buoyancy tendency equation that represents a positive external feedback process acting on the buoyancy fluctuations. Linear stability analyses and nonlinear integrations are presented for the case of constant heat flux boundary conditions. Although the large wavenumber modes grow the fastest from a state of rest, the nonlinear solutions progressively evolve to cells of small wavenumber. Applications to mesoscale cellular convection in the atmosphere are discussed.     

Models of solar differential rotation

Abstract

Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect.

In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms^?1) or a slow (～ 1 ms^?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions.

The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded.     

Convection in a rotating annulus having sidewalls of height-dependent thermal conductance

Abstract

Thermal convection in a vertically-mounted, rotating annulus of a particular design proposed by Davies and Walin (1977) is investigated. The annulus used in the present study differs from the conventional type in some important aspects: the sidewalls are finitely conducting, and the thermal conductance of the sidewalls is height-dependent. The theoretical model due to Davies and Walin is briefly recounted. The present study aims to verify the theoretical model; we have acquired numerical solutions to the governing Navier-Stokes equations. The numerical results are supportive of the theoretical contentions. The near-linear dependence of the isothermal slope on the parameter D, which is a function of Ω and ΔT, is corroborated within reasonable limits. New data on the vertical and radial structures of the meridional and azimuthal flows are presented. The numerical results also confirm that the shape of the sidewall thickness has a substantial influence on the meridional flow patterns. In the bulk of the interior flow field, the dominant azimuthal flow field and the temperature field are linked by the thermal wind relation.     

A review of: “Problems of Solar and Stellar Oscillations. Iau Colloquium No. 66. Edited by D. O. Gough. D. Reidel Publishing Company, 493 pp. Df1195.00/U.S. $84.50 (ISBN 90-277-1554-8). 1983.”

Abstract

The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed.     

Convection driven by centrifugal bouyancy in a rotating annulus

Abstract

Drift rates and amplitudes of convection columns driven by centrifugal bouyancy in a cylindrical fluid annulus rotating about a vertical axis have been measured by thermistor probes. Conical top and bottom boundaries of the annular fluid region are responsible for the prograde Rossby wave like dynamics of the convection columns. A constant positive temperature difference between the outer and the inner cylindrical boundaries is generated by the circulation of thermostatically controled water. Mercury and water have been used as converting fluids. The measurements extend the earlier visual observations of Busse and Carrigan (1974) and provide quantitative data for an eventual comparison with nonlinear theories of thermal Rossby waves. The measured drift frequencies are in general agreement with linear theory. Of particular interest is the decline of the amplitude of convection with increasing Rayleigh number in a region beyond the onset of convection.     

Bénard convection in a rotating fluid

Abstract

The steady nonlinear regime of Bénard convection in a uniformly rotating fluid is treated using a two-dimensional primitive-equation numerical model with rigid boundaries. Quantitative comparisons with laboratory heat transport data for water are made in the parameter ranges for which the experimental flows are approximately two-dimensional and steady. When an experimentally realistic spatial periodicity is imposed upon the numerical solution, the model simulates the experimental determinations of Nusselt number fairly accurately. In particular, it predicts the observed non-monotonic dependence on Taylor number. When spatial periodicities corresponding to those of the linear stability problem are specified, however, the accuracy of the simulation is less and the Taylor number dependence is monotonic.     

A note on the heat transfer by the axisymmetric thermal convection in a rotating fluid annulus

Abstract

Some new measurements are presented of the axisymmetric heat transport in a differentially heated rotating fluid annulus. Both rigid and free upper surface cases are studied, for Prandtl numbers of 7 and 45, from low to high rotation rates. The rigid lid case is extended to high rotation rates by suppressing the baroclinic waves, that would normally develop at some intermediate rotation rate, with the use of sloping endwalls.

A parameter P is defined as the square of the ratio of the (non-rotating) thermal sidewall layer thickness to the Ekman layer thickness. For small P the heat transport remains unaffected by the rotation, but as P increases to order unity the Ekman layer becomes thin enough to inhibit the radial mass transport, and hence the heat flux. No explicit Prandtl number dependence is observed. Also this scaling allows the identification of the region in which the azimuthal velocity reaches its maximum. Direct comparisons are drawn with previous experimental and numerical results, which show what can be interpreted as an inhibiting effect of increasing curvature on the heat transport.     

Linear Stability of Thermal Convection in Rotating Systems with Fixed Heat Flux Boundaries

Linear stability of rotating thermal convection in a horizontal layer of Boussinesq fluid under the fixed heat flux boundary condition is examined by the use of a vertically truncated system up to wavenumber one. When the rotation axis is in the vertical direction, the asymptotic behavior of the critical convection for large rotation rates is almost the same as that under the fixed temperature boundary condition. However, when the rotation axis is horizontal and the lateral boundaries are inclined, the mode with zero horizontal wavenumber remains as the critical mode regardless of the rotation rate. The neutral curve has another local minimum at a nonzero horizontal wavenumber, whose asymptotic behavior coincides with the critical mode under the fixed temperature condition. The difference of the critical horizontal wavenumber between those two geometries is qualitatively understood by the difference of wave characteristics; inertial waves and Rossby waves, respectively.     

Convection with heat flux prescribed on the boundaries of the system. I. The effect of temperature dependence of material properties

Abstract

We study the bifurcation to steady two-dimensional convection with the heat flux prescribed on the fluid boundaries. The fluid is weakly non-Boussinesq on account of a slight temperature dependence of its material properties. Using expansions in the spirit of shallow water theory based on the preference for large horizontal scales in fixed flux convection, we derive an evolution equation for the horizontal structure of convective cells. In the steady state, this reduces to a simple nonlinear ordinary differential equation. When the horizontal scales of the cells exceed a certain critical size, the bifurcation to steady convection is subcritical and the degree of subcriticality increases with increasing cell size.