Α2 Dynamo Research Articles

Budget equations and astrophysical non-linear mean-field dynamos

ABSTRACTSolar, stellar and galactic large-scale magnetic fields are originated due to a combined action of non-uniform (differential) rotation and helical motions of plasma via mean-field dynamos. Usually, non-linear mean-field dynamo theories take into account algebraic and dynamic quenching of alpha effect and algebraic quenching of turbulent magnetic diffusivity. However, the theories of the algebraic quenching do not take into account the effect of modification of the source of turbulence by the growing large-scale magnetic field. This phenomenon is due to the dissipation of the strong large-scale magnetic field resulting in an increase of the total turbulent energy. This effect has been studied using the budget equation for the total turbulent energy (which takes into account the feedback of the generated large-scale magnetic field on the background turbulence) for (i) a forced turbulence, (ii) a shear-produced turbulence, and (iii) a convective turbulence. As the result of this effect, a non-linear dynamo number decreases with increase of the large-scale magnetic field, so that that the mean-field αΩ, α2, and α2Ω dynamo instabilities are always saturated by the strong large-scale magnetic field.

Magnetic fields of low-mass main sequences stars: non-linear dynamo theory and mean-field numerical simulations

ABSTRACTOur theoretical and numerical analysis have suggested that for low-mass main sequences stars (of the spectral classes from M5 to G0) rotating much faster than the Sun, the generated large-scale magnetic field is caused by the mean-field α2Ω dynamo, whereby the α2 dynamo is modified by a weak differential rotation. Even for a weak differential rotation, the behaviour of the magnetic activity is changed drastically from aperiodic regime to non-linear oscillations and appearance of a chaotic behaviour with increase of the differential rotation. Periods of the magnetic cycles decrease with increase of the differential rotation, and they vary from tens to thousand years. This long-term behaviour of the magnetic cycles may be related to the characteristic time of the evolution of the magnetic helicity density of the small-scale field. The performed analysis is based on the mean-field simulations (MFS) of the α2Ω and α2 dynamos and a developed non-linear theory of α2 dynamo. The applied MFS model was calibrated using turbulent parameters typical for the solar convective zone.

On the helicity characteristics and induced emf of magnetic-Coriolis wave packets

SUMMARY In a rapidly rotating Boussinesq fluid, buoyant anomalies radiate low-frequency inertial wave packets that disperse along the rotation axis. The wave packets lead to axially elongated vortices, which propagate negative (positive) kinetic helicity upwards (downwards) with respect to the rotation vector. The kinetic helicity carried by the inertial wave packets is near-maximal relative to the velocity and vorticity fields. In classical mean-field theory, kinetic helicity is often associated with the α-effect, which is thought to be an important ingredient for planetary dynamos. The modification of inertial wave packets in the presence of a transverse uniform magnetic field is investigated here, motivated by small-scale dynamics in planetary cores, where a large-scale magnetic field affects fluid motions. We study numerically the dispersion of wave packets from an isolated buoyant source and from a random layer of buoyant anomalies, while varying the Lehnert number Le—the ratio of the frequencies of Alfvén and inertial waves. We find that for Le < 0.1, the vortices are columnar and continue to segregate kinetic helicity so that it is negative (positive) above (below) the buoyant source. Importantly, the wave packets induce an α-effect, which remains strong and coherent for Earth-like values of the Lehnert number (Le < 0.1). The interaction of wave packets emitted by multiple neighbouring buoyant sources results in an α-effect that is stronger than the α-effect induced by wave packets launched from an isolated buoyant source, and we provide an analytical explanation for this. The coherence of the α-effect induced by the wave packets, for Earth-like values of the Lehnert number, lends support to the α2 dynamo model driven by helical waves.

Large-scale dynamo action of magnetized Taylor–Couette flows

ABSTRACT A conducting Taylor–Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler–Spruit type. The flows are unstable against non-axisymmetric perturbations which form electromotive forces defining α effect and eddy diffusivity. If both degenerated modes with m = ±1 are excited with the same power then the global α effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite α effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an α2 dynamo. Moreover, an αΩ dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of Rm ≃ 2000 to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, ηT ≃ η. By basically lower or larger eddy diffusivities the dynamo instability is suppressed.

Magnetic helicity and fluxes in an inhomogeneous α2 dynamo

