Buoyancy-driven perturbations in a rapidly rotating,electrically conducting fluid: part II – dynamo action

The helicity, electromotive force and α-effect produced in a homogeneous, rapidly rotating, electrically conducting fluid by an isolated source of buoyancy at small Elsasser number are calculated, visualized and analyzed. Due to physical symmetries of the system, the integrals of helicity and electromotive force over all space are zero. However, each has a significant non-zero value when integrated over the cross section of the Taylor column. The local α-effect is found to be significantly anisotropic; it is strongest when the applied magnetic field is toroidal and the resulting EMF is parallel to the applied field.     

Small-scale magnetic buoyancy and magnetic pumping effects in a turbulent convection

We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused by three kinds of the inhomogeneities; i.e., inhomogeneous turbulence, the nonuniform fluid density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy and magnetic pumping effects which are due to the effect of the mean magnetic field on the large-scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic field, while they do not affect the turbulent magnetic diffusion of the poloidal field.     

Mean-field concept and direct numerical simulations of rotating magnetoconvection and the geodynamo

Mean-field theory describes magnetohydrodynamic processes leading to large-scale magnetic fields in various cosmic objects. In this study magnetoconvection and dynamo processes in a rotating spherical shell are considered. Mean fields are defined by azimuthal averaging. In the framework of mean-field theory, the coefficients which determine the traditional representation of the mean electromotive force, including derivatives of the mean magnetic field up to the first order, are crucial for analyzing and simulating dynamo action. Two methods are developed to extract mean-field coefficients from direct numerical simulations of the mentioned processes. While the first method does not use intrinsic approximations, the second one is based on the second-order correlation approximation. There is satisfying agreement of the results of both methods for sufficiently slow fluid motions. Both methods are applied to simulations of rotating magnetoconvection and a quasi-stationary geodynamo. The mean-field induction effects described by these coefficients, e.g., the α-effect, are highly anisotropic in both examples. An α²-mechanism is suggested along with a strong γ-effect operating outside the inner core tangent cylinder. The turbulent diffusivity exceeds the molecular one by at least one order of magnitude in the geodynamo example. With the aim to compare mean-field simulations with corresponding direct numerical simulations, a two-dimensional mean-field model involving all previously determined mean-field coefficients was constructed. Various tests with different sets of mean-field coefficients reveal their action and significance. In the magnetoconvection and geodynamo examples considered here, the match between direct numerical simulations and mean-field simulations is only satisfying if a large number of mean-field coefficients are involved. In the magnetoconvection example, the azimuthally averaged magnetic field resulting from the numerical simulation is in good agreement with its counterpart in the mean-field model. However, this match is not completely satisfactory in the geodynamo case anymore. Here the traditional representation of the mean electromotive force ignoring higher than first-order spatial derivatives of the mean magnetic field is no longer a good approximation.     

On the role of magnetostrophic waves in geodynamo

It is shown that magnetostrophic waves which are generated in the equatorial plane of the Earth’s core due to the instability of the equatorial jet and which propagate almost transversely to the rotational axis off the tangent cylinder, have a negative helicity in the northern hemisphere and positive helicity in the southern hemisphere. When the wave trains propagate through the regions with a constant azimuthal magnetic field caused by the Ω-effect, this helicity distribution induces an electromotive force (emf) (due to the α-effect), which may lead to the maintenance of the initial dipole field by the scenario of the α-Ω dynamo.     

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.     

On hydromagnetic instabilities and the mean electromotive force in a non-uniformly stratified Earth’s core affected by viscosity

Linear magnetoconvection in a model of a non-uniformly stratified horizontal rotating fluid layer with a toroidal magnetic field is investigated for no-slip and finitely electrically conductive boundaries and with very thin stably stratified upper sublayer. The basic parabolic temperature profile is determined by the temperature difference between the boundaries and by the homogeneous heat source distribution in the layer. This results in a density pattern, in which a stably stratified upper sublayer is present. The developed diffusive perturbations (modes) are strongly affected by the complicated coupling of viscous, thermal and magnetic diffusive processes. The calculations were performed for various values of Roberts number (q ≪ 1 and q = O(1)). The mean electromotive force produced by the developed hydromagnetic instabilities is investigated to find the modes, which can be appropriate for creating the α-effect. It was found that the azimuthal part of the EMF is dominant for westward modes when the Elsasser number Λ ≲ O(1).     

