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.     

Magnetic helicity generation in the frame of Kazantsev model

Using a magnetic dynamo model, suggested by Kazantsev (J. Exp. Theor. Phys. 1968, vol. 26, p. 1031), we study the small-scale helicity generation in a turbulent electrically conducting fluid. We obtain the asymptotic dependencies of dynamo growth rate and magnetic correlation functions on magnetic Reynolds numbers. Special attention is devoted to the comparison of a longitudinal correlation function and a function of magnetic helicity for various conditions of asymmetric turbulent flows. We compare the analytical solutions on small scales with numerical results, calculated by an iterative algorithm on non-uniform grids. We show that the exponential growth of current helicity is simultaneous with the magnetic energy for Reynolds numbers larger than some critical value and estimate this value for various types of asymmetry.     

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.     

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.     

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.     

Scaling Law for Galactic Magnetic Fields

A nonlinear mean field dynamo in turbulent disks and spherical shells is discussed. We use a nonlinearity in the dynamo which includes the effect of delayed back-reaction of the mean magnetic field on the magnetic part of the — effect. This effect is determined by an evolutionary equation. The axisymmetric case is considered. An analytical expression (in a single-mode approximation) is derived which gives the magnitude of the mean magnetic field as a function of rotation and the parameters for turbulent disks. The value obtained for the mean magnetic field is in agreement with observations for galaxies.     

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.     

Galactic Dynamo and Spiral Arms - 3D MHD Simulations

In the present project we investigate the evolution of a three-dimensional (3D), large-scale galactic magnetic field under the influence of gas flows in spiral arms and in the presence of dynamo action. Our principal goal is to check how the dynamical evolution of gaseous spiral arms affects the global magnetic field structure and to what extent our models could explain the observed spiral patterns of polarization B-vectors in nearby galaxies. A two-step scheme is used: the N-body simulations of a two-component, self-gravitating disk provide the time-dependent velocity fields which are then used as the input to solve the mean-field dynamo equations. We found that the magnetic field is directly influenced by large-scale non-axisymmetric density wave flows yielding the magnetic field locally well-aligned with gaseous spiral arms in a manner similar to that discussed already by Otmianowska-Mazur et al. 1997. However, an additional field amplification, introduced by a non-zero -term in the dynamo equations, is required to cause a systematic increase of magnetic energy density against the diffusive losses. Our simulated magnetic fields are also used to construct the models of a high-frequency (Faraday rotation-free) polarized radio emission accounting for effects of projection and limited resolution, thus suitable for direct comparisons with observations.     

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.     

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.     

Galactic Spiral Arms and Dynamo Control Parameters

We discuss the effects of galactic spiral arms on the -coefficient, turbulent diffusivity and turbulent energy density of the interstellar turbulence. We argue that the -coefficient and the dynamo number are larger in the interarm regions, whereas the kinetic energy density of turbulence is larger in the arms; the turbulent magnetic diffusivity can be only weakly affected by the spiral pattern.     

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.     

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.     

Turbulent “Dynamo” in the Passive Scalar Problem

A turbulent magnetic dynamo can be considered as the evolution of a vector field in a turbulent fluid flow. The problem of evolution of scalar fields (e.g., number density of small particles) in a turbulent fluid flow is similar to the turbulent magnetic dynamo. The dynamo instability results in generation of magnetic field. The most important effect which can cause a generation of mean magnetic field in a turbulent fluid flow is the -effect: = – (1/3) u · ( × u), where u is the turbulent velocity field with the correlation time . A similar instability in the passive scalar problem results in formation of large-scale inhomogeneous structures in a spatial distribution of particles due to the -effect: = u_p ( · u_p), where u _p is the random velocity field of the particles which they acquire in a turbulent fluid velocity field. The effect is caused by inertia of particles which results in divergent velocity field of the particles. This results in additional turbulent nondiffusive flux of particles. The mean-field dynamics of inertial particles are studied by considering the stability of the equilibrium solution of the derived evolution equation for the mean number density of the particles in the limit of large Péclet numbers. The resulting equation is reduced to an eigenvalue problem for a Schrödinger equation with a variable mass, and a modified Rayleigh-Ritz variational method is used to estimate the lowest eigenvalue (corresponding to the growth rate of the instability). This estimate is in good agreement with obtained numerical solution of the Schrödinger equation. Similar effects arise during turbulent transport of gaseous admixtures (or light noninertial particles) in a low-Mach-number compressible fluid flow. The discussed effects are important in planetary and atmospheric physics (cloud formation, pollutant dynamics, preferential concentration of particles in protoplanetary disks and also planetesimals in them).     

Kinematic dynamo action in a sphere: Effects of periodic time-dependent flows on solutions with axial dipole symmetry

