Magnetic Helicity Flux Research Articles

Magnetic helicity fluxes in an α 2 dynamo embedded in a halo

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical α quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with the earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that “softer” (halo) boundary conditions give a lower energy of the saturated mean magnetic field.

Can catastrophic quenching be alleviated by separating shear and α effect?

The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic α operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the α layer. We perform a systematic study of a simple axisymmetric two-layer αΩ dynamo in a spherical shell for Reynolds numbers in the range 1 ≤ R m ≤ 105. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary.

Trend of photospheric magnetic helicity flux in active regions generating halo coronal mass ejections

[NASA ADS]

Shear-driven and diffusive helicity fluxes in αΩ dynamos

We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux.

Magnetic helicity fluxes in αΩ dynamos

AbstractIn turbulent dynamos the production of large-scale magnetic fields is accompanied by a separation of magnetic helicity in scale. The large- and small-scale parts increase in magnitude. The small-scale part can eventually work against the dynamo and quench it, especially at high magnetic Reynolds numbers. A one-dimensional mean-field model of a dynamo is presented where diffusive magnetic helicity fluxes within the domain are important. It turns out that this effect helps to alleviate the quenching. Here we show that internal magnetic helicity fluxes, even within one hemisphere, can be important for alleviating catastrophic quenching.

TIME EVOLUTION OF CORONAL MAGNETIC HELICITY IN THE FLARING ACTIVE REGION NOAA 10930

To study the three-dimensional (3D) magnetic field topology and its long-term evolution associated with the X3.4 flare of 2006 December 13, we investigate the coronal relative magnetic helicity in the flaring active region (AR) NOAA 10930 during the time period of December 8-14. The coronal helicity is calculated based on the 3D nonlinear force-free magnetic fields reconstructed by the weighted optimization method of Wiegelmann, and is compared with the amount of helicity injected through the photospheric surface of the AR. The helicity injection is determined from the magnetic helicity flux density proposed by Pariat et al. using Solar and Heliospheric Observatory/Michelson Doppler Imager magnetograms. The major findings of this study are the following. (1) The time profile of the coronal helicity shows a good correlation with that of the helicity accumulation by injection through the surface. (2) The coronal helicity of the AR is estimated to be -4.3$\times$10$^{43}$ Mx$^{2}$ just before the X3.4 flare. (3) This flare is preceded not only by a large increase of negative helicity, -3.2$\times$10$^{43}$ Mx$^{2}$, in the corona over ~1.5 days but also by noticeable injections of positive helicity though the photospheric surface around the flaring magnetic polarity inversion line during the time period of the channel structure development. We conjecture that the occurrence of the X3.4 flare is involved with the positive helicity injection into an existing system of negative helicity.

Open and closed boundaries in large-scale convective dynamos

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect-conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime increases as Rm^(1/2). Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

THE PHOTOSPHERIC ENERGY AND HELICITY BUDGETS OF THE FLUX-INJECTION HYPOTHESIS

The flux-injection hypothesis for driving coronal mass ejections (CMEs) requires the transport of substantial magnetic energy and helicity flux through the photosphere concomitant with the eruption. Under the magnetohydrodynamics approximation, these fluxes are produced by twisting magnetic field and/or flux emergence in the photosphere. A CME trajectory, observed 2000 September 12 and fitted with a flux-rope model constrains energy and helicity budgets for testing the flux-injection hypothesis. Optimal velocity profiles for several driving scenarios are estimated by minimizing the photospheric plasma velocities for a cylindrically symmetric flux-rope magnetic field subject to the flux budgets required by the flux-rope model. Ideal flux injection, involving only flux emergence, requires hypersonic upflows in excess of the solar escape velocity 617 km/s over an area of 6x10^8 km^2 to satisfy the energy and helicity budgets of the flux-rope model. These estimates are compared with magnetic field and Doppler measurements from Solar Heliospheric Observatory/Michelson Doppler Imager on 2000 September 12 at the footpoints of the CME. The observed Doppler signatures are insufficient to account for the required energy and helicity budgets of the flux-injection hypothesis.

