Search NASA⌕ Search

Engineering topics

Chen, Li-Jen

Publications and source records attributed to Chen, Li-Jen.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Collaborative Research: Vlasov–Maxwell Simulations to resolve electron heating and dissipation, in quasiperpendicular shocks (Final Technical Report)

The project is to study electron heating at quasi-perpendicular shocks. We have analyzed measurements from the Magnetosphere Multiscale (MMS) spacecraft and particle-in-cell simulations to address the question. One key discovery is that waves and are excited in the shock transition region, which further develop into complicated structures and thin current layers, which contribute to accelerate and heat electrons. We further investigated the excitation and evolution of these waves and structures, and their association with ion and electron dynamics. In addition, the analysis technique of Vlasov equation terms is further applied to study electron dynamics at the dayside magnetopause, to understand the effect of density, velocity and temperature variations on electron distributions and kinetics.

58 GEOSCIENCES↗
DOE OSTI2024

Anomalous Resistivity and Electron Heating by Lower Hybrid Drift Waves during Magnetic Reconnection with a Guide Field

The lower hybrid drift wave (LHDW) has been a candidate for anomalous resistivity and electron heating inside the electron diffusion region of magnetic reconnection. In a laboratory reconnection layer with a finite guide field, quasielectrostatic LHDW (ES-LHDW) propagating along the direction nearly perpendicular to the local magnetic field is excited in the electron diffusion region. ES-LHDW generates large density fluctuations ($δn_{\text e}$, about 25% of the mean density) that are correlated with fluctuations in the out-of-plane electric field ($δE_Y$, about twice larger than the mean reconnection electric field). With a small phase difference (~ 30°) between two fluctuating quantities, the anomalous resistivity associated with the observed ES-LHDW is twice larger than the classical resistivity and accounts for 20% of the mean reconnection electric field. After we verify the linear relationship between $δn_{\text e}$ and $δE_Y$, anomalous electron heating by LHDW is estimated by a quasilinear analysis. The estimated electron heating is about 2.6 ± 0.3 MW / m 3 , which exceeds the classical Ohmic heating of about 2.0 ± 0.2 MW / m 3 . This LHDW-driven heating is consistent with the observed trend of higher electron temperatures when the wave amplitude is larger. Finally, presented results provide the first direct estimate of anomalous resistivity and electron heating power by LHDW, which demonstrates the importance of wave-particle interactions in magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Oblique propagation and temperature effects on the resonant right-hand ion beam instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2020

Electron Energy Partition across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Electron Energy Partition across Interplanetary Shocks. I. Methodology and Data Product

Analyses of 15,314 electron velocity distribution functions (VDFs) within ±2 hr of 52 interplanetary (IP) shocks observed by the Wind spacecraft near 1 au are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar VDF, while both the halo and beam/strahl components were best fit to bi-kappa VDF. This is the first statistical study to show that the core electron distribution is better fit to a self-similar VDF than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The ranges of values defined by the lower and upper quartiles for the kappa exponents are κ ec ~ 5.40–10.2 for the core, κ eh ~ 3.58–5.34 for the halo, and κ eb ~ 3.40–5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents is s ec ~ 2.00–2.04, and those of asymmetric bi-self-similar core exponents are p ec ~ 2.20–4.00 for the parallel exponent and q ec ~ 2.00–2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

MMS Measurements of the Vlasov Equation: Probing the Electron Pressure Divergence Within Thin Current Sheets

We investigate the kinetic structure of electron-scale current sheets found in the vicinity of the magnetopause and embedded in the magnetosheath within the reconnection exhaust. A new technique for computing terms of the Vlasov equation using Magnetospheric Multiscale (MMS) measurements is presented and applied to study phase space density gradients and the kinetic origins of the electron pressure divergence found within these current sheets. Crescent-shaped structures in ∇(⟂2)f(e) give rise to bipolar and quadrupolar signatures in v · ∇f(e) measured near the maximum ∇ · P(e) inside the current layers. The current density perpendicular to the magnetic field is strong (J⟂ ∼2 μA/sq.m), and the thickness of the current layers ranges from 3 to 5 electron inertial lengths. The electron flows supporting the current layers mainly result from the combination of E × B and diamagnetic drifts. We find nonzero J · E′ within the current sheets even though they are observed apart from typical diffusion region signatures.

Shuster, J. R.↗
DOE OSTI2019

Whistler wave generation by electron temperature anisotropy during magnetic reconnection at the magnetopause