Much work on turbulent three‐dimensional dynamos has been done using triply periodic domains, in which there are no magnetic helicity fluxes. Here, we present simulations where the turbulent intensity is still nearly homogeneous, but now there is a perfect conductor boundary condition on one end and a vertical field or pseudovacuum conditions on the other. This leads to migratory dynamo waves. Good agreement with a corresponding analytically solvable α2 dynamo is found. Magnetic helicity fluxes are studied in both types of models. It is found that at moderate magnetic Reynolds numbers, most of the magnetic helicity losses occur on large scales. Whether this changes at even larger magnetic Reynolds numbers, as required for alleviating the catastrophic dynamo quenching problem, remains still unclear.

Generation of a Large-scale Magnetic Field in a Convective Full-sphere Cross-helicity Dynamo

Abstract We study the effects of the cross-helicity in the full-sphere large-scale mean-field dynamo models of a 0.3 M ⊙ star rotating with a period of 10 days. In exploring several dynamo scenarios that stem from magnetic field generation by the cross-helicity effect, we found that the cross-helicity provides the natural generation mechanisms for the large-scale scale axisymmetric and nonaxisymmetric magnetic field. Therefore, the rotating stars with convective envelopes can produce a large-scale magnetic field generated solely due to the turbulent cross-helicity effect (we call it γ 2-dynamo). Using mean-field models we compare the properties of the large-scale magnetic field organization that stems from dynamo mechanisms based on the kinetic helicity (associated with the α 2 dynamos) and cross-helicity. For the fully convective stars, both generation mechanisms can maintain large-scale dynamos even for the solid body rotation law inside the star. The nonaxisymmetric magnetic configurations become preferable when the cross-helicity and the α-effect operate independently of each other. This corresponds to situations with purely γ 2 or α 2 dynamos. The combination of these scenarios, i.e., the γ 2 α 2 dynamo, can generate preferably axisymmetric, dipole-like magnetic fields at strengths of several kGs. Thus, we found a new dynamo scenario that is able to generate an axisymmetric magnetic field even in the case of a solid body rotation of the star. We discuss the possible applications of our findings to stellar observations.

Parametric study and optimization trends for the Von-Kármán-sodium dynamo experiment

We present magneto-hydrodynamic simulations of liquid sodium flow performed with the PLUTO compressible MHD code. We investigated the influence of the remanent magnetic field orientation and intensity, the impinging velocity field due to Ekman pumping as well as the impeller dimensions on the magnetic field collimation by helical flows in-between the impeller blades. For a simplified Cartesian geometry, we model the flow dynamics of a multi-blade impeller inspired by the Von-Kármán-Sodium experiment. This study shows that a remanent magnetic field oriented in the toroidal direction is the less efficient configuration to collimate the magnetic field, although if the radial or vertical components are not negligible, the collimation is significantly improved. As the intensity of the remanent magnetic field increases, the system magnetic energy becomes larger, but the magnetic field collimation efficiency remains the same, so the gain of magnetic energy is smaller as the remanent magnetic field intensity increases. The magnetic field collimation is modified if the impinging velocity field changes: the collimation is weaker if the impinging velocity increases from Γ = 0.8 to 0.9 and slightly larger if the impinging velocity decreases from Γ = 0.8 to 0.7. The analysis of the impeller dimensions points out that the most efficient configuration to collimate the magnetic field requires a ratio between the impeller blade height and the base longitude between 0.375 and 0.5. The largest enhancement of the hypothetical α2 dynamo loop, compared to the hypothetical Ω-α dynamo loop, is observed for the model that mimics the TM 73 impeller configuration rotating in the unscooping direction with a remanent magnetic field of 10–3 T orientated in the radial or vertical direction. The optimization trends obtained in the parametric analysis are also confirmed by simulations with a higher resolution and turbulence degree.

Internally driven inertial waves in geodynamo simulations

© The Authors 2017. Inertial waves are oscillations in a rotating fluid, such as the Earth's outer core, which result from the restoring action of the Coriolis force. In an earlier work, it was argued by Davidson that inertial waves launched near the equatorial regions could be important for the α2 dynamo mechanism, as they can maintain a helicity distribution which is negative (positive) in the north (south). Here,we identify such internally driven inertial waves, triggered by buoyant anomalies in the equatorial regions in a strongly forced geodynamo simulation. Using the time derivative of vertical velocity, ∂uz/∂t, as a diagnostic for traveling wave fronts, we find that the horizontal movement in the buoyancy field near the equator is well correlated with a corresponding movement of the fluid far from the equator. Moreover, the azimuthally averaged spectrum of ∂uz/∂t lies in the inertial wave frequency range. We also test the dispersion properties of the waves by computing the spectral energy as a function of frequency, ω, and the dispersion angle, θ. Our results suggest that the columnar flow in the rotation-dominated core, which is an important ingredient for the maintenance of a dipolar magnetic field, is maintained despite the chaotic evolution of the buoyancy field on a fast timescale by internally driven inertial waves.