Book Review of Accretion Power in Astrophysics (3rd Edition) by Juhan Frank,Andrew King,Derek Raine Cambridge University Press, 2002, XIV+384 pp., £27.95, hardback (ISBN 0 521 62053 8), paperback (ISBN 0 521 62957 8).

The velocity, pressure, perturbation magnetic field, helicity and electromotive force driven by an isolated buoyant parcel in an unbounded, rapidly rotating, electrically conducting fluid in the limit of small Elsasser number and very small Ekman number are calculated, visualized and analyzed. On the scale of the parcel, the solution is identical to that obtained in the limit of small Ekman number and zero Elsasser number. On the scale of the Taylor-column, it is elongated in the direction of the applied magnetic field and compressed in the direction perpendicular to it. The α-effect calculated by averaging the electromotive force on planes normal to rotation is strongly anisotropic: near the parcel and in the inner part of the Taylor-column it is strongest when the applied magnetic field is perpendicular to rotation and gravity; in the outer part of the Taylor-column it is strongest when the applied magnetic field is in the same plane as rotation and gravity.     

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.     

Tuning of the mean-field geodynamo model

Parker’s two-dimensional (2D) dynamo model with an algebraic form of nonlinearity for the α-effect is considered. The model uses geostrophic distributions for the α-effect and differential rotation, which are derived from the three-dimensional (3D) convection models. The resulting configurations of the magnetic field in the liquid core are close to the solutions in Braginsky’s Z-model. The implications of the degree of geostrophy observed in the 3D dynamo models for the behavior of the mean magnetic field are explored. It is shown that the reduction in geostrophy leads to magnetic field reversals accompanied by the relative growth of the nondipole component of the field on the surface of the liquid core. The simulations with a random α-effect which causes turbulent pulsations are carried out. The approach is capable of producing realistic sequences of magnetic reversals.     

Geomagnetic polarity reversals in a turbulent core

The behavior of the main magnetic field components during a polarity transition is investigated using the α²-dynamo model for magnetic field generation in a turbulent core. It is shown that rapid reversals of the dipole field occur when the helicity, a measure of correlation between turbulent velocity and vorticity, changes sign. Two classes of polarity transitions are possible. Within the first class, termed component reversals, the dipole field reverses but the toroidal field does not. Within the second class, termed full reversals, both dipole and toroidal fields reverse. Component reversals result from long term fluctuations in core helicity; full reversals result from short term fluctuations. A set of time-evolution equations are derived which govern the dipole field behavior during an idealized transition. Solutions to these equations exhibit transitions in which the dipole remains axial while its intensity decays rapidly toward zero, and is regenerated with reversed polarity. Assuming an electrical conductivity of 3 × 10⁵ mho m^?1 for the fluid core, the time interval required to complete the reversal process can be as short as 7500 years. This time scale is consistent with paleomagnetic observations of the duration of reversals. A possible explanation of the cause of reversals is proposed, in which the core's net helicity fluctuates in response to fluctuations in the level of turbulence produced by two competing energy sources—thermal convection and segregation of the inner core. Symmetry considerations indicate that, in each hemisphere, helicity generated by heat loss at the core-mantle boundary may have the opposite sign of helicity generated by energy release at the inner core boundary. Random variations in rates of energy release can cause the net helicity and the α-effect to change sign occasionally, provoking a field reversal. In this model, energy release by inner core formation tends to destabilize stationary dynamo action, causing polarity reversals.     

湍流输送的非线性热力学性质