Scale Locality of Magnetohydrodynamic Turbulence

We investigate the scale locality of cascades of conserved invariants at high kinetic and magnetic Reynold's numbers in the "inertial-inductive range" of magnetohydrodynamic (MHD) turbulence, where velocity and magnetic field increments exhibit suitable power-law scaling. We prove that fluxes of total energy and cross helicity-or, equivalently, fluxes of Elsässer energies-are dominated by the contributions of local triads. Flux of magnetic helicity may be dominated by nonlocal triads. The magnetic stretching term may also be dominated by nonlocal triads, but we prove that it can convert energy only between velocity and magnetic modes at comparable scales. We explain the disagreement with numerical studies that have claimed conversion nonlocally between disparate scales. We present supporting data from a 1024{3} simulation of forced MHD turbulence.

Equatorial magnetic helicity flux in simulations with different gauges

AbstractWe 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 &amp; Co. KGaA, Weinheim)

The critical role of magnetic helicity in astrophysical large-scale dynamos

The role of magnetic helicity in astrophysical large-scale dynamos is reviewed and compared with cases where there is no energy supply and an initial magnetic field can only decay. In both cases magnetic energy tends to get redistributed to larger scales. Depending on the efficiency of magnetic helicity fluxes the decay of a helical field can speed up. Likewise, the saturation of a helical dynamo can speed up through magnetic helicity fluxes. The astrophysical importance of these processes is reviewed in the context of the solar dynamo and an estimated upper limit for the magnetic helicity flux of 1046 Mx2/cycle is given.

Small-scale magnetic helicity losses from a mean-field dynamo

Using mean-field models with a dynamical quenching formalism, we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity density or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong for both large- and small-scale fields. In the presence of shear, there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region, there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behaviour emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain.

THE SATURATION LIMIT OF THE MAGNETOROTATIONAL INSTABILITY

Simulations of the magnetorotational instability (MRI) in a homogeneous shearing box have shown that the asymptotic strength of the magnetic field declines steeply with increasing resolution. Here I model the MRI driven dynamo as a large scale dynamo driven by the vertical magnetic helicity flux. This growth is balanced by large scale mixing driven by a secondary instability. The saturated magnetic energy density depends almost linearly on the vertical height of the typical eddies. The MRI can drive eddies with arbitrarily large vertical wavenumber, so the eddy thickness is either set by diffusive effects, by the magnetic tension of a large scale vertical field component, or by magnetic buoyancy effects. In homogeneous, zero magnetic flux, simulations only the first effect applies and the saturated limit of the dynamo is determined by explicit or numerical diffusion. The exact result depends on the numerical details, but is consistent with previous work, including the claim that the saturated field energy scales as the gas pressure to the one quarter power (which we interpret as an artifact of numerical dissipation). The magnetic energy density in a homogeneous shearing box will tend to zero as the resolution of the simulation increases, but this has no consequences for the dynamo or for angular momentum transport in real accretion disks. The claim that the saturated state depends on the magnetic Prandtl number may also be an artifact of simulations in which microphysical transport coefficients set the MRI eddy thickness. Finally, the efficiency of the MRI dynamo is a function of the ratio of the Alfv\'en velocity to the product of the pressure scale height and the local shear. As this approaches unity from below the dynamo reaches maximum efficiency.

MAGNETIC HELICITY EXCHANGE BETWEEN NEIGHBORING ACTIVE REGIONS

Addressing the long-lasting problem of the magnetic helicity distribution in the solar corona, a proof for magnetic helicity exchange between two neighboring emerging active regions (ARs) was found: when AR 9188 emerged, it first started to accumulate positive helicity, while the later neighboring emerging AR 9192 accumulated negative helicity. At a later time, after the bright connecting loops became visible between the two ARs, AR 9188 also suddenly started to gain negative helicity. At the same time AR 9192 started to loose negative helicity. It was found that the magnetic helicity fluxes of the two ARs changed simultaneously by almost the same amount. At one instant it was even possible to determine that the connecting loop between the two ARs carried negative helicity. We excluded the possibility that magnetic flux emergence was causing the observed variation of the magnetic helicity. Hence, magnetic helicity was indeed transferred from the late emerging AR 9192 to AR 9188 via an unbalanced magnetic torque along the loop. Such kinds of helicity transfer might be a common mechanism of redistribution of magnetic helicity in the solar atmosphere, which has just not been widely observed yet.

