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.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε² being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε^8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k _c . For waves, the largest growth rate is of the order ε^4/3. In the case of Eckhaus instability the growth rate is of the order of α².     

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.     

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.     

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.     

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.     

The onset of magnetoconvection at large Prandtl number in a rotating layer I. Finite magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress—Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of order unity, as is the ratio of thermal to magnetic diffusivity. Attention is focused on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection. The case of main interest is the layer confined between electrically-insulating no-slip walls, but the analysis is guided by a parallel study based on illustrative boundary conditions that are mathematically simpler.     

An experimental investigation of convection in a rotating sphere subject to time varying thermal boundary conditions

Abstract

A series of experiments has been undertaken to investigate the onset of convection in a rapidly rotating fluid filled sphere. The boundary is subjected to a time varying temperature allowing the simulation of radial temperature profiles associated with internal heating. The system is similar to that treated theoretically by Roberts (1968), Busse (1970) and Soward (1977). It is found that Busse's modification of Roberts' linear analysis, taking into account velocity perturbations which are antisymmetric about the equatorial plane, provides a good estimate of the temperature gradient required to initiate convection. As observed in the experiments of Carrigan and Busse (1983) and predicted by linear theory, convection appears in the form of rolls or columns, aligned parallel to the rotation axis. As in earlier experiments, observed azimuthal wavenumbers are consistently smaller than predicted which we postulate to be a consequence of nonlinear effects. Owing to the presence of a centrifugally driven thermal wind, the predicted azimuthal drift of the rolls has not been observed.     

Thermal convection in a twisted horizontal magnetic field

The onset of convection in a layer of an electrically conducting fluid heated from below is considered in the case when the layer is permeated by a horizontal magnetic field of strength B ₀ the orientation of which varies sinusoidally with height. The critical value of the Rayleigh number for the onset of convection is derived as a function of the Chandrasekhar number Q. With increasing Q the height of the convection rolls decreases, while their horizontal wavelength slowly increases. Potential applications to the penumbral filaments of sunspots are briefly discussed.     

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.     

Finite prandtl number convection in spherical shells

Abstract

Finite amplitude convection in spherical shells with spherically symmetric gravity and heat source distribution is considered. The nonlinear problem of three-dimensional convection in shells with stress-free and isothermal boundaries is solved by expanding the dependent variables in terms of powers of the amplitude of convection. The preferred mode of convection is determined by a stability analysis in which arbitrary infinitesimal disturbances are superimposed on the steady solutions. The shell is assumed to be thick and only shells for which the ratio ζ of outer radius to inner radius is 2 or 3 are considered. Three cases, two of which lead to a self adjoint problem, are treated in this paper. The stable solutions are found to be l=2 modes for ζ=3 where l is the degree of the spherical harmonics and an l=3 non-axisymmetric mode which exhibits the symmetry of a tetrahedron for ζ=2. These stable solutions transport the maximum amount of heat. The Prandtl number dependence of the heat transport is computed for the various solutions analyzed in the paper.     

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.     

Effect of a uniform magnetic field on nonlinear magnetocenvection in a rotating fluid spherical shell

Abstract

Dynamic interaction between magnetic field and fluid motion is studied through a numerical experiment of nonlinear three-dimensional magnetoconvection in a rapidly rotating spherical fluid shell to which a uniform magnetic field parallel to its spin axis is applied. The fluid shell is heated by internal heat sources to maintain thermal convection. The mean value of the magnetic Reynolds number in the fluid shell is 22.4 and 10 pairs of axially aligned vortex rolls are stably developed. We found that confinement of magnetic flux into anti-cyclonic vortex rolls was crucial on an abrupt change of the mode of magnetoconvection which occurred at Δ = 1 ～ 2, where A is the Elsasser number. After the mode change, the fluid shell can store a large amount of magnetic flux in itself by changing its convection style, and the magnetostrophic balance among the Coriolis, Lorentz and pressure forces is established. Furthermore, the toroidal/poloidal ratio of the induced magnetic energy becomes less than unity, and the magnetized anti-cyclones are enlarged due to the effect of the magnetic force. Using these key ideas, we investigated the causes of the mode change of magnetoconvection. Considering relatively large magnetic Reynolds number and a rapid rotation rate of this model, we believe that these basic ideas used to interpret the present numerical experiment can be applied to the dynamics in the Earth's and other planetary cores.     

The onset of magnetoconvection at large Prandtl number in a rotating layer II. Small magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress-Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of unit order. In Part I of this series, it was also supposed that the ratio thermal diffusivity diffusivity/magnetic diffusivity is O(1), but here we suppose that this ratio is large. The character of the solution is changed in this limit. In the case of main interest, when the layer is confined between electrically-insulating no-slip walls, the solution is significantly different from the solution when the mathematically simpler, illustrative boundary conditions also considered in Part I are employed. As in Part I, attention is focussed on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection.     

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.     

Convection in a rotating magnetic system and Taylor's constraint

Abstract

A system is considered in which electrically conducting fluid is contained between two rigid horizontal planes and bounded laterally by a circular cylinder. The fluid is permeated by a strong azimuthal magnetic field. The strength of the field increases linearly with distance from the vertical axis of the cylinder, about which the entire system rotates rapidly. An unstable temperature gradient is maintained by heating the fluid from below and cooling from above. When viscosity and inertia are neglected, an arbitrary geostrophic velocity, which is aligned with the applied azimuthal magnetic field and independent of the axial coordinate, can be superimposed on the basic axisymmetric state. In this inviscid limit, the geostrophic velocity which occurs at the onset of convection is such that the net torque on geostrophic cylinders vanishes (Taylor's condition). The mathematical problem which describes the ensuing marginal convection is nonlinear, and was discussed previously for the planar case by Soward (1986). Here new features are isolated which result from the cylindrical geometry. New asymptotic solutions are derived valid when Taylor's condition is relaxed to include viscous effects.     

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.     

Planform of convection with strongly temperature-dependent viscosity

Abstract

This paper experimentally investigates the convective planform near critical in a fluid layer whose temperature-dependent viscosity varies from top to bottom by up to a factor of 1500. Convection occurs in three different planforms: rolls, hexagons and squares. The square planform, which appears only for fluids with viscosity variation greater than about 50, replaces the hexagonal convection pattern as the Rayleigh number increases much above critical. The large amplitude of hexagonal convection with strong viscosity variation precludes studying the hexagon-square transition with perturbation methods of the type used to study the hexagon-roll transitions at smaller viscosity variations.     

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.     

Shear flow instabilities in a rotating system

Abstract

The stability of a plane parallel shear flow with the profile U(z) = tanh z is considered in a rotating system with the axis of rotation in the z-direction. The establishment of the basic flow requires a baroclinic state, but baroclinic effects are suppressed in the stability analysis by assuming a limit of high thermal conductivity. It is shown that the strongest growing disturbance changes from a purely transverse form in the limit of vanishing rotation rate to a nearly longitudinal form as the angular velocity of rotation increases. An analytical solution of the stability equation is obtained for vanishing growth rates of the transverse form of the instability. But, in general, the solution of the problem requires numerical integrations which demonstrate that the preferred direction of the wave vector of the instability is towards the left of the direction of the mean flow.