We investigate the nonlinear evolution of the $m = 1$ resistive internal kink mode in a two-dimensional (2D) configuration containing a central region of negative current density, also known as the "current hole" setup. The finite-element code FINMHD is used to solve a reduced set of incompressible Magnetohydrodynamic (MHD) equations with a current-vorticity formulation. First, the kink instability linearly develops in agreement with the general theory of resistive internal kink mode, and it subsequently leads to the formation of a current sheet. At relatively low Lundquist number $S$, a magnetic reconnection process proceeds with a rate predicted by the Sweet-Parker regime. Conversely, when $S$ exceeds a critical value that is $S_c \simeq 10^4$, the current sheet is disrupted by the formation of plasmoids on a slightly sub-Alfv\'enic time scale. In the latter case, a stochastic reconnection regime exhibiting Petshek-type features enriched by plasmoids is reached. A relatively fast normalized reconnection rate value of order $0.02$ is also measured. Finally, we compare our results with those obtained in similar 2D previous studies using ideal MHD instabilities to initiate the process, and discuss their relevance for the general theory of plasmoid chains formation and associated fast reconnection regime.

The process of fast magnetic reconnection supported by the formation of plasmoid chains in the high Lundquist number ($S$) regime is investigated using a recently developed adaptive finite-element magnetohydrodynamic (MHD) code. We employ a two-dimensional incompressible model with a set of reduced visco-resistive MHD equations. The tilt instability setup is chosen to provide a three-step mechanism, where two curved current sheets initially form on an Alfv\'enic time scale followed by a second phase of super-Alfv\'enic growth of plasmoid chains for $S \ge S_c$ (Baty 2020). A third phase is reached where an ensuing stochastic time-dependent reconnection regime with a fast time-averaged rate independent of $S$ is obtained. We reveal the multi-scale current structures during magnetic reconnection, where merging events of plasmoids give rise to monster plasmoids with shocks bounding the outflow regions. At high enough $S$ values (typically for $S \sim 100 S_c$), a dynamical Petschek-type reconnection is achieved with pairs of slow-mode shocks emanating from a small central region containing a few plasmoids. Finally, we briefly discuss the relevance of our results to explain the flaring activity in solar corona and internal disruptions in tokamaks.

We investigate by means of two-dimensional incompressible magnetohydrodynamic (MHD) numerical simulations, the onset phase of the fast collisional magnetic reconnection regime that is supported by the formation of plasmoid chains when the Lundquist number S exceeds a critical value. The present study extends previous results obtained at magnetic Prandtl number Pm = 1 (Baty 2020) to a range of different Pm values. We use FINMHD code where a set of reduced visco-resistive MHD equations is employed to form two quasi-singular current layers as a consequence of the tilt instability. The results reinforce the conclusion that, a phase of sudden super-Alfvenic growth (when Pm is not too high) of plasmoid chains is obtained, following a previous quiescent phase during current sheet formation on a slower Alfvenic time scale. We compare our results with predictions from the general theory of the plasmoid instability. We also discuss the importance of this onset phase to reach the ensuing stochastic time-dependent reconnection regime, where a fast time-averaged rate independent of S is obtained. Finally, we briefly discuss the relevance of our results to explain the flaring activity in solar corona and internal disruptions in tokamaks.

We investigate, by means of two-dimensional incompressible magnetohydrodynamic (MHD) numerical simulations, the fast collisional magnetic reconnection regime that is supported by the formation of plasmoid chains when the Lundquist number $S$ exceeds a critical value (at magnetic Prandtl number, $P_m = 1$). A recently developed characteristic-Galerkin finite-element code, FINMHD, that is specifically designed for this aim in a reduced visco-resistive MHD framework, is employed. Contrary to previous studies, a different initial setup of two repelling current channels is chosen in order to form two quasi-singular current layers on an Alfv\'enic time scale as a consequence of the tilt instability. If $S < 5 \times 10^3$, a subsequent stationary reconnection process is obtained with a rate scaling as $S^{-1/2}$ as predicted by the classical Sweet-Parker model. Otherwise, a stochastic time-dependent reconnection regime occurs, with a fast time-averaged rate independent of $S$ and having a normalized value of $0.014$. The latter regime is triggered by the formation of two chains of plasmoids disrupting the current sheets with a sudden super-Alfv\'enic growth following a quiescent phase, in agreement with the general theory of the plasmoid instability proposed by Comisso et al. [Phys. Plasmas 23, 100702 (2016)]. Moreover, the non-monotonic dependence of the plasmoid growth rate with $S$ following an asymptotically decreasing logarithmic law in the infinite $S$-limit is confirmed. We also closely compare our results to those obtained during the development of the coalescence instability setup in order to assess the generality of the mechanism. Finally, we briefly discuss the relevance of our results to explain the flaring activity in solar corona and internal disruptions in tokamaks.

