Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma 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 109 records · Page 6

Simulations of plasmas and fluids using anti-symmetric models

ALMA (Anti-symmetric, Large-Moment, Accelerated) is a fast, flexible, and scalable toolkit designed to solve hyperbolic conservation law systems in hybrid supercomputers. Here this manuscript describes the theoretical background and implementation of ALMA, which uses the anti-symmetric formulation of fluids to obtain simple, robust, and easily paralellizable code. Practical GPU acceleration is realized on entire applications with an overall gain factor of 2 to 4. ALMA also provides a parallel, GPU accelerated sparse solver based on geometric multigrid, capable of diagonalizing linear systems with 239 unknowns. Here we demonstrate ALMA's scaling and performance in petascale supercomputers and use standard fluid models to verify the overall approach with canonical benchmark problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MRT 7365 Power flow physics and key physics phenomena: EMPIRE verification suite

This milestone work baselines electromagnetic particle-in-cell capability of the EMPIRE plasma simulation code to model key processes germane to the physics of electrode plasmas arising in magnetically-insulated transmission lines operating at or near 20 MA. This evaluation is done so through the provision of benchmark verification problems designed to exercise the individual and combined physics models on a small-scale surrogate geometry for the final-feed-to-load region of the Z accelerator under representative operating conditions. In this report, we overview our test designs, and present a portfolio of simulation results along with performance assessments which altogether establish state-of-the-art. In particular, two main verification categories are covered this report: (1) Z-relevant desorption physics (Temkin isotherm), and (2) two approaches to simulate electrode plasma creation and dynamics (automatic creation versus self-consistent creation through direct simulation Monte Carlo collisions).

43 PARTICLE ACCELERATORS↗

Achieving Higher Fidelity Pulsed-Power Simulations Using Advanced Gap Closure Modeling

This project has progressed several physics models in the EMPIRE plasma simulation code to achieve higher fidelity simulations of high-current diode operation in pulsed-power accelerators. In this report, we present details for the following major work products: (1) a set of verification problems covering all key processes involved in gap closure physics has been designed; this suite has facilitated feature vetting and overall model maturation, (2) a new EMPIRE exemplar has been developed: the Radiographic Integrated Test Stand 6 (RITS-6) diode, and (3) An exemplar for the Saturn accelerator exemplar was enabled by the models matured under this work to self-consistently simulate further into the diode pulse than previously possible (bipolar flow regime). These developments have lead to the highest confidence EMPIRE power flow predictions of the Saturn accelerator to date. Additionally, three modeling approaches for simulating electrode plasmas have been investigated. We report on these results and provide recommendations.

43 PARTICLE ACCELERATORS↗

Interpretive modelling of boron transport in the boundary plasma of WEST experiments with the impurity powder dropper

Abstract Boron (B) powder injection is a potential alternative to glow discharge boronization as a wall conditioning method for tokamaks. This technique is currently being studied in WEST experiments, during which B powder is injected by an Impurity Powder Dropper developed by PPPL. In order to interpret and analyse experimental trends, and to help develop future experiments, a modelling workflow using a boundary plasma simulation (SOLEDGE-EIRENE) and powder ablation simulation (Dust Injection Simulator) was developed and tested. The effect of adding a B neutral source to simulated deuterium + oxygen (D + O) plasmas was compared to experimental data from the WEST C5 campaign, where B powder was injected in a dedicated experiment. While the impact of B injection on radiated power P rad measurements at the upper divertor was similar, there were significant differences in measurements of P rad , outer strike point electron temperature T e OSP and O-II line intensity at the lower divertor between experiment and simulation. This discrepancy suggests that those parameters were affected by phenomena not present in the simulations, with the most likely candidates being reduced D recycling and a reduced O sourcing from the divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of spatiotemporal smoothing on the two-plasmon–decay instability

Higher levels of hot electrons from the two-plasmon–decay instability are observed when smoothing by spectral dispersion (SSD) is turned off in directly driven inertial confinement fusion experiments at the Omega Laser Facility. This finding is explained using a hot-spot model based on speckle statistics and simulation results from the laser–plasma simulation environment. The model accurately reproduces the relative increase in hot-electron activity at two different drive intensities, although it slightly overestimates the absolute number of hot electrons in all cases. Extrapolating from the current ≈ 360-GHz system while adhering to the logic of the hot-spot model suggests that larger SSD bandwidth should significantly mitigate hot-electron generation, and legacy 1-THz OMEGA experiments appear to support this conclusion. These results demonstrate that it is essential to account for laser speckles and spatiotemporal smoothing to obtain quantitative agreement with experiments. A compilation of hot-electron data from the past two decades reveals several other important points: (1) many prior experiments are more easily understood using recent results from multibeam absolute instability theory; and (2) experiments with ignition-scale conditions produce less hot electrons compared to OMEGA spherical experiments for a given vacuum overlapped intensity, which is a promising result for validating performance predictions based on hydrodynamic scaling relations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ionization dynamics of intense laser-produced argon plasmas revealed by NLTE modeling

The ionization dynamics and transient behavior of under-dense plasma irradiated by an intense laser are investigated. We report two significant effects in the ionization behavior: (1) a surprisingly large delay in ionization response and (2) a stepwise ionization process which involves collisional and laser-driven photoionization (LDP) processes. Ionization induced by intense lasers can exhibit delayed responses due to rapid changes in conditions, particularly when atomic transition processes occur more slowly than the relevant time scales. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by LDP—plays an important role in this ionization delay, with collisional excitation acting as the bottleneck. Even low-energy photons (∼3.5 eV) can predominantly ionize plasmas, challenging the conventional belief that such energies are insufficient to overcome the binding energy of bound electrons. These findings underscore the necessity of including such processes into plasma simulations for various laser-plasma experiments.

Collisional excitation↗

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hermes-3

Hermes-3 is a plasma simulation model built on BOUT++, developed as part of research studies of the edge and divertor regions of magnetically confined fusion plasmas such as tokamaks. The main purpose of this model is multi-species simulation of plasmas containing a mixture of deuterium, tritium, helium and other plasma and neutral gas species.

Dudson, BenjaminD↗

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