湍流输送是一种热力学不可逆过程,本文利用非线性热力学研究了湍流输送的特征. 将热力学流对热力学力以平衡态作为参考态进行Taylor展开,可以得到湍流输送系数是系统宏观参量梯度的Taylor级数. 线性湍流输送系数是Reynolds湍流闭合方案的K闭合湍流输送系数;而湍流输送系数非线性项则是系统偏离热力学平衡态所造成的热力学非线性效应. 湍流输送系数这一热力学性质提供了一种热力学湍流闭合方案. 线性湍流输送系数是正定的,湍流输送只能使系统宏观参量均匀化;而在远离平衡态的热力学非线性区,可能导致湍流输送系数负黏性现象. 在最小熵产生态的条件下,热力学流对热力学力Taylor展开的各级系数间存在一种递推关系. 利用这种递推关系大大减少了由实验确定的Taylor级数的系数个数.     

Buoyancy-driven perturbations in a rapidly rotating,electrically conducting fluid: part I – flow and magnetic field

The steady velocity, perturbation pressure and perturbation magnetic field, driven by an isolated buoyant parcel of Gaussian shape in a rapidly rotating, unconfined, incompressible electrically conducting fluid in the presence of an imposed uniform magnetic field, are obtained by means of the Fourier transform in the limit of small Ekman number. Lorentz and inertial forces are neglected. The solution requires at most evaluation of a single integral and is found in closed form in some spatial regions. The solution has structure on two disparate scales: on the scale of the buoyant parcel and on the scale of the Taylor column, which is elongated in the direction of the rotation axis. The detailed structures of the flow and pressure depend linearly on the relative orientation of gravity and rotation, with the solution for arbitrary orientation being a linear combination of two limiting cases in which these vectors are colinear (polar case) and perpendicular (equatorial case). The perturbation magnetic field depends additionally on the relative orientation of the imposed magnetic field, and three limiting cases of interest are presented in which gravity and rotation are colinear (polar–toroidal case), gravity and imposed field are colinear (equatorial–radial case) and all three are mutually perpendicular (equatorial–toroidal case). Visualization and analysis of the velocity and perturbation magnetic field vectors are facilitated by dividing these vector fields into geostrophic and ageostrophic protions. In all cases, the geostrophic and ageostrophic portions have different structure on the Taylor-column scale. The buoyancy force is balanced by a pressure force in the polar case and by a flux of momentum in the equatorial case. The pressure force and momentum flux do not decay in strength with increasing axial distance. Far from the parcel, the axial mass flux varies as the inverse one-third power of distance from the parcel. The velocity has a single geostrophic vortex in the polar case and two vortices in the equatorial case. The perturbation magnetic field has two, four and one geostrophic vortices in the polar–toroidal, equatorial–radial and equatorial–toroidal cases, respectively. To facilitate comparison of the present results with numerical simulations carried out in a finite domain, a set of boundary conditions are developed, with may be applied at a finite distance from the parcel.     

Problem of the Rotating Magnetoconvection in Variously Stratified Fluid Layer Revisited

The linear magnetoconvection in the rotating uniformly as well as non-uniformly stratified horizontal layer with azimuthal magnetic field is investigated for the various mechanical and electrical boundary conditions and especially, for various values of Roberts number. The developed diffusive perturbations (modes) are strongly influenced not only by the mentioned properties of boundaries but also by complicated coupling of viscous, thermal and magnetic diffusive processes. The mean electromotive force produced by developed hydromagnetic instabilities is also investigated to determine the hydromagnetic processes which are appropriate for -effect. The presented paper is an unification of hitherto published results of the authors and gives a short survey of many developments of corresponding model by Soward (1979).     

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.     

A comparison of numerical schemes to solve the magnetic induction eigenvalue problem in a spherical geometry

We compare various methods of solving the magnetic induction eigenvalue problem in a sphere, each using toroidal–poloidal decomposition and spherical harmonics, but with a different radial discretisation. In the case of quiescent flow where only diffusion acts upon the magnetic field, we benchmark numerical convergence against the analytic decay rates, and find that a Galerkin scheme based on Chebyshev polynomials with an associated projection chosen such that the diffusion operator is self-adjoint, exhibits the fastest convergence of the schemes described. The importance of the speed of convergence becomes heightened with the introduction of a non-quiescent flow because of the reduction in the magnetic field length scales. We find that sufficiently converged solutions are generally difficult to locate unless we use the optimal Galerkin scheme.     

