Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

UCLA parallel PIC framework

The UCLA Parallel PIC Framework (UPIC) has been developed to provide trusted components for the rapid construction of new, parallel Particle-in-Cell (PIC) codes. The Framework uses object-based ideas in Fortran95, and is designed to provide support for various kinds of PIC codes on various kinds of hardware. The focus is on student programmers. The Framework supports multiple numerical methods, different physics approximations, different numerical optimizations and implementations for different hardware. It is designed with "defensive" programming in mind, meaning that it contains many error checks and debugging helps. Above all, it is designed to hide the complexity of parallel processing. It is currently being used in a number of new Parallel PIC codes.

Norton, Charles D.↗
NASA NTRSDate not supplied

Plasmaspheric Plume Turbulence: Signature of an Electrostatic Corotation Convection Shear-Layer Instability

As demonstrated by IMAGE-EUV observations, plasmaspheric drainage plumes are a common feature emanating from the plasmasphere during periods of enhanced convection. Concurrent IMAGE-EUV and LANL-MPA observations reveal that plasmaspheric plumes are not a uniform density structure that extends out to geosynchronous orbit and beyond. Instead, in–situ observations reveal a large spatial-scale structure populated by small-scale density structures suggestive of the presence of a turbulent process. Here, we investigate the occurrence and nature of plasmaspheric drainage plumes as observed by concurrent EUV-MPA observations, as well as explore the possibility that the in–situ observed small-scale density structures are the signature of an instability produced by the sheared velocity found within the plasmaspheric layers separating the corotation of the main plasmasphere from the convection-driven flow generating the plume. Particle-in-cell simulations indicate that the generation of shear-flow driven instabilities is possible under plasmaspheric plasma conditions. The initial results of these simulations indicate that the plasmaspheric instabilities generated are electrostatic in nature, consist with in–situ plume observations.

Mark L Adrian↗
NASA NTRS2022

Modeling the Lunar Wake Response to a CME Using a Hybrid PIC Model

In the solar wind, a low-density wake region forms downstream of the nightside lunar surface. In this study, we use a series of 3D hybrid particle-in-cell simulations to model the response of the lunar wake to a passing coronal mass ejection (CME). Average plasma parameters are derived from the Wind spacecraft located at 1 au during three distinct phases of a passing halo (Earth-directed) CME on 2015 June 22. Each set of plasma parameters, representing the shock/plasma sheath, a magnetic cloud, and plasma conditions we call the mid-CME phase, are used as the time-static upstream boundary conditions for three separate simulations. These simulation results are then compared with results that use nominal solar wind conditions. Results show a shortened plasma void compared to nominal conditions and a distinctive rarefaction cone originating from the terminator during the CME's plasma sheath phase, while a highly elongated plasma void reforms during the magnetic cloud and mid-CME phases. Developments of electric and magnetic field intensification are also observed during the plasma sheath phase along the central wake, while electrostatic turbulence dominates along the plasma void boundaries and 2–3 lunar radii RM downstream in the central wake during the magnetic cloud and mid-CME phases. The simulations demonstrate that the lunar wake responds in a dynamic way with the changes in the upstream solar wind during a CME.

Anthony P. Rasca↗
NASA NTRS2022

The Effect of Solar Wind Turbulence on Parallel and Oblique Firehose Instabilities

We consider the firehose instability coexisting with the omnipresent ambient solar wind turbulence. The characteristic temporal and spatial scales of the turbulence are comparable to those of the instability. Therefore, turbulence may violate the common assumption of a uniform and stationary background used to describe instabilities and make the properties of the instabilities different. To investigate this effect, we perform threedimensional hybrid simulations with particle-in-cell ions and a quasi-neutralizing electron fluid. We find that the turbulence significantly reduces the growth rates and saturation levels of both instabilities. Comparing the cases with and without turbulence, the former results in a higher temperature anisotropy in the asymptotic marginally stable state at large times. In the former case, the distribution function averaged over the simulation box is also closer to the initial one.

S A Markovskii↗
NASA NTRS2022

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.

magnetic reconnection↗
NASA NTRS2022

Dayside Transient Phenomena and Their Impact on the Magnetosphere and Ionosphere