Two magnetopause reconnection events of the Magnetospheric Multiscale mission with whistler wave activity are presented. The whistler mode around half of the electron cyclotron frequency is excited near the magnetospheric separatrix. In both events, there are positive correlations between the whistler wave and the lower hybrid drift instability (LHDI) activities, suggesting a possible role of LHDI in the whistler wave generation. A sudden change in the electron pitch angle distribution (PAD) function of energetic electrons is observed right after intense LHDI activity. This change in the PAD leads to temperature anisotropy in energetic electrons which is responsible for the whistler wave excitation. The measured dispersion relation demonstrates that the whistler wave propagates toward the X line nearly parallel to the magnetic field line. Furthermore, a linear analysis with the measured distribution function verifies that the whistler mode is excited by the temperature anisotropy in energetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Ion Behaviors in the Reconnection Diffusion Region of a Corrugated Magnetotail Current Sheet

The current sheet structure and ion behaviors in a magnetotail reconnection diffusion regionare investigated. The multispacecraft analysis suggests a corrugated current sheet structure, interpretedas due to a apping motion that propagates along geocentric solar magnetospheric along the +y direction inthe Geocentric Solar Magnetospheric (GSM) coordinate. The electric eld (E) and ion distributions havesimilarities with those in a planar current sheet. Energetic ions move along the current direction, suggestingthe acceleration by the observed reconnection E during the meandering motion. Counterstreaming ionsalong the current sheet normal suggest the acceleration by the Hall E that is observed to be the dominantcomponent. However, at certain locations, E and counterstreaming ions signicantly deviate from the localnormal direction, and more than one pair of counterstreaming populations exist, possibly because thecorrugated current sheet enables ions entering the current sheet at different locations with differentvelocities to mix together

Wang, Shan↗
DOE OSTI2018

Observational Evidence of Magnetic Reconnection in the Terrestrial Bow Shock Transition Region

Here, we report evidence of magnetic reconnection in the transition region of the Earth's bow shock when the angle between the shock normal and the immediate upstream magnetic field is 65°. An ion-skin-depth-scale current sheet exhibits the Hall current and field pattern, electron outflow jet, and enhanced energy conversion rate through the nonideal electric field, all consistent with a reconnection diffusion region close to the X-line. In the diffusion region, electrons are modulated by electromagnetic waves. An ion exhaust with energized field-aligned ions and electron parallel heating are observed in the same shock transition region. The energized ions are more separated from the inflowing ions in velocity above the current sheet than below, possibly due to the shear flow between the two inflow regions. The observation suggests that magnetic reconnection may contribute to shock energy dissipation.

58 GEOSCIENCES↗
NASA NTRS2018

Energy Conversion and Partition in the Asymmetric Reconnection Diffusion Region

We investigate the energy conversion and partition in the asymmetric reconnection diffusion region using two-dimensional particle-in-cell simulations and Magnetosphere Multiscale (MMS) mission observations. Under an upstream condition with equal temperatures in the two inflow regions, the simulation analysis indicates that the energy partition between ions and electrons depends on the distance from the X-line. Within the central electron diffusion region (EDR), nearly all dissipated electromagnetic field energies are converted to electrons. From the EDR to the ion diffusion region (IDR) scales, the rate of the electron energy gain decreases to be lower than that of ions. A magnetopause reconnection event inside the IDR observed by MMS shows comparable ion and electron energy gains, consistent with the simulation result in the transition region from EDR to IDR. At the EDR scale, the electron energization is mainly by the reconnection electric field (E(sub r)); in-plane electric fields (E(sub xz)) provide additional positive contributions near the X-line and do negative work on electrons beyond the EDR. The guide field reduces the electron energization by both E(sub r) and E(sub xz) in the EDR. For ion energization, E(sub r) and E(sub xz) have comparable contributions near the time of the peak reconnection rate, while E(sub xz) dominants at later time. At the IDR scale, the guide field causes asymmetry in the amount of the energy gain and energization mechanisms between two exhausts but does not have significant effects on energy partition. Our study advances understanding of ion and electron energization in asymmetric reconnect IDRs.

Wang, Shan↗
DOE OSTI2018

Whistler wave generation by anisotropic tail electrons during asymmetric magnetic reconnection in space and laboratory

