Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasmoid Instability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Plasmoid Instability in Forming Current Sheets

The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most previous studies have focused on Sweet–Parker current sheets, which are unattainable in typical astrophysical systems. Here we derive a general set of scaling laws for the plasmoid instability in resistive and visco-resistive current sheets that evolve over time. Our method relies on a principle of least time that enables us to determine the properties of the reconnecting current sheet (aspect ratio and elapsed time) and the plasmoid instability (growth rate, wavenumber, inner layer width) at the end of the linear phase. After this phase the reconnecting current sheet is disrupted and fast reconnection can occur. The scaling laws of the plasmoid instability are not simple power laws, and they depend on the Lundquist number (S), the magnetic Prandtl number (P m ), the noise of the system (${\psi }_{0}$), the characteristic rate of current sheet evolution ($1/\tau $), and the thinning process. We also demonstrate that previous scalings are inapplicable to the vast majority of astrophysical systems. Furthermore, we explore the implications of the new scaling relations in astrophysical systems such as the solar corona and the interstellar medium. In both of these systems, we show that our scaling laws yield values for the growth rate, wavenumber, and aspect ratio that are much smaller than the Sweet–Parker–based scalings.

Sun: flares↗
DOE OSTI2018

Role of the Plasmoid Instability in Magnetohydrodynamic Turbulence

The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At a sufficiently large magnetic Reynolds number ( R m = 10 6 ), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index –2.2 as proposed by recent analytic theories. In conclusion, we also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Scalings pertaining to current sheet disruption mediated by the plasmoid instability

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. As a result, the analytical scaling relations are validated with numerical solutions of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Plasmoid instability in the semi-collisional regime

We investigate analytically and numerically the semi-collisional regime of the plasmoid instability, defined by the inequality δSP ≫ ρs ≫ δin, where δSP is the width of a Sweet-Parker current sheet, ρs is the ion sound Larmor radius, and δin is width of boundary layer that arises in the plasmoid instability analysis. Theoretically, this regime is predicted to exist if the Lundquist number S and the length of the current sheet L are such that (L/ρs ) 14/9 < S < (L/ρs ) 2 (for a sinusoidal-like magnetic configuration; for a Harris-type sheet the lower bound is replaced with (L/ρs ) 8/5 ). These bounds are validated numerically by means of simulations using a reduced gyrokinetic model (Zocco & Schekochihin, Physics of Plasmas 18, 2011) conducted with the code Viriato. Importantly, this regime is conjectured to allow for plasmoid formation at relatively low, experimentally accessible, values of the Lundquist number. Our simulations obtain plasmoid instability at values of S as low as ∼ 250. The simulations do not prescribe a Sweet-Parker sheet; rather, one is formed self-consistenly during the nonlinear evolution of the initial tearing mode configuration. This proves that this regime of the plasmoid instability is realizable, at least at the relatively low values of the Lundquist number that are accessible to current dedicated experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Plasmoid instability in the semi-collisional regime

Here, we investigate analytically and numerically the semi-collisional regime of the plasmoid instability, defined by the inequality $δ_{SP} \gg ρ_s \gg δ_{in}$, where $δ_{SP}$ is the width of a Sweet–Parker current sheet, $ρ_s$ is the ion sound Larmor radius and $δ_{in}$ is the width of the boundary layer that arises in the plasmoid instability analysis. Theoretically, this regime is predicted to exist if the Lundquist number S and the length of the current sheet L are such that ($L/ρ_s$) 14/9 < S < ($L/ρ_s$) 2 (for a sinusoidal-like magnetic configuration; for a Harris-type sheet the lower bound is replaced with ($L/ρ_s$) 8/5 ). These bounds are validated numerically by means of simulations using a reduced gyrokinetic model conducted with the code Viriato. Importantly, this regime is conjectured to allow for plasmoid formation at relatively low, experimentally accessible, values of the Lundquist number. Our simulations obtain plasmoid instability at values of S as low as ~250. The simulations do not prescribe a Sweet–Parker sheet; rather, one is formed self-consistently during the nonlinear evolution of the initial tearing mode configuration. This proves that this regime of the plasmoid instability is realizable, at least at the relatively low values of the Lundquist number that are accessible to current dedicated experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Plasmoid Instability in Evolving Current Sheets and Onset of Fast Reconnection

The scaling of the plasmoid instability maximum linear growth rate with respect to the Lundquist number S in a Sweet–Parker current sheet, ${\gamma }_{\max }\sim {S}^{1/4}$, indicates that at high S, the current sheet will break apart before it approaches the Sweet–Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We carry out a series of two-dimensional simulations and develop diagnostics to separate fluctuations from an evolving background. It is found that the fluctuation amplitude starts to grow only when the linear growth rate is sufficiently high $({\gamma }_{\max }{\tau }_{A}\gt O(1))$ to overcome advection loss and the stretching effect due to the outflow. The linear growth rate continues to rise until the sizes of plasmoids become comparable to the inner layer width of the tearing mode. At this point, the current sheet is disrupted and the instability enters the early nonlinear regime. The growth rate suddenly decreases, but the reconnection rate starts to rise rapidly, indicating that current sheet disruption triggers the onset of fast reconnection. We identify important timescales of the instability development, as well as scalings for the linear growth rate, current sheet width, and dominant wavenumber at disruption. These scalings depend not only on the Lundquist number, but also on the noise amplitude. A phenomenological model that reproduces scalings from simulation results is proposed. The model incorporates the effect of reconnection outflow, which is crucial for yielding a critical Lundquist number S c below which disruption does not occur. As a result, the critical Lundquist number S c is not a constant value, but has a weak dependence on the noise amplitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2017

