From stable dipolar towards reversing numerical dynamos

We are using a three-dimensional convection-driven numerical dynamo model without hyperdiffusivity to study the characteristic structure and time variability of the magnetic field in dependence of the Rayleigh number (Ra) for values up to 40 times supercritical. We also compare a variety of ways to drive the convection and basically find two dynamo regimes. At low Ra, the magnetic field at the surface of the model is dominated by the non-reversing axial dipole component. At high Ra, the dipole part becomes small in comparison to higher multipole components. At transitional values of Ra, the dynamo vacillates between the dipole-dominated and the multipolar regime, which includes excursions and reversals of the dipole axis. We discuss, in particular, one model of chemically driven convection, where for a suitable value of Ra, the mean dipole moment and the temporal evolution of the magnetic field resemble the known properties of the Earth’s field from paleomagnetic data.     

The initial temporal evolution of a feedback dynamo for Mercury

Various possibilities are currently under discussion to explain the observed weakness of the intrinsic magnetic field of planet Mercury. One of the possible dynamo scenarios is a dynamo with feedback from the magnetosphere. Due to its weak magnetic field, Mercury exhibits a small magnetosphere whose subsolar magnetopause distance is only about 1.7 Hermean radii. We consider the magnetic field due to magnetopause currents in the dynamo region. Since the external field of magnetospheric origin is antiparallel to the dipole component of the dynamo field, a negative feedback results. For an αΩ-dynamo, two stationary solutions of such a feedback dynamo emerge: one with a weak and the other with a strong magnetic field. The question, however, is how these solutions can be realized. To address this problem, we discuss various scenarios for a simple dynamo model and the conditions under which a steady weak magnetic field can be reached. We find that the feedback mechanism quenches the overall field to a low value of about 100–150 nT if the dynamo is not driven too strongly.     

The geodynamo as a bistable oscillator

Abstract

Our intent is to provide a simple and quantitative understanding of the variability of the axial dipole component of the geomagnetic field on both short and long time scales. To this end we study the statistical properties of a prototype nonlinear mean field model. An azimuthal average is employed, so that (1) we address only the axisymmetric component of the field, and (2) the dynamo parameters have a random component that fluctuates on the (fast) eddy turnover time scale. Numerical solutions with a rapidly fluctuating α reproduce several features of the geomagnetic field: (1) a variable, dominantly dipolar field with additional fine structure due to excited overtones, and sudden reversals during which the field becomes almost quadrupolar, (2) aborted reversals and excursions, (3) intervals between reversals having a Poisson distribution. These properties are robust, and appear regardless of the type of nonlinearity and the model parameters. A technique is presented for analysing the statistical properties of dynamo models of this type. The Fokker-Planck equation for the amplitude a of the fundamental dipole mode shows that a behaves as the position of a heavily damped particle in a bistable potential ∝(1 ? a ²)², subject to random forcing. The dipole amplitude oscillates near the bottom of one well and makes occasional jumps to the other. These reversals are induced solely by the overtones. Theoretical expressions are derived for the statistical distribution of the dipole amplitude, the variance of the dipole amplitude between reversals, and the mean reversal rate. The model explains why the reversal rate increases with increasing secular variation, as observed. Moreover, the present reversal rate of the geodynamo, once per (2?3) × 10⁵ year, is shown to imply a secular variation of the axial dipole moment of ～ 15% (about the current value). The theoretical dipole amplitude distribution agrees well with the Sint-800 data.     

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.     

Periodic and aperiodic dynamo waves

Abstract

In order to show that aperiodic magnetic cycles, with Maunder minima, can occur naturally in nonlinear hydromagnetic dynamos, we have investigated a simple nonlinear model of an oscillatory stellar dynamo. The parametrized mean field equations in plane geometry have a Hopf bifurcation when the dynamo number D=1, leading to Parker's dynamo waves. Including the nonlinear interaction between the magnetic field and the velocity shear results in a system of seven coupled nonlinear differential equations. For D>1 there is an exact nonlinear solution, corresponding to periodic dynamo waves. In the regime described by a fifth order system of equations this solution remains stable for all D and the velocity shear is progressively reduced by the Lorentz force. In a regime described by a sixth order system, the solution becomes unstable and successive transitions lead to chaotic behaviour. Oscillations are aperiodic and modulated to give episodes of reduced activity.     

On dynamo action in the high-conductivity limit

Abstract

