On the onset of thermal convection in slowly rotating fluid shells

Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Y^l _l (θ,φ) is preferred.     

Anisotropic turbulent thermal diffusion and thermal convection in a rapidly rotating fluid sphere

Estimates of the molecular values of magnetic, viscous and thermal diffusion suggest that the state of the Earth’s core is turbulent and that complete numerical simulation of the geodynamo is not realizable at present. Large eddy simulation of the geodynamo with modelling of the sub-grid scale turbulence must be used. Current geodynamo models effectively model the sub-grid scale turbulence with isotropic diffusivities larger than the molecular values appropriate for the core. In the Braginsky and Meytlis (1990) picture of core turbulence the thermal and viscous diffusivities are enhanced up to the molecular magnetic diffusivity in the directions of the rotation axis and mean magnetic field. We neglect the mean magnetic field herein to isolate the effects of anisotropic thermal diffusion, enhanced or diminished along the rotation axis, and explore the instability of a steady conductive basic state with zero mean flow in the Boussinesq approximation. This state is found to be more stable (less stable) as the thermal diffusion parallel to the rotation axis is increased (decreased), if the transverse thermal diffusion is fixed. To examine the effect of simultaneously varying the diffusion along and transverse to the rotation axis, the Frobenius norm is used to control for the total thermal diffusion. When the Frobenius norm of the thermal diffusion tensor is fixed, it is found that increasing the thermal diffusion parallel to the rotation axis is destabilising. This result suggests that, for a fixed total thermal diffusion, geodynamo codes with anisotropic thermal diffusion may operate at lower modified Rayleigh numbers.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

On the onset of convection in rotating spherical shells

Abstract

The linear problem of the onset of convection in rotating spherical shells is analysed numerically in dependence on the Prandtl number. The radius ratio η=r _i/r _o of the inner and outer radii is generally assumed to be 0.4. But other values of η are also considered. The goal of the analysis has been the clarification of the transition between modes drifting in the retrograde azimuthal direction in the low Taylor number regime and modes traveling in the prograde direction at high Taylor numbers. It is shown that for a given value m of the azimuthal wavenumber a single mode describes the onset of convection of fluids of moderate or high Prandtl number. At low Prandtl numbers, however, three different modes for a given m may describe the onset of convection in dependence on the Taylor number. The characteristic properties of the modes are described and the singularities leading to the separation with decreasing Prandtl number are elucidated. Related results for the problem of finite amplitude convection are also reported.     

Non-linear marginal convection in a rotating magnetic system

Abstract

An inviscid, electrically conducting fluid is contained between two rigid horizontal planes and bounded laterally by two vertical walls. The fluid is permeated by a strong uniform horizontal magnetic field aligned with the side wall boundaries and the entire system rotates rapidly about a vertical axis. The ratio of the magnitudes of the Lorentz and Coriolis forces is characterized by the Elsasser number, A, and the ratio of the thermal and magnetic diffusivities, q. By heating the fluid from below and cooling from above the system becomes unstable to small perturbations when the adverse density gradient as measured by the Rayleigh number, R, is sufficiently large.

With the viscosity ignored the geostrophic velocity, U, which is aligned with the applied magnetic field, is independent of the coordinate parallel to the rotation axis but is an arbitrary function of the horizontal cross-stream coordinate. At the onset of instability the value of U taken ensures that Taylor's condition is met. Specifically the Lorentz force, which results from marginal convection must not cause any acceleration of the geostrophic flow. It is found that the critical Rayleigh number characterising the onset of instability is generally close to the corresponding value for the usual linear problem, in which Taylor's condition is ignored and U is chosen to vanish. Significant differences can occur when q is small owing to a complicated flow structure. There is a central interior region in which the local magnetic Reynolds number, R_m , based on U is small of order q and on exterior region in which R_m is of order unity.     

Generation of magnetic fields by convection in a rotating sphere,I

The magnetohydrodynamic dynamo problem is solved for an electrically conducting spherical fluid shell with spherically symmetric distributions of gravity and heat sources. The dynamics of motions generated by thermal buoyancy are dominated by the effects of rotation of the fluid shell. Dynamos are found for low and intermediate values of the Taylor number, T ? 10⁵, if the scale of the nonaxisymmetric component of the velocity field is sufficiently small. The generation of magnetic fields of quadrupolar symmetry is preferred at Rayleigh numbers close to the critical value R_c for onset of convection. As the Rayleigh number increases, the generation of dipolar magnetic fields becomes preferred.     

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.     

Baroclinically-driven flows and dynamo action in rotating spherical fluid shells

