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.     

Kinematic dynamo action with helical symmetry in an unbounded fluid conductor: Part 2. Further development of an explicit solution for the prototype case of lortz

Abstract

An explicit example of a steady prototype Lortz dynamo is elaborated in terms of a previously derived illustrative, exact, closed form solution to the nonlinear dynamo equations. The eigenvalue character of the dynamo problem is now introduced which simplifies the solution. The magnetic field lines, which lie on circular cylinders, and velocity streamline pattern are then displayed and discussed. Analysis of the magnetic energy balance by way of the Poynting flux reveals the existence of a finite critical cylinder across which zero net magnetic energy flows, thereby proving that the material inside is a self-excited dynamo, despite the fact that the total magnetic energy is unbounded.     

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.     

Models for convectively driven hydromagnetic dynamos

Abstract

Models of a convectively driven hydromagnetic dynamo are constructed using a truncated modal expansion. The resulting nonlinear partial differential equations are integrated numerically. The results confirm that rotation is a necessary condition for effective dynamo action, and suggest that equipartition of kinetic and magnetic energies is qualitatively valid, and that toroidal field energies can be much larger than poloidal.     

Dynamo in Asymmetric Square Convection

Linear and nonlinear dynamo action is investigated for square patterns in nonrotating and weakly rotating Boussinesq Rayleigh-Bénard convection in a plane horizontal layer. The square-pattern solutions may or may not be symmetric to up-down reflections. Vertically symmetric solutions correspond to checkerboard patterns. They do not possess a net kinetic helicity and are found to be incapable of kinematic dynamo action at least up to magnetic Reynolds numbers of , 12 000. There also exist vertically asymmetric squares, characterized by rising (descending) motion in the centers and descending (rising) motion near the boundaries, among them such that possess full horizontal square symmetry and others lacking also this symmetry. The flows lacking both the vertical and horizontal symmetries possess kinetic helicity and show kinematic dynamo action even without rotation. The generated magnetic fields are concentrated in vertically oriented filamentary structures. Without rotation these dynamos are, however, always only kinematic, not nonlinear dynamos since the back-reaction of the magnetic field then forces the solution into the basin of attraction of a roll pattern incapable of dynamo action. But with rotation added parameter regions are found where stationary asymmetric squares are also nonlinear dynamos. These nonlinear dynamos are characterized by a subtle balance between the Coriolis and Lorentz forces. In some parameter regions also nonlinear dynamos with flows in the form of oscillating squares or stationary modulated rolls are found.     

Large-scale convective dynamos in a stratified rotating plane layer

We study the effect of stratification on large-scale dynamo action in convecting fluids in the presence of background rotation. The fluid is confined between two horizontal planes and both boundaries are impermeable, stress-free and perfectly conducting. An asymptotic analysis is performed in the limit of rapid rotation (τ???1 where τ is the Taylor number). We analyse asymptotic magnetic dynamo solutions in rapidly rotating systems generalising the results of Soward [A convection-driven dynamo I. The weak field case. Philos. Trans. R. Soc. Lond. A 1974, 275, 611–651] to include the effects of compressibility. We find that in general the presence of stratification delays the efficiency of large-scale dynamo action in this regime, leading to a reduction of the onset of dynamo action and in the nonlinear regime a diminution of the large-scale magnetic energy for flows with the same kinetic energy.     

Convective dynamos with intermediate and strong fields

Abstract

The weak-field Benard-type dynamo treated by Soward is considered here at higher levels of the induced magnetic field. Two sources of instability are found to occur in the intermediate field regime M ～ T ^1/12, where M and T are the Hartmann and Taylor numbers. On the time scale of magnetic diffusion, solutions may blow up in finite time owing to destabilization of the convection by the magnetic field. On a faster time scale a dynamic instability related to MAC-wave instability can also occur. It is therefore concluded that the asymptotic structure of this dynamo is unstable to virtual increases in the magnetic field energy.

In an attempt to model stabilization of the dynamo in a strong-field regime we consider two approximations. In the first, a truncated expansion in three-dimensional plane waves is studied numerically. A second approach utilizes an ad hoc set of ordinary differential equations which contains many of the features of convection dynamos at all field energies. Both of these models exhibit temporal intermittency of the dynamo effect.     

Fast dynamo action in the Ponomarenko dynamo

Abstract

An analysis of small-scale magnetic fields shows that the Ponomarenko dynamo is a fast dynamo; the maximum growth rate remains of order unity in the limit of large magnetic Reynolds number. Magnetic fields are regenerated by a “stretch-diffuse” mechanism. General smooth axisymmetric velocity fields are also analysed; these give slow dynamo action by the same mechanism.     

Kinematic dynamos associated with large scale fluid motions

Abstract

