Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle methods”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
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

Core-edge modeling of gyrokinetic turbulence by coupling the delta-f and total-f models in the XGC code

Whole device gyrokinetic turbulence is modeled using a delta-f model in the core coupled to a total-f model in the edge. This new model is implemented in the gyrokinetic code X-Point Included Gyrokinetic Code (XGC) and works within a single simulation. It employs a canonical Maxwellian background in the core and the usual local Maxwellian in the edge. These different core and edge backgrounds are coupled over an overlap region. A thermal-bath operator is used to keep δf small in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic study of shock formation and particle acceleration in laser-driven quasi-parallel magnetized collisionless shocks

Quasi-parallel magnetized collisionless shocks are believed to be one of the most efficient accelerators in the universe. Compared to quasi-perpendicular shocks, quasi-parallel shocks are more difficult to form in the laboratory and to simulate because of their large spatial scales and long formation times. Our two-dimensional particle-in-cell simulations show that the early stages of quasi-parallel shock formation are achievable in experiments planned for the National Ignition Facility and that particles accelerated by diffusive shock acceleration (DSA) are expected to be observable in the experiment. Repetitive ion acceleration by crossings of the shock front, a key feature of DSA, is seen in the simulations. Other characteristic features of quasi-parallel shocks such as upstream wave excitation by energetic ions are also observed, and energy partition between the ions and the electrons in the downstream of the shock is briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient backward x-ray emission in a finite-length plasma irradiated by a laser pulse of picosecond duration

Motivated by experiments employing picosecond-long, kilojoule laser pulses, we examined x-ray emission in a finite-length underdense plasma irradiated by such a pulse using two-dimensional particle-in-cell simulations. We found that, in addition to the expected forward emission, the plasma also efficiently emits in the backward direction. Our simulations reveal that the backward emission occurs when the laser exits the plasma. The longitudinal plasma electric field generated by the laser at the density down-ramp turns around some of the laser-accelerated electrons and re-accelerates them in the backward direction. As the electrons collide with the laser, they emit hard x rays. The energy conversion efficiency is comparable to that for the forward emission, but the effective source size is smaller. We show that the picosecond laser duration is required for achieving a spatial overlap between the laser and the backward energetic electrons. At peak laser intensity of 1.4×1020 W/cm2, backward-emitted photons (energies above 100 keV and 10° divergence angle) account for 2×10−5 of the incident laser energy. This conversion efficiency is three times higher than that for similarly selected forward-emitted photons. The source size of the backward photons (5 μm) is three times smaller than the source size of the forward photons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

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

Plasma gradient effect on direct laser acceleration

The transfer of a high-intensity laser pulse energy to a high-energy electron beam via the direct laser acceleration mechanism is shown to be significantly enhanced through control of the plasma density gradient. Experiments performed using the OMEGA EP facility's high-intensity beams altered the plasma density and gradients by changing the Mach number and the angle of the gas-jet nozzle to the laser axis. When a long density gradient at the rear of the target is used, the total high-energy electron number measured was enhanced by 4.5 times compared to a shorter rear gradient. Complementary two-dimensional simulations, which follow the laser field evolution and the corresponding electron dynamics, strongly support the key trends observed in the experiment. The effect is twofold, the long density gradient provides the longest acceleration distance while it minimizes the formation of the sheath field as the electron beam exits into the vacuum. This study shows the importance of tailoring the plasma density.

Laser plasma interactions↗
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 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 OSTI2025

Nonlocal, diamagnetic electromagnetic effects in magnetically insulated transmission lines

We identify the time-dependent physics responsible for the critical reduction of current losses in magnetically insulated transmission lines (MITLs) due to uninsulated space charge-limited currents of electrons emitted by field stress. A drive current of sufficiently short pulse length introduces a strong enough time dependence that steady-state results alone become inadequate for the complete understanding of current losses. The time-dependent physics can be described as a nonlocal, diamagnetic electromagnetic response of space charge limited currents. As the pulse length is increased or equivalently, the MITL length reduced, these time-dependent effects diminish and current losses converge to those predicted by the well-known Child–Langmuir law in the external (vacuum) fields. We present a simple one-dimensional (1D) model that encapsulates the essence of this physics. We find excellent agreement with 2D particle-in-cell simulations for two MITL geometries, Cartesian parallel plate and azimuthally symmetric straight coaxial. Based on the 1D model, we explore various scaling dependencies of MITL losses with relevant parameters, e.g., peak current, pulse length, geometrical dimensions, etc. We propose an improved physics model of magnetic insulation in the form of a Hull curve, which could also help improve predictions of current losses by common circuit element codes, such as BERTHA. Finally, we describe how to calculate the temperature rise due to electron impact within the 1D model.

Computer simulation↗
DOE OSTI2025

Increased electron, positron, and x-ray production from high intensity laser interactions using micro-wire targets

We report increases in energetic electrons, positrons, and x-rays emitted from high-intensity laser interactions (10 18−20 W/cm 2 ) with structured silicon micro-wires on the surface of a 1 mm gold converter target using a 10 ps laser pulse. A total of four different wire configurations are tested, where the gaps (7–28 μm) between the wires and the thicknesses (3–6 μm) of the wires are varied, while the height remains constant (⁠ ~25 μm). We observe the largest enhancement in electrons, positrons, and x-rays with the sparsest wire configurations. The electron temperature (T e ≈6 MeV) remains consistent across all shots, regardless of whether wires or planar targets are used. This suggests that the observed enhancement is due to increased laser light absorption by the accelerated electrons over a long scale length. Two-dimensional particle-in-cell simulations confirm that absorption is significantly enhanced with the wire target. Additionally, specific simulations examining laser pointing on different parts of the wire structure reveal that, while the final electron spectrum remains largely insensitive, the angular distribution is highly sensitive to these variations.

Bremsstrahlung↗
DOE OSTI2025

On the role of self-generated magnetic fields in focused ion beams driven by intense picosecond laser pulses

We investigate the physics of ion beam focusing in the interaction of intense picosecond laser pulses with concave foils, using two-dimensional particle-in-cell simulations. Curved targets have been observed to deliver brighter multi-MeV ion beams than flat targets, and they could be used to ignite pre-compressed ICF fusion capsules in a scheme called ion fast ignition. Focusing is achieved by shaping the target and relying on an ion acceleration mechanism in which the ions are accelerated normal to the target surface; subsequently, the plasma expansion grows magnetic fields that advect with the plasma and deflect ions away from the target symmetry axis, causing the focal spot to shift downstream from its geometric location. The drop in both electric and magnetic field strengths in the expanding plasma makes the ion beam envelope asymmetric with respect to the focusing plane.

Focused ion beam↗
DOE OSTI2025

Multimodal azimuthal oscillations in electron beam generated E×B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E×B) fields are promising for low-damage material processing. However, these plasmas can be subject to the formation of azimuthally propagating structures that enhance the radial transport of energetic charged species, which can harm the gentle processing capability of the plasma. In this work, we investigate the azimuthal structure formation in an e-beam generated E×B plasma using experimental diagnostics and 2D3V particle-in-cell simulations. Our findings demonstrate the formation of multiple simultaneously occurring azimuthally propagating modes that exhibit nontrivial radial dependence. It is suggested that the multimodal azimuthal spectrum is caused in part by the radial outflow of ions in the plasma.

electron beams↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