Search NASA⌕ Search

SEARCH · Search NASA

Results for “PICS”

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

Characterization of sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements with 3D hybrid particle-in-cell code

Hard x-rays produced by intense laser-produced fast electrons interacting with solids are a vital source for producing radiographs of high-density objects and implosion cores for inertial confinement fusion. Accurate calculation of hard x-ray sources requires a three-dimensional (3D) simulation geometry that fully models the electron transport dynamics, including electron recirculation and the generation of absolute photon yields. To date, 3D simulations of laser-produced bremsstrahlung photons over tens of picoseconds and code benchmarking have not been performed definitively. In this study, we characterize sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements for refluxing and non-refluxing targets using the 3D hybrid particle-in-cell (PIC), Large Scale Plasma code. Bremsstrahlung radiation and escaped electron data were obtained by focusing a 50-TW Leopard laser (15 J, 0.35 ps, 2 × 10 19 W/cm 2 ) on a 100-μm-thick Cu foil and a Cu with a large plastic backing (Cu–CH target). Data for both the Cu and Cu–CH targets were reproduced for simulations with a given set of electron parameters. Comparison of the simulations revealed that the hard x-ray emission from the Cu target was significantly longer in duration than that from the Cu–CH target. Here, the benchmarked hybrid PIC code could prove to be a powerful tool in the design and optimization of time- and angular-dependent bremsstrahlung sources for flash x-ray and gamma-ray radiography.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

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

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗
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 OSTI2024

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

We present two-dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein–Greene–Kruskal modes, which are exact nonlinear steady-state solutions of the Vlasov–Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 20482 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Physics↗
DOE OSTI2024

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (V b ), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (V d ) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (V s ) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterizing W sources in the all-W wall, all-RF WEST tokamak environment

In this work, experimental data, together with interpretive modeling tools, are examined to study trends in the tungsten (W) source in the all-W environment of the WEST tokamak, both from the divertor and from the main chamber. In particular, a poloidal limiter protecting an ion cyclotron resonance heating (ICRH) antenna is used as proxy for main chamber sourcing. The key study is carried out by stepping up lower hybrid current drive (LHCD) power, as the only auxiliary power source. Limiter and divertor W sources exhibit a qualitatively similar proportionality to the total power crossing the separatrix, P SEP , although the main chamber source remains substantially lower than the divertor source, for the range of P SEP accessible in the experiments. Intepretive modeling of the limiter source is carried out with a particle-in-cell (PIC) sheath model coupled to a surface sputtering model. Oxygen is used as a proxy for all light impurity species allowing for characterization of the critical W erosion regions. To get a good quantitative match to the data, it is necessary to assume that the oxygen arrives at the surface mostly at high ionization stages (4+ and above). A separate simulation with SOLEDGE-EIRENE, constrained to measured upstream scrape-off-layer plasma profiles, gives oxygen fractional abundances that are compatible with the PIC simulation result. This is understood to arise from transport processes that dominate over recombination. Substituting the LHCD by ICRH, in an equivalent experiment, the local W source exhibits a 3× enhancement. This can be matched by the simulation, by assuming local RF electric field rectification, based on ~100 eV peak-to-peak, near-antennna electric field. This work has highlighted the particular importance of understanding the ion charge state balance of light impurities as these are most likely the dominant sputtering species in fusion devices with high-Z walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

2D radial-azimuthal particle-in-cell benchmark for E × B discharges

In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Parametric studies of stream instability-induced higher harmonics in plasma ionization breakdown near an emissive dielectric surface

Here, in this work, we comprehensively investigate the generation of higher harmonic (HH) electric fields normal to the applied rf electric field in multipactor-coexisting plasma breakdown by fully kinetic particle-in-cell (PIC) simulations and a theoretical model. Firstly, a base case at driving frequency 1 GHz, transverse rf electric field amplitude 3 MV m –1 , and background gas pressure 0.2 Torr, is studied in detail. The enhanced harmonic frequency observed is around ten times the fundamental rf frequency, significantly lower than the Langmuir frequency. A theoretical model reveals that the fundamental mechanism of HHs generation is stream–plasma instability, which originates from stream-like secondary electron emission interacting with plasma. The resulting HH frequency and the growth rate of its oscillating amplitude from the theoretical model, agree well with the PIC simulations. With increasing pressure, the HH oscillations are found to be significantly reduced. This is because at higher pressure the gas ionization rate is higher, which causes a more rapidly increasing plasma density, leaving less time for the growth of instability. Furthermore, the parameter space in terms of background gas pressure and rf field amplitude within which the HHs appear is revealed. Finally, the effect of the driving rf frequency on HHs is also investigated, and it is found that the instability-induced oscillating HHs field remains when the driving frequency is increased to 2 GHz, however, it is highly reduced at higher driving frequency of 5 GHz, as oscillations at the fundamental frequency start playing a more important role.

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