一种磁张量探测系统载体的磁张量补偿方法

针对磁张量系统载体产生的磁张量值对系统测量精度产生很大影响的问题,以及现有磁补偿模型存在非线性、分体式和参数多的问题,提出一种磁张量系统载体的一体化线性磁张量补偿方法.分析了载体硬磁材料产生固有磁张量值和软磁材料产生感应磁张量值的微观机理,并推导了相应的数学表达式,结合固有磁场影响和感应磁场影响建立了载体磁张量补偿模型.模型中含有20个载体磁张量补偿系数,对模型求解得到补偿系数,结合三分量磁场测量值即可达到对载体磁张量的补偿.实测实验表明,磁张量补偿方法计算得到的载体磁张量值与载体实际产生的磁张量值仅差32nT/m,可以有效完成对磁张量系统的载体磁张量补偿.     

Non-local effects in the mean-field disc dynamo I. An asymptotic expansion

Abstract

We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion.     

Processing of magnetotelluric data

The processing of magnetotelluric data involves concepts from electromagnetic theory, time series analysis and linear systems theory for reducing natural electric and magnetic field variations recorded at the earth's surface to forms suitable for studying the electrical properties of the earth's interior.The electromagnetic field relations lead to either a scalar transfer impedance which couples an electric component to an orthogonal magnetic component at the surface of a plane-layered earth, or a tensor transfer impedance which couples each electric component to both magnetic components in the vicinity of a lateral inhomogeneity.A number of time series spectral analysis methods can be used for estimating the complex spectral coefficients of the various field quantities. These in turn are used for estimating the nature of the transfer function or tensor impedance. For two dimensional situations, the tensor impedance can be rotated to determine the principal directions of the electrical structure.In general for real data, estimates of the apparent resistivity are more stable when calculated from the tensor elements rather than from simple orthogonal field ratios (Cagniard estimates), even when the fields are measured in the principal coordinates.     

Solar polar magnetic field

The solar polar magnetic field has attracted the attention of researchers since the polar magnetic field reversal was revealed in the middle of the last century (Babcock and Livingston, 1958). The polar magnetic field has regularly reversed because the magnetic flux is transported from the sunspot formation zone owing to differential rotation, meridional circulation, and turbulent diffusion. However, modeling of these processes leads to ambiguous conclusions, as a result of which it is sometimes unclear whether a transport model is actual. Thus, according to the last Hinode data, the problem of a standard transport model (Shiota et al., 2012) consists in that a decrease in the polar magnetic flux in the Southern Hemisphere lags behind such a decrease in the flux in the Northern Hemisphere (from 2008 to June 2012). On the other hand, Svalgaard and Kamide (2012) consider that the asymmetry in the sign reversal simply results from the asymmetry in the emerging flux in the sunspot formation region. A detailed study of the polar magnetic flux evolution according to the Solar Dynamics Observatory (SDO) data for May 2010–December 2012 is illustrated in the present work. Helioseismic & Magnetic Imager (HMI) magnetic data in the form of a magnetic field component along the line of sight (the time resolution is 720 s) are used here. The magnetic fluxes in sunspot formation regions and at high latitudes have been compared.     

Effects of spatial diffusion in the inner plasma sheet and the structure of the diffusion tensor

A standard pair of equations is used to describe the behaviour of a single monoenergetic particle (proton or electron) population on a geomagnetic flux tube drifting in the magnetosphere. When particle losses from the drifting flux tube into the ionosphere are neglected, this behaviour is adiabatic in a thermodynamic sense. For a population of particles with an isotropic pitch-angle distribution, the generalization of that system of equations is obtained by adding radial and azimuthal spatial diffusion terms. The magnetic field is taken to be dipolar in the inner magnetosphere. The potential electric field is assumed to consist of magnetospheric convection and corotation components. Experimental data are used to estimate the radial equatorial profiles of the plasma sheet pressure. Assuming that the local time and L-shell variations are separable and supposing steady-state conditions, the expressions for the diffusion tensor components are evaluated. The influence of spatial diffusion on the radial and azimuthal profiles of the plasma pressure in the inner plasma sheet is also discussed.