Nonlinear dynamical modeling of solar cycles using dynamo formulation with turbulent magnetic helicity

Abstract

We describe a sequence of two-dimensional numerical simulations of inflection point instability in a stably stratified shear flow near the ground. The fastest growing Kelvin-Helmholtz modes are studied in detail; in particular we investigate the growth inhibiting effect of the ground which is predicted by linear theory and the Reynolds number dependence of the process of growth to finite amplitude. We consider flows which are both above and below the critical Reynolds number (Re = 300) which has been reported by Woods (1969) to mark the boundary between flows which have turbulent final states and those which do not. A global energy budget reveals a fundamental difference in character of the finite amplitude billows in these two Reynolds number regimes. However, for relatively high Reynolds numbers (Re = 10³) we do not find any explicit evidence for secondary instability. Above the transition Reynolds number the modified mean flow induced by wave growth is characterized by a splitting of the original shear layer and of the in version in which it is embedded.     

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.     

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.     

The DRESDYN project: liquid metal experiments on dynamo action and magnetorotational instability

ABSTRACT

Magnetic fields of planets, stars and galaxies are generated by self-excitation in moving electrically conducting fluids. Once produced, magnetic fields can play an active role in cosmic structure formation by destabilising rotational flows that would be otherwise hydrodynamically stable. For a long time, both hydromagnetic dynamo action as well as magnetically triggered flow instabilities had been the subject of purely theoretical research. Meanwhile, however, the dynamo effect has been observed in large-scale liquid sodium experiments in Riga, Karlsruhe and Cadarache. In this paper, we summarise the results of liquid metal experiments devoted to the dynamo effect and various magnetic instabilities such as the helical and the azimuthal magnetorotational instability and the Tayler instability. We discuss in detail our plans for a precession-driven dynamo experiment and a large-scale Tayler–Couette experiment using liquid sodium, and on the prospects to observe magnetically triggered instabilities of flows with positive shear.     

An overshoot solar dynamo with a strong return meridional flow

The meridional circulation plays an essential role in determining the basic mechanism of the dynamo action in the case of a low eddy diffusivity. Flux-transport dynamos with strong return flow and a deep stagnation point are discussed in the case of a positive α-effect located in the overshoot layer and a rotation law consistent with helioseismology. By means of a linear dynamo model, it will be shown that the migration of the toroidal belts at lower latitudes and the periods of the activity cycles are consistent with the observations. Moreover, at variance with previous investigations, the typical critical dynamo numbers of dipolar solutions are significantly smaller than those of quadrupolar solutions even in the regime of strong flow.     

Energy requirement for a working dynamo

There has for many years been interest in finding necessary conditions for dynamo action. These are usually expressed in terms of bounds on integrated properties of the flow. The bounds can clearly be improved when the flow structure can be taken into account. Recent research presents techniques for finding optimised dynamos (that is with the lowest dynamo threshold) subject to constraints, (e.g. with fixed mean square vorticity). It is natural to ask if such an optimum solution can exist when the mean square velocity is fixed. The aim of this note is to show that this is not the case and in fact that a steady or periodic dynamo can exist in a bounded conductor with an arbitrarily small value of the kinetic energy.     

A generalized two-disk dynamo model

Abstract

A generalized two-disk dynamo model is considered that includes mechanical friction; this model is intended to simulate in its broad character the behavior of the geodynamo. Fixed points, limit cycles and chaotic attractors are located for different input parameters of the model. The chaotic regimes are of several kinds as are the “routes to chaos”. Several approximate models, helpful for studying the dynamo are discussed. A number of essential differences from the well-known Rikitake dynamo are demonstrated.     

On magnetic boundary conditions for non-spectral dynamo simulations

We address issues associated with non-local magnetic boundary conditions for non-spectral dynamo simulations. We introduce an integro-differential formulation for a domain bounded by an insulating outer domain. We show how to combine the flexibility of a local discretisation with a rigorous formulation of magnetic boundary conditions in arbitrary geometries. This formulation substantiates from mathematical point of view a new method for numerical solution of magnetohydrodynamic problems with non-local boundary conditions based on coupling finite volumes and boundary elements. Finally, we discuss practical efficiency of this new method.     

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.     

Broken ergodicity,magnetic helicity,and the MHD dynamo

