Engineering Papers⌕ Search

Engineering topics

Drake, J. F.

Publications and source records attributed to Drake, J. F..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

On the Short-scale Spatial Variability of Electron Inflows in Electron-only Magnetic Reconnection in the Turbulent Magnetosheath Observed by MMS

We investigate the detailed properties of electron inflow in an electron-only reconnection event observed by the four Magnetospheric Multiscale (MMS) spacecraft in the Earth's turbulent magnetosheath downstream of the quasi-parallel bow shock. The lack of ion coupling was attributed to the small-scale sizes of the current sheets, and the observed bidirectional super-Alfvénic electron jets indicate that the MMS spacecraft crossed the reconnecting current sheet on both sides of an active X-line. Remarkably, the MMS spacecraft observed the presence of large asymmetries in the two electron inflows, with the inflows (normal to the current sheet) on the two sides of the reconnecting current layer differing by as much as a factor of four. Furthermore, even though the four MMS spacecraft were separated by less than seven electron skin depths, the degree of inflow asymmetry was significantly different at the different spacecraft. The asymmetry in the inflow speeds was larger with increasing distances downstream from the reconnection site, and the asymmetry was opposite on the two sides of the X-line. We compare the MMS observations with a 2D kinetic particle-in-cell (PIC) simulation and find that the asymmetry in the inflow speeds stems from in-plane currents generated via the combination of reconnection-mediated inflows and parallel flows along the magnetic separatrices in the presence of a large guide field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Faster Form of Electron Magnetic Reconnection with a Finite Length X-Line

Observations in Earth’s turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Energy Flux Densities near the Electron Dissipation Region in Asymmetric Magnetopause Reconnection

Magnetic reconnection is of fundamental importance to plasmas because of its role in releasing and repartitioning stored magnetic energy. Previous results suggest that this energy is predominantly released as ion enthalpy flux along the reconnection outflow. Using Magnetospheric Multiscale data we find the existence of very significant electron energy flux densities in the vicinity of the magnetopause electron dissipation region, orthogonal to the ion energy outflow. Finally, these may significantly impact models of electron transport, wave generation, and particle acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Magnetic Reconnection in Three Dimensions: Modeling and Analysis of Electromagnetic Drift Waves in the Adjacent Current Sheet

Magnetic reconnection at the subsolar magnetopause is persistently accompanied by strong fluctuations of the magnetic field (B), plasma density (n), and all components of the electric field (E) and current (J). The strongest fluctuations are at frequencies below the lower hybrid frequency and observed in a thin, intense current sheet adjacent to the electron diffusion region. In this current sheet, the background magnitudes of B and n are changing considerably, creating an inhomogeneous plasma environment. In this paper, we show that the fluctuations in B and n are consistent with an oscillatory displacement of the current sheet in the surface normal direction. The displacement is propagating parallel to the magnetic reconnection X-line. Wavelengths are on the order of or longer than the thickness of the current sheet to which they are confined, so we label these waves electromagnetic drift waves. E and J fluctuations are more complex than a simple displacement. They have significant variations in the component along B, which suggest that the drift waves also may be confined along B. The current sheet is supported by an electron drift driven by normal electric field, which, in turn, is balanced by an ion pressure gradient. We suggest that wave growth comes from an instability related to the drift between the electrons and ions. We discuss the possibility that drift waves can displace or penetrate into the electron diffusion region giving magnetic reconnection three-dimensional structure. Drift waves may corrugate the X-line, possibly breaking the X-line and generating turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Instabilities and turbulence in low-$β$ guide field reconnection exhausts with kinetic Riemann simulations

The role of turbulence in low-$β$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D, the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress the electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD is controlled by the ratio of the electron thermal speed to Alfvén speed, with the lower thermal speed being more unstable. In the 3D simulations, an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfvén speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with 1 $\gg$ $β >$ $m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

Here using 2.5 dimensional kinetic particle-in-cell simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasisteady reconnection rates. Note that in these simulations, the ion Larmor radius is comparable to the ion inertial length. Moreover, the transition to a more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasisteady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with the exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfvén speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments, we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Particle heating and energy partition in low- β guide field reconnection with kinetic Riemann simulations