Choosing a simple class of flows, with characteristics that may be present in the Earth's core, we study the ability to generate a magnetic field when the flow is permitted to oscillate periodically in time. The flow characteristics are parameterised by D, representing a differential rotation, M, a meridional circulation, and C, a component characterising convective rolls. The dynamo action of all solutions with fixed parameters (steady flows) is known from earlier studies. Dynamo action is sensitive to these flow parameters and fails spectacularly for much of the parameter space where magnetic flux is concentrated into small regions, leading to high diffusion. In addition, steady flows generate only steady or regularly reversing oscillatory fields and cannot therefore reproduce irregular geomagnetic-type reversal behaviour. Oscillations of the flow are introduced by varying the flow parameters in time, defining a closed orbit in the space ( D,?M ). When the frequency of the oscillation is small, the net growth rate of the magnetic field over one period approaches the average of the growth rates for steady flows along the orbit. At increased frequency time-dependence appears to smooth out flux concentrations, often enhancing dynamo action. Dynamo action can be impaired, however, when flux concentrations of opposite signs occur close together as smoothing destroys the flux by cancellation. It is possible to produce geomagnetic-type reversals by making the orbit stray into a region where the steady flows generate oscillatory fields. In this case, however, dynamo action was not found to be enhanced by the time-dependence. A novel approach is being taken to solve the time-dependent eigenvalue problem where, by combining Floquet theory with a matrix-free Krylov-subspace method, we can avoid large memory requirements for storing the matrix required by the standard approach.     

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.     

On the modified Taylor constraint

A dynamo driven by motions unaffected by viscous forces is termed magnetostrophic. Although such a model might describe magnetic field generation in Earth’s core well, a magnetostrophic dynamo has not yet been found even though Taylor [Proc. R. Soc. Lond. A 1963, 274, 274–283] devised an apparently viable method of finding one. His method for determining the fluid velocity from the magnetic field and the energy source involved only the evaluation of integrals along lines parallel to the Earth’s axis of rotation and the solution of a second-order ordinary differential equation. It is demonstrated below that an approximate solution of this equation for a broad family of magnetic fields is immediate. Furthermore inertia, which was neglected in Taylor’s theory, is restored here, so that the modified theory includes torsional waves, whose existence in the Earth’s core has been inferred from observations of the length of day. Their theory is reconsidered.     

Broken ergodicity in magnetohydrodynamic turbulence

Turbulent magnetofluids appear in various geophysical and astrophysical contexts, in phenomena associated with planets, stars, galaxies and the universe itself. In many cases, large-scale magnetic fields are observed, though a better knowledge of magnetofluid turbulence is needed to more fully understand the dynamo processes that produce them. One approach is to develop the statistical mechanics of ideal (i.e. non-dissipative), incompressible, homogeneous magnetohydrodynamic (MHD) turbulence, known as “absolute equilibrium ensemble” theory, as far as possible by studying model systems with the goal of finding those aspects that survive the introduction of viscosity and resistivity. Here, we review the progress that has been made in this direction. We examine both three-dimensional (3-D) and two-dimensional (2-D) model systems based on discrete Fourier representations. The basic equations are those of incompressible MHD and may include the effects of rotation and/or a mean magnetic field B _o. Statistical predictions are that Fourier coefficients of the velocity and magnetic field are zero-mean random variables. However, this is not the case, in general, for we observe non-ergodic behavior in very long time computer simulations of ideal turbulence: low wavenumber Fourier modes that have relatively large means and small standard deviations, i.e. coherent structure. In particular, ergodicity appears strongly broken when B _o?=?0 and weakly broken when B _o?≠?0. Broken ergodicity in MHD turbulence is explained by an eigenanalysis of modal covariance matrices. This produces a set of modal eigenvalues inversely proportional to the expected energy of their associated eigenvariables. A large disparity in eigenvalues within the same mode (identified by wavevector k ) can occur at low values of wavenumber k?=?| k |, especially when B _o?=?0. This disparity breaks the ergodicity of eigenvariables with smallest eigenvalues (largest energies). This leads to coherent structure in models of ideal homogeneous MHD turbulence, which can occur at lowest values of wavenumber k for 3-D cases, and at either lowest or highest k for ideal 2-D magnetofluids. These ideal results appear relevant for unforced, decaying MHD turbulence, so that broken ergodicity effects in MHD turbulence survive dissipation. In comparison, we will also examine ideal hydrodynamic (HD) turbulence, which, in the 3-D case, will be seen to differ fundamentally from ideal MHD turbulence in that coherent structure due to broken ergodicity can only occur at maximum k in numerical simulations. However, a nonzero viscosity eliminates this ideal 3-D HD structure, so that unforced, decaying 3-D HD turbulence is expected to be ergodic. In summary, broken ergodicity in MHD turbulence leads to energetic, large-scale, quasistationary magnetic fields (coherent structures) in numerical models of bounded, turbulent magnetofluids. Thus, broken ergodicity provides a large-scale dynamo mechanism within computer models of homogeneous MHD turbulence. These results may help us to better understand the origin of global magnetic fields in astrophysical and geophysical objects.     

Spherical dynamos with anisotropic α-effect

Abstract

The mean field induction equation of Steenbeck, Krause and Rädler (1966) is solved for anisotropic α _ik -tensors of varying anisotropy. The attention is restricted to α²-dynamos in spheres with steady axisymmetric magnetic fields. The ratio of the electrical conductivity outside and inside the sphere is varied, but in all cases it is found that a steady dynamo does not exist when the anisotropy of the α _ik -tensor exceeds a critical value. Such a critical value does not exist in the exceptional case of the Fermi boundary condition. The results emphasize the important effect of boundaries on the existence of solutions of the dynamo problem.