Strongly stratified hydromagnetic turbulence has previously been found to produce magnetic flux concentrations if the domain is large enough compared with the size of turbulent eddies. Mean-field simulations (MFS) using parameterizations of the Reyno
lds and Maxwell stresses show a negative effective magnetic pressure instability and have been able to reproduce many aspects of direct numerical simulations (DNS) regarding the growth rate of this large-scale instability, shape of the resulting magnetic structures, and their height as a function of magnetic field strength. Unlike the case of an imposed horizontal field, for a vertical one, magnetic flux concentrations of equipartition strength with the turbulence can be reached. This results in magnetic spots that are reminiscent of sunspots. Here we want to find out under what conditions magnetic flux concentrations with vertical field occur and what their internal structure is. We use a combination of MFS, DNS, and implicit large-eddy simulations to characterize the resulting magnetic flux concentrations in forced isothermal turbulence with an imposed vertical magnetic field. We confirm earlier results that in the kinematic stage of the large-scale instability the horizontal wavelength of structures is about 10 times the density scale height. At later times, even larger structures are being produced in a fashion similar to inverse spectral transfer in helically driven turbulence. Using turbulence simulations, we find that magnetic flux concentrations occur for different values of the Mach number between 0.1 and 0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are being confined by downflows along the tubes and corresponding inflow from the sides, which keep the field concentrated.
Energetic nonthermal particles (cosmic rays, CRs) are accelerated in supernova remnants, relativistic jets and other astrophysical objects. The CR energy density is typically comparable with that of the thermal components and magnetic fields. In this
review we discuss mechanisms of magnetic field amplification due to instabilities induced by CRs. We derive CR kinetic and magnetohydrodynamic equations that govern cosmic plasma systems comprising the thermal background plasma, comic rays and fluctuating magnetic fields to study CR-driven instabilities. Both resonant and non-resonant instabilities are reviewed, including the Bell short-wavelength instability, and the firehose instability. Special attention is paid to the longwavelength instabilities driven by the CR current and pressure gradient. The helicity production by the CR current-driven instabilities is discussed in connection with the dynamo mechanisms of cosmic magnetic field amplification.
The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbu
lent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.
When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spat
ial and temporal scale separations are poor. The approximate form of the kernel 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 alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.
We determine the timescales associated with turbulent diffusion and isotropization in closure models using anisotropically forced and freely decaying turbulence simulations and to study the applicability of these models. We compare the results from a
nisotropically forced three-dimensional numerical simulations with the predictions of the closure models and obtain the turbulent timescales mentioned above as functions of the Reynolds number. In a second set of simulations, turning the forcing off enables us to study the validity of the closures in freely decaying turbulence. Both types of experiments suggest that the timescale of turbulent diffusion converges to a constant value at higher Reynolds numbers. Furthermore, the relative importance of isotropization is found to be about 2.5 times larger at higher Reynolds numbers than in the more viscous regime.
We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Suns surface and following their evolution for a pre
scribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.
Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using bo
th the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both forces are chosen to have Roberts flow-like geometry. Results: Examples with an MHD background are studied where the previously used quasi-kinematic test-field method breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the same results as the imposed-field method, where the field-aligned component of the actual electromotive force from the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, $alpha$ is found to be quenched proportional to the fourth power of the field strength, regardless of the type of background studied.
(abridged) Aims: To study turbulent transport coefficients that describe the evolution of large-scale magnetic fields in turbulent convection. Methods: We use the test field method together with 3D numerical simulations of turbulent convection with s
hear and rotation to compute turbulent transport coefficients describing the evolution of large-scale magnetic fields in mean-field theory in the kinematic regime. 1D mean-field models are used with the derived turbulent transport coefficients to compare with direct simulations. Results: The alpha-effect increases monotonically as rotation increases. Turbulent diffusivity, eta_t, is proportional to the square of the turbulent vertical velocity. Whereas eta_t decreases approximately inversely proportional to the wavenumber of the field, the alpha-effect and turbulent pumping show a more complex behaviour. In the presence of shear and no rotation a small alpha-effect is induced which does not seem to show any consistent trend as a function of shear. If the shear is large enough, this small alpha is able to excite a dynamo in the mean-field model. The coefficient responsible for driving the shear-current effect shows several sign changes as a function of depth but is also able to contribute to dynamo action in the mean-field model. The growth rates in these cases are well below those in direct simulations suggesting that an incoherent alpha-shear dynamo may also act in them. If both rotation and shear are present, the alpha-effect is more pronounced. The combination of the shear-current and Omega x J-effects is also stronger than in the case of shear only, but subdominant to the alpha-shear dynamo. The results of direct simulations are consistent with mean-field models where all of these effects are taken into account without the need to invoke incoherent effects.
Dynamo action owing to helically forced turbulence and large-scale shear is studied using direct numerical simulations. The resulting magnetic field displays propagating wave-like behavior. This behavior can be modelled in terms of an alphaOmega dyna
mo. In most cases super-equipartition fields are generated. By varying the fraction of helicity of the turbulence the regeneration of poloidal fields via the helicity effect (corresponding to the alpha-effect) is regulated. The saturation level of the magnetic field in the numerical models is consistent with a linear dependence on the ratio of the fractional helicities of the small and large-scale fields, as predicted by a simple nonlinear mean-field model. As the magnetic Reynolds number (Rm) based on the wavenumber of the energy-carrying eddies is increased from 1 to 180, the cycle frequency of the large-scale field is found to decrease by a factor of about 6 in cases where the turbulence is fully helical. This is interpreted in terms of the turbulent magnetic diffusivity, which is found to be only weakly dependent on Rm.
The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesi
an domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.