Kinetic Riemann simulations have been completed to investigate particle heating during guide field reconnection in the low-β environment of the inner heliosphere and the solar corona. The reconnection exhaust is bounded by two rotational discontinuities (RDs), and two slow shocks (SSs) form within the exhaust as in magnetohydrodynamic (MHD) models. At the RDs, ions are accelerated by the magnetic field tension to drive the reconnection outflow as well as flows in the out-of-plane direction. The out-of-plane flows stream toward the midplane and meet to drive the SSs. The SSs differ greatly from those in the MHD model. The turbulence at the shock fronts and both upstream and downstream is weak, and so the shocks are laminar and produce little dissipation. Downstream of the SSs, the counterstreaming ion beams lead to higher density, which leads to a positive potential between the SSs which acts to confine the downstream electrons to maintain charge neutrality. The potential accelerates electrons from upstream of the SSs to the downstream region and traps a small fraction but only modestly increases the downstream electron temperature above the upstream value. In the low-β limit, the released magnetic energy is split between bulk flow and ion heating with little energy going to electrons. That the model produces neither strong electron heating nor an energetic electron component implicates that other mechanisms, such as multiple x-line reconnection, are required to explain energetic electron production in large flares. The model can be tested with the expected data from the Parker Solar Probe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A computational model for exploring particle acceleration during reconnection in macroscale systems

A new computational model suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems is presented. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfvén waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Wave Generation and Heat Flux Suppression in Astrophysical Plasma Systems

Heat flux suppression in collisionless plasmas for a large range of plasma β is explored using two-dimensional particle-in-cell simulations with a strong, sustained thermal gradient. We find that a transition takes place between whistler-dominated (high-β) and double-layer-dominated (low-β) heat flux suppression. Whistlers saturate at small amplitude in the low beta limit and are unable to effectively suppress the heat flux. Electrostatic double layers (DLs) suppress the heat flux to a mostly constant factor of the free-streaming value once this transition happens. The DL physics is an example of ion–electron coupling and occurs on a scale of roughly the electron Debye length. The scaling of ion heating associated with the various heat flux driven instabilities is explored over the full range of β explored. The range of plasma-βs studied in this work makes it relevant to the dynamics of a large variety of astrophysical plasmas, including the intracluster medium of galaxy clusters, hot accretion flows, stellar and accretion disk coronae, and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Characterizing Ion Flows Across a Magnetotail Dipolarization Jet: Characterizing Ion Flows Across a Magnetotail Dipolarization Jet

The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the magnetic field and Earth-directed ion flow structure depend on the dawn-dusk width. A reversal in the usual Hall magnetic field near the center of the current sheet on the duskside of larger jets is observed. For small widths, the maximum velocity of the earthward flow is significantly reduced below the theoretical limit of the upstream Alfvén speed. However, the ion flow speed approaches this limit once the width exceeds the ion Larmor radius based on the normal magnetic field, $B^z$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Localized and Intense Energy Conversion in the Diffusion Region of Asymmetric Magnetic Reconnection

We assess a high-resolution simulation of magnetopause reconnection observed by the Magnetospheric Multiscale mission and explain the occurrence of strongly localized dissipation with an amplitude more than an order of magnitude larger than expected. Unlike symmetric reconnection, wherein reconnection of the ambient reversed magnetic field drives the dissipation, we discover that the annihilation of the self-generated, out-of-plane (Hall) magnetic field plays the dominant role. Electrons flow along the magnetosheath separatrices, converge in the diffusion region, and jet past the X-point into the magnetosphere. The resulting accumulation of negative charge generates intense parallel electric fields that eject electrons along the magnetospheric separatrices and produce field-aligned beams. Many of these features match Magnetospheric Multiscale observations.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2018

Electron magnetic reconnection without ion coupling in Earth’s turbulent magnetosheath

Magnetic reconnection in current sheets is a magnetic-to-particle energy conversion process that is fundamental to many space and laboratory plasma systems. In the standard model of reconnection, this process occurs in a minuscule electron-scale diffusion region. On larger scales, ions couple to the newly reconnected magnetic-field lines and are ejected away from the diffusion region in the form of bi-directional ion jets at the ion Alfvén speed. Much of the energy conversion occurs in spatially extended ion exhausts downstream of the diffusion region. In turbulent plasmas, which contain a large number of small-scale current sheets, reconnection has long been suggested to have a major role in the dissipation of turbulent energy at kinetic scales. However, evidence for reconnection plasma jetting in small-scale turbulent plasmas has so far been lacking. Here we report observations made in Earth’s turbulent magnetosheath region (downstream of the bow shock) of an electron-scale current sheet in which diverging bi-directional super-ion-Alfvénic electron jets, parallel electric fields and enhanced magnetic-to-particle energy conversion were detected. Contrary to the standard model of reconnection, the thin reconnecting current sheet was not embedded in a wider ion-scale current layer and no ion jets were detected. Observations of this and other similar, but unidirectional, electron jet events without signatures of ion reconnection reveal a form of reconnection that can drive turbulent energy transfer and dissipation in electron-scale current sheets without ion coupling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Nonlinear Electrostatic Steepening of Whistler Waves: The Guiding Factors and Dynamics in Inhomogeneous Systems

