Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell (PIC)”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗
NASA NTRS1989

A general concurrent algorithm for plasma particle-in-cell simulation codes

The general concurrent particle-in-cell (GCPIC) algorithm has been used to implement an electrostatic particle-in-cell code on a 32-node hypercube parallel computer. The GCPIC algorithm decomposes the PIC code by dividing the particle simulation physical domain into subdomains that are equal in number to the number of processors; all subdomains will accordingly possess approximately equal numbers of particles. The portion of the code which updates particle positions and velocities is nearly 100 percent efficient when the number of particles increases linearly with that of hypercube processors.

Liewer, Paulett C.↗
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↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Relativistic Particle-In-Cell Simulation Studies of Prompt and Early Afterglows from GRBs

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities {e.g., the Weibel (filamentation) 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.-/

Nishikawa, Ken-Ichi↗
NASA NTRS2008

New Relativistic Particle-In-Cell Simulation Studies of Prompt and Early Afterglows from GRBs

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electro-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the collisionless relativistic shock particle acceleration is due to plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel (filamentation) 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.

Nishikawa, Ken-ichi↗
NASA NTRS2007

Relativistic Particle-in-Cell Simulation Studies of Prompt and Early Afterglows Observed by GLAST

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet. 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.

Mizuno, Y.↗
DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗
NASA NTRS2006

Particle Demagnetization in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS2006

How Does Collisionless Magnetic Reconnection Work in the Presence of a Guide Magnetic Field?

The dissipation mechanism of guide field magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as 2.5 and three-dimensional PIC simulations of guide-field magnetic reconnection. We will show that diffusion region scale sizes in moderate and large guide field cases are determined by electron Larmor radii, and that analytical estimates of diffusion region dimensions need to include description of the heat flux tensor. The dominant electron dissipation process appears to be based on thermal electron inertia, expressed through nongyrotropic electron pressure tensors. We will argue that this process remains viable in three dimensions by means of a detailed comparison of high resolution particle-in-cell simulations.

Hesse, Michael↗
NASA NTRS2014

Open Boundary Particle-in-Cell Simulation of Dipolarization Front Propagation

First results are presented from an ongoing open boundary 2-1/2D particle-in-cell simulation study of dipolarization front (DF) propagation in Earth's magnetotail. At this stage, this study is focused on the compression, or pileup, region preceding the DF current sheet. We find that the earthward acceleration of the plasma in this region is in general agreement with a recent DF force balance model. A gyrophase bunched reflected ion population at the leading edge of the pileup region is reflected by a normal electric field in the pileup region itself, rather than through an interaction with the current sheet. We discuss plasma wave activity at the leading edge of the pileup region that may be driven by gradients, or by reflected ions, or both; the mode has not been identified. The waves oscillate near but above the ion cyclotron frequency with wavelength several ion inertial lengths. We show that the waves oscillate primarily in the perpendicular magnetic field components, do not propagate along the background magnetic field, are right handed elliptically (close to circularly) polarized, exist in a region of high electron and ion beta, and are stationary in the plasma frame moving earthward. We discuss the possibility that the waves are present in plasma sheet data, but have not, thus far, been discovered.

dipolarization front↗
DOE OSTI2025

Azimuthal and axial structures in 3D particle-in-cell simulation of Penning discharge

We report the results of a 3D particle-in-cell simulation of cylindrical Penning discharge with the axial magnetic field in the so-called reflex configuration, where the cathode and anti-cathode are biased to the same negative potential. The discharge is supported by thermal electron emission from the cathode. Electron and ion collisions, including ionization, are fully accounted for. The emphasis is on a specific regime in which the plasma potential at the center of the discharge is positive with respect to the chamber walls, serving as an anode. Spatial and temporal scales of the observed azimuthal and axial fluctuations and structures are characterized. It is suggested that azimuthal structures are caused by the dissipative gradient-drift instability. We find that the axial fluctuations related to the plasma-beam instabilities are weakly correlated with the azimuthal perturbations of the density, so that the azimuthal modes rotate as a whole and do not show any axial shear. The behavior of the electric potential is more involved, demonstrating intermittent standing wave and propagating structures in the axial direction that modulate the electron transport, producing the standing wave pattern (along the z-direction) in the radial electron flux.

3D PIC simulation of discharge↗
DOE OSTI2026

Proton beam enhancements from ultrafast laser interactions with compound parabolic concentrators

Target normal sheath acceleration (TNSA) is among the most widely studied laser–plasma ion acceleration mechanisms. In this Letter, we report on studies of proton acceleration from flat Cu targets with cone-like focusing structures called compound parabolic concentrators (CPCs) coupled to their front surface. The CPC acts as a non-imaging focusing optic that enhances the laser intensity at the target's front surface by ∼2×, from 9 × 10 18 to 2 × 10 19 W/cm 2 . This effect drives a reduction in the particle source size from the typical laser spot size of 120 μm to the CPC tip size of 65 μm and significant enhancements in the peak energy and temperature (2.6 × increase) of the resulting TNSA proton beam. We also observe increased opening angle of the beam. 2D PIC simulations have been conducted and replicate the experimental behavior. These increases have implications for long focal length facilities that require higher energy proton beams, which now can be achieved without significant infrastructure changes.

3D printing↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
NASA NTRS2008

Production of Magnetic Turbulence by Cosmic Rays Drifting Upstream of Supernova Remnant Shocks

I will present results of our recent two- and three-dimensional Particle-In-Cell simulations of magnetic-turbulence production by cosmic-ray ions drifting upstream of supernova remnant shocks. These studies' aim is twofold: test recent predictions of strong amplification in short wavelength, non-resonant wave modes, and study the subsequent evolution of the magnetic turbulence, including its backreaction on cosmic-ray trajectories. We confirm that the drifting cosmic rays give rise to a turbulent magnetic field, but show that an oblique filamentary mode grows more rapidly than the non-resonant parallel modes found in analytical theory. The field perturbations grow more slowly than estimated using a quasi-linear analytical approach for the parallel plane-wave mode, and saturate in amplitude at deltaB/B approximately equal to 1. The backreaction of the magnetic turbulence on the particles leads to an alignment of the bulk-flow velocities of the cosmic rays and the background medium. This is an essential characteristic of cosmic ray-modified shocks: the upstream flow speed is continuously changed by the cosmic rays. The reduction of relative drift between cosmic rays and background medium accounts for the saturation of the instability at only moderate magnetic-field amplitudes. It is possible that the prolonged magnetic field growth observed in recent MHD simulations results from a cosmic-ray current assumed to be constant and thus immune to the backreaction from the turbulent field. We speculate that the parallel plane-wave mode found in analytical treatments very quickly leads co filamentation, which we observe in our PIC modeling and is also apparent in the MHD simulations.

Stroman, Thomas↗