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.     

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.     

Necessary conditions for the magnetohydrodynamic dynamo

Abstract

A magnetohydrodynamic, dynamo driven by convection in a rotating spherical shell is supposed to have averages that are independent of time. Two cases are considered, one driven by a fixed temperature difference R and the other by a given internal heating rate Q. It is found that when q, the ratio of thermal conductivity to magnetic diffusivity, is small, R must be of order q ^?4/3 and Q of order q ^?2 for dynamo action to be possible; q is small in the Earth's core, so it is hoped that the criteria will prove useful in practical as well as theoretical studies of dynamic dynamos. The criteria can be further strengthened when the ohmic dissipation of the field is significant in the energy balance. The development includes the derivation of two necessary conditions for dynamo action, both based on the viscous dissipation rate of the velocity field that drives the dynamo.     

The weak taylor state in an αω-dynamo

Abstract

A spherical αω-dynamo is studied for small values of the viscous coupling parameter ε ～ v^1/2, paying attention particularly to large dynamo numbers. The present study is a follow-up of the work by Hollerbach et al. (1992) with their choice of α-effect and Archimedean wind including also the constraint of magnetic field symmetry (or antisymmetry) due to equatorial plane. The magnetic field scaled by ε^1/2 is independent of ε in the solutions for dynamo numbers smaller than a certain value of D _b (the Ekman state) which are represented by dynamo waves running from pole to equator or vice-versa. However, for dynamo numbers larger than D _b the solution bifurcates and subsequently becomes dependent on ε. The bifurcation is a consequence of a crucial role of the meridional convection in the mechanism of magnetic field generation. Calculations suggest that the bifurcation appears near dynamo number about 33500 and the solutions for larger dynamo numbers and ε = 0 become unstable and fail, while the solutions for small but non-zero ε are characterized by cylindrical layers of local maximum of magnetic field and sharp changes of geostrophic velocity. Our theoretical analysis allows us to conclude that our solution does not take the form of the usual Taylor state, where the Taylor constraint should be satisfied due to the special structure of magnetic field. We rather obtained the solution in the form of a “weak” Taylor state, where the Taylor constraint is satisfied partly due to the amplitude of the magnetic field and partly due to its structure. Calculations suggest that the roles of amplitude and structure are roughly fifty-fifty in our “weak” Taylor state solution and thus they can be called a Semi-Taylor state. Simple estimates show that also Ekman state solutions can be applicable in the geodynamo context.     

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.     

A nonlinear dynamo driven by rapidly rotating convection

In Kim et al. (Kim, E., Hughes, D.W. and Soward, A.M., “An investigation into high conductivity dynamo action driven by rotating convection”, Geophys. Astrophys. Fluid Dynam. 91, 303–332 ().) we investigated kinematic dynamo action driven by rapidly rotating convection in a cylindrical annulus. Here we extend this work to consider self-consistent nonlinear dynamo action in which the back-reaction of the Lorentz force on the flow is taken into account. In particular, we investigate, as a function of magnetic Prandtl number, the evolution of an initially weak magnetic field in two different types of convective flow – one chaotic and the other integrable. On saturation, the latter shows a systematic dependence on the magnetic Prandtl number whereas the former appears not to. In addition, we show how, in keeping with the findings of Cattaneo et al. (Cattaneo, F., Hughes, D.W. and Kim, E., “Suppression of chaos in a simplified nonlinear dynamo model”, Phys. Rev. Lett. 76, 2057–2060 ().), saturation of the growth of the magnetic field is brought about, for the originally chaotic flow, by a strong suppression of chaos.     

Penetrative convection in rotating fluids: A model for the base of stellar convection zones

Abstract

To model penetrative convection at the base of a stellar convection zone we consider two plane parallel, co-rotating Boussinesq layers coupled at their fluid interface. The system is such that the upper layer is unstable to convection while the lower is stable. Following the method of Kondo and Unno (1982, 1983) we calculate critical Rayleigh numbers R_c for a wide class of parameters. Here, R_c is typically much less than in the case of a single layer, although the scaling R_c～T^2/3 as T → ∞ still holds, where T is the usual Taylor number. With parameters relevant to the Sun the helicity profile is discontinuous at the interface, and dominated by a large peak in a thin boundary layer beneath the convecting region. In reality the distribution is continuous, but the sharp transition associated with a rapid decline in the effective viscosity in the overshoot region is approximated by a discontinuity here. This source of helicity and its relation to an alpha effect in a mean-field dynamo is especially relevant since it is a generally held view that the overshoot region is the location of magnetic field generation in the Sun.     

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.     

A Route to Magnetic Field Reversals: An Example of an ABC-Forced Nonlinear Dynamo