Whistler mode chorus waves are particularly important in outer radiation belt dynamics due to their key role in controlling the acceleration and scattering of electrons over a very wide energy range. The efficiency of wave-particle resonant interactions is defined by whistler wave properties which have been described by the approximation of plane linear waves propagating through the cold plasma of the inner magnetosphere. However, recent observations of extremely high-amplitude whistlers suggest the importance of nonlinear wave-particle interactions for the dynamics of the outer radiation belt. Oblique chorus waves observed in the inner magnetosphere often exhibit drastically nonsinusoidal (with significant power in the higher harmonics) waveforms of the parallel electric field, presumably due to the feedback from hot resonant electrons. We have considered the nature and properties of such nonlinear whistler waves observed by the Van Allen Probes and Time History of Events and Macroscale Interactions define during Substorms in the inner magnetosphere, and we show that the significant enhancement of the wave electrostatic component can result from whistler wave coupling with the beam-driven electrostatic mode through the resonant interaction with hot electron beams. Being modulated by a whistler wave, the electron beam generates a driven electrostatic mode significantly enhancing the parallel electric field of the initial whistler wave. Finally, we confirm this mechanism using a self-consistent particle-in-cell simulation. The nonlinear electrostatic component manifests properties of the beam-driven electron acoustic mode and can be responsible for effective electron acceleration in the inhomogeneous magnetic field.

58 GEOSCIENCES↗
DOE OSTI2018

Suppression of Electron Thermal Conduction by Whistler Turbulence in a Sustained Thermal Gradient

The dynamics of weakly magnetized collisionless plasmas in the presence of an imposed temperature gradient along an ambient magnetic field is explored with particle-in-cell simulations and modeling. Two thermal reservoirs at different temperatures drive an electron heat flux that destabilizes off-angle whistler-type modes. The whistlers grow to large amplitude, δB/B 0 ≃ 1, and resonantly scatter the electrons, significantly reducing the heat flux. Interestingly, the resulting steady-state heat flux is largely independent of the thermal gradient. The rate of thermal conduction is instead controlled by the finite propagation speed of the whistlers, which act as mobile scattering centers that convect the thermal energy of the hot reservoir. The findings are relevant to thermal transport in high-β astrophysical plasmas such as hot accretion flows and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Super-Alfvénic Propagation and Damping of Reconnection Onset Signatures

The quadrupolar out-of-plane Hall magnetic leld generated during collisionless reconnection propagates away from the x line as a kinetic Alfvén wave (KAW). While it has been shown that this KAW carries substantial Poynting flux and propagates super-Alfvenically, how this KAW damps as it propagates away from the x line is not well understood. In this study, this damping is examined using kinetic particle-in-cell simulations of antiparallel symmetric magnetic reconnection in a one-dimensional current sheet equilibrium. In the reconnection simulations, the KAW wave vector has a typical magnitude comparable to an inverse fluid Larmor radius (effectively an inverse ion Larmor radius) and a direction of 85–89°relative to the local magnetic field. We find that the damping of the reconnection KAW is consistent with linear Landau damping results from a numerical Vlasov dispersion solver. This knowledge allows us to generalize our damping predictions to regions in the magnetotail and solar corona where the magnetic geometry can be approximated as a current sheet. For the magneto tail, the KAW from reconnection will not damp away before propagating the approximately 20 Earth radii associated with global magneto tail distances. For the solar corona, on the other hand, these KAWs will completely damp before reaching the distances comparable to the flare loop length.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Effect of a Guide Field on Local Energy Conversion During Asymmetric Magnetic Reconnection: Particle-in-Cell Simulations

We use theory and simulations to study how the out-of-plane (guide) magnetic field strength modifies the location where the energy conversion rate between the electric field and the plasma is appreciable during asymmetric magnetic reconnection, motivated by observations (Genestreti et al., 2017). For weak guide fields, energy conversion is maximum on the magnetospheric side of the X line, midway between the X line and electron stagnation point. As the guide field increases, the electron stagnation point gets closer to the X line, and energy conversion occurs closer to the electron stagnation point. We motivate one possible nonrigorous approach to extend the theory of the stagnation point location to include a guidefield. The predictions are compared to two-dimensional particle-in-cell (PIC) simulations with vastly different guide fields. The simulations have upstream parameters corresponding to three events observed with Magnetospheric Multiscale (MMS). The predictions agree reasonably well with the simulation results,capturing trends with the guide field. The theory correctly predicts that the X line and stagnation points approach each other as the guide field increases. The results are compared to MMS observations, Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) observations of each event,and a global resistive-magnetohydrodynamics simulation of the 16 October 2015 event. The PIC simulation results agree well with the global observations and simulation but differ in the strong electric fields and energy conversion rates found in MMS observations. The observational, theoretical, and numerical resultssuggest that the strong electric fields observed by MMS do not represent a steady global reconnection rate.

79 ASTRONOMY AND ASTROPHYSICS↗