Linear Magnetoconvection in Rotating Fluid Spheres Permeated by a Uniform Axial Magnetic Field

A linear analysis of thermally driven magnetoconvection is carried out with emphasis on its application to convection in the Earth's core. We consider a rotating and self-gravitating fluid sphere (or spherical shell) permeated by a uniform magnetic field parallel to the spin axis. In rapidly rotating cases, we find that five different convective modes appear as the uniform field is increased; namely, geostrophic, polar convective, magneto-geostrophic, fast magnetostrophic and slow magnetostrophic modes. The polar convective (P) and magneto-geostrophic (E) modes seem to be of geophysical interest. The P mode is characterized by such an axisymmetric meridional circulation that the fluid penetrates the equatorial plane, suggesting that generation of quadrapole from dipole fields could be explained by a linear process. The E mode is characterized by a few axially aligned columnar rolls which are almost two-dimensional due to a modified Proudman-Taylor theorem.     

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.     

Magnetoconvection

Abstract

Cowling investigated the effect of an imposed magnetic field on convection in order to explain the origin of sunspots. After summarizing the classical linear theory of Boussinesq magnetoconvection, this review proceeds to more recent nonlinear results. Weakly nonlinear theory is used to establish the relevant bifurcation structure, which involves steady, oscillatory and chaotic solutions. Behaviour found in numerical experiments can then be related to these analytical results. Thereafter, attention is focused on the astrophysically relevant problem of fully compressible magnetoconvection. Steady two-dimensional nonlinear solutions show two important effects: stratification introduces an asymmetry between rising and falling fluid, while compressibility leads to evacuated magnetic flux sheets. Time-dependent behaviour includes transitions between standing waves and travelling waves, as well as changes in horizontal scale, leading to the development of more complicated spatial structures. Work on three-dimensional models, which is now in progress, will lead to a better understanding of the structure of a sunspot.     

Sunspots

Abstract

Some examples of research on structure and formation of sunspots are briefly recollected in historical sequence. They relate to many facets of sunspots, first: magnetic inhibition of convection, the conjecture of a fiat penumbra, the stratification beneath the umbra, the observable magnetic profile, the Evershed effect as syphon flow, the concept of a magnetopause; next: cooling by Alfven waves, evolution and stability, the “bright ring”, the observed change of umbra brightness with the phase of the sunspot cycle, the hypothetical cluster of separate flux strands underneath the umbra, the profile of the magnetopause, the structure of the penumbra and the inclination of its field and finally: the concept of a deep penumbra with volume currents, exchange convection and the concept of a second current sheet separating umbra and penumbra.

Of course, the rigorous theoretical modeling of local magnetoconvection is an essential tool for our understanding of all these processes. I do not deal with it here, but the reader has a fascinating review of magnetoconvection already in his hands (Weiss, 1991).     

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.     

Formation of magnetic flux tubes in cylindrical wedge geometry

Three-dimensional (3D) MHD numerical simulations have not been able to demonstrate convincingly the spontaneous formation of large vertical flux tubes. Two-dimensional (2D) magnetoconvection in axisymmetric cylinders forms a central magnetic flux tube surrounded by annular convection rings. To study the robustness of this type of solution in three dimensions, the nonlinear resistive MHD equations are solved numerically in a 3D cylindrical wedge from an initially uniform vertical magnetic field. It is shown that the 2D result is retrieved for small domain radii. However, for larger radii the central axis loses its importance and in this case many convection cells form in the numerical domain. Magnetic flux is captured between cells where flow converges and the reduced amount of flux that congregates at the central axis is eroded by the surrounding convection.     

Turbulent magnetic fields in the Sun