This paper builds on a speculation by Moffatt (1979) on an apparent conflict between two results of dynamo theory in the high conductivity limit. Firstly, the finding by Bondi and Gold (1950) on the boundedness of the magnetic dipole moment of a perfectly conducting fluid body is, for a sphere, extended to all magnetic multipole moments. Secondly, a refined version is considered of the simple spherical mean-field dynamo model proposed by Krause and Steenbeck (1967). Some constraints on the mean electromotive force near the boundary of the conducting body are taken into account, which have not been recognized up to now. In the framework of the second order correlation approximation it is shown that it is just these constraints that ensure the boundedness of the magnetic multipole moments in the high conductivity limit. Thus the apparent conflict is resolved. In this context another possible source of error in mean-field dynamo models is pointed out. The present theory also adds insight into dynamo process in cosmical objects, in a way that is briefly discussed.     

Three-dimensional domino model in geodynamo

The Lagrangian formalism is applied to consider temporal evolution of the ensemble of interacting magnetohydrodynamical cyclones governed by Langevin-type equations in a rotating medium. This problem is relevant for fast-rotating convective objects such as the cores of planets and a number of stars, where the Rossby numbers are far below unity and the geostrophic balance of the forces takes place. The paper presents the results of modeling for both the two-dimensional (2D) case when the cyclones can rotate relative to the rotation axis of the whole system in the vertical plane, and for the case of spatial rotation by two angles. It is shown that variations in the heat flux on the outer boundary of the spherical shell modulate the frequency of the reversals of the mean dipole magnetic field, which agrees with the three-dimensional (3D) modeling of the planetary dynamo. Applications of the model for giant planets are discussed, and an explanation for some episodes in the history of the geomagnetic field in the past is suggested.     

Dipole moment scaling for convection-driven planetary dynamos

Scaling laws are derived for the time-average magnetic dipole moment in rotating convection-driven numerical dynamo models. Results from 145 dynamo models with a variety of boundary conditions and heating modes, covering a wide section of parameter space, show that the time-average dipole moment depends on the convective buoyancy flux F. Two distinct regimes are found above the critical magnetic Reynolds number for onset of dynamo action. In the first regime the external magnetic field is dipole-dominant, whereas for larger buoyancy flux or slower rotation the external field is dominated by higher multipoles and the dipole moment is reduced by a factor of 10 or more relative to the dipolar regime. For dynamos driven by basal heating, the dipole moment M increases like M  F^1/3 in the dipolar regime. Reversing dipolar dynamos tend to cluster near the multipolar transition, which is shown to depend on a local Rossby number parameter. The geodynamo lies close to this transition, suggesting an explanation for polarity reversals and the possibility of a weaker dipole earlier in Earth history. Internally heated dynamos generate smaller dipole moments overall and show a gradual transition from dipolar to multipolar states. Our scaling yields order of magnitude agreement with the dipole moments of Earth, Jupiter, Saturn, Uranus, Neptune, and Ganymede, and predicts a multipolar-type dynamo for Mercury.     

Kinematic turbulent dynamo in a shear flow

Abstract

We consider the turbulent dynamo action in a differentially rotating flow by making use of a kinematic approach when the effect of a generated magnetic field on turbulent motions is neglected. The mean electromotive force is calculated in a quasilinear approximation. Differential rotation can stretch turbulent magnetic field lines and break the symmetry of turbulence in such a way that turbulent motions become suitable for the generation of a large scale magnetic field. The presence of shear changes the type of an equation governing the mean magnetic field. Due to shear stresses the mean magnetic field can be generated by a turbulent dynamo action even in a uniform turbulence. The growth rate depends on the length scale of the mean field being faster for the field with a smaller length scale.     

Nonlinear dynamics of a stellar dynamo in a spherical shell

Abstract

A simple mean-field model of a nonlinear stellar dynamo is considered, in which dynamo action is supposed to occur in a spherical shell, and where the only nonlinearity retained is the influence of the Lorentz forces on the zonal flow field. The equations are simplified by truncating in the radial direction, while full latitudinal dependence is retained. The resulting nonlinear p.d.e.'s in latitude and time are solved numerically, and it is found that while regular dynamo wave type solutions are stable when the dynamo number D is sufficiently close to its critical value, there is a wide variety of stable solutions at larger values of D. Furthermore, two different types of dynamo can coexist at the same parameter values. Implications for fields in late-type stars are discussed.     

On the impossibility of mean-field dynamos with some spherical symmetry of the motions

