No Arabic abstract
Plasma relaxation in the presence of an initially braided magnetic field can lead to self-organization into relaxed states that retain non-trivial magnetic structure. These relaxed states may be in conflict with the linear force-free fields predicted by the classical Taylor theory, and remain to be fully understood. Here, we study how the individual field line helicities evolve during such a relaxation, and show that they provide new insights into the relaxation process. The line helicities are computed for numerical resistive-magnetohydrodynamic simulations of a relaxing braided magnetic field with line-tied boundary conditions, where the relaxed state is known to be non-Taylor. Firstly, our computations confirm recent analytical predictions that line helicity will be predominantly redistributed within the domain, rather than annihilated. Secondly, we show that self-organization into a relaxed state with two discrete flux tubes may be predicted from the initial line helicity distribution. Thirdly, for this set of line-tied simulations we observe that the sub-structure within each of the final tubes is a state of uniform line helicity. This uniformization of line helicity is consistent with Taylor theory applied to each tube individually. However, it is striking that the line helicity becomes significantly more uniform than the force-free parameter.
Models for astrophysical plasmas often have magnetic field lines that leave the boundary rather than closing within the computational domain. Thus, the relative magnetic helicity is frequently used in place of the usual magnetic helicity, so as to restore gauge invariance. We show how to decompose the relative helicity into a relative field-line helicity that is an ideal-magnetohydrodynamic invariant for each individual magnetic field line, and vanishes along any field line where the original field matches the reference field. Physically, this relative field-line helicity is a magnetic flux, whose specific definition depends on the gauge of the reference vector potential on the boundary. We propose a particular `minimal gauge that depends only on the reference field and minimises this boundary contribution, so as to reveal topological information about the original magnetic field. We illustrate the effect of different gauge choices using the Low-Lou and Titov-Demoulin models of solar active regions. Our numerical code to compute appropriate vector potentials and relative field-line helicity in Cartesian domains is open source and freely available.
One of the greatest challenges in solar physics is understanding the heating of the Suns corona. Most theories for coronal heating postulate that free energy in the form of magnetic twist/stress is injected by the photosphere into the corona where the free energy is converted into heat either through reconnection or wave dissipation. The magnetic helicity associated with the twist/stress, however, is expected to be conserved and appear in the corona. In previous work we showed that helicity associated with the small-scale twists undergoes an inverse cascade via stochastic reconnection in the corona, and ends up as the observed large-scale shear of filament channels. Our ``helicity condensation model accounts for both the formation of filament channels and the observed smooth, laminar structure of coronal loops. In this paper, we demonstrate, using helicity- and energy-conserving numerical simulations of a coronal system driven by photospheric motions, that the model also provides a natural mechanism for heating the corona. We show that the heat generated by the reconnection responsible for the helicity condensation process is sufficient to account for the observed coronal heating. We study the role that helicity injection plays in determining coronal heating and find that, crucially, the heating rate is only weakly dependent on the net helicity preference of the photospheric driving. Our calculations demonstrate that motions with 100% helicity preference are least efficient at heating the corona; those with 0% preference are most efficient. We discuss the physical origins of this result and its implications for the observed corona.
Solar and stellar dynamos shed small-scale and large-scale magnetic helicity of opposite signs. However, solar wind observations and simulations have shown that some distance above the dynamo both the small-scale and large-scale magnetic helicities have reversed signs. With realistic simulations of the solar corona above an active region now being available, we have access to the magnetic field and current density along coronal loops. We show that a sign reversal in the horizontal averages of the magnetic helicity occurs when the local maximum of the plasma beta drops below unity and the field becomes nearly fully force free. Hence, this reversal is expected to occur well within the solar corona and would not directly be accessible to in-situ measurements with the Parker Solar Probe or SolarOrbiter. We also show that the reversal is associated with subtle changes in the relative dominance of structures with positive and negative magnetic helicity.
We propose a novel approach to reconstruct the surface magnetic helicity density on the Sun or sun-like stars. The magnetic vector potential is determined via decomposition of vector magnetic field measurements into toroidal and poloidal components. The method is verified using data from a non-axisymmetric dynamo model. We apply the method to vector field synoptic maps from Helioseismic and Magnetic Imager (HMI) onboard of Solar Dynamics Observatory (SDO) to study evolution of the magnetic helicity density during solar cycle 24. It is found that the mean helicity density of the non-axisymmetric magnetic field of the Sun evolves in a way which is similar to that reported for the current helicity density of the solar active regions. It has predominantly the negative sign in the northern hemisphere, and it is positive in the southern hemisphere. Also, the hemispheric helicity rule for the non-axisymmetric magnetic field showed the sign inversion at the end of cycle 24. Evolution of magnetic helicity density of large-scale axisymmetric magnetic field is different from that expected in dynamo theory. On one hand, the mean large- and small-scale components of magnetic helicity density display the hemispheric helicity rule of opposite sign at the beginning of cycle 24. However, later in the cycle, the two helicities exhibit the same sign in contrast with the theoretical expectations.
Magnetic helicity fluxes in turbulently driven alpha^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 10^5. 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 $alphaOmega$ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.