Multiple-Scale Physics During Magnetic Reconnection

Magnetic reconnection is a key fundamental process in magnetized plasmas wherein the global magnetic topology is modified and stored energy is transferred from fields to particles. Reconnection is an inherently local process, and mechanisms to couple global-scale dynamics are not well understood. This dissertation explores two different mechanisms for cross-scale coupling during magnetic reconnection. As one example, we theoretically examine reconnection in a collisionless plasma using particle-in-cell simulations and demonstrate that large scale reconnection physics can couple to and drive microscopic instabilities, even in two-dimensional systems if significant scale separation exists between the Debye length and the electron skin depth. The physics underlying these instabilities is explained using simple theoretical models, and their potential connection to existing discrepancies between laboratory experiments and numerical simulations is explored. In three-dimensional systems, these instabilities are shown to generate anomalous resistivity that balances a substantial fraction of the electric field. In contrast, we also use experiments to investigate cross-scale couplings during reconnection in a collisional plasma. A leading candidate for coupling global and local scales is the hierarchical breakdown of elongated, reconnecting current sheets into numerous smaller current sheets -– the plasmoid instability. In the Magnetic Reconnection Experiment (MRX), recent hardware improvements have extended the accessible parameter space allowing for the study of long-lived, elongated current sheets. Moreover, by using Argon, reproducible and collisional plasmas are produced, which allow for a detailed statistical study of collisional reconnection. As a result, we have conclusively measured the onset of sub-ion-scale plasmoids during resistive, anti-parallel reconnection for the first time. The current sheet thickness is intermediate between ion and electron kinetic scales such that the plasma is in the Hall-MHD regime. Surprisingly, plasmoids are observed at Lundquist numbers < 100 well below theoretical predictions (> 10,000). The number of plasmoids scales with both Lundquist number and current sheet aspect ratio. The Hall quadrupolar fields are shown to suppress plasmoids. Finally, plasmoids are shown to couple local and global physics by enhancing the reconnection rate. These results are compared with prior studies of tearing and plasmoid instability, and implications for astrophysical plasmas, laboratory experiments, and theoretical studies of reconnection are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Influence of 3D plasmoid dynamics on the transition from collisional to kinetic reconnection

Within the resistive magnetohydrodynamic model, high-Lundquist number reconnection layers are unstable to the plasmoid instability, leading to a turbulent evolution where the reconnection rate can be independent of the underlying resistivity. However, the physical relevance of these results remains questionable for many applications. First, the reconnection electric field is often well above the runaway limit, implying that collisional resistivity is invalid. Furthermore, both theory and simulations suggest that plasmoid formation may rapidly induce a transition to kinetic scales, due to the formation of thin current sheets. Here, this problem is studied for the first time using a first-principles kinetic simulation with a Fokker-Planck collision operator in 3D. The low-β reconnecting current layer thins rapidly due to Joule heating before the onset of the oblique plasmoid instability. Linear growth rates for standard (k y = 0) tearing modes agree with semicollisional boundary layer theory, but the angular spectrum of oblique (|k y |>0) modes is significantly narrower than predicted. In the nonlinear regime, flux-ropes formed by the instability undergo complex interactions as they are advected and rotated by the reconnection outflow jets, leading to a turbulent state with stochastic magnetic field. In a manner similar to previous 2D results, super-Dreicer fields induce a transition to kinetic reconnection in thin current layers that form between flux-ropes. These results may be testable within new laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Relativistic resistive magnetohydrodynamic reconnection and plasmoid formation in merging flux tubes

Herein we apply the general relativistic resistive magnetohydrodynamics code BHAC to perform a 2D study of the formation and evolution of a reconnection layer in between two merging magnetic flux tubes in Minkowski space–time. Small-scale effects in the regime of low resistivity most relevant for dilute astrophysical plasmas are resolved with very high accuracy due to the extreme resolutions obtained with adaptive mesh refinement. Numerical convergence in the highly non-linear plasmoid-dominated regime is confirmed for a sweep of resolutions. We employ both uniform resistivity and non-uniform resistivity based on the local, instantaneous current density. For uniform resistivity we find Sweet–Parker reconnection, from $η$ = 10-2 down to $η$ = 10 -4 , for a reference case of magnetization $σ$ = 3.33 and plasma-$β$ = 0.1. For uniform resistivity $η$ = 5 × 10 -5 the tearing mode is recovered, resulting in the formation of secondary plasmoids. The plasmoid instability enhances the reconnection rate to $v_{rec}$ ~ 0.03$c$ compared to $v_{rec}$ ~ 0.01$c$ for $η$ = 10 -4 . For non-uniform resistivity with a base level $η_0$ = 10 -4 and an enhanced current-dependent resistivity in the current sheet, we find an increased reconnection rate of $v_{rec}$ ~ 0.1$c$. The influence of the magnetization $σ$ and the plasma-$β$ is analysed for cases with uniform resistivity $η$ = 5 × 10 -5 and $η$ = 10 -4 in a range 0.5 ≤ $σ$ ≤ 10 and 0.01 ≤ $β$ ≤ 1 in regimes that are applicable for black hole accretion discs and jets. The plasmoid instability is triggered for Lundquist numbers larger than a critical value of $S_c$ ≈ 8000.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Blob Formation and Ejection in Coronal Jets due to the Plasmoid and Kelvin–Helmholtz Instabilities

