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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stability of the modulator in a plasma-modulated plasma accelerator

We explore the regime of operation of the modulator stage of a recently proposed laser-plasma accelerator scheme [Phys. Rev. Lett. 127, 184801 (2021)], dubbed the plasma-modulated plasma accelerator (P-MoPA). The P-MoPA scheme offers a potential route to high-repetition-rate, GeV-scale plasma accelerators driven by picosecond-duration laser pulses from, for example, kilohertz thin-disk lasers. The first stage of the P-MoPA scheme is a plasma modulator in which a long, high-energy “drive” pulse is spectrally modulated by copropagating in a plasma channel with the low-amplitude plasma wave driven by a short, low-energy “seed” pulse. The spectrally modulated drive pulse is converted to a train of short pulses, by introducing dispersion, which can resonantly drive a large wakefield in a subsequent accelerator stage with the same on-axis plasma density as the modulator. In this paper we derive the 3D analytic theory for the evolution of the drive pulse in the plasma modulator and show that the spectral modulation is independent of transverse coordinate, which is ideal for compression into a pulse train. We then identify a transverse mode instability (TMI), similar to the TMI observed in optical fiber lasers, which sets limits on the energy of the drive pulse for a given set of laser-plasma parameters. We compare this analytic theory with particle-in-cell (PIC) simulations and find that even higher energy drive pulses can be modulated than those demonstrated in the original proposal.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic K α emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity ( ∼ 10 20 W cm − 2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement ( 1.6 × ) of K α emission from microstructured targets compared to flat foils, but with a similar intensity and profile of He α and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the K α enhancement. The reasons for the lower-than-expected enhancement of K α emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of K α emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of K α emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Expansion-Driven Self-Magnetization of High-Energy-Density Plasmas

Understanding plasma self-magnetization is one of the fundamental challenges in both laboratory and astrophysical plasmas. Self-magnetization can modify plasma transport properties, altering the dynamical evolution of plasmas. Multiple high-energy-density (HED) experiments have observed the formation of ion-scale magnetic filaments of megagauss strength, though their origin remains debated. Here, in this study, we conduct 2D collisional particle-in-cell (PIC) simulations with a laser ray-tracing module for a fully self-consistent simulation of the plasma ablation, expansion, and magnetization. The simulations use a planar geometry, effectively suppressing the Biermann magnetic fields, to focus on anisotropy-driven instabilities. The laser intensity is varied between 10 13 and 10 14 W/cm 2 , which is relevant to HED and inertial fusion experiments where collisions must be considered. We find that, above a critical intensity, the plasma rapidly self-magnetizes via an expansion-driven Weibel process, producing a plasma beta of 100 (𝛽 = 8⁢𝜋⁢𝑘 𝐵 ⁢𝑛 𝑒 ⁢𝑇 𝑒 /𝐵 2 ) and Hall parameter 𝜔 ce ⁢𝜏 𝑒 >1 within the first few hundred picoseconds. The magnetic field is sufficiently strong to modify plasma heat transport, and simulations with an artificially suppressed magnetic field show noticeably different temperature profiles.

Lezhnin, K. V. [Princeton Plasma Physics Laborator↗
DOE OSTI2024

Theoretical Assessment of the Transition Between Electron Emission Mechanisms for Nonplanar Diodes

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Conformal mapping↗
DOE OSTI2021

First Experiments and Radiographs on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) Dense Plasma Focus

We report a dense plasma focus (DPF) is a relatively compact coaxial plasma gun, which completes its discharge as a Z-pinch. These devices are designed to operate at a variety of scales to produce short (<100 ns) pulses of ions, X-rays, and/or neutrons. LLNL recently constructed and brought into operation a new device, the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF, which is designed for radiography and high-yield operations. This device has been commissioned and has achieved neutron yields of up to 3.8E11 neutrons/pulse at 2.5-MA peak current while operating at up to 1 MJ of stored energy in its original pulsed power configuration. MJOLNIR is equipped with a wide range of diagnostics, including nuclear activation detectors, neutron time-of-flight (nToF) detectors, a fast-framing camera, optical light gates, and a time-gated neutron and X-ray imager. LLNL also runs unique particle-in-cell (PIC) simulations of DPF discharges in the Chicago code and has gained significant insight into the various physical factors that influence neutron yield. MJOLNIR is one of the first DPFs whose design and continual upgrades are heavily influenced by model predictions. In this article, we describe insights from modeling, device operation, and recent results. Comparisons between modeling predictions and measurements, as well as X-ray and neutron images are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Direct Measurement of Backstreaming Ions from the Axis 2 Septum Dump and RF Generation

The DARHT-II accelerator produces a 2-kA, 17-MeV beam in a 1600-ns pulse. After exiting the accelerator, the pulse is sliced into four short pulses by a kicker and quadrupole septum and then transported for several meters to a high Z material target for conversion to x-rays for radiography. The un-kicked beam is diverted to a graphite dump. The interaction of this beam with the dump produces ions that propagate back to the target beamline and affect the beam properties. Coupling of the electron beam and ions to the beam dump vacuum chamber can also produce significant amounts of radiofrequency fields (RF) that disrupt the beam transport to the target. This requires a change in the nominal tune in order to minimize the heating of the carbon beam dump and suppress desorption of ions. Initially, ions are generated through stimulated desorption and ionization of the monolayers of water vapor stuck to the beam dump surface. These ions propagate upstream at high speeds and act as an electrostatic focusing lens by reducing the beam space charge. As the beam size on the dump decreases, the dump is sufficiently heated so that a surface plasma is formed and much greater quantities of ions are generated through temperature limited thermionic emission. The interaction of these ions with the beam creates RF which causes significant disruption of the beam. The results of a 3D PIC simulation show that ions can propagate back to the target beam line in about 500 ns. A beam position monitor (BPM) located just upstream of the septum was modified to include four capacitive pickups in addition to four of the usual inductive pickups. The inductive and capacitive pickups provide a measure of the beam current. The capacitive pickups are very sensitive to low velocity ions. As ions from the beam dump pass through the BPM the inductive and capacitive pickup will record two different measures of the beam current. Initial results from this BPM will be discussed along with other measurements which indicate the presence of backstreaming ions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Plasma Photonics for the Control of Kilojoule Lasers

The main goal of this project was to develop plasma photonic devices (plasma lenses or ellipsoidal plasma mirrors (EPMs)) for controlling the focus of high-energy (kJ-class) laser systems. In order to defocus a kJclass laser beam, a plasma optic must produce a very specific plasma density profile. For a plasma lens, this profile must be peaked on axis and fall off parabolicly with radius. We developed four plasma-lens target designs that were predicted to yield the ideal plasma profile required for defocusing kJ-class laser beams. Particle-in-cell (PIC) simulations using the code OSIRIS were used to explore the design. The resulting plasma lens must also be underdense, which sets an upper limit on the solid density of the plasma-lens targets to ~ 1 mg/cm 3 or less. This low of a density is highly challenging to fabricate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic Kα emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity (~10 20 W cm -2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement (1.6×) of Kα emission from microstructured targets compared to flat foils, but with a similar intensity and profile of Heα and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the Kα enhancement. The reasons for the lower-than-expected enhancement of Kα emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of Kα emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of Kα emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

New model for the ion collection by cylindrical probes over a wide range of collisionality

Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2021

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
NASA NTRS1991

Computer simulation of astrophysical plasmas

The role of sophisticated numerical models and simulations in the field of plasma astrophysics is discussed. The need for an iteration between microphysics and macrophysics in order for astrophysical plasma physics to produce quantitative results that can be related to astronomical data is stressed. A discussion on computational requirements for simulations of astrophysical plasmas contrasts microscopic plasma simulations with macroscopic system models. An overview of particle-in-cell simulations (PICS) is given and two examples of PICS of astrophysical plasma are discussed including particle acceleration by collisionless shocks in relativistic plasmas and magnetic field reconnection in astrophysical plasmas.

Max, Claire E.↗
NASA NTRS1991

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1990

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