Search NASA⌕ Search

SEARCH · Search NASA

Results for “Current sheet thinning”

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 OSTI2023

Applying Magnetic Curvature to MMS Data to Identify Thin Current Sheets Relative to Tail Reconnection

Thin current sheets (TCSs) have been postulated to be a necessary precondition for reconnection onset. Magnetic reconnection X-lines in the magnetotail have been observed to be more common duskward of midnight. We take advantage of the MMS tetrahedral formation during the 2017–2020 MMS tail seasons to calculate the thickness of the cross-tail neutral sheet relative to ion gyroradius. While a similar technique was applied to Cluster data, current sheet thickness over a broader range of radial distances has not been robustly explored before this study. We compare our analysis to recent theories regarding mechanisms of tail current sheet thinning and to recent simulations. We find MMS spent more than twice as long in ion-scale TCSs in the pre-midnight sector than post-midnight, despite nearly even plasma sheet dwell time. The dawn-dusk asymmetry in the distribution of Ion Diffusion Regions, as previously reported in relation to regions of TCSs, is also analyzed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Model of a Thin Current Sheet in the Earth’s Magnetotail with a Kinetic Description of Magnetized Electrons

Steady-state two-dimensional thin current sheets (TCSs) in collisionless space plasma, similar to the current sheet of the near-Earth magnetotail, are considered. The magnetic field of the sheet is orthogonal to the current and has a nonzero normal component, the electrons are magnetized, and the ions are unmagnetized. To solve the problem of kinetic description of electrons in such current sheets in a general form, they are described by the Vlasov equation in the drift approximation, the general solution to which is found in the form of a function of three independent integrals of the system of drift motion equations. An important case is considered in which the electron guiding centers obey a Maxwell–Boltzmann distribution in a stationary electromagnetic field. The obtained results make it possible to create one- and two-dimensional numerical− analytical models of current sheets in which unmagnetized ions are described by the Vlasov equation, which should be solved numerically, whereas the contribution of magnetized electrons is taken into account analytically. To numerically solve the time-independent Vlasov equation, a new method is proposed that allows one to perform the bulk of computations by using graphic processors. On the basis of the new theory, a one-dimensional numerical−analytical model of a steady-state TCS in the near-Earth magnetotail is presented, TCS configurations are calculated, and the role of electrostatic effects and electron pressure anisotropy is analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Drift Instabilities in Thin Current Sheets Using a Two–Fluid Model With Pressure Tensor Effects

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two–fluid 10–moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10–moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. As a result, comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Three-Dimensional Magnetic Reconnection With a Spatially Confined X-Line Extent: Implications for Dipolarizing Flux Bundles and the Dawn-Dusk Asymmetry

Using 3-D particle-in-cell simulations, we study magnetic reconnection with the X-line being spatially confined in the current direction.We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial ($d_i$) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the current direction is $≲ O$(10$d_i$). When the thin current sheet extent is long enough, we find that it consists of two distinct regions; a suppressed reconnecting region (on the ion-drifting side) exists adjacent to the active region where reconnection proceeds normally as in a 2-D case with a typical fast rate value ≃ 0.1. The extent of this suppression region is $≃ O$(10$d_i$), and it suppresses reconnection when the thin current sheet extent is comparable or shorter. The time scale of current sheet thinning toward fast reconnection can be translated into the spatial scale of this suppression region, because electron drifts inside the ion diffusion region transport the reconnected magnetic flux, which drives outflows and furthers the current sheet thinning, away from this region. This is a consequence of the Hall effect in 3-D. While the existence of this suppression region may explain the shortest possible azimuthal extent of dipolarizing flux bundles at Earth, it may also explain the dawn-dusk asymmetry observed at the magnetotail of Mercury, which has a global dawn-dusk extent much shorter than that of Earth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Ion-scale Current Structures in Short Large-amplitude Magnetic Structures

Magnetic reconnection is an explosive phenomenon widely happen in plasma environments. Recent studies show that reconnection could occur in the shock transition region. We study thin current layers in the region in front of the shock and find evidence of reconnection. The analysis also shows the possible generation mechanism of the thin current sheet that further enables reconnection. A large-scale wave is further generated, and the thin current sheets are generated through compression or secondary waves. The thin current structures also lead to electron heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reconnection and particle acceleration in three-dimensional current sheet evolution in moderately magnetized astrophysical pair plasma

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Physics↗
DOE OSTI2021

Magnetic nulls in interacting dipolar fields

The prominence of nulls in reconnection theory is due to the expected singular current density and the indeterminacy of field lines at a magnetic null. Electron inertia changes the implications of both features. Magnetic field lines are distinguishable only when their distance of closest approach exceeds a distance Δ d . Electron inertia ensures Δ d ≳ c/ω pe . The lines that lie within a magnetic flux tube of radius Δ d at the place where the field strength B is strongest are fundamentally indistinguishable. If the tube, somewhere along its length, encloses a point where B = 0 vanishes, then distinguishable lines come no closer to the null than ≈ (a 2 c/ω pe ) 1/3 , where a is a characteristic spatial scale of the magnetic field. Here, the behaviour of the magnetic field lines in the presence of nulls is studied for a dipole embedded in a spatially constant magnetic field. In addition to the implications of distinguishability, a constraint on the current density at a null is obtained, and the time required for thin current sheets to arise is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron and Proton Energization in 3D Reconnecting Current Sheets in Semirelativistic Plasma with Guide Magnetic Field