Ribbons characterize magnetohydrodynamic magnetic fields better than lines: a lesson from dynamo theory

Blackman and Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here, we present such a schematic ribbon diagram for the α2 dynamo that tracks magnetic helicity and provides distinct scales of large-scale magnetic helicity, small-scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new ‘2.5 scale’ minimalist generalization of the helicity-evolving equations for the α2 dynamo that separately allows for these three distinct length-scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large-scale field first grows at a rate independent of the magnetic Reynolds number RM before quenching to an RM-dependent regime. But we also show that the larger the ratio of the wavenumber where the small-scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large-scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in magnetohydrodynamic are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines.

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 α2 Ω type, which is generally expected to be relevant for the Sun, whereas the latter case refers to α2 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.

Growth rate of magnetic field during starburst phase in a galaxy

AbstractDuring an active star formation phase, it is conceivable that the turbulent level of a galaxy increases and consequently an increase in magnetic field by turbulent dynamo process. We point out that the thickness of the dynamo region is sensitive to the turbulent level. We examine the linear growth rate of three different dynamos, namely, αω, α2ω and α2 dynamos. We find that the dependence of the growth rate on turbulent level is quite different for different dynamos. Moreover, the growth rate is not necessary monotonic to turbulence level even in the linear regime.

CATASTROPHIC QUENCHING IN αΩ DYNAMOS REVISITED

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astrophysical large-scale dynamos. The recognition of this fact led to the development of the dynamical α quenching formalism, which predicts catastrophically low mean fields in open systems. Here, we show that in oscillatory αΩ dynamos this formalism predicts an unphysical magnetic helicity transfer between scales. An alternative technique is proposed where this artifact is removed by using the evolution equation for the magnetic helicity of the total field in the shearing advective gauge. In the traditional dynamical α quenching formalism, this can be described by an additional magnetic helicity flux of small-scale fields that does not appear in homogeneous α2 dynamos. In αΩ dynamos, the alternative formalism is shown to lead to larger saturation fields than what has been obtained in some earlier models with the traditional formalism. We have compared the predictions of the two formalisms to results of direct numerical simulations, finding that the alternative formulation provides a better fit. This suggests that worries about catastrophic dynamo behavior in the limit of large magnetic Reynolds number are unfounded.

Scaling and intermittency in incoherent α-shear dynamo

We consider mean-field dynamo models with fluctuating α effect, both with and without large-scale shear. The α effect is chosen to be Gaussian white noise with zero mean and a given covariance. In the presence of shear, we show analytically that (in infinitely large domains) the mean-squared magnetic field shows exponential growth. The growth rate of the fastest growing mode is proportional to the shear rate. This result agrees with earlier numerical results of Yousef et al. and the recent analytical treatment by Heinemann, McWilliams & Schekochihin who use a method different from ours. In the absence of shear, an incoherent α2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third- and fourth-order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian.

Modeling spatio‐temporal nonlocality in mean‐field dynamos

AbstractWhen scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α2 dynamos with α linearly depending on position, which makes this dynamo oscillatory, too. In both cases, the critical values of the dynamo number is lowered due to spatio‐temporal nonlocality.When scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α2 dynamos

Magnetic helicity fluxes and their effect on stellar dynamos

AbstractMagnetic helicity fluxes in turbulently driven α2 dynamos are studied to demonstrate their ability to alleviate catastrophic quenching. A one-dimensional mean-field formalism is used to achieve magnetic Reynolds numbers of the order of 105. We study both diffusive magnetic helicity fluxes through the mid-plane as well as those resulting from the recently proposed alternate dynamic quenching formalism. By adding shear we make a parameter scan for the critical values of the shear and forcing parameters for which dynamo action occurs. For this αΩ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.

A SHELL MODEL TURBULENT DYNAMO

Turbulent dynamo phenomena, observed almost everywhere in astrophysical objects and also in the laboratory in the recent VKS2 experiment, are investigated using a shell model technique to describe magnetohydrodynamic turbulence. Detailed numerical simulations at very high Rossby numbers (α2 dynamo) show that as the magnetic Reynolds number increases, the dynamo action starts working and different regimes are observed. The model, which displays different large-scale coherent behaviors corresponding to different regimes, is able to reproduce the magnetic field reversals observed both in a geomagnetic dynamo and in the VKS2 experiment. While rough quantitative estimates of typical times associated with the reversal phenomenon are consistent with paleomagnetic data, the analysis of the transition from oscillating intermittent through reversal and finally to stationary behavior shows that the nature of the reversals we observe is typical of α2 dynamos and completely different from VKS2 reversals. Finally, the model shows that coherent behaviors can also be naturally generated inside the many-mode dynamical chaotic model, which reproduces the complexity of fluid turbulence, as described by the shell technique.