Dayside transients, such as hot flow anomalies, foreshock bubbles, magnetosheath jets, flux transfer events, and surface waves, are frequently observed upstream from the bow shock, in the magnetosheath, and at the magnetopause. They play a significant role in the solar wind-magnetosphere-ionosphere coupling. Foreshock transient phenomena, associated with variations in the solar wind dynamic pressure, deform the magnetopause, and in turn generates field-aligned currents (FACs) connected to the auroral ionosphere. Solar wind dynamic pressure variations and transient phenomena at the dayside magnetopause drive magnetospheric ultra low frequency (ULF) waves, which can play an important role in the dynamics of Earth’s radiation belts. These transient phenomena and their geoeffects have been investigated using coordinated in-situ spacecraft observations, spacecraft-borne imagers, ground-based observations, and numerical simulations. Cluster, THEMIS, Geotail, and MMS multi-mission observations allow us to track the motion and time evolution of transient phenomena at different spatial and temporal scales in detail, whereas ground-based experiments can observe the ionospheric projections of transient magnetopause phenomena such as waves on the magnetopause driven by hot flow anomalies or flux transfer events produced by bursty reconnection across their full longitudinal and latitudinal extent. Magnetohydrodynamics (MHD), hybrid, and particle-in-cell (PIC) simulations are powerful tools to simulate the dayside transient phenomena. This paper provides a comprehensive review of the present understanding of dayside transient phenomena at Earth and other planets, their geoeffects, and outstanding questions.

Bow shock↗
NASA NTRS2021

Lower-Hybrid Drift Waves and Their Interaction With Plasmas in A 3D Symmetric Reconnection Simulation With Zero Guide Field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where is the ion (electron) gyroradius. The short-wavelength mode with develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

Shan Wang↗
NASA NTRS2021

Reconnection and particle acceleration in 3D 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 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in 3D, we conduct a comprehensive particle-in-cell simulation study comparing 2D and 3D 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 2D, 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 3D simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3D current-sheet evolution is not always fundamentally classical reconnection with perturbing 3D 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 2D. Intriguingly, nonthermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in 2D and 3D. Though the variety of underlying current-sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Gregory R. Werner↗
NASA NTRS2021

On the Presence and Thermalization of Cold Ions in the Exhaust of Antiparallel Symmetric Reconnection

Using fully kinetic 2.5 dimensional particle-in-cell simulations of anti-parallel symmetric magnetic reconnection, we investigate how initially cold ions are captured by the reconnection process, and how they evolve and behave in the exhaust. We find that initially cold ions can remain cold deep inside the exhaust. Cold ions that enter the exhaust downstream of active separatrices, closer to the dipolarization front, appear as cold counter-streaming beams behind the front. In the off-equatorial region, these cold ions generate ion-acoustic waves that aid in the thermalization both of the incoming and outgoing populations. Closest to the front, due to the stronger magnetization, the ions can remain relatively cold during the neutral plane crossing. In the intermediate exhaust, the weaker magnetization leads to enhanced pitch angle scattering and reflection. Cold ions that enter the exhaust closer to the X line, at active separatrices, evolve into a thermalized exhaust. Here, the cold populations are heated through a combination of thermalization at the separatrices and pitch angle scattering in the curved magnetic field around the neutral plane. Depending on where the ions enter the exhaust, and how long time they have spent there, they are accelerated to different energies. The superposition of separately thermalized ion populations that have been accelerated to different energies form the hot exhaust population.

Cecilia Norgren↗
NASA NTRS2021

Magnetospheric Multiscale Observations of an Expanding Oxygen Wave in Magnetic Reconnection

Heavier plasma species such as oxygen ions can have a large impact on the magnetic reconnection process. It has been hypothesized that the acceleration of demagnetized oxygen ions by the Hall electric field will lead to the formation of an oxygen wave that expands into the exhaust. By comparing data from NASA's Magnetospheric Multiscale mission to a fully kinetic particle-in-cell simulation, we can for the first time provide observational evidence of such an expanding oxygen wave. The wave is characterized by an oxygen jet consisting of cold ions directed toward the neutral sheet associated with a density cavity. This density cavity forms as the EO are subject to collective acceleration by the Hall electric field leaving behind a region of low-density oxygen ions. Our results are important for the understanding of the role and effect of oxygen ions in magnetic reconnection.

Hakon Midthun Kolsto↗
NASA NTRS2021

The Micro‐Macro Coupling of Mass‐Loading in Symmetric Magnetic Reconnection With Cold Ions

We investigate how magnetic reconnection is influenced by an inflow of a dense cold ion population. We compare two 2.5D Particle-In-Cell simulations, one containing the cold population and one without. We find that the cold population influences the reconnection process on both global and kinetic scales, and that the dominant contribution can be explained through mass-loading. We provide an analysis of how these multiscale changes are related through kinetic processes in the ion diffusion region, the so-called micro-macro coupling of mass-loading. The inertia of the cold ion population is found to be the significant link that connects the changes on different scales. The cold and warm populations exhibit counter streaming behavior when and after the ion diffusion region reorganizes itself in response to the arrival of the cold population. This signature of the cold population should be observable by spacecraft observatories such as MMS.

