Koen Kemel

DETECTION OF NEGATIVE EFFECTIVE MAGNETIC PRESSURE INSTABILITY IN TURBULENCE SIMULATIONS

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

THE NEGATIVE EFFECTIVE MAGNETIC PRESSURE IN STRATIFIED FORCED TURBULENCE

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to converge for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models, the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough.

Active Region Formation through the Negative Effective Magnetic Pressure Instability

The negative effective magnetic-pressure instability operates on scales encompassing many turbulent eddies, which correspond to convection cells in the Sun. This instability is discussed here in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details, for example the onset of the instability occurs at the same depth. This depth increases with increasing field strength, such that the growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.

Spontaneous Formation of Magnetic Flux Concentrations in Stratified Turbulence

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNSs) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in the formation of large-scale inhomogeneous magnetic structures from an initially uniform magnetic field. Earlier investigations of this instability in DNSs of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma β are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L2/ηt, where ηt is the turbulent diffusivity and L is the scale of the domain. DNSs are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.

Properties of the negative effective magnetic pressure instability

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large-scale instability. In this study, ...

Mean-field transport in stratified and/or rotating turbulence

Context. The large-scale magnetic fields of stars and galaxies are often described in the framework of mean-field dynamo theory. At moderate magnetic Reynolds numbers, the transport coefficients de ...

Model of driven and decaying magnetic turbulence in a cylinder

Using mean-field theory, we compute the evolution of the magnetic field in a cylinder with outer perfectly conducting boundaries and imposed axial magnetic and electric fields. The thus injected magnetic helicity in the system can be redistributed by magnetic helicity fluxes down the gradient of the local current helicity of the small-scale magnetic field. A weak reversal of the axial magnetic field is found to be a consequence of the magnetic helicity flux in the system. Such fluxes are known to alleviate so-called catastrophic quenching of the α effect in astrophysical applications. A stronger field reversal can be obtained if there is also a significant kinetic α effect. Application to the reversed field pinch in plasma confinement devices is discussed.

Non-uniformity effects in the negative effective magnetic pressure instability

In direct numerical simulations of strongly stratified turbulence we have previously studied the development of large scale magnetic structures starting from a uniform background field. This is caused by an instability resulting from a negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure, and was qualitatively reproduced in mean-field simulations (MFS) where this pressure reduction was modeled as a function of the mean magnetic field normalized by the equipartition field. We now investigate the effect of mean current density on the turbulent pressure reduction. In our MFS, such currents are associated with sharp gradients of the growing structures. We find that an enhanced mean current density increases the suppression of the turbulent pressure.

The negative magnetic pressure effect in stratified turbulence

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux tubes at the bottom of the convection zone with field strengths two orders of magnitude above equipartition, remain numerically unverified at best. As such, in recent years the idea of a formation of sunspots near the top of the convection zone has generated some interest. The presence of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective transport coefficient formalism in the mean-field formulation. The question is what happens to these coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive behaviour. In the present work we investigate the formation of large-scale magnetic structures as a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in the presence of stratification. In direct numerical simulations of forced turbulence in a stratified box, we verify the existence of this effect. This phenomenon can cause formation of large-scale magnetic structures even from initially uniform large-scale magnetic field.

Turbulent magnetic pressure instability in stratified turbulence

A reduction of total mean turbulent pressure due to the presence of magnetic fields was previously shown to be a measurable effect in direct numerical simulations. However, in the studied parameter regime the formation of large-scale structures, as anticipated from earlier mean-field simulations, was not found. An analysis of the relevant mean-field parameter dependency and the parameter domain of interest is conducted in order to clarify this apparent discrepancy.