Simon Candelaresi

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.

SUPERFLARE OCCURRENCE AND ENERGIES ON G, K AND M TYPE DWARFS

Kepler data from G-, K-, and M-type stars are used to study conditions that lead to superflares with energies above 10(34) erg. From the 117,661 stars included, 380 show superflares with a total of ...

Decay of helical and nonhelical magnetic knots.

We present calculations of the relaxation of magnetic field structures that have the shape of particular knots and links. A set of helical magnetic flux configurations is considered, which we call n-foil knots of which the trefoil knot is the most primitive member. We also consider two nonhelical knots; namely, the Borromean rings as well as a single interlocked flux rope that also serves as the logo of the Inter-University Centre for Astronomy and Astrophysics in Pune, India. The field decay characteristics of both configurations is investigated and compared with previous calculations of helical and nonhelical triple-ring configurations. Unlike earlier nonhelical configurations, the present ones cannot trivially be reduced via flux annihilation to a single ring. For the n-foil knots the decay is described by power laws that range form t(-2/3) to t(-1/3), which can be as slow as the t(-1/3) behavior for helical triple-ring structures that were seen in earlier work. The two nonhelical configurations decay like t(-1), which is somewhat slower than the previously obtained t(-3/2) behavior in the decay of interlocked rings with zero magnetic helicity. We attribute the difference to the creation of local structures that contain magnetic helicity which inhibits the field decay due to the existence of a lower bound imposed by the realizability condition. We show that net magnetic helicity can be produced resistively as a result of a slight imbalance between mutually canceling helical pieces as they are being driven apart. We speculate that higher order topological invariants beyond magnetic helicity may also be responsible for slowing down the decay of the two more complicated nonhelical structures mentioned above.

Magnetic helicity transport in the advective gauge family

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e., when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

Magnetic-field decay of three interlocked flux rings with zero linking number.

The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic helicity need to be considered to characterize the decay of the field.

Self-Organizing Knotted Magnetic Structures in Plasma

We perform full-magnetohydrodynamics simulations on various initially helical configurations and show that they reconfigure into a state where the magnetic field lines span nested toroidal surfaces. This relaxed configuration is not a Taylor state, as is often assumed for relaxing plasma, but a state where the Lorentz force is balanced by the hydrostatic pressure, which is lowest on the central ring of the nested tori. Furthermore, the structure is characterized by a spatially slowly varying rotational transform, which leads to the formation of a few magnetic islands at rational surfaces. We then obtain analytic expressions that approximate the global structure of the quasistable linked and knotted plasma configurations that emerge, using maps from S^{3} to S^{2} of which the Hopf fibration is a special case. The knotted plasma configurations have a highly localized magnetic energy density and retain their structure on time scales much longer than the Alfvénic time scale.

Mimetic Methods for Lagrangian Relaxation of Magnetic Fields

We present a new code that performs a relaxation of a magnetic field toward a force-free state (Beltrami field) using a Lagrangian numerical scheme. Beltrami fields are of interest for the dynamics of many technical and astrophysical plasmas as they are the lowest energy states that the magnetic field can reach. The numerical method strictly preserves the magnetic flux and the topology of magnetic field lines. In contrast to other implementations we use mimetic operators for the spatial derivatives in order to improve accuracy for high distortions of the grid. Compared with schemes using direct derivatives we find that the final state of the simulation approximates a force-free state with a significantly higher accuracy. We implement the scheme in a code which runs on graphical processing units, which leads to an enhanced computing speed compared to previous relaxation codes.

Kinetic helicity needed to drive large-scale dynamos

Magnetic field generation on scales that are large compared with the scale of the turbulent eddies is known to be possible via the so-called α effect when the turbulence is helical and if the domain is large enough for the α effect to dominate over turbulent diffusion. Using three-dimensional turbulence simulations, we show that the energy of the resulting mean magnetic field of the saturated state increases linearly with the product of normalized helicity and the ratio of domain scale to eddy scale, provided this product exceeds a critical value of around unity. This implies that large-scale dynamo action commences when the normalized helicity is larger than the inverse scale ratio. Our results show that the emergence of small-scale dynamo action does not have any noticeable effect on the large-scale dynamo. Recent findings by Pietarila Graham et al. [Phys. Rev. E 85, 066406 (2012)] of a smaller minimal helicity may be an artifact due to the onset of small-scale dynamo action at large magnetic Reynolds numbers. However, the onset of large-scale dynamo action is difficult to establish when the kinetic helicity is small. Instead of random forcing, they used an ABC flow with time-dependent phases. We show that such dynamos saturate prematurely in a way that is reminiscent of inhomogeneous dynamos with internal magnetic helicity fluxes. Furthermore, even for very low fractional helicities, such dynamos display large-scale fields that change direction, which is uncharacteristic of turbulent dynamos.

MAGNETIC FIELD RELAXATION AND CURRENT SHEETS IN AN IDEAL PLASMA

We investigate the existence of magnetohydrostatic equilibria for topologically complex magnetic fields. The approach employed is to perform ideal numerical relaxation experiments. We use a newly developed Lagrangian relaxation scheme that exactly preserves the magnetic field topology during the relaxation. Our configurations include both twisted and sheared fields, of which some fall into the category for which Parker predicted no force-free equilibrium. The first class of field considered contains no magnetic null points, and field lines connect between two perfectly conducting plates. In these cases, we observe only resolved current layers of finite thickness. In further numerical experiments, we confirm that magnetic null points are loci of singular currents.

Braided magnetic fields: equilibria, relaxation and heating

We examine the dynamics of magnetic flux tubes containing non-trivial field line braiding (or linkage), using mathematical and computational modelling, in the context of testable predictions for the laboratory and their significance for solar coronal heating. We investigate the existence of braided force-free equilibria, and demonstrate that for a field anchored at perfectly-conducting plates, these equilibria exist and contain current sheets whose thickness scales inversely with the braid complexity—as measured for example by the topological entropy. By contrast, for a periodic domain braided exact equilibria typically do not exist, while approximate equilibria contain thin current sheets. In the presence of resistivity, reconnection is triggered at the current sheets and a turbulent relaxation ensues. We finish by discussing the properties of the turbulent relaxation and the existence of constraints that may mean that the final state is not the linear force-free field predicted by Taylors hypothesis.