Fast electron transport dynamics and energy deposition in magnetized, imploded cylindrical plasma

Inertial confinement fusion approaches involve the creation of high-energy-density states through compression. High gain scenarios may be enabled by the beneficial heating from fast electrons produced with an intense laser and by energy containment with a high-strength magnetic field. Here, we report experimental measurements from a configuration integrating a magnetized, imploded cylindrical plasma and intense laser-driven electrons as well as multi-stage simulations that show fast electrons transport pathways at different times during the implosion and quantify their energy deposition contribution. The experiment consisted of a CH foam cylinder, inside an external coaxial magnetic field of 5 T, that was imploded using 36 OMEGA laser beams. Two-dimensional (2D) hydrodynamic modelling predicts the CH density reaches 9.0 g cm − 3 , the temperature reaches 920 eV and the external B-field is amplified at maximum compression to 580 T. At pre-determined times during the compression, the intense OMEGA EP laser irradiated one end of the cylinder to accelerate relativistic electrons into the dense imploded plasma providing additional heating. The relativistic electron beam generation was simulated using a 2D particle-in-cell (PIC) code. Finally, three-dimensional hybrid-PIC simulations calculated the electron propagation and energy deposition inside the target and revealed the roles the compressed and self-generated B-fields play in transport. During a time window before the maximum compression time, the self-generated B-field on the compression front confines the injected electrons inside the target, increasing the temperature through Joule heating. For a stronger B-field seed of 20 T, the electrons are predicted to be guided into the compressed target and provide additional collisional heating. This article is part of a discussion meeting issue ‘Prospects for high gain inertial fusion energy (part 2)’.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electron acceleration at oblique angles via stimulated Raman scattering at laser irradiance >10 16 W cm –2 μm 2

The generation of hot, directional electrons via laser-driven stimulated Raman scattering (SRS) is a topic of great importance in inertial confinement fusion (ICF) schemes. Little recent research has been dedicated to this process at high laser intensity, in which back, side, and forward scatter simultaneously occur in high energy density plasmas, of relevance to, for example, shock ignition ICF. In this paper, we present an experimental and particle-in-cell (PIC) investigation of hot electron production from SRS in the forward and near-forward directions from a single speckle laser of wavelength λ 0 = 1.053μm, peak laser intensities in the range I 0 = 0.2–1.0×10 17 Wcm –2 and target electron densities between ne = 0.3–1.6% n c , where n c is the plasma critical density. As the intensity and density are increased, the hot electron spectrum changes from a sharp cutoff to an extended spectrum with a slope temperature T = 34 ± 1 keV and maximum measured energy of 350 keV experimentally. Multidimensional PIC simulations indicate that the high energy electrons are primarily generated from SRS-driven electron plasma wave phase fronts with k vectors angled ~50° with respect to the laser axis. These results are consistent with analytical arguments that the spatial gain is maximized at an angle which balances the tendency for the growth rate to be larger for larger scattered light wave angles until the kinetic damping of the plasma wave becomes important. The efficiency of generated high energy electrons drops significantly with a reduction in either laser intensity or target electron density, which is a result of the rapid drop in growth rate of Raman scattering at angles in the forward direction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Energy-preserving coupling of explicit particle-in-cell with Monte Carlo collisions

The particle-in-cell (PIC) and Monte Carlo collisions (MCC) methods are workhorses of many numerical simulations of physical systems. Recently, it was pointed out that, while the two methods can be exactly—or nearly—energy-conserving independently, combining the two leads to anomalous numerical heating. This paper reviews the standard explicit PIC-MCC algorithm, elucidates the origins of the anomalous numerical heating, and explains how to couple the two methods such that the anomalous numerical heating is avoided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