Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ion density waves driving the formation of filamentary dust structures

The PlasmaKristall-4 (PK-4) experiment on the International Space Station allows for the study of the 3-dimensional interaction between plasma and dust particles. Previous simulations of the PK-4 environment have discovered fast moving ionization waves in the dc discharge. These ionization waves vary the plasma parameters by up to an order of magnitude, which may affect the mechanisms responsible for the self-organization of chains seen in the PK-4 experiment. Here, we adapt a molecular dynamics simulation to employ temporally varying plasma conditions in order to investigate the effect on the dust charging and electrostatic potential. In order to describe the differences between the average of the plasma conditions and the time-varying plasma condition, we present a model to reproduce the interaction that takes into account the negative potential from the dust grain and the positive potential from the ion wake.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling and verification of dynamic field ionization for laser-target interactions

Integrating field ionization models into kinetic plasma simulations is required for a variety of applications, especially when field strengths vary from low to high regimes, such as in laser-target interactions. The introduction of new physics models into kinetic codes requires a rigorous verification of their accuracy through well-defined verification problems. In this work, the field ionization model that has been included in the kinetic plasma code VPIC is presented, along with the detailed approach adopted for its integration. This model includes a comprehensive range of field ionization mechanisms: multiphoton ionization, tunneling ionization, and barrier suppression ionization. New verification problems employed to evaluate the ionization model's fidelity are outlined, and the simulation parameters that affect the accuracy of simulation results are explored. Additionally, this work addresses the impact of field ionization on computational performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Monte-Carlo method for Coulomb collisions in moderately coupled plasmas

Here, we introduce a new approach for coupling cumulative and single-scattering distributions, extending Monte Carlo methods for Coulomb collisions in weakly coupled plasmas into the moderately coupled regime. Using time steps typical of cumulative-scattering methods, we show that this approach not only recovers leading-order corrections proportional to the inverse Coulomb logarithm but also solves the full Boltzmann collision integral. A key feature of the method is its simplicity: it can be integrated into cumulative-scattering methods for weakly coupled plasmas with only ~ 10 additional lines of code.

Boltzmann collision integral↗
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 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 OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗
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 OSTI2026

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗
DOE OSTI2025

Measuring the growth of Alfvén wave parametric decay instability using counter-propagating waves: Theory and simulations

The parametric decay instability (PDI) of Alfvén waves—where a pump Alfvén wave decays into a backward-propagating child Alfvén wave and a forward ion acoustic wave—is a fundamental nonlinear wave-wave interaction and holds significant implications for space and laboratory plasmas. However, to date there has been no direct experimental measurement of PDI. Here, we propose a novel and experimentally viable scheme to quantify the growth of Alfvén wave PDI on a linear device using a large pump Alfvén wave and a small counter-propagating seed Alfvén wave, with the seed-wave frequency tuned to match the backward Alfvén wave generated by standard PDI. Using hybrid simulations, we show that energy transfer from the pump to the seed reduces the latter's spatial damping. By comparing seed-wave amplitudes with and without the pump wave, this damping reduction can be used as a direct and reliable proxy for PDI growth. The method is validated in our simulations across a range of plasma and wave parameters and agrees well with theoretical predictions. Notably, the scheme exhibits no threshold for PDI excitation and is, in principle, readily implementable under current laboratory conditions. Finally, this scheme is a critical step toward solving the challenge of experimentally accessing Alfvén wave PDI and provides an elegant method that may be used to validate fundamental theories of parametric instabilities in controlled laboratory settings.

heliosphere↗
DOE OSTI2026

Flying focus with arbitrary directionality for spatiotemporal control of laser intensity

Flying focus techniques produce laser pulses whose focal points travel at arbitrary, controllable velocities. While this flexibility can enhance a broad range of laser-based applications, existing techniques constrain the motion of the focal point to the propagation direction of the pulse. Here, we introduce a flying focus configuration that decouples the motion of the focus from the propagation direction. Here, a chirped laser pulse focused and diffracted by a diffractive lens and grating creates a focal point that can move both along and transverse to the propagation direction. The focal length of the lens, grating period, and chirp can be tuned to control the direction and velocity of the focus. Simulations demonstrate this control for a holographic configuration suited to high-power pulses, in which two off-axis pump beams with different focal lengths encode the equivalent phase of a chromatic lens and grating in a gas or plasma. For low-power pulses, conventional solid-state or adaptive optics can be used instead. Multi-dimensional control over the focal trajectory enables new configurations for applications, including laser wakefield acceleration of ions, nonlinear Thomson scattering, and surface-plasmon emission of THz radiation.

Classical optics↗
DOE OSTI2025

Axion pulsarscope

Electromagnetic fields surrounding pulsars may source coherent ultralight axion signals at the known rotational frequencies of the neutron stars, which can be detected by laboratory experiments (e.g., pulsarscopes). As a promising case study, we model axion emission from the well-studied Crab pulsar, which would yield a prominent signal at 𝑓 ≈ 29.6 Hz regardless of whether the axion contributes to the dark matter abundance. We estimate the relevant sensitivity of future axion dark matter detection experiments such as DMRadio-GUT, Dark SRF, and CASPEr, assuming different magnetosphere models to bracket the uncertainty in astrophysical modeling. For example, depending on final experimental parameters, the Dark SRF experiment could probe axions with any mass 𝑚 𝑎 ≪ 10 −13 eV down to 𝑔 𝑎⁢𝛾⁢𝛾 ∼3 × 10 −13 GeV −1 with one year of data and assuming the vacuum magnetosphere model. These projected sensitivities may be degraded depending on the extent to which the magnetosphere is screened by charge-filled plasma. The promise of pulsar-sourced axions as a clean target for direct detection experiments motivates dedicated simulations of axion production in pulsar magnetospheres.

axions↗
DOE OSTI2025

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Single-Shot Reconstruction of Electron Beam Longitudinal Phase Space in a Laser Wakefield Accelerator

We report on a single-shot longitudinal phase-space reconstruction diagnostic for electron beams in a laser wakefield accelerator via the experimental observation of distinct periodic modulations in the angularly resolved spectra. Such modulated angular spectra arise as a result of the direct interaction between the ultrarelativistic electron beam and the laser driver in the presence of the wakefield. A constrained theoretical model for the coupled oscillator, assisted by a genetic algorithm, can recreate the experimental electron spectra and, thus, fully reconstructs the longitudinal phase-space distribution of the electron beam with a temporal resolution of approximately 1.3 fs. In particular, it reveals the slice energy spread of the electron beam, which is important to measure for applications such as x-ray free electron lasers. In our experiment, the root-mean-square energy spread retrieved is bounded at 9.9 MeV, corresponding to a 0.9%–3.0% relative spread, despite the overall GeV energy beam having approximately 100% relative energy spread.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces

This project focused on the development and integration of high-performance simulation tools to predict the operating behavior of Plasma-Facing Components in magnetic confinement fusion systems. A key objective at Illinois was to assess the impact of the dynamic interplay between the evolving material surface and the magnetized plasma sheath, and characterize the impact of tungsten-based PFCs on plasma contamination, including phenomena such as surface erosion, dynamic recycling of fuel species, and tritium retention, which are critical for the success of future magnetic fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