A standard approach to the kinematic dynamo problem is that pioneered by Bullard and Gellman (1954), which utilizes the toroidal-poloidal separation and spherical harmonic expansion of the magnetic and velocity fields. In these studies, the velocity field is given as a combination of small number of toroidal and poloidal harmonics, with their radial dependences prescribed by some physical considerations. Starting from the original paper of Bullard and Gellman (1954), a number of authors repeated such analyses on different combination of velocity fields, including the most recent and comprehensive effort by Dudley and James (1989). In this paper, we re-examine the previous kinematic dynamo models, using the computer algebra approach initiated by Kono (1990). This method is particularly suited to this kind of research since different velocity fields can be treated by a single program. We used the distribution of magnetic energies in various harmonics to infer the convergence of the results.

The numerical results obtained in this study for the models of Bullard and Gellman (1954), Lilley (1970), Gubbins (1973), Pekeris et al. (1973), Kumar and Roberts (1975), and Dudley and James (1989) are consistent with the previously reported results, in particular, with the extensive calculation of Dudley and James. In addition, we found that the combination of velocities used by Lilley can support the dynamo action if the radial dependence of the velocity is modified.

We also examined the helicity distributions in these dynamo models, to see if there is any correlation between the helicity and the efficiency of dynamo action. A successful dynamo can result both from the cases in which the helicity distributions are symmetric or antisymmetric with respect to the equator. In both cases, it appears that the dynamo action is efficient if the volume integral of helicity over a hemisphere is large.     

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.     

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.     

The galactic dynamo: Axisymmetric and non-axisymmetric modes

Abstract

We discuss recent developments in the theory of large-scale magnetic structures in spiral galaxies. In addition to a review of galactic dynamo models developed for axisymmetric disks of variable thickness, we consider the possibility of dominance of non-axisymmetric magnetic modes in disks with weak deviations from axial symmetry. Difficulties of straightforward numerical simulation of galactic dynamos are discussed and asymptotic solutions of the dynamo equations relevant for galactic conditions are considered. Theoretical results are compared with observational data.     

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.     

Evolution of aligned states within nonlinear dynamos

ABSTRACT

The Archontis dynamo is a rare example of an MHD dynamo within which forcing drives a dynamo where the flow and magnetic fields are almost perfectly aligned and the energies are approximately equal. In this paper, I expand upon our knowledge of the dynamo by showing that the intermediate steady states of the kinetic and magnetic energies observed by Cameron and Galloway are not a necessary feature of aligned dynamos. Furthermore, I show that the steady state into which the flow and magnetic fields eventually evolve is remarkably robust to the addition of time dependence and asymmetry to the forcing.     

Dynamo models with magnetic diffusivity-quenching

Abstract

The magnetic influence on a turbulent plasma also produces a complicated structure of the eddy diffusivity tensor rather than a simple and traditional quenching of the eddy diffusivity. Dynamo models in plane (galaxy) and spherical (star) geometries with differential relation are developed here to answer the question whether the dynamo mechanism is still yielding stable configurations. Magnetic saturation of the dynamos is always found—at magnetic energies exceeding the energy-equipartition value.

We find that the effect of magnetic back-reaction on the turbulent diffusivity depends highly on whether the dynamo is oscillatory or not. The steady modes are extremely influenced. They saturate at field strengths strongly exceeding its turbulence-equipartition value. Subcritical excitation is even found for strong seed fields. The oscillating dynamos, however, only provide a small effect. Hence, the strong over-equipartition of the internal solar magnetic fields revealed by studies of flux-tube dynamics cannot be explained with the theory presented. Also the run of the cycle frequency with the dynamo number is too smooth in order to explain observations of stellar chromospheric activity.     

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.     

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.     

Mechanisms for magnetic field generation in precessing cubes

ABSTRACT

It is shown that flows in precessing cubes develop at certain parameters large axisymmetric components in the velocity field which are large enough to either generate magnetic fields by themselves, or to contribute to the dynamo effect if inertial modes are already excited and acting as a dynamo. This effect disappears at small Ekman numbers. The critical magnetic Reynolds number also increases at low Ekman numbers because of turbulence and small-scale structures.     

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.     

The Dynamical Effects of Hyperviscosity on Numerical Geodynamo Models

The most fundamental difficulty in the construction of an Earth-like dynamo model is associated with the constraint caused by the rapid rotation of the Earth. To stabilise numerical codes, many workers have introduced hyperviscosity into the governing equations. One of the major effects introduced by hyperviscosity is to offset the rotational constraint, and, consequently, to alter the key dynamics of an Earth-like dynamo. In this paper, an Earth-like convection model with or without the presence of an imposed magnetic field is investigated with or without the effect of hyperviscosity. A nonlinear dynamo model with the mean field approximation is also used to examine the dynamical effect of hyperviscosity. The results suggest that great care should be taken when hyperviscosity is employed in geodynamo models.