Nigel Weiss recounts his Presidential Address 2001, given to the RAS A&G Ordinary Meeting on 9 February 2001.
Recent high-resolution observations, from the ground and from space, have revealed the fine structure of magnetic features at the surface of the Sun. At the same time, advances in computing power have at last made it possible to develop models of turbulent magnetoconvection that can be related to these observations. The key features of flux emergence and annihilation, as observed by the MDI experiment on SOHO, are reproduced in kinematic calculations, while three-dimensional numerical experiments reveal the dynamical processes that are involved. The pattern of convection depends on the strength of the magnetic field: as the mean field decreases, slender rising plumes give way to a regime where magnetic flux is separated from the motion and then to one where locally intense magnetic fields nestle between broad and vigorously convecting plumes. Moreover, turbulent convection is itself able to act as a small-scale dynamo, generating disordered fields near the solar surface.     

Thermal Generation of Alfvén Waves in Oscillatory Magnetoconvection: Diffusively Modified Modes

Thermal instabilities in the form of oscillatory magnetoconvection representing diffusively modified Alfvén waves in an electrically-conducting Bénard fluid layer of rigid walls in the presence of a vertical magnetic field are investigated. Emphasis of the article is on the transition from a nearly undamped Alfvén wave to diffusively modified Alfvén waves, and on the effect of physically realisable magnetic field boundary conditions on magnetoconvection. It is found that the extra magnetic dissipation in the magnetic Hartmann boundary layers can enhance oscillatory magnetoconvection in the form of strongly modified Alfvén waves. Oscillatory magnetoconvection produced solely by the Alfvén wave mechanism can be the most unstable mode even in the presence of a strong viscous effect. This article also represents the first study on the effect of an electrically conducting wall on magnetoconvection which is associated with a nonlinear eigenvalue problem. We find that the electrically perfectly conducting condition does not yield a good approximation for magnetoconvection with an electrically highly conducting wall. The size of oscillation frequency with an electrically highly conducting wall can be more than a factor of 2 larger than that obtained using the perfectly conducting condition.     

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.     

Mean-field concept and direct numerical simulations of rotating magnetoconvection and the geodynamo

Mean-field theory describes magnetohydrodynamic processes leading to large-scale magnetic fields in various cosmic objects. In this study magnetoconvection and dynamo processes in a rotating spherical shell are considered. Mean fields are defined by azimuthal averaging. In the framework of mean-field theory, the coefficients which determine the traditional representation of the mean electromotive force, including derivatives of the mean magnetic field up to the first order, are crucial for analyzing and simulating dynamo action. Two methods are developed to extract mean-field coefficients from direct numerical simulations of the mentioned processes. While the first method does not use intrinsic approximations, the second one is based on the second-order correlation approximation. There is satisfying agreement of the results of both methods for sufficiently slow fluid motions. Both methods are applied to simulations of rotating magnetoconvection and a quasi-stationary geodynamo. The mean-field induction effects described by these coefficients, e.g., the α-effect, are highly anisotropic in both examples. An α²-mechanism is suggested along with a strong γ-effect operating outside the inner core tangent cylinder. The turbulent diffusivity exceeds the molecular one by at least one order of magnitude in the geodynamo example. With the aim to compare mean-field simulations with corresponding direct numerical simulations, a two-dimensional mean-field model involving all previously determined mean-field coefficients was constructed. Various tests with different sets of mean-field coefficients reveal their action and significance. In the magnetoconvection and geodynamo examples considered here, the match between direct numerical simulations and mean-field simulations is only satisfying if a large number of mean-field coefficients are involved. In the magnetoconvection example, the azimuthally averaged magnetic field resulting from the numerical simulation is in good agreement with its counterpart in the mean-field model. However, this match is not completely satisfactory in the geodynamo case anymore. Here the traditional representation of the mean electromotive force ignoring higher than first-order spatial derivatives of the mean magnetic field is no longer a good approximation.     

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.     

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.     

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.     

Reducing the non-axisymmetry of a planetary dynamo and an application to saturn

Abstract