The dynamics of stably stratified stellar radiative zones is of considerable interest due to the availability of increasingly detailed observations of Solar and stellar interiors. This article reports the first non-axisymmetric and time-dependent simulations of flows of anelastic fluids driven by baroclinic torques in stably stratified rotating spherical shells – a system serving as an elemental model of a stellar radiative zone. With increasing baroclinicity a sequence of bifurcations from simpler to more complex flows is found in which some of the available symmetries of the problem are broken subsequently. The poloidal component of the flow grows relative to the dominant toroidal component with increasing baroclinicity. The possibility of magnetic field generation thus arises and this paper proceeds to provide some indications for self-sustained dynamo action in baroclinically-driven flows. We speculate that magnetic fields in stably stratified stellar interiors are thus not necessarily of fossil origin as it is often assumed.     

Mean-field equations for weakly nonlinear two-scale perturbations of forced hydromagnetic convection in a rotating layer

We consider stability of regimes of hydromagnetic thermal convection in a rotating horizontal layer with free electrically-conducting boundaries, to perturbations involving large spatial and temporal scales. Equations governing the evolution of weakly nonlinear mean perturbations are derived under the assumption that the α-effect is insignificant in the leading-order (e.g. due to a symmetry of the system). The mean-field equations generalise the standard equations of hydromagnetic convection: New terms emerge – a second-order linear operator representing the combined eddy diffusivity and quadratic terms associated with the eddy advection. If the perturbed CHM regime is nonsteady and insignificance of the α-effect in the system does not rely on the presence of a spatial symmetry, the combined eddy diffusivity operator also involves a nonlocal pseudodifferential operator. If the perturbed CHM state is almost symmetric, α-effect terms appear in the mean-field equations as well. Near a point of a symmetry-breaking bifurcation, cubic nonlinearity emerges in the equations. All the new terms are in general anisotropic. A method for evaluation of their coefficients is presented; it requires solution of a significantly smaller number of auxiliary problems than in a straightforward approach.     

Finite-Amplitude thermal convection within a self-gravitating fluid sphere

Abstract

Finite-difference calculations have been carried out to determine the structure of finite-amplitude thermal convection within a self-gravitating fluid sphere with uniform heat release. For a fixed-surface boundary condition single-cell convection breaks up into double-cell convection at a Rayleigh number of 3 × 10⁴, at a Rayleigh number of 5 × 10⁵ four-cell convection is observed. With a free-surface boundary condition only single cell convection is obtained up to a Rayleigh number of 5 × 10⁶.     

Specific features of development of two-dimensional convection in a rotating compressible conducting fluid

An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E _M = 3 × 10⁻³. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.     

Finite amplitude convection and magnetic field generation in a rotating spherical shell

Abstract

Finite amplitude solutions for convection in a rotating spherical fluid shell with a radius ratio of η=0.4 are obtained numerically by the Galerkin method. The case of the azimuthal wavenumber m=2 is emphasized, but solutions with m=4 are also considered. The pronounced distinction between different modes at low Prandtl numbers found in a preceding linear analysis (Zhang and Busse, 1987) is also found with respect to nonlinear properties. Only the positive-ω-mode exhibits subcritical finite amplitude convection. The stability of the stationary drifting solutions with respect to hydrodynamic disturbances is analyzed and regions of stability are presented. A major part of the paper is concerned with the growth of magnetic disturbances. The critical magnetic Prandtl number for the onset of dynamo action has been determined as function of the Rayleigh and Taylor numbers for the Prandtl numbers P=0.1 and P=1.0. Stationary and oscillatory dynamos with both, dipolar and quadrupolar, symmetries are close competitors in the parameter space of the problem.     

Three-dimensional convection in spherical shells

Abstract

In this study, the equations of the three-dimensional convective motion of an infinite Prandtl number fluid are solved in spherical geometry, for Rayleigh numbers up to 15 times the critical number. An iterative method is used to find stationary solutions. The spherical parts of the operators are treated using a Galerkin collocation method while the radial and time dependences are expressed using finite difference methods. A systematic search for stationary solutions has led to eight different stream patterns for a low Rayleigh number (1.28 times the critical number). They can be classified as:

I) Axisymmetrical solutions, analogous to rolls in plane geometry.

II) Solutions which have several ascending plumes within a large area of ascending current, and also several descending plumes within an area of descending current. This type of flow is analogous to bimodal circulation in plane geometry.

III) Solutions characterized by isolated ascending (or descending) plumes separated from each other by a closed polyhedral network of descending (or ascending) currents. This type of circulation is called ‘polygonal’ in analogy with hexagonal circulation in plane geometry.

The behaviour of each of the eight solutions has been studied by increasing the Rayleigh number up to 15 times the critical number. A trend towards transitions from type (I) and type (II) solutions to type (III) solutions is observed. It is inferred that only the “polygonal” solutions are stable for a Rayleigh number greater than 15 times the critical number.     

Thermal and magnetic instabilities in a rapidly rotating fluid sphere

Abstract