Abstract

A spherical mean-field dynamo model is considered in which both the mean motion and the mean electro-motive force due to fluctuating motions show some spherical symmetry. It is shown that under some reasonable assumptions the magnetic field is bound to decay to zero.     

A shell model for a large-scale turbulent dynamo

This article addresses the interesting and important problem of large-scale magnetic field generation in turbulent flows, using a self-consistent dynamo model recently developed. The main idea of this model is to consider the induction equation for the large-scale magnetic field, integrated consistently with the turbulent dynamics at smaller scales described by a magnetohydrodynamic shell model. The questions of dynamo action threshold, magnetic field saturation, magnetic field reversals, nature of the dynamo transition and the changes of small-scale turbulence as a consequence of the dynamo onset are discussed. In particular, the stability curve obtained by the model integration is shown in a very wide range of values of the magnetic Prandtl number not yet accessible by direct numerical simulation but more realistic for natural dynamos. Moreover, from our analysis it is shown that the large-scale dynamo transition displays a hysteretic behaviour and therefore a subcritical nature. The model successfully reproduces magnetic polarity reversals, showing the capability to generate persistence times which are increasing for decreasing magnetic diffusivity. Moreover, when the system reaches a statistically stationary dynamo state, where the large-scale magnetic field can abruptly reverse its polarity (magnetic reversal state) or not, keeping the same polarity (steady state), it shows an unmistakable tendency towards the energy equipartition for the turbulence at small scale.     

A unified approach to a class of slow dynamos

Abstract

Dynamo action in a highly conducting fluid with small magnetic diffusivity η is particularly sensitive to the topology of the flow. The sites of rapid magnetic field regeneration, when they occur, appear to be located at the stagnation points or in regions where the particle paths are chaotic. Elsewhere only slow dynamo action is to be expected. Two such examples are the nearly axially symmetric dynamo of Braginsky and the generalisation to smooth velocity fields of the Ponomarenko dynamo. Here a method of solution is developed, which applies to both these examples and is applicable to other situations, where magnetic field lines are close to either closed or spatially periodic contours. Particular attention is given to field generation in the neighbourhood of resonant surfaces where growth rates may be intermediate between the slow diffusive and fast convective time scales. The method is applied to the case of the two-dimensional ABC-flows, where it is shown that such intermediate dynamo action can occur on resonant surfaces.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

Torus dynamo in the outer rings of galaxies

ABSTRACT

The magnetic fields in the inner parts of some spiral galaxies are understood quite well. Their generation is connected with the dynamo mechanism that is based on the joint action of turbulent diffusion and the α-effect. Usually the galactic dynamo is described with the so-called no-z approximation which takes into account that the galaxy disc is quite thin, with the implication that some spatial derivatives may be replaced by algebraic expressions. Some galaxies have outer rings that are situated at some distance from the galactic centre. The magnetic field can be described there also using the no-z model. As the thickness of such objects is comparable with their width, it is necessary to take into account the z-dependence of the field. We have studied the magnetic field evolution using the no-z approximation and torus dynamo model for the torus with rectangular cross-section in the axisymmetric case.     

Residual fields from extinct dynamos

Abstract

The paper explores some of the many facets of the problem of the generation of magnetic fields in convective zones of declining vigor and/or thickness. The ultimate goal of such work is the explanation of the magnetic fields observed in A-stars. The present inquiry is restricted to kinematical dynamos, to show some of the many possibilities, depending on the assumed conditions of decline of the convection. The examples serve to illustrate in what quantitative detail it will be necessary to describe the convection in order to extract any firm conclusions concerning specific stars.

The first illustrative example treats the basic problem of diffusion from a layer of declining thickness. The second adds a buoyant rise to the field in the layer. The third treats plane dynamo waves in a region with declining eddy diffusivity, dynamo coefficient, and large-scale shear. The dynamo number may increase or decrease with declining convection, with an increase expected if the large-scale shear does not decline as rapidly as the eddy diffusivity. It is shown that one of the components of the field may increase without bound even in the case that the dynamo number declines to zero.     

The flux ejection dynamo effect

Abstract

