Search NASA⌕ Search

SEARCH · Search NASA

Results for “PIC simulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of the Wavelength of the Plasma Waves on Cross-Field Electron Transport in Partially Magnetized Plasmas

In this article, the effects of increasing plasma density on the cross-field electron transport driven by kinetic instabilities are investigated. A 2-D collisionless particle-in-cell (PIC) simulation is performed on cases with only singly charged ions and with a mixture of singly and doubly charged ions, at different levels of plasma density. It is observed using the PIC simulations that the level of plasma density affects the plasma waves that are driven by the microinstabilities, including the electron cyclotron drift instability (ECDI) and ion–ion two-stream instability (IITSI). Here, the simulation results suggest that the enhancement of cross-field electron transport is correlated with the excitation of smaller wavelength modes and the increased amplitude of the plasma waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Kinetic Simulations of Current-Sheet Formation and Reconnection at a Magnetic X Line

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.

Black, C.↗
NASA NTRS2006

3-D RPIC simulations of relativistic jets: Particle acceleration, magnetic field generation, and emission

Nonthermal radiation observed from astrophysical systems containing (relativistic) jets and shocks, e.g., supernova remnants, active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Fermi acceleration is the mechanism usually assumed for the acceleration of particles in astrophysical environments. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet, rather than by the scattering of particles back and forth across the shock as in Fermi acceleration. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the .shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. We will review recent PIC simulations which show particle acceleration in jets.

Nishikawa, K.-I.↗
DOE OSTI2024

Proton Cyclotron Waves and Pickup Ion Ring Distribution Instabilities Upstream of Mars

Proton cyclotron waves (PCWs) upstream of Mars are thought to be associated with the instabilities of pickup ions. The instabilities of pickup ions of ring distributions have larger maximum growth rates compared with beam distributions. However, observations revealed a notably‐reduced occurrence rate of PCWs for pickup ions of ring distributions. Linear instability analysis and a corresponding two‐dimensional particle‐in‐cell (PIC) simulation are performed to investigate the instabilities of pickup ion ring distributions upstream of Mars. Linear instability analysis indicates that a pickup ion ring distribution is unstable to the ion cyclotron, ion Bernstein, and mirror instabilities. The corresponding PIC simulation confirms the linear analysis results and further demonstrates that the pickup ions are scattered toward an isotropic shell distribution by the waves excited. Interestingly, the saturation energy of the waves is much lower than that driven by a corresponding pickup ion beam distribution, and mirror waves eventually dominate the system.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2010

Hybrid-PIC Computer Simulation of the Plasma and Erosion Processes in Hall Thrusters

HPHall software simulates and tracks the time-dependent evolution of the plasma and erosion processes in the discharge chamber and near-field plume of Hall thrusters. HPHall is an axisymmetric solver that employs a hybrid fluid/particle-in-cell (Hybrid-PIC) numerical approach. HPHall, originally developed by MIT in 1998, was upgraded to HPHall-2 by the Polytechnic University of Madrid in 2006. The Jet Propulsion Laboratory has continued the development of HPHall-2 through upgrades to the physical models employed in the code, and the addition of entirely new ones. Primary among these are the inclusion of a three-region electron mobility model that more accurately depicts the cross-field electron transport, and the development of an erosion sub-model that allows for the tracking of the erosion of the discharge chamber wall. The code is being developed to provide NASA science missions with a predictive tool of Hall thruster performance and lifetime that can be used to validate Hall thrusters for missions.

Hofer, Richard R.↗
DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗
NASA NTRS2011

Current-Sheet Formation and Reconnection at a Magnetic X Line in Particle-in-Cell Simulations

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.

Black, C.↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Raman gap and collisional absorption

One of the long-standing puzzles observed in many laser-plasma experiments is the gap in the Raman backscattering spectrum. This gap is characterized by the absence of backscattered light between some critical wavelength and twice the incident laser wavelength. The latter is associated with the absolute Raman instability from the quarter-critical density surface. Supported by particle-in-cell (PIC) simulations, it is suggested that the gap can result from the collisional damping of the backscattered light. A linear analysis of the competition between the Raman growth rate and the damping rate in a non-homogenous plasma predicts the gap's existence and width as a function of the system's parameters. The theory is compared with the PIC simulations and past experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