Modelling and observations of photospheric magnetic helicity

Mounting observational evidence of the emergence of twisted magnetic flux tubes through the photosphere have now been published. Such flux tubes, formed by the solar dynamo and transported through the convection zone, eventually reach the solar atmosphere. Their accumulation in the solar corona leads to flares and coronal mass ejections. Since reconnections occur during the evolution of the flux tubes, the concepts of twist and magnetic stress become inappropriate. Magnetic helicity, as a well preserved quantity, in particular in plasma with high magnetic Reynolds number, is a more suitable physical quantity to use, even if reconnection is involved.Only recently, it has been realized that the flux of magnetic helicity can be derived from magnetogram time series. This paper reviews the advances made in measuring the helicity injection rate at the photospheric level, mostly in active regions. It relates the observations to our present theoretical understanding of the emergence process. Most of the helicity injection is found during magnetic flux emergence, whereas the effect of differential rotation is small, and the long-term evolution of active regions is still puzzling. The photospheric maps of the injection of magnetic helicity provide new spatial information about the basic properties of the link between the solar activity and its sub-photospheric roots. Finally, the newest techniques to measure photospheric flows are reviewed.

MAGNETIC RELAXATION IN THE SOLAR CORONA

This is a mathematical study of the long-lived hydromagnetic structures produced in the tenuous solar corona by the turbulent, resistive relaxation of a magnetic field under the condition of extremely high electrical conductivity. The relaxation theory of Taylor, originally developed for a laboratory device, is extended to treat the open atmosphere where the relaxing field must interact with its surrounding fields. A boundary-value problem is posed for a two-dimensional model that idealizes the corona as the half Cartesian plane filled with a potential field (1) that is anchored to a rigid, perfectly conducting base and (2) that embeds a force-free magnetic field in the form of a flux-rope oriented horizontally and perpendicular to the Cartesian plane. The flux-rope has a free boundary, which is an unknown in the construction of a solution for this atmosphere. Pairs of magnetostatic solutions are constructed to represent the initial and final states of a flux-rope relaxation that conserve both the total magnetic helicity and total axial magnetic flux, using a numerical iterative method specially developed for this study. The collection of numerical solutions found provides an insight into the interplay among several hydromagnetic properties in the formation of long-lived coronal structures. In particular, the study shows (1) that the outward spread of reconnection between a relaxing flux-rope and its external field may be arrested at some outer magnetic flux surface within which a constant-α force-free field emerges as the minimum-energy state and (2) that this outward spread is complicated by an inward, partial collapse of the relaxing flux-rope produced by a loss of internal magnetic pressure.

Paradigm shifts in solar dynamo modeling

AbstractSelected topics in solar dynamo theory are being highlighted. The possible relevance of the near-surface shear layer is discussed. The role of turbulent downward pumping is mentioned in connection with earlier concerns that a dynamo-generated magnetic field would be rapidly lost from the convection zone by magnetic buoyancy. It is argued that shear-mediated small-scale magnetic helicity fluxes are responsible for the success of some of the recent large-scale dynamo simulations. These fluxes help in disposing of excess small-scale magnetic helicity. This small-scale magnetic helicity, in turn, is generated in response to the production of an overall tilt in each Parker loop. Some preliminary calculations of this helicity flux are presented for a system with uniform shear. In the Sun the effects of magnetic helicity fluxes may be seen in coronal mass ejections shedding large amounts of magnetic helicity.

Large-scale dynamos in turbulent convection with shear