The mean-field effects of cyclonic convection become increasingly complex when the cyclonic rotation exceeds ½-π. Net helicity is not required, with negative turbulent diffusion, for instance, appearing in mirror symmetric turbulence. This paper points out a new dynamo effect arising in convective cells with strong asymmetry in the rotation of updrafts as against downdrafts. The creation of new magnetic flux arises from the ejection of reserve flux through the open boundary of the dynamo region. It is unlike the familiar α-effect in that individual components of the field may be amplified independently. Several formal examples are provided to illustrate the effect. Occurrence in nature depends upon the existence of fluid rotations of the order of π in the convective updrafts. The flux ejection dynamo may possibly contribute to the generation of field in the convective core of Earth and in the convective zone of the sun and other stars.     

Geomagnetic polarity reversals in a turbulent core

The behavior of the main magnetic field components during a polarity transition is investigated using the α²-dynamo model for magnetic field generation in a turbulent core. It is shown that rapid reversals of the dipole field occur when the helicity, a measure of correlation between turbulent velocity and vorticity, changes sign. Two classes of polarity transitions are possible. Within the first class, termed component reversals, the dipole field reverses but the toroidal field does not. Within the second class, termed full reversals, both dipole and toroidal fields reverse. Component reversals result from long term fluctuations in core helicity; full reversals result from short term fluctuations. A set of time-evolution equations are derived which govern the dipole field behavior during an idealized transition. Solutions to these equations exhibit transitions in which the dipole remains axial while its intensity decays rapidly toward zero, and is regenerated with reversed polarity. Assuming an electrical conductivity of 3 × 10⁵ mho m^?1 for the fluid core, the time interval required to complete the reversal process can be as short as 7500 years. This time scale is consistent with paleomagnetic observations of the duration of reversals. A possible explanation of the cause of reversals is proposed, in which the core's net helicity fluctuates in response to fluctuations in the level of turbulence produced by two competing energy sources—thermal convection and segregation of the inner core. Symmetry considerations indicate that, in each hemisphere, helicity generated by heat loss at the core-mantle boundary may have the opposite sign of helicity generated by energy release at the inner core boundary. Random variations in rates of energy release can cause the net helicity and the α-effect to change sign occasionally, provoking a field reversal. In this model, energy release by inner core formation tends to destabilize stationary dynamo action, causing polarity reversals.     

Self-consistent dynamo models driven by hydromagnetic instabilities

The dynamics of the Earth's core are dominated by a balance between Lorentz and Coriolis forces. Previous studies of possible (magnetostrophic) hydromagnetic instabilities in this regime have been confined to geophysically unrealistic flows and fields. In recent papers we have treated rather general fields and flows in a spherical geometry and in a computationally simple plane-layer model. These studies have highlighted the importance of differential rotation in determining the spatial structure of the instability. Here we have proceeded to use these results to construct a self-consistent dynamo model of the geomagnetic field. An iterative procedure is employed in which an α-effect is calculated from the form of the instability and is then used in a mean field dynamo model. The mean zonal field calculated there is then input back into the hydromagnetic stability problem and a new α-effect calculated. The whole procedure is repeated until the input and output zonal fields are the same to some tolerance.     

Cross helicity and related dynamo

The turbulent cross helicity is directly related to the coupling coefficients for the mean vorticity in the electromotive force and for the mean magnetic-field strain in the Reynolds stress tensor. This suggests that the cross-helicity effects are important in the cases where global inhomogeneous flow and magnetic-field structures are present. Since such large-scale structures are ubiquitous in geo/astrophysical phenomena, the cross-helicity effect is expected to play an important role in geo/astrophysical flows. In the presence of turbulent cross helicity, the mean vortical motion contributes to the turbulent electromotive force. Magnetic-field generation due to this effect is called the cross-helicity dynamo. Several features of the cross-helicity dynamo are introduced. Alignment of the mean electric-current density J with the mean vorticity Ω , as well as the alignment between the mean magnetic field B and velocity U , is supposed to be one of the characteristic features of the dynamo. Unlike the case in the helicity or α effect, where J is aligned with B in the turbulent electromotive force, we in general have a finite mean-field Lorentz force J ?×? B in the cross-helicity dynamo. This gives a distinguished feature of the cross-helicity effect. By considering the effects of cross helicity in the momentum equation, we see several interesting consequences of the effect. Turbulent cross helicity coupled with the mean magnetic shear reduces the effect of turbulent or eddy viscosity. Flow induction is an important consequence of this effect. One key issue in the cross-helicity dynamo is to examine how and how much cross helicity can be present in turbulence. On the basis of the cross-helicity transport equation, its production mechanisms are discussed. Some recent developments in numerical validation of the basic notion of the cross-helicity dynamo are also presented.