Whistler wave generation near the magnetospheric separatrix during reconnection at the dayside magnetopause is studied with data from the Magnetospheric Multiscale (MMS) mission. The dispersion relation of the whistler mode is measured for the first time near the reconnection region in space, which shows that whistler waves propagate nearly parallel to the magnetic field line. A linear analysis indicates that the whistler waves are generated by temperature anisotropy in the electron tail population. This is caused by loss of electrons with a high velocity parallel to the magnetic field to the exhaust region. There is a positive correlation between activities of whistler waves and the lower-hybrid drift instability (LHDI) both in laboratory and space, indicating the enhanced transport by LHDI may be responsible for the loss of electrons with a high parallel velocity.

Dispersion relation↗
DOE OSTI2017

Electron heating and energy inventory during asymmetric reconnection in a laboratory plasma

Electron heating and the energy inventory during asymmetric reconnection are studied in the laboratory plasma with a density ratio of about 8 across the current sheet. Features of asymmetric reconnection such as the large density gradients near the low-density-side separatrices, asymmetric in-plane electric field, and bipolar out-of-plane magnetic field are observed. Unlike the symmetric case, electrons are also heated near the low-density-side separatrices. The measured parallel electric field may explain the observed electron heating. Although large fluctuations driven by lower-hybrid drift instabilities are also observed near the low-density-side separatrices, laboratory measurements and numerical simulations reported here suggest that they do not play a major role in electron energization. The average electron temperature increase in the exhaust region is proportional to the incoming magnetic energy per an electron/ion pair but exceeds scalings of the previous space observations. This discrepancy is explained by differences in the boundary condition and system size. The profile of electron energy gain from the electric field shows that there is additional electron energy gain associated with the electron diamagnetic current besides a large energy gain near the X-line. This additional energy gain increases electron enthalpy, not the electron temperature. Finally, a quantitative analysis of the energy inventory during asymmetric reconnection is conducted. Unlike the symmetric case where the ion energy gain is about twice more than the electron energy gain, electrons and ions obtain a similar amount of energy during asymmetric reconnection.

Energy inventory↗
DOE OSTI2017

Parallel electron heating in the magnetospheric inflow region: Electron Heating in Reconnection

We analyzed 20 reconnection events observed by the Magnetosphere MultiScale mission, finding that in a subregion of the magnetospheric inflow, the electron temperature parallel to the magnetic field (T e|| ) decreases toward the separatrix with increasing densities. Such Te|| variation is not consistent with the heating through the parallel potential mechanism for magnetospheric electrons. Associated with the change in T e|| , low-energy gyrotropic magnetosheath-like electrons are observed. These electrons may come from the magnetosheath, penetrate to the magnetospheric side possibly due to waves in the lower hybrid frequency range, and do not experience as much energization as those from deeper in the magnetosphere. In addition, as T e|| decreases, the high-energy field-aligned phase-space densities decrease or fluctuate significantly over time. Wave-particle interaction might produce energy conversion and cause the high-energy phase-space density variations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Two-Scale Ion Meandering Caused by the Polarization Electric Field During Asymmetric Reconnection

Ion velocity distribution functions (VDFs) from a particle-in-cell simulation of asymmetric reconnection are investigated to reveal a two-scale structure of the ion diffusion region (IDR). Ions bouncing in the inner IDR are trapped mainly by the electric field normal to the current sheet (N direction), while those reaching the outer IDR are turned back mainly by the magnetic force. The resulting inner layer VDFs have counter-streaming populations along N with decreasing counter-streaming speeds away from the midplane while maintaining the out-of-plane speed, and the outer layer VDFs exhibit crescent shapes toward the out-of-plane direction. Observations of the above VDF features and the normal electric fields provide evidence for the two-scale meandering motion.

Wang, Shan↗
NASA NTRS2016

Finite Gyroradius Effects in the Electron Outflow of Asymmetric Magnetic Reconnection

We present observations of asymmetric magnetic reconnection showing evidence of electron demagnetization in the electron outflow. The observations were made at the magnetopause by the four Magnetospheric Multiscale (MMS) spacecraft, separated by approximately 15 km. The reconnecting current sheet has negligible guide field, and all four spacecraft likely pass close to the electron diffusion region just south of the X line. In the electron outflow near the X line, all four spacecraft observe highly structured electron distributions in a region comparable to a few electron gyroradii. The distributions consist of a core with T(sub parallel) greater than T(sub perpendicular) and a nongyrotropic crescent perpendicular to the magnetic field. The crescents are associated with finite gyroradius effects of partly demagnetized electrons. These observations clearly demonstrate the manifestation of finite gyroradius effects in an electron-scale reconnection current sheet.

Norgren, C.↗