If a conducting fluid shell is undergoing spin-axisymmetric differential rotation and overlies the dynamo generating region of a planet then it is capable of greatly reducing the non-spin-axisymmetric components of the generated field, provided the appropriate magnetic Reynolds number is large. The model, closely related to the electromagnetic skin effect, is quantified and applied to Saturn. The observed small dipole tilt (～ 1°) of Saturn's magnetic field can be explained because of the presence of a stably stratified conducting layer overlying the dynamo region. This layer is a predicted consequence of the thermal evolution, arises because of the limited solubility of helium in metallic hydrogen (Stevenson, 1980), and appears to be required by the Voyager infrared observations indicating depletion of helium from Saturn's atmosphere. The much larger dipole tilt angles of Jupiter and the Earth indicate the absence of any such stable, differentially rotating layer with a large magnetic Reynolds number.     

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.     

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.     

Magnetic tracers,a probe of the solar convective layers

Abstract

A meridional circulation of sunspots has been measured through the digital analysis of the Meudon spectroheliograms from 1978 to 1983. Old and young sunspots follow a zonal meridional circulation, in several bands of latitude, in which two adjacent bands have opposite motions. This meridional circulation pattern is time-dependent. Using the H _α filaments as magnetic field tracers, a large-scale magnetic pattern has been found that was also obtained independently by direct measurement of the magnetic field (Hoeksema, 1988).

The coincidence of a large-scale magnetic pattern with a zonal meridional circulation suggests the existence of azimuthal rolls below the surface, and these azimuthal rolls can explain a number of properties of the solar cycle. New rolls occur with increasing proximity to the Equator, thereby indicating the direction of propagation of the dynamo wave. The occurrence of rolls is very favorable to the emergence of the magnetic regions. The rolls also influence the magnetic complexity of the active regions. They modulate the surface rotation through the Coriolis force, which accelerates or decelerates the fluid particles. They therefore offer a plausible explanation of the torsional oscillation pattern.

There are a number of problems raised by such an unexpected circulation pattern: for example, the coexistence of axisymmeric rolls with hypothetical giant cells, the location of the dynamo source below or within the convective zone, and the coupling of the radiative interior and the convective layers. To resolve these important issues, continuous observational studies are needed of the manifestation of solar activity, as well as of radius and luminosity variations. So, we have aimed our paper at an audience of theoreticians in the hope that they take up the challenges we describe.     

Dynamos with weakly convecting outer layers: implications for core-mantle boundary interaction

Convection in the Earth's core is driven much harder at the bottom than the top. This is partly because the adiabatic gradient steepens towards the top, partly because the spherical geometry means the area involved increases towards the top, and partly because compositional convection is driven by light material released at the lower boundary and remixed uniformly throughout the outer core, providing a volumetric sink of buoyancy. We have therefore investigated dynamo action of thermal convection in a Boussinesq fluid contained within a rotating spherical shell driven by a combination of bottom and internal heating or cooling. We first apply a homogeneous temperature on the outer boundary in order to explore the effects of heat sinks on dynamo action; we then impose an inhomogeneous temperature proportional to a single spherical harmonic Y ₂² in order to explore core-mantle interactions. With homogeneous boundary conditions and moderate Rayleigh numbers, a heat sink reduces the generated magnetic field appreciably; the magnetic Reynolds number remains high because the dominant toroidal component of flow is not reduced significantly. The dipolar structure of the field becomes more pronounced as found by other authors. Increasing the Rayleigh number yields a regime in which convection inside the tangent cylinder is strongly affected by the magnetic field. With inhomogeneous boundary conditions, a heat sink promotes boundary effects and locking of the magnetic field to boundary anomalies. We show that boundary locking is inhibited by advection of heat in the outer regions. With uniform heating, the boundary effects are only significant at low Rayleigh numbers, when dynamo action is only possible for artificially low magnetic diffusivity. With heat sinks, the boundary effects remain significant at higher Rayleigh numbers provided the convection remains weak or the fluid is stably stratified at the top. Dynamo action is driven by vigorous convection at depth while boundary thermal anomalies dominate in the upper regions. This is a likely regime for the Earth's core.     

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.