A linear analysis is used to study the stability of a rapidly rotating, electrically-conducting, self-gravitating fluid sphere of radius r ₀, containing a uniform distribution of heat sources and under the influence of an azimuthal magnetic field whose strength is proportional to the distance from the rotation axis. The Lorentz force is of a magnitude comparable with that of the Coriolis force and so convective motions are fully three-dimensional, filling the entire sphere. We are primarily interested in the limit where the ratio q of the thermal diffusivity κ to the magnetic diffusivity η is much smaller than unity since this is possibly of the greatest geophysical relevance.

Thermal convection sets in when the temperature gradient exceeds some critical value as measured by the modified Rayleigh number R_c. The critical temperature gradient is smallest (R_c reaches a minimum) when the magnetic field strength parameter Λ ? 1. [R_c and Λ are defined in (2.3).] The instability takes the form of a very slow wave with frequency of order κ/r ² ₀ and its direction of propagation changes from eastward to westward as Λ increases through Λ _c ? 4.

When the fluid is sufficiently stably stratified and when Λ > Λ_m ? 22 a new mode of instability sets in. It is magnetically driven but requires some stratification before the energy stored in the magnetic field can be released. The instability takes the form of an eastward propagating wave with azimuthal wavenumber m = 1.     

The onset of convection in a rotating layer of viscous fluid with an imposed magnetic field: the dependence on the prandlt numbers

The onset of Boussinesq convection in a horizontal layer of an electrically conducting incompressible fluid is considered. The layer rotating about a vertical axis is heated from below; a vertical magnetic field is imposed. Rigid electrically insulating boundaries are assumed. The loss of stability of the trivial steady state, which occurs as the Rayleigh numbers increase, can be accompanied by the development of a monotonic or an oscillatory instability, depending on the parameter values of the problem at hand (the Taylor number, the Chandrasekhar number, the kinematic and the magnetic Prandtl numbers). When the instability is monotonic, the emerging convective rolls themselves are also unstable if the Taylor number is sufficiently large (the so-called Küppers-Lortz instability takes place). In the present work it is studied how the critical value of the Rayleigh number, the type of the trivial steady state instability, and the critical value of the Taylor number for the Küppers-Lortz instability depend on the kinematic and the magnetic Prandtl numbers. We consider the values of the Prandtl number not exceeding 1, which is typical for the outer core of the Earth.     

Numerical simulations of thermal convection in a rapidly rotating spherical shell cooled inhomogeneously from above

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell with and without inhomogeneous temperature anomalies on the top boundary have been carried out using a three-dimensional, time-dependent, spectral-transform code. The spherical shell of Boussinesq fluid has inner and outer radii the same as those of the Earth's liquid outer core. The Taylor number is 10⁷, the Prandtl number is 1, and the Rayleigh number R is 5R_c (R_c is the critical value of R for the onset of convection when the top boundary is isothermal and R is based on the spherically averaged temperature difference across the shell). The shell is heated from below and cooled from above; there is no internal heating. The lower boundary of the shell is isothermal and both boundaries are rigid and impermeable. Three cases are considered. In one, the upper boundary is isothermal while in the others, temperature anomalies with (l,m) = (3,2) and (6,4) are imposed on the top boundary. The spherically averaged temperature difference across the shell is the same in all three cases. The amplitudes of the imposed temperature anomalies are equal to one-half of the spherically averaged temperature difference across the shell. Convective structures are strongly controlled by both rotation and the imposed temperature anomalies suggesting that thermal inhomogeneities imposed by the mantle on the core have a significant influence on the motions inside the core. The imposed temperature anomaly locks the thermal perturbation structure in the outer part of the spherical shell onto the upper boundary and significantly modifies the velocity structure in the same region. However, the radial velocity structure in the outer part of the shell is different from the temperature perturbation structure. The influence of the imposed temperature anomaly decreases with depth in the shell. Thermal structure and velocity structure are similar and convective rolls are more columnar in the inner part of the shell where the effects of rotation are most dominant.     

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.     

Implementation of a high-order combined compact difference scheme in problems of thermally driven convection and dynamo in rotating spherical shells

We present an improved solution method for modeling thermally driven convection and dynamo in a rotating spherical shell. In this method, we introduce a high-order three-point combined compact difference scheme (CCDS) on non-uniform grid points in radius, while spherical harmonic expansion is conventionally performed in the angular direction. The governing equations in the spectral form are time-stepped together with the implicit CCDS up to the second derivative. To improve stability of the scheme, a boundary closure scheme is developed on non-uniform mesh. Numerical comparison with a published benchmark solution at moderate Ekman and Rayleigh numbers demonstrates that accuracy and convergence of the CCDS is fairly good and superior to the existing finite difference scheme using more stencil. With this scheme, we could more accurately solve problems of convection and also dynamo action in planetary core with less grid points.