Abstract Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to plasma energy. We compare fully 3D simulations with 2D in two different planes to isolate reconnection and DKI effects. In zero guide field, these processes yield distinct energy conversion signatures: ions gain more energy than electrons in 2Dxy(reconnection), while the opposite is true in 2Dyz(DKI), and the 3D result falls in between. The flux-rope instability, which occurs only in 3D, allows more magnetic energy to be released than in 2D, but the rate of energy conversion in 3D tends to be lower. Increasing the guide magnetic field strongly suppresses DKI, and in all cases slows and reduces the overall amount of energy conversion; it also favors electron energization through a process by which energy is first stored in the motional electric field of flux ropes before energizing particles. Understanding the evolution of the energy partition thus provides insight into the role of various plasma processes, and is important for modeling radiation from astrophysical sources such as accreting black holes and their jets.

Astronomy & Astrophysics↗
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 OSTI2020

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Here, we use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number S limit (above ~104), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate ~0.01–0.13 is weakly dependent on S. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a superdiffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.

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 OSTI2018

Magnetosheath Reconnection Before Magnetopause Reconnection Driven by Interplanetary Tangential Discontinuity: A Three‐Dimensional Global Hybrid Simulation With Oblique Interplanetary Magnetic Field

Abstract Terrestrial dayside dynamics associated with a southward turning, oblique interplanetary magnetic field (IMF) carried by an interplanetary tangential discontinuity (TD) is investigated by performing a three‐dimensional global‐scale hybrid simulation systematically for cases in which the incoming solar wind TD possesses various magnetic field rotation angles ΔΦ = 90° to 180° and half widths w = 2 d i 0 to w = 30 d i 0 , where d i 0 is the ion inertial length in the solar wind. Overall, the TD is compressed while being transmitted into the magnetosheath, with different compression processes downstream of the quasi‐parallel (Q‐||) and quasi‐perpendicular (Q‐ ⊥ ) shocks. It is found that magnetosheath reconnection may take place downstream of both the Q‐|| and Q‐ ⊥ shocks due to interaction of the directional TD with the bow shock and magnetopause, but the existence of magnetosheath reconnection depends on w and ΔΦ . Magnetosheath flux ropes are formed through three‐dimensional patchy reconnection in the thinned current sheet, with a longer rope lengths under a larger ΔΦ . There exists a dawn‐dusk asymmetry in the spatial extent of the flux ropes, which becomes more significant as ΔΦ decreases. When ΔΦ decreases to 90° , no reconnection flux ropes are found. Magnetopause reconnection is initiated when the magnetic fluxes with a southward turning IMF (on the sunward side) reach the magnetopause, and the magnetopause flux ropes can be mixed with the magnetosheath ones. Our simulation demonstrates that the effects of a southward turning of the IMF may not be a simple field direction change that leads to reconnection only at the magnetopause.

Guo, Zhifang↗
DOE OSTI2023

Interplay of Three-Dimensional Instabilities and Magnetic Reconnection in the Explosive Onset of Magnetospheric Substorms

Magnetospheric substorms are preceded by a slow growth phase of magnetic flux loading and current sheet thinning in the tail. Extensive data sets have provided evidence of the triggering of instabilities at substorm onset, including magnetic reconnection and ballooning instabilities. Using an exact kinetic magnetotail equilibrium we present particle-in-cell simulations which capture the explosive nature of substorms through a disruption of the dipolarization front by the ballooning instability. We use self-consistent particle tracking to determine the nonthermal particle acceleration mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Effects of plasmoid formation on sawtooth process in a tokamak

For realistic values of the Lundquist number in tokamak plasmas, 1/1 magnetic islands lead to the formation of a secondary thin current sheet, which breaks up into a chain of small magnetic islands, called plasmoids. The role of plasmoid dynamics during the sawtooth reconnection process in fusion plasmas remains an unresolved issue. In this study, systematic simulations are performed to investigate the resistive internal kink mode using the full resistive magnetohydrodynamics equations implemented in the NIMROD code in a simplified tokamak geometry. For Lundquist number S≥1.6×10 7 , the secondary current sheet is found to be unstable to plasmoids during the nonlinear resistive kink mode evolution with a critical aspect ratio of the current sheet of ~70. The merging of small plasmoids leads to the formation of a monster plasmoid that can significantly affect the primary island evolution. This may provide an explanation for the partial reconnection observed in sawtooth experiments.

Computer simulation↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗
DOE OSTI2021

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Spectral Properties and Energy Injection in Mercury's Magnetotail Current Sheet

Mercury's magnetotail hosts a thin and highly dynamic current sheet (CS), where magnetic reconnection and strong fluctuations frequently occur. Here, we statistically analyze magnetic field power spectra across 370 magnetotail CSs observed by MESSENGER. About 20% of the events are quasi-laminar, showing single power-law spectra, whereas ∼80% are turbulent, exhibiting a spectral break separating inertial and kinetic ranges. A dawn–dusk asymmetry is identified: inertial-range slopes are systematically shallower on the dawnside, whereas kinetic-range slopes are steeper, indicating more developed turbulence there, consistent with the higher occurrence of reconnection-related processes on the dawnside. Component analysis shows that the transverse components, orthogonal to the tail-aligned principal field ( B X ), display shallow slopes near −1 in the inertial range, suggesting energy injection at ion scales rather than a classical inertial range. These results demonstrate that Mercury's unique plasma environment fundamentally reshapes the initiation of turbulence and the redistribution of energy in the magnetotail.

Mercury↗
DOE OSTI2021

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