No Arabic abstract
Magnetic reconnection plays a critical role in many astrophysical processes where high energy emission is observed, e.g. particle acceleration, relativistic accretion powered outflows, pulsar winds and probably in dissipation of Poynting flux in GRBs. The magnetic field acts as a reservoir of energy and can dissipate its energy to thermal and kinetic energy via the tearing mode instability. We have performed 3d nonlinear MHD simulations of the tearing mode instability in a current sheet. Results from a temporal stability analysis in both the linear regime and weakly nonlinear (Rutherford) regime are compared to the numerical simulations. We observe magnetic island formation, island merging and oscillation once the instability has saturated. The growth in the linear regime is exponential in agreement with linear theory. In the second, Rutherford regime the island width grows linearly with time. We find that thermal energy produced in the current sheet strongly dominates the kinetic energy. Finally preliminary analysis indicates a P(k) 4.8 power law for the power spectral density which suggests that the tearing mode vortices play a role in setting up an energy cascade.
This paper studies the growth rate of reconnection instabilities in thin current sheets in the presence of both resistivity and viscosity. In a previous paper, Pucci and Velli (2014), it was argued that at sufficiently high Lundquist number S it is impossible to form current sheets with aspect ratios L/a which scale as $L/asim S^alpha$ with $alpha > 1/3$ because the growth rate of the tearing mode would then diverge in the ideal limit $Srightarrowinfty$. Here we extend their analysis to include the effects of viscosity, (always present in numerical simulations along with resistivity) and which may play a role in the solar corona and other astrophysical environments. A finite Prandtl number allows current sheets to reach larger aspect ratios before becoming rapidly unstable in pile-up type regimes. Scalings with Lundquist and Prandtl numbers are discussed as well as the transition to kinetic reconnection
We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $Stoinfty$. Our simulations confirm that the tearing instability growth rate can be as fast as $gammaapprox 0.6,{tau_A}^{-1}$, where $tau_A$ is the ideal Alfvenic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvenic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.
The recent realization that Sweet-Parker current sheets are violently unstable to the secondary tearing (plasmoid) instability implies that such current sheets cannot occur in real systems. This suggests that, in order to understand the onset of magnetic reconnection, one needs to consider the growth of the tearing instability in a current layer as it is being formed. Such an analysis is performed here in the context of nonlinear resistive MHD for a generic time-dependent equilibrium representing a gradually forming current sheet. It is shown that two onset regimes, single-island and multi-island, are possible, depending on the rate of current sheet formation. A simple model is used to compute the criterion for transition between these two regimes, as well as the reconnection onset time and the current sheet parameters at that moment. For typical solar corona parameters this model yields results consistent with observations.
Using fully kinetic simulations, we study the suppression of asymmetric reconnection in the limit where the diamagnetic drift speed >> Alfven speed and the magnetic shear angle is moderate. We demonstrate that the slippage between electrons and the magnetic flux facilitates reconnection, and can even result in fast reconnection that lacks one of the outflow jets. Through comparing a case where the diamagnetic drift is supported by the temperature gradient with a companion case that has a density gradient instead, we identify a robust suppression mechanism. The drift of the x-line is slowed down locally by the asymmetric nature of the current sheet and the resulting tearing modes, then the x-line is run over and swallowed by the faster-moving following flux.
Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.