FOS: Physical sciences↗
DOE OSTI2021

Relativistically transparent magnetic filaments: scaling laws, initial results and prospects for strong-field QED studies

Abstract Relativistic transparency enables volumetric laser interaction with overdense plasmas and direct laser acceleration of electrons to relativistic velocities. The dense electron current generates a magnetic filament with field strength of the order of the laser amplitude (>10 5 T). The magnetic filament traps the electrons radially, enabling efficient acceleration and conversion of laser energy into MeV photons by electron oscillations in the filament. The use of microstructured targets stabilizes the hosing instabilities associated with relativistically transparent interactions, resulting in robust and repeatable production of this phenomenon. Analytical scaling laws are derived to describe the radiated photon spectrum and energy from the magnetic filament phenomenon in terms of the laser intensity, focal radius, pulse duration, and the plasma density. These scaling laws are compared to 3D particle-in-cell (PIC) simulations, demonstrating agreement over two regimes of focal radius. Preliminary experiments to study this phenomenon at moderate intensity ( a 0 ∼ 30) were performed on the Texas Petawatt Laser. Experimental signatures of the magnetic filament phenomenon are observed in the electron and photon spectra recorded in a subset of these experiments that is consistent with the experimental design, analytical scaling and 3D PIC simulations. Implications for future experimental campaigns are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of a guide field on the turbulence like properties of magnetic reconnection

The effect of an external guide field on the turbulence-like properties of magnetic reconnection is studied using five different 2.5D kinetic particle-in-cell (PIC) simulations. The magnetic energy spectrum is found to exhibit a slope of approximately −5/3 in the inertial range, independent of the guide field. On the contrary, the electric field spectrum in the inertial range steepens more with the guide field and approaches a slope of −5/3. In addition, spectral analysis of the different terms of the generalized Ohm's law is performed and found to be consistent with PIC simulations of turbulence and MMS observations. Finally, the guide field effect on the energy transfer behavior is examined using the von Kármán–Howarth (vKH) equation based on incompressible Hall-MHD. The general characteristics of the vKH equation with constant rate of energy transfer in the inertial range are consistent in all the simulations. This suggests that the qualitative behavior of energy spectrum and energy transfer in reconnection are similar to that of turbulence, indicating that reconnection fundamentally involves an energy cascade.

Physics↗
DOE OSTI2023

Synchrotron Firehose Instability

We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Particle Acceleration, Magnetic Field Generation and Emission from Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), supernova remnants, and Galactic microquasar systems usually have power-law emission spectra. Fermi acceleration is the mechanism usually assumed for the acceleration of particles in astrophysical environments. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that particle acceleration occurs within the downstream jet, rather than by the scattering of particles back and forth across the shock as in Fermi acceleration. Shock acceleration' is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different spectral properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. We will review recent PIC simulations of relativistic jets and try to make a connection with observations.

Nishikawa, K.-I.↗
DOE OSTI2025

Fully plasma-based electron injector for a linear collider or XFEL

We demonstrate through high-fidelity particle-in-cell (PIC) simulations a simple approach for efficiently generating 20 + GeV electron beams with the necessary charge, energy spread, and emittance for use as an injector in a future linear collider or a next generation XFEL. A high quality injected bunch is generated by self-focusing an unmatched electron driver in a nonlinear plasma wakefield. Over pump depletion distances, the drive beam dynamics and self-loading effects lead to high energy, low-energy spread output beams. For plasma densities of 10 18 c⁢m −3 , PIC simulation results indicate that self-injected beams with 0.52 n⁢C charge can be accelerated to 20 GeV with projected core energy spreads of ≲ 1%, normalized slice emittances of 110n⁢m, peak normalized brightness of ≳ 10 19 A/m 2 /rad 2 , and transfer efficiencies of ≳ 44%.

Particle acceleration in plasmas↗
DOE OSTI2023

Effective Resistivity in Relativistic Collisionless Reconnection

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing nonthermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations, we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without a guide field. In X-points, where diverging flows result in a nondiagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain, we identify the mechanisms driving the nonideal electric field using a full Ohm law by means of a statistical analysis based on our PIC simulations. We show that the nonideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated nonuniform effective resistivity model that is negligible on global scales and becomes significant only locally in X-points. It captures the properties of collisionless reconnection with the aim of mimicking its essentials in nonideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