Solving the problem of fast eruptive events in magnetically dominated astrophysical plasmas requires the use of particularly well adapted numerical tools. Indeed, the central mechanism based on magnetic reconnection is determined by a complex behavior with quasi-singular forming current layers enriched by their associated small scale magnetic islands called plasmoids. A new code for the solution of two dimensional dissipative magnetohydrodynamics (MHD) equations in cartesian geometry specifically developed to this end is thus presented. A current-vorticity formulation representative of an incompressible model is chosen in order to follow the formation of the current sheets and the ensuing magnetic reconnection process. A finite element discretization using triangles with quadratic basis functions on an unstructured grid is employed, and implemented via a highly adaptive characteristic- Galerkin scheme. The adaptivity of the code is illustrated on simplified test equations and finally for magnetic reconnection associated to the non linear development of the tilt instability between two repelling current channels. Varying the Lundquist number S, has allowed to study the transition between the steady-state Sweet-Parker reconnection regime (for S < 10^4) and plasmoids dominated reconnection one (for S > 10^5). The implications for the understanding of the mechanism explaining the fast conversion of free magnetic energy in astrophysical environments such as in solar corona are briefly discussed.

We study the magneto-hydrodynamic tearing instability occurring in a double current sheet configuration when a guide field is present. This is investigated by means of resistive relativistic magneto-hydrodynamic (RRMHD) simulations. Following the dynamics of the double tearing mode (DTM), we are able to compute synthetic synchrotron spectra in the explosive reconnection phase. The pulsar striped wind model represents a site where such current sheets are formed, including a guide field. The variability of the Crab nebula/pulsar system, seen as flares, can be therefore naturally explained by the DTM explosive phase in the striped wind. Our results indicate that the Crab GeV flare can be explained by the double tearing mode in the striped wind region if the magnetization parameter $\sigma$ is around $10^5$.

Magnetic reconnection associated to the double tearing mode (DTM) is investigated by means of resistive relativistic magnetohydrodynamic (RRMHD) simulations. A linearly unstable double current sheet system in two dimensional cartesian geometry is considered. For initial perturbations of large enough longitudinal wavelengths, a fast reconnection event is triggered by a secondary instability that is structurally driven by the nonlinear evolution of the magnetic islands. The latter reconnection phase and time scale appear to weakly depend on the plasma resistivity and magnetization parameter. We discuss the possible role of such explosive reconnection dynamics to explain the MeV flares observed in the Crab pulsar nebula. Indeed the time scale and the critical minimum wavelength give constraints on the Lorentz factor of the striped wind and on the location of the emission region respectively.

The long term magnetohydrodynamic stability of magnetized transonic/supersonic jets is numerically investigated using a spatial approach. We focus on two-dimensional linearly-unstable slab configurations where the jet is embedded in a flow-aligned uniform magnetic field of weak amplitude. We compare our results with previous studies using a temporal approach where longitudinally periodic domains were adopted. The finite-volume based versatile advection code (VAC) is used to solve the full set of ideal compressible MHD equations. We follow the development of Kelvin-Helmholtz modes that are driven by a white noise perturbation continuously introduced at the jet inlet. No noticeable difference is observed in spatial simulations versus analogous temporal ones during the linear and early non-linear evolution of the configuration. However, in the case of transonic flows, a different long-term scenario occurs in our spatial runs. Indeed, after the large-scale disruption of the flow, a sheath region of enhanced magnetic field encompassing the jet core forms along the whole flow. This provides a partial stabilization mechanism leading to enhanced stability for later times, which is almost independent of the initial magnitude of the magnetic field. The implication of this mechanism for the stability of astrophysical jets is discussed.

The Kelvin-Helmholtz (KH) instability occurring in a single shear flow configuration that is embedded in a uniform flow-aligned magnetic field, is revisited by means of high resolution two-dimensional (2D) magnetohydrodynamic (MHD) simulations. First, the calculations extend previous studies of magnetized shear flows to a higher compressibility regime. The nonlinear evolution of an isolated KH billow emerging from the fastest growing linear mode for a convective sonic Mach number $M_{cs}=0.7$ layer is in many respects similar to its less compressible counterpart (Mach $M_{cs}=0.5$). In particular, the disruptive regime where locally amplified, initially weak magnetic fields, control the nonlinear saturation process is found for Alfv\'en Mach numbers $4\simlt M_A \simlt 30$. The most notable difference between $M_{cs}=0.7$ versus $M_{cs}=0.5$ layers is that higher density contrasts and fast magnetosonic shocklet structures are observed. Second, the use of adaptive mesh refinement allows to parametrically explore much larger computational domains, including up to 22 wavelengths of the linearly dominant mode. A strong process of large-scale coalescence is found, whatever the magnetic field regime. It proceeds through continuous pairing/merging events between adjacent vortices up to the point where the final large-scale vortical structure reaches the domain dimensions. This pairing/merging process is attributed to the growth of subharmonic modes and is mainly controlled by relative phase differences between them. These grid-adaptive simulations demonstrate that even in very weak magnetic field regimes ($M_A \simeq 30$), the large-scale KH coalescence process can trigger tearing-type reconnection events previously identified in cospatial current-vortex sheets.