We perform 2D resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible to the formation of blobs that are ejected in the jet. Our simulation with low plasma β (Case I) shows that magnetic islands form easily and propagate upward in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 Å and 211 Å) and the cool channel (304 Å) in observations by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory . The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s{sup −1}. These parameters are similar to the properties of blobs observed in extreme-ultraviolet (EUV) jets by AIA. The Kelvin–Helmholtz instability develops in our simulation with moderately high plasma β (Case II) and leads to the formation of bright vortex-like blobs above the multiple high magnetosonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magnetosonic Mach number regions disappear. In the lower plasma β case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current-sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin–Helmholtz instability along the jet are both possible causes of the formation of blobs observed at EUV wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Observations and Modeling of the Onset of Fast Reconnection in the Solar Transition Region

Magnetic reconnection is a fundamental plasma process that plays a critical role not only in energy release in the solar atmosphere, but also in fusion, astrophysics, and other space plasma environments. One of the challenges in explaining solar observations in which reconnection is thought to play a critical role is to account for the transition of the dynamics from a slow quasi-continuous phase to a fast and impulsive energetic burst of much shorter duration. Despite the theoretical progress in identifying mechanisms that might lead to rapid onset, a lack of observations of this transition has left models poorly constrained. High-resolution spectroscopic observations from NASA's Interface Region Imaging Spectrograph now reveal tell-tale signatures of the abrupt transition of reconnection from a slow phase to a fast, impulsive phase during UV bursts or explosive events in the Sun's atmosphere. Our observations are consistent with numerical simulations of the plasmoid instability, and provide evidence for the onset of fast reconnection mediated by plasmoids and new opportunities for remote-sensing diagnostics of reconnection mechanisms on the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Regimes of magnetic reconnection in colliding laser-produced magnetized plasma bubbles

Here, we conduct a multiparametric study of driven magnetic reconnection relevant to recent experiments on colliding magnetized laser-produced plasmas using particle-in-cell simulations. Varying the background plasma density, plasma resistivity, and plasma bubble geometry, the 2D simulations demonstrate a rich variety of reconnection behaviors and show the coupling between magnetic reconnection and the global hydrodynamical evolution of the system. We consider both the collision between two radially expanding bubbles where reconnection is seeded by the pre-existing X-point and the collision between two flows in a quasi-1D geometry with initially anti-parallel fields where reconnection must be initiated by the tearing instability. At a baseline case of low-collisionality and low background density, the current sheet is strongly compressed to below scale of the ion-skin-depth scale, and rapid, multi-plasmoid reconnection results. Increasing the plasma resistivity, we observe a collisional slow-down of reconnection and stabilization of plasmoid instability for Lundquist numbers less than approximately S ~10 3 . Second, increasing the background plasma density modifies the compressibility of the plasma and can also slow down or even prevent reconnection, even in completely collisionless regimes, by preventing the current sheet from thinning down to the scale of the ion-skin depth. These results have implications for understanding recent and future experiments, and signatures for these processes for proton-radiography diagnostics of these experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Collaborative Research: Unveiling the Interplay Between Magnetic Reconnection and Turbulence — Theory and Modeling

This is the final technical report for the project "Collaborative Research: Unveiling the Interplay Between Magnetic Reconnection and Turbulence — Theory and Modeling." This project investigated the interplay between two fundamental plasma processes — turbulence and magnetic reconnection — and the associated heating and particle acceleration. In particular, we focused on a new regime where the plasmoid instability mediates the turbulent energy cascade. In this new regime, reconnecting current sheets are disrupted by the growth of plasmoids/flux ropes on time scales shorter than the typical turnover times of eddies, speeding up the energy cascade and steepening the turbulence energy spectrum. This project addressed the following outstanding open questions: (1) What are the fundamental differences between plasmoid‐mediated turbulence cascade in three dimensions (3D) and two dimensions (2D)? (2) What are the essential features of plasmoid‐mediated turbulence in weakly collisional or collisionless regimes beyond resistive MHD? (3) How does the interplay of turbulence and reconnection regulate the energy release and dissipation in space and astrophysical plasmas? To address these questions, we have carried out three investigations: (1) Three-dimensional plasmoid-mediated turbulence energy cascade; (2) plasmoid-mediated reconnection and turbulence in three-dimensional Hall MHD; (3) reconnection and heating in coronal loops. This technical report summarizes the key findings of these investigations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