Magnetic reconnection↗
NASA NTRS2023

Equilibrium Selection via Current Sheet Relaxation and Guide Field Amplification

Although there is a continuous spectrum of current sheet equilibria, how a particular equilibrium is selected by a given system remains a mystery. Yet, only a limited number of equilibrium solutions are used for analyses of magnetized plasma phenomena. Here we present the exact process of equilibrium selection, by analyzing the relaxation process of a disequilibrated current sheet under a finite guide field. It is shown via phase-space analyses and particle-in-cell simulations that the current sheet relaxes in such a way that the guide field is locally amplified, yielding a mixed equilibrium from the spectrum. Comparisons to spacecraft observations and solar wind current sheet statistics demonstrate that such mixed equilibria are ubiquitous and exist as underlying local structures in various physical environments.

Young-Dae Yoon↗
NASA NTRS2023

Electron Acceleration and Heating during Magnetic Reconnection in the Earth's Quasi-parallel Bow Shock

We perform a 2.5-dimensional particle-in-cell simulation of a quasi-parallel shock, using parameters for the Earth's bow shock, to examine electron acceleration and heating due to magnetic reconnection. The shock transition region evolves from the ion-coupled reconnection dominant stage to the electron-only reconnection dominant stage, as time elapses. The electron temperature enhances locally in each reconnection site, and ion-scale magnetic islands generated by ion-coupled reconnection show the most significant enhancement of the electron temperature. The electron energy spectrum shows a power law, with a power-law index around 6. We perform electron trajectory tracing to understand how they are energized. Some electrons interact with multiple electron-only reconnection sties, and Fermi acceleration occurs during multiple reflections. Electrons trapped in ion-scale magnetic islands can be accelerated in another mechanism. Islands move in the shock transition region, and electrons can obtain larger energy from the in-plane electric field than the electric potential in those islands. These newly found energization mechanisms in magnetic islands in the shock can accelerate electrons to energies larger than the achievable energies by the conventional energization due to the parallel electric field and shock drift acceleration. This study based on the selected particle analysis indicates that the maximum energy in the nonthermal electrons is achieved through acceleration in ion-scale islands, and electron-only reconnection accounts for no more than half of the maximum energy, as the lifetime of sub-ion-scale islands produced by electron-only reconnection is several times shorter than that of ion-scale islands.

Solar magnetic reconnection↗
NASA NTRS2023

Direct Observation of Electron Temperature Anisotropy Localized to One Separatrix during Electron-Only Magnetic Reconnection in a Laboratory Plasma

Anisotropic electron heating, Te∥/Te⊥ > 1 (relative to the local magnetic field) during electron-only magnetic reconnection with a large guide field is directly measured in a laboratory plasma through multi-dimensional incoherent Thomson scattering measurements of the electron velocity distribution function. The preferentially parallel electron heating is localized to one separatrix in the reconnection plane and anisotropies of 1.5 are observed. The localization of the heating to one separatrix and the anisotropy are reproduced with a 2D particle-in-cell simulation. The characteristics of the anisotropic heating are consistent with predictions for electron energization by the parallel reconnection electric field under strong guide field. The effective electron temperature is found to increase throughout the outflow region, a possible indication of the effects of collisions and the fully 3D nature of magnetic reconnection in the experiment.

Peiyun Shi↗
NASA NTRS2024

Enablement or Suppression of Collisionless Magnetic- Reconnection By the Background Plasma Beta and Guide Field

How magnetic reconnection is triggered or suppressed is an important outstanding problem. By considering pinching of a current sheet that has formed at non-equilibrium, we show that the background plasma beta is a major controlling factor in the onset and nature of magnetic reconnection. A high plasma beta inhibits a current sheet from pinching down to kinetic scales required for collisionless reconnection, while a low beta facilitates it. A simple adiabatic model provides a good prediction for the reconnection-enabled regions in thickness versus peak plasma beta space, which are confirmed by a series of particle-in-cell simulations with varying initial parameters. A strong dependency of the peak reconnection rate on the plasma beta is clearly predicted with reconnection being favored in low beta conditions. A finite guide field is an additional source of reconnection suppression, consistent with previous observations that reconnection requires a large enough magnetic shear angle for high-beta situations.

Young-Dae Yoon↗
NASA NTRS2003

UCLA Parallel PIC Framework

The UCLA Parallel PIC Framework (UPIC) has been developed to provide trusted components for the rapid construction of new, parallel Particl-in-Cell (PIC) codes.

particle-in-cell↗