We consider an unforced, incompressible, turbulent magnetofluid constrained by concentric inner and outer spherical surfaces. We define a model system in which normal components of the velocity, magnetic field, vorticity, and electric current are zero on the boundaries. This choice allows us to find a set of Galerkin expansion functions that are common to both velocity and magnetic field, as well as vorticity and current. The model dynamical system represents magnetohydrodynamic (MHD) turbulence in a spherical domain and is analyzed by the methods similar to those applied to homogeneous MHD turbulence. We find a statistical theory of ideal (i.e. no dissipation) MHD turbulence analogous to that found in the homogeneous case, including the prediction of coherent structure in the form of a large-scale quasistationary magnetic field. This MHD dynamo depends on broken ergodicity, an effect that is enhanced when total magnetic helicity is increased relative to total energy. When dissipation is added and large scales are only weakly damped, quasiequilibrium may occur for long periods of time, so that the ideal theory is still pertinent on a global scale. Over longer periods of time, the selective decay of energy over magnetic helicity further enhances the effects of broken ergodicity. Thus, broken ergodicity is an essential mechanism and relative magnetic helicity is a critical parameter in this model MHD dynamo theory.     

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.     

Reversal sequence in a multiple scale dynamo mechanism

We investigate the temporal evolution of the magnetic dipole field intensity of the Earth through a multiscale dynamo mechanism. On a large range of spatio-temporal scales, the helical motions of the fluid flow are given by a schematic model of a fully developed turbulence. The system construction is symmetric with respect to left-handed and right-handed cyclones. The multiscale cyclonic turbulence coupled with a differential rotation (schematic ω dynamo) or alone (schematic ² dynamo) is the ingredient of the loop through which poloidal and toroidal fields are built from one another. Two kinds of reaction of the magnetic field on the flow are considered: the presence of a magnetic field first favours a larger-scale organization of the flow, and, second, impedes this flow by the effect of growing Lorentz forces. We obtain the general features of the geomagnetic field intensity observed over geological times and describe a general mechanism for reversals, excursions and secular variation. The mechanism happens to keep a memory during the chrons and loose it during the events (excursions and inversions).     

Magnetic Helicity in fast dynamos

The behaviour of magnetic helicity in kinematic dynamos at large magnetic Reynolds number is considered. Hughes, et al . [ Phys. Lett. A 223 , 167-172 (1996)] observe that the relative helicity tends to zero in the limit of large magnetic Reynolds number. This paper gives upper bounds on the helicity, by relating the helicity spectrum to the energy spectrum. These bounds are confirmed by numerical simulation and the distribution of helicity over scales is considered. Although it is found that the total helicity becomes small in the limit of high conductivity, there can remain significant, but cancelling, helicity at large and small scales of the field. This is illustrated by considering the evolution of helicity in the stretch-twist-fold dynamo picture.     

On the inverse dynamo problem in a simplified mean-field model

Inverse dynamo theory seeks to gain information about the motion of a liquid conductor from measurements of the magnetic field in the surrounding vacuum. We consider here a highly simplified model problem, namely a steady α²-dynamo in plane geometry with an α-field varying only in the z-direction normal to the conductor–vacuum interface. Based on perturbation theory about constant-α solutions, we find as many integral conditions on α(z) as modes are present in the vacuum field. This result is corroborated by the complete solution of a special case.     

A fully implicit model for simulating dynamo action in a Cartesian domain

We present a fully implicit numerical method to solve the incompressible MHD equations in a strongly rotating Cartesian domain. The equations are solved in a primitive variable formulation using a finite volume discretization. In order to use massively parallel computers, we applied a domain decomposition approach in space. The performance of this model is compared with an earlier model, which treated the convective terms of the equations in an explicit manner. Our results indicate that although the fully implicit method needs about three times the memory of the implicit–explicit method, it is superior in terms of computational efficiency. As an application of this model, we investigated the influence of the Prandtl number in the range of 0.01–1000 on the dynamics of the dynamo.     

The solar magnetic flux mid-term periodicities and the solar dynamo

We investigate here the fluctuations in the total, open and closed solar magnetic flux (SMF) for the period 1971–1999 by means of the maximum entropy method in the frequency range 5×10⁻⁹–10⁻⁷ Hz (6 yr to 120 days). We use monthly data for the total, open and closed magnetic solar fluxes. Periodicities found in the series are similar showing that there is some relationship between the fluxes. The most important finding of this work is the existence of fluctuations at around 1.3 and 1.7 yr in the SMF with alternating importance during consecutive even and odd solar cycles. These fluctuations are directly related with variations present in cosmic rays, solar wind parameters and geomagnetic activity indexes. A quasi-triennial periodicity previously found in sunspots and other solar phenomena is also of importance. The SMF is generated by the action of the solar dynamo; therefore, it is through the magnetic flux that the solar dynamo influences several heliospheric phenomena.