Homotopic Arnold tongues deformation of the MHD α2‐dynamo

AbstractWe consider a mean–field α2–dynamo with helical turbulence parameter α(r)=α0+γΔα(r) and a boundary homotopy with parameter β∈[0,1] interpolating between Dirichlet (idealized, β=0) and Robin (physically realistic, β=1) boundary conditions. It is shown that the zones of oscillatory solutions at β=1 end up at the diabolical points for β=0 under the homotopic deformation. The underlying network of the diabolical points for β=0 substantially determines the choreography of eigenvalues and thus the character of the dynamo instability for β=1. Using perturbation theory we derive the first–order approximations to the resonance (Arnold's) tongues in the α0βγ‐space, which turn out to be cones in the vicinity of the diabolical points, selected by the Fourier coefficients of Δα(r). The space orientation of the 3D tongues is determined by the Krein signature of the modes involved in the diabolical crossings at the apexes of the cones. The Krein space induced geometry of the resonance zones explains the subtleties in finding α‐profiles leading to oscillatory dynamos, and it explicitly predicts the locations of the spectral exceptional points, which are important ingredients in the recent theories of polarity reversals of the geomagnetic field. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Asymptotic methods for spherically symmetric MHD α2‐dynamos

AbstractWe consider two models of spherically‐symmetric MHD α2–dynamos; one with idealized boundary conditions (BCs); and one with physically realistic BCs. As it has been shown in our previous work, the eigenvalues λ of a model with idealized BCs and constant α–profile α0 are linear functions of α0 and form a mesh in the (α0, λ)–plane. The nodes of the spectral mesh correspond to double‐degenerate eigenvalues of algebraic and geometric multiplicity 2 (diabolical points). It was found that perturbations of the constant α –profile lead to a resonant unfolding of the diabolical points with selection rules of the resonant unfolding defined by the Fourier coefficients of the perturbations. In the present contribution we present new exact results on the spectrum of the model with physically realistic BCs and constant α. For non‐degenerate (simple) eigenvalues perturbation gradients are found at any particular α0. We briefly discuss the spectral behavior of the α2–dynamo operator over a family of homotopic deformations of the BCs between idealized ones and physically realistic ones. Furthermore, we demonstrate that although the spectral singularities are lifted, a memory about their locations remains deeply imprinted in the homotopic family of spectral deformations due to a hidden underlying invariance. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

On Krein space related perturbation theory for MHD α2‐dynamos

AbstractThe spectrum of the spherically symmetric α2–dynamo is studied in the case of idealized boundary conditions. Starting from the exact analytical solutions of models with constant α –profiles a perturbation theory and a Galerkin technique are developed in a Krein‐space approach. With the help of these tools a very pronounced α –resonance pattern is found in the deformations of the spectral mesh as well as in the unfolding of the diabolical points located at the nodes of this mesh. Non‐oscillatory as well as oscillatory dynamo regimes are obtained. An estimation technique is developed for obtaining the critical α –profiles at which the eigenvalues enter the right spectral half‐plane with non‐vanishing imaginary components (at which overcritical oscillatory dynamo regimes form). (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

A multi-layered kinematic dynamo model: implications of a stratified upper layer in the Earth's core

It has been suggested that there exists a stably stratified electrically conducting layer at the top of the Earth's outer fluid core and that lateral temperature gradients in the lower mantle is capable of a driving thermal-wind-type flow near the core–mantle boundary. We investigate how such a flow in a stable layer could influence the geomagnetic field and the geodynamo using a very simple two-dimensional kinematic dynamo model in Cartesian geometry. The dynamo has four layers representing the inner core, convecting lower outer core, stable upper core, and insulating mantle. An α2 dynamo operates in the convecting outer core and a horizontal shear flow is imposed in the stable layer. Exact dynamo solutions are obtained for a range of parameters, including different conductivities for the stable layer and inner core. This allows us to connect our solutions with known, simpler solutions of a single-layer α2 dynamo, and thereby assess the effects of the extra layers. We confirm earlier results that a stable, static layer can enhance dynamo action. We find that shear flows produce dynamo wave solutions with a different spatial structure from the steady α2 dynamos solutions. The stable layer controls the behavior of the dynamo system through the interface conditions, providing a new means whereby lateral variations on the boundary can influence the geomagnetic field.