We are investigating numerically the nonlinear behaviour of a space-periodic MHD system with ABC forcing. Most computations are performed for magnetic Reynolds numbers increasing from 0 to 60 and a fixed kinematic Reynolds number, small enough for the trivial solution with a zero magnetic field to be stable to velocity perturbations. At the critical magnetic Reynolds number for the onset of instability of the trivial solution the dominant eigenvalue of the kinematic dynamo problem is real. In agreement with the bifurcation theory new steady states with non-vanishing magnetic field appear in this bifurcation. Subsequent bifurcations are investigated. A regime is detected, where chaotic variations of the magnetic field orientation (analogous to magnetic field reversals) are observed in the temporal evolution of the system.     

Oscillatory large-scale dynamos from Cartesian convection simulations

We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo is expected to be of α² Ω type, which is generally expected to be relevant for the Sun, whereas the latter case refers to α² dynamos that are more likely to occur in more rapidly rotating stars whose differential rotation is small. We perform a parameter study where the shear flow and the rotational influence are varied to probe the relative importance of both types of dynamos. Oscillatory solutions are preferred both in the kinematic and saturated regimes when the negative ratio of shear to rotation rates, q?≡??S/Ω, is between 1.5 and 2, i.e. when shear and rotation are of comparable strengths. Other regions of oscillatory solutions are found with small values of q, i.e. when shear is weak in comparison to rotation, and in the regime of large negative qs, when shear is very strong in comparison to rotation. However, exceptions to these rules also appear so that for a given ratio of shear to rotation, solutions are non-oscillatory for small and large shear, but oscillatory in the intermediate range. Changing the boundary conditions from vertical field to perfect conductor ones changes the dynamo mode from oscillatory to quasi-steady. Furthermore, in many cases an oscillatory solution exists only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. However, the cases with rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic boundary conditions. The strengths of total and large-scale components of the magnetic field in the saturated state, however, are sensitive to the chosen boundary conditions.     

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.     

Towards the interpretation of Uranus's eccentric magnetic field

Abstract

Coriolis forces stimulate dynamo action in a rapidly-rotating fluid by promoting complexities in the pattern of fluid motions, notably departures from symmetry about the axis of rotation. This pattern and its time variations determine the instantaneous form and temporal behaviour of the magnetic field so produced. Instantaneous magnetic fields will usually exhibit in their broad-scale features approximate alignment with the rotation axis. This is borne out by observations of the magnetic fields of the Earth, Jupiter and Saturn, and it is likely on general grounds that Neptune will be found to have an aligned magnetic field. But, as is shown by laboratory and theoretical studies of thermal convection in rapidly-rotating fluids, for some ranges of rotation speed, rate of heating, etc. certain patterns can occur which in electrically-conducting fluids would produce magnetic fields exhibiting departures from alignment with the rotation axis, which instantaneously could be quite pronounced but would average out to very small values over sufficiently long periods of time. These findings indicate obvious strategies for theoretical studies towards the interpretation of Uranus's eccentric magnetic field (which need not invoke departures from axial symmetry in the thermal, mechanical or electrical boundary conditions of the dynamo region within the planet) and for further observational studies.     

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.     

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.     

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.     

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.     

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.     

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.     

Maximally-efficient-generation approach in the dynamo theory

Abstract

We propose a method of derivation of global asymptotic solutions of the hydromagnetic dynamo problem at large magnetic Reynolds number. The procedure reduces to matching the local asymptotic forms for the magnetic field generated near individual extrema of generation strength. The basis of the proposed method, named here the Maximally-Efficient-Generation Approach (MEGA), is the assertion that properties of global asymptotic solutions of the kinematic dynamo are determined by the distribution of the generation strength near its leading extrema and by the number and distribution of the extrema.

The general method is illustrated by the global asymptotic solution of the α²-dynamo problem in a slab. The nature of oscillatory solutions revealed earlier in numerical simulations and the reasons for the dominance of even magnetic modes in slab geometry are clarified.

Applicability of the asymptotic solutions at moderate values of the asymptotic parameter is also discussed. We confirm this applicability using comparisons with complementary asymptotic expansions and numerical simulations. In particular, this justifies application of the MEGA solutions to estimation of the generation threshold.     

Properties of nonlinear dynamo waves

Abstract

Dynamo theory offers the most promising explanation of the generation of the sun's magnetic cycle. Mean field electrodynamics has provided the platform for linear and nonlinear models of solar dynamos. However the nonlinearities included arc (necessarily) arbitrarily imposed in these models. This paper conducts a systematic survey of the role of nonlinearities in the dynamo process, by cousidering the behaviour of dynamo waves in the nonlinear regime. It is demonstrated that only by considering realistic nonlinearities that are non-local in space and time can modulation of the basic dynamo wave be achieved. Moreover this modulation is greatest when there is a large separation of timescales provided by including a low magnetic Prandtl number in the equation for the velocity perturbations.