Aims. To study the existence of large-scale convective dynamos under the influence of shear and rotation. Methods. Three-dimensional numerical simulations of penetrative compressible convection with uniform horizontal shear are used to study dynamo action and the generation of large-scale magnetic fields. We consider cases where the magnetic Reynolds number is either marginal or moderately supercritical with respect to small-scale dynamo action in the absence of shear and rotation. Our magnetic Reynolds number is based on the wavenumber of the depth of the convectively unstable layer. The effects of magnetic helicity fluxes are studied by comparing results for the magnetic field with open and closed boundaries. Results. Without shear no large-scale dynamos are found even if the ingredients necessary for the α-effect (rotation and stratification) are present in the system. When uniform horizontal shear is added, a large-scale magnetic field develops, provided the boundaries are open. In this case the mean magnetic field contains a significant fraction of the total field. For those runs where the magnetic Reynolds number is between 60 and 250, an additional small-scale dynamo is expected to be excited, but the field distribution is found to be similar to cases with smaller magnetic Reynolds number where the small-scale dynamo is not excited. In the case of closed (perfectly conducting) boundaries, magnetic helicity fluxes are suppressed and no large-scale fields are found. Similarly, poor large-scale field development is seen when vertical shear is used in combination with periodic boundary conditions in the horizontal directions. If, however, open (normal-field) boundary conditions are used in the x-direction, a large-scale field develops. These results support the notion that shear not only helps to generate the field, but it also plays a crucial role in driving magnetic helicity fluxes out of the system along the isocontours of shear, thereby allowing efficient dynamo action.

Survey of Magnetic Helicity Injection in Regions Producing X‐Class Flares

Virtually all X-class flares produce a coronal mass ejection (CME), and each CME carries magnetic helicity into the heliosphere. Using magnetograms from the Michelson Doppler Imager on the Solar and Heliospheric Observatory, we surveyed magnetic helicity injection into 48 X-flare-producing active regions recorded by the MDI between 1996 July and 2005 July. Magnetic helicity flux was calculated according to the method of Chae for the 48 X-flaring regions and for 345 non-X-flaring regions. Our survey revealed that a necessary condition for the occurrence of an X-flare is that the peak helicity flux has a magnitude >6 × 1036 Mx2 s-1. X-flaring regions also consistently had a higher net helicity change during the ~6 day measurement intervals than nonflaring regions. We find that the weak hemispherical preference of helicity injection, positive in the south and negative in the north, is caused by the solar differential rotation, but it tends to be obscured by the intrinsic helicity injection, which is more disorganized and tends to be of opposite sign. An empirical fit to the data shows that the injected helicity over the range 1039-10 43 Mx2 s-1 is proportional to magnetic flux squared. Similarly, over a range of 0.3-3000 days, the time required to generate the helicity in a CME is inversely proportional to the magnetic flux squared. Most of the X-flare regions generated the helicity needed for a CME in a few days to a few hours.

Tests and Comparisons of Velocity‐Inversion Techniques

Recently, several methods that measure the velocity of magnetized plasma from time series of photospheric vector magnetograms have been developed. Velocity fields derived using such techniques can be used both to determine the fluxes of magnetic energy and helicity into the corona, which have important consequences for understanding solar flares, coronal mass ejections, and the solar dynamo, and to drive time-dependent numerical models of coronal magnetic fields. To date, these methods have not been rigorously tested against realistic, simulated data sets, in which the magnetic field evolution and velocities are known. Here we present the results of such tests using several velocity-inversion techniques applied to synthetic magnetogram data sets, generated from anelastic MHD simulations of the upper convection zone with the ANMHD code, in which the velocity field is fully known. Broadly speaking, the MEF, DAVE, FLCT, IM, and ILCT algorithms performed comparably in many categories. While DAVE estimated the magnitude and direction of velocities slightly more accurately than the other methods, MEF's estimates of the fluxes of magnetic energy and helicity were far more accurate than any other method's. Overall, therefore, the MEF algorithm performed best in tests using the ANMHD data set. We note that ANMHD data simulate fully relaxed convection in a high-β plasma, and therefore do not realistically model photospheric evolution.