Nonhelical mean‐field dynamos in a sheared turbulence

Mechanisms of nonhelical large‐scale dynamos (shear‐current dynamo and effect of homogeneous kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large‐scale shear are discussed. We have found that the shearcurrent dynamo can act even in random flows with small Reynolds numbers. However, in this case mean‐field dynamo requires small magnetic Prandtl numbers (i.e., when Pm < Pm^cr < 1). The threshold in the magnetic Prandtl number, Pm^cr = 0.24, is determined using second order correlation approximation (or first‐order smoothing approximation) for a background random flow with a scale‐dependent viscous correlation time τ_c = (νk 2)^–1 (where ν is the kinematic viscosity of the fluid and k is the wave number). For turbulent flows with large Reynolds numbers shear‐current dynamo occurs for arbitrary magnetic Prandtl numbers. This dynamo effect represents a very generic mechanism for generating large‐scale magnetic fields in a broad class of astrophysical turbulent systems with large‐scale shear. On the other hand, mean‐field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not realistic for a broad class of astrophysical systems because it requires a very specific random forcing of kinetic helicity fluctuations that contains, e.g., low‐frequency oscillations. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Do mean‐field dynamos in nonrotating turbulent shear‐flows exist?

A plane‐shear flow in a fluid with forced turbulence is considered. If the fluid is electrically‐conducting then a mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic tensor. It is checked whether in this case self‐excitation of a large‐scale magnetic field is possible (so‐called W̄ × J̄ ‐dynamo) and the answer is NO. The calculations reveal the cross‐stream components of the EMF perpendicular to the mean current having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic Prandtl number of about unity have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather strong. If, on the other hand, the turbulence field is strati.ed in the vertical direction then a dynamo‐active α ‐effect is produced. The critical magnetic Reynolds number for such a self‐excitation in a simple shear flow is slightly above 10 like for the other – but much more complicated – flow patterns used in existing dynamo experiments with liquid sodium or gallium. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dimensionless measures of turbulent magnetohydrodynamic dissipation rates

The magnetic Reynolds number, R _M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R _M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number,   R _M,e  , as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range   R _M,e∼ O (1) ≪ R _M  . We also distinguish   R _M,e  from     where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That   R _M,e  and     approach unity even if   R _M≫ 1  highlights that, just as in hydrodynamic turbulence, energy dissipation of large-scale structures in turbulent flows via a cascade can be much faster than the dissipation of large-scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.     

Suppression of the large‐scale Lorentz force by turbulence

The components of the total stress tensor (Reynolds stress plus Maxwell stress) are computed within the quasilinear approximation for a driven turbulence influenced by a large‐scale magnetic background field. The conducting fluid has an arbitrary magnetic Prandtl number and the turbulence without the background field is assumed as homogeneous and isotropic with a free Strouhal number St. The total large‐scale magnetic tension is always reduced by the turbulence with the possibility of a ‘catastrophic quenching’ for large magnetic Reynolds number Rm so that even its sign is reversed. The total magnetic pressure is enhanced by turbulence in the high‐conductivity limit but it is reduced in the low‐conductivity limit. Also in this case the sign of the total pressure may reverse but only for special turbulences with sufficiently large St > 1. The turbulence‐induced terms of the stress tensor are suppressed by strong magnetic fields. For the tension term this quenching grows with the square of the Hartmann number of the magnetic field. For microscopic (i.e. small) diffusivity values the magnetic tension term becomes thus highly quenched even for field amplitudes much smaller than their equipartition value. In the opposite case of large‐eddy simulations the magnetic quenching is only mild but then also the turbulence‐induced Maxwell tensor components for weak fields remain rather small (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the interaction between differential rotation and magnetic fields in the Sun

We have performed 3-D numerical simulations of compressible convection under the influence of rotation and magnetic fields in spherical shells. They aim at understanding the subtle coupling between convection, rotation and magnetic fields in the solar convection zone. We show that as the magnetic Reynolds number is increased in the simulations, the magnetic energy saturates via nonlinear dynamo action, to a value smaller but comparable to the kinetic energy contained in the shell, leading to increasingly strong Maxwell stresses that tend to weaken the differential rotation driven by the convection. These simulations also indicate that the mean toroidal and poloidal magnetic fields are small compared to their fluctuating counterparts, most of the magnetic energy being contained in the non-axisymmetric fields. The intermittent nature of the magnetic fields generated by such a turbulent convective dynamo confirms that in the Sun the large-scale ordered dynamo responsible for the 22-year cycle of activity can hardly be located in the solar convective envelope.     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the role of magnetic fields in star formation

Magnetic fields are observed in star forming regions. However simulations of the late stages of star formation that do not include magnetic fields provide a good fit to the properties of young stars including the initial mass function (IMF) and the multiplicity. We argue here that the simulations that do include magnetic fields are unable to capture the correct physics, in particular the high value of the magnetic Prandtl number, and the low value of the magnetic diffusivity. The artificially high (numerical and uncontrolled) magnetic diffusivity leads to a large magnetic flux pervading the star forming region. We argue further that in reality the dynamics of high magnetic Prandtl number turbulence may lead to local regions of magnetic energy dissipation through reconnection, meaning that the regions of molecular clouds which are forming stars might be essentially free of magnetic fields. Thus the simulations that ignore magnetic fields on the scales on which the properties of stellar masses, stellar multiplicities and planet-forming discs are determined, may be closer to reality than those which include magnetic fields, but can only do so in an unrealistic parameter regime.     

Strong mean field dynamos require supercritical helicity fluxes

Several one and two dimensional mean field models are analyzed where the effects of current helicity fluxes and boundaries are included within the framework of the dynamical quenching model. In contrast to the case with periodic boundary conditions, the final saturation energy of the mean field decreases inversely proportional to the magnetic Reynolds number. If a nondimensional scaling factor in the current helicity flux exceeds a certain critical value, the dynamo can operate even without kinetic helicity, i.e. it is based only on shear and current helicity fluxes, as first suggested by Vishniac & Cho (2001, ApJ 550, 752). Only above this threshold is the current helicity flux also able to alleviate catastrophic quenching. The fact that certain turbulence simulations have now shown apparently non‐resistively limited mean field saturation amplitudes may be suggestive of the current helicity flux having exceeded this critical value. Even below this critical value the field still reaches appreciable strength at the end of the kinematic phase, which is in qualitative agreement with dynamos in periodic domains. However, for large magnetic Reynolds numbers the field undergoes subsequent variations on a resistive time scale when, for long periods, the field can be extremely weak. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Axisymmetry vs. nonaxisymmetry of a Taylor‐Couette flow with azimuthal magnetic fields

The instability of a supercritical Taylor‐Couette flow of a conducting fluid with resting outer cylinder under the influence of a uniform axial electric current is investigated for magnetic Prandtl number Pm = 1. In the linear theory the critical Reynolds number for axisymmetric perturbations is not influenced by the current‐induced axisymmetric magnetic field but all axisymmetric magnetic perturbations decay. The nonaxisymmetric perturbations with m = 1 are excited even without rotation for large enough Hartmann numbers (“Tayler instability”). For slow rotation their growth rates scale with the Alfvén frequency of the magnetic field but for fast rotation they scale with the rotation rate of the inner cylinder. In the nonlinear regime the ratio of the energy of the magnetic m = 1 modes and the toroidal background field is very low for the non‐rotating Tayler instability but it strongly grows if differential rotation is present. For super‐Alfv´enic rotation the energies in the m = 1 modes of flow and field do not depend on the molecular viscosity, they are almost in equipartition and contain only 1.5 % of the centrifugal energy of the inner cylinder. The geometry of the excited magnetic field pattern is strictly nonaxisymmetric for slow rotation but it is of the mixed‐mode type for fast rotation – contrary to the situation which has been observed at the surface of Ap stars. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear simulations of magnetic Taylor‐Couette flow with currentfree helical magnetic fields

Themagnetorotational instability (MRI) in cylindrical Taylor‐Couette flow with external helical magnetic field is simulated for infinite and finite aspect ratios. We solve the MHD equations in their small Prandtl number limit and confirm with timedependent nonlinear simulations that the additional toroidal component of the magnetic field reduces the critical Reynolds number from O (10⁶) (axial field only) to O (10³) for liquid metals with their small magnetic Prandtl number. Computing the saturated state we obtain velocity amplitudes which help designing proper experimental setups. Experiments with liquid gallium require axial field ∼50 Gauss and axial current ∼4 kA for the toroidal field. It is sufficient that the vertical velocity u_z of the flow can be measured with a precision of 0.1 mm/s.We also show that the endplates enclosing the cylinders do not destroy the traveling wave instability which can be observed as presented in earlier studies. For TC containers without and with endplates the angular momentum transport of the MRI instability is shown as to be outwards. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stability of the solar tachocline with magnetic fields

The stability of magnetic fields in the solar tachocline is investigated. We present stability limits for higher azimuthal wave numbers and results on the dependence of the stability on the location of toroidal magnetic fields in latitude. While the dependence of the wave number with the largest growth rate on the magnetic field strength and the magnetic Prandtl number is small, the dependence on the magnetic Reynolds number Rm indicates that lowest azimuthal modes are excited for very high Rm. Upon varying the latitudinal position of the magnetic field belts, we find slightly lower stability limits for high latitudes, and very large stability limits at latitudes below 10°, with little dependence on latitude in between. An increase of the maximum possible field was achieved by adding a poloidal field. The upper limit for the toroidal field which can be stored in the radiative tachocline is then 1000 G, compared to about 100 G for a purely toroidal field as was found in an earlier work. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar Magnetic Carpet II: Coronal Interactions of Small-Scale Magnetic Fields

This paper is the second in a series of studies working towards constructing a realistic, evolving, non-potential coronal model for the solar magnetic carpet. In the present study, the interaction of two magnetic elements is considered. Our objectives are to study magnetic energy build-up, storage and dissipation as a result of emergence, cancellation, and flyby of these magnetic elements. In the future these interactions will be the basic building blocks of more complicated simulations involving hundreds of elements. Each interaction is simulated in the presence of an overlying uniform magnetic field, which lies at various orientations with respect to the evolving magnetic elements. For these three small-scale interactions, the free energy stored in the field at the end of the simulation ranges from 0.2 – 2.1×10²⁶ ergs, whilst the total energy dissipated ranges from 1.3 – 6.3×10²⁶ ergs. For all cases, a stronger overlying field results in higher energy storage and dissipation. For the cancellation and emergence simulations, motion perpendicular to the overlying field results in the highest values. For the flyby simulations, motion parallel to the overlying field gives the highest values. In all cases, the free energy built up is sufficient to explain small-scale phenomena such as X-ray bright points or nanoflares. In addition, if scaled for the correct number of magnetic elements for the volume considered, the energy continually dissipated provides a significant fraction of the quiet Sun coronal heating budget.     

The solar wind-magnetosphere dynamo and the magnetospheric substorm

In a quiet condition, the solar wind kinetic energy is converted into electrical energy. A small part of this energy is dissipated as heat energy in the polar ionosphere. We identify at least three types of magnetospheric disturbances which are not associated with an increase of the heat production and call them reversible disturbances, while the magnetospheric substorm is an irreversible disturbance which is associated with a large increase of the heat production.The magnetosphere appears to have an inherent internal instability by which a large amount of heat energy is sporadically produced in the polar upper atmosphere at the expense of the magnetic energy in the magnetotail. A positive feed-back process may be responsible for the growth of the instability and for the expansive phase, while the recovery phase sets in when some process begins to suppress the positive feed-back process.     

The dual role of shear in large‐scale dynamos

The role of shear in alleviating catastrophic quenching by shedding small‐scale magnetic helicity through fluxes along contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force‐free degeneracy of Beltrami fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition field strength those estimates result in a 5‐fold decrease of the magnetic diffusivity as the magnetic Reynolds number based on the wavenumber of the energy‐carrying eddies is increased from 2 to 600. Finally, the role of shear in driving turbulence and large‐scale fields by the magneto‐rotational instability is emphasized. New simulations are presented and the 3π /4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a useful diagnostic tool in identifying mean‐field dynamo action in simulations and to distinguish this from other scenarios invoking magnetic buoyancy as a means to explain migration away from the midplane. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

耀斑大尺度电流片的湍流耗散和能谱分析

磁重联在宇宙的许多动力学现象中都是非常核心的过程.磁流体动力学(MHD)数值模拟是研究磁重联过程以及相应物理图像的一种很有效的手段.通过不同的参数组合,来研究MHD数值模拟中磁雷诺数和空间分辨率对磁重联率、数值耗散和能谱分布的影响.对得到的数据进行分析后,发现磁雷诺数对磁重联率和能谱分布有一定的影响.磁雷诺数越大,磁重联过程进入非线性阶段所需的特征时间越短,磁重联率就越早发生跃升.磁雷诺数R_m对耗散开始发挥作用的Kolmogorov微观尺度l_ko有明显影响:R_m越大,l_ko就越小.研究了磁重联过程中包括数值耗散在内的额外耗散对重联过程的影响.结果表明,撕裂模不稳定性开始之前的额外耗散以纯数值耗散为主,撕裂模不稳定性出现之后,额外耗散出现同步跃升,说明不稳定性导致的湍流明显增强了耗散的效果,相当于在局部湍流区引入了超电阻.能谱分析进一步表明,大尺度电流片的l_ko完全可能出现在宏观的MHD尺度上.     

The Tayler instability of toroidal magnetic fields in a columnar gallium experiment

The nonaxisymmetric Tayler instability of toroidal magnetic fields due to axial electric currents is studied for conducting incompressible fluids between two coaxial cylinders without endplates. The inner cylinder is considered as so thin that the limit of R_in → 0 can be computed. The magnetic Prandtl number is varied over many orders of magnitudes but the azimuthal mode number of the perturbations is fixed to m = 1. In the linear approximation the critical magnetic field amplitudes and the growth rates of the instability are determined for both resting and rotating cylinders. Without rotation the critical Hartmann numbers do not depend on the magnetic Prandtl number but this is not true for the corresponding growth rates. For given product of viscosity and magnetic diffusivity the growth rates for small and large magnetic Prandtl number are much smaller than those for Pm = 1. For gallium under the influence of a magnetic field at the outer cylinder of 1 kG the resulting growth time is 5 s. The minimum electric current through a container of 10 cm diameter to excite the instability is 3.20 kA. For a rotating container both the critical magnetic field and the related growth times are larger than for the resting column (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar Magnetic Carpet III: Coronal Modelling of Synthetic Magnetograms

This article is the third in a series working towards the construction of a realistic, evolving, non-linear force-free coronal-field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48-hour time series of synthetic magnetograms produced from the model of Meyer et al. (Solar Phys. 272, 29, 2011). Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5?–?0.8 Mm. The energy dissipated is currently found to be too small to account for the heating of the entire quiet-Sun corona. However, the form and location of energy-dissipation regions qualitatively agree with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release.     

Coronal activity from dynamos in astrophysical rotators

We show that a steady mean-field dynamo in astrophysical rotators leads to an outflow of relative magnetic helicity and thus magnetic energy available for particle and wind acceleration in a corona. The connection between energy and magnetic helicity arises because mean-field generation is linked to an inverse cascade of magnetic helicity. To maintain a steady state in large magnetic Reynolds number rotators, there must then be an escape of relative magnetic helicity associated with the mean field, accompanied by an equal and opposite contribution from the fluctuating field. From the helicity flow, a lower limit on the magnetic energy deposited in the corona can be estimated. Steady coronal activity including the dissipation of magnetic energy, and formation of multi-scale helical structures therefore necessarily accompanies an internal dynamo. This highlights the importance of boundary conditions which allow this to occur for non-linear astrophysical dynamo simulations. Our theoretical estimate of the power delivered by a mean-field dynamo is consistent with that inferred from observations to be delivered to the solar corona, the Galactic corona, and Seyfert 1 AGN coronae.     

An accretion disc model with a magnetic wind and turbulent viscosity

A model is presented for an accretion disc with turbulent viscosity and a magnetically influenced wind. The magnetic field is generated by a dynamo in the disc, involving the turbulence and radial shear. Disc-wind solutions are found for which the wind mass flux is sufficient to play a major part in driving an imposed steady inflow, but small enough for most material to be accreted on to the central object. Constraints arise for the magnetic Reynolds and Prandtl numbers in terms of the turbulent Mach number and vertical length-scale of the disc's horizontal magnetic field. It is shown that the imposition of a stellar boundary condition enhances the wind mass flux in the very inner region of the disc and may result in jet formation.     

Simulations of the anisotropic kinetic and magnetic alpha effects

Using simulations of isotropically forced helical turbulence the contributions to kinetic and magnetic alpha effects are computed. It is shown that for the parameter regimes considered in an earlier publication (Brandenburg & Subramanian 2005), the expressions for isotropic and anisotropic alpha effects give quantitatively similar results. Both kinetic and magnetic alpha effects are proportional to a relaxation time whose value, in units of the turnover time, is shown to be approximately unity and independent of the magnetic Reynolds number. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)