Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma models”

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 559 records · Page 31

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗

Fast ion transport by sawtooth instability in the presence of ICRF-NBI synergy in JET plasmas

JET experiments have shown that the three-ion scenarios using waves in the ion cyclotron range of frequencies (ICRF) is an efficient way to build fast ion population through beam ion acceleration by radio frequency (RF) waves. Here, such a heating scheme is applied to plasmas with at least two thermal ion species. Analysis of mixed discharges with complex heating schemes requires a workflow that allows to model thermal and fast ion transport consistently. This paper is dedicated to modelling of a mixed plasma discharge with significant fraction of fast ions and contributes to development of fast ion transport models. For interpretive analysis with the TRANSP code a JET hydrogen-deuterium (HD) plasma discharge with neutral beam injection (NBI) and ICRF heating has been chosen. The task is complicated by NBI-ICRF synergy and plasma magnetohydrodynamic activity, like sawtooth crashes. D beam ions accelerated by RF waves form a high energy tail in fast ion distribution. Significant difference between the neutron rate computed by TRANSP and measured one is observed if the same diffusivity for electrons and ions is assumed. Sensitivity studies show that uncertainties in input plasma parameters and thermal ion transport models are crucial for modelling mixed plasma discharges and increased D transport is required to reach the plasma composition consistent with diagnostic measurements at the plasma edge. Fast ion redistribution by a sawtooth instability is characterised by non-resonant transport due to reconnection of magnetic field lines and resonant transport caused by resonance interaction between the instability and fast ions. With ORBIT simulations it has been shown that resonant interaction strongly affects fast ions of high energies, like beam ions accelerated by RF waves and fusion products. For the considered case, fast ion profiles simulated by ORBIT remain peaked after the sawtooth crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

CHEQUP v0.1

CHEQUP (Castro-based Hofi Expansion with QUasineutral Plasma) is a simulation code for modeling the formation of hydrodynamic optical-field-ionized (HOFI) plasma channels, which are used as waveguides in laser-plasma acceleration experiments. This includes experiments performed at LBNL's BELLA facility as well as other laser facilities across the world. CHEQUP extends the open-source Castro hydrodynamics framework with physics modules tailored for modeling HOFI plasma channels -- including multi-species ionization and three-body recombination for mixtures of hydrogen, nitrogen, helium, and argon ; a two-temperature model tracking electron and heavy-species temperatures separately ; and coupling with other codes of the BLAST ecosystem (https://blast.lbl.gov/) such as WarpX, via the openPMD standard. CHEQUP inherits from Castro the ability to run on modern GPU architectures (NVIDIA CUDA, AMD HIP) and supports adaptive mesh refinement (AMR) for efficient multi-scale resolution. Compared to existing tools, CHEQUP would be, to our knowledge, the first open-source code implementing the full HOFI channel formation physics, and the first implementation capable of running on GPUs. This enables significantly faster, large-scale parameter scans critical for the design of next-generation LPA-based accelerators and light sources.

Lehe, Remi [Lawrence Berkeley National Laboratory ↗

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dataset, Code, and Models for Training Deep Learning Potentials for Low Temperature Plasma-Surface Interactions

This repository contains datasets, training scripts, and finished models, and test simulations used in the development of DeepREBO— a machine-learned interatomic potential trained to emulate the REBO2 empirical potential. The data was generated to study deep potential development for simulations of plasma-surface interactions. It uses an active learning framework, starting from a minimal dataset and iteratively expanding it. Included are those generated datasets, the trained models, and simulations used to evaluate the performance of the training process. This resource supports reproducibility and provides a reference framework for training deep potentials in plasma-surface interaction studies.

active learning↗

Real-time capable modeling of ICRF heating on NSTX and WEST via machine learning approaches

Abstract A real-time capable core Ion Cyclotron Range of Frequencies (ICRF) heating model on NSTX and WEST is developed. The model is based on two nonlinear regression algorithms, the random forest ensemble of decision trees and the multilayer perceptron neural network. The algorithms are trained on TORIC ICRF spectrum solver simulations of the expected flat-top operation scenarios in NSTX and WEST assuming Maxwellian plasmas. The surrogate models are shown to successfully capture the multi-species core ICRF power absorption predicted by the original model for the high harmonic fast wave and the ion cyclotron minority heating schemes while reducing the computational time by six orders of magnitude. Although these models can be expanded, the achieved regression scoring, computational efficiency and increased model robustness suggest these strategies can be implemented into integrated modeling frameworks for real-time control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic mesh generation and field line reconstruction for scrape-off layer and divertor modeling in stellarators

The design of divertor targets and baffles for optimal heat and particle exhaust from magnetically confined fusion plasmas requires a combination of fast, low-fidelity models (such as EMC3-Lite (Feng 2022 Plasma Phys. Control. Fusion 64 125012)) for scoping studies and high-fidelity ones (such as EMC3-EIRENE (Feng et al 2014 Contrib. Plasma Phys. 54 426–31)) for verification. Both of those approaches benefit from a magnetic flux tube mesh for fast interpolation and mapping of field line segments (Feng et al 2005 Phys. Plasmas 12 052505). A new automated mesh generator for unstructured quadrilateral flux tubes with adaptive refinement is presented and integrated into FLARE (Frerichs 2024 Nucl. Fusion 64 106034). For HSX with an extended first wall, it is found that several layers of flux tubes can span the entire half field period before splitting is required. This is an advantage over the traditional setup of the EMC3-EIRENE mesh where careful construction of several sub- domains is required already for the much tighter present first wall. In particular, there is no longer the need to manually construct a suitable outer boundary for the mesh. The divide and conquer paradigm with unstructured mesh layout offers a powerful alternative for fast head load approximation that is suitable for integration into optimization workflows. Further examples for W7-X and CTH demonstrate the versatile application range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Reference Governor for Plasma-Shape Control in Tokamaks

Tokamaks are torus-shaped devices where a high-temperature, ionized gas (i.e. a plasma) is confined by means of helical magnetic fields. The ultimate goal of tokamaks is producing net electrical energy by means of nuclear-fusion reactions within the plasma. The shape of the plasma inside the tokamak is closely related with its confinement characteristics and, therefore, with the produced fusion power. In order to achieve and maintain plasma shapes that maximize the fusion yield, feedback controllers have been successfully developed during the years. However, machine-protection and plasma-stability requirements impose hard constraints on the plasma-shape targets that can be safely achieved by a feedback controller. In this work, a reference governor is proposed to provide, in real time, plasma-shape targets that fulfill the aforementioned safety requirements. The reference governor calculates such plasma-shape targets in response to changes in the state of the closed-loop system, as well as disturbances in the plasma energy and current. Here, the reference-governor design is based on a dynamical model of the plasma electromagnetic response. Simulation tests of the reference governor have been carried out for different plasma scenarios in the DIII-D tokamak to emulate realistic cases where safety may be compromised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗

The impact of non-local parallel electron transport on plasma-impurity reaction rates in tokamak scrape-off layer plasmas

Abstract Plasma-impurity reaction rates are a crucial part of modelling tokamak scrape-off layer (SOL) plasmas. To avoid calculating the full set of rates for the large number of important processes involved, a set of effective rates are typically derived which assume Maxwellian electrons. However, non-local parallel electron transport may result in non-Maxwellian electrons, particularly close to divertor targets. Here, the validity of using Maxwellian-averaged rates in this context is investigated by computing the full set of rate equations for a fixed plasma background from kinetic and fluid SOL simulations. We consider the effect of the electron distribution as well as the impact of the electron transport model on plasma profiles. Results are presented for lithium, beryllium, carbon, nitrogen, neon and argon. It is found that electron distributions with enhanced high-energy tails can result in significant modifications to the ionisation balance and radiative power loss rates from excitation, on the order of 50%–75% for the latter. Fluid electron models with Spitzer-Härm or flux-limited Spitzer-Härm thermal conductivity, combined with Maxwellian electrons for rate calculations, can increase or decrease this error, depending on the impurity species and plasma conditions. Based on these results, we also discuss some approaches to experimentally observing non-local electron transport in SOL plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of an implicit electromagnetic capability for a hybrid gyrokinetic ion-fluid electron model

Here we report on the development and implementation of a hybrid kinetic ion–fluid electron model for electromagnetic COGENT simulations of edge plasmas. COGENT is a finite-volume gyrokinetic code that employs a locally field-aligned coordinate system combined with a mapped multi-block grid technology to handle strongly anisotropic edge plasma turbulence. The simulation model involves the long-wavelength limit of the ion gyrokinetic equation coupled to the vorticity and Ohm's law equations for the electromagnetic field perturbations. In order to handle the fast Alfvén wave time scales, an implicit-explicit time integration approach with a physics-based preconditioner is used. The model is successfully applied to the simulations of ion-scale resistive-drift ballooning turbulence in a toroidal annulus geometry. Substantial speed-up over a fully explicit time integration approach is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of non-equilibrium ionization and excitation on radiation absorption in plasma plumes induced by ablation of metal targets with nanosecond laser pulses

Ionization and radiation absorption in nanosecond laser-induced plasma plumes are routinely modeled using the Saha–Boltzmann equilibrium ionization model (EQM). However, the equilibrium assumption can be inaccurate during the laser pulse when non-equilibrium effects significantly impact radiation absorption. In the present work, the EQM and non-equilibrium collisional-radiative model (CRM) are compared to reveal the effect of plasma non-equilibrium on radiation absorption in non-homogeneous plumes and degree of plasma shielding. Simulations of plume expansion induced by irradiation of a copper target in 1 atm argon background gas with a 10 ns Gaussian pulse at a fluence from 8 Jcm −2 to 14 Jcm −2 are performed with a hybrid computational model that couples a lumped particle direct simulation Monte Carlo method with either CRM or EQM. The simulations show that the EQM strongly underestimates the effects of ionization and radiation absorption compared to CRM and, contrary to the CRM, predict strong ionization of the background gas. The differences between the models are explained by the qualitatively different coupling between plume expansion and dynamics of ionization and exitation processes in the CRM and EQM under conditions when the characteristic times for most radiation- and electron-induced processes are longer than the pulse duration. The CRM-based predictions are also found to agree much better with available experimental data. In conclusion, these results indicate that the model of Saha–Boltzmann equilibrium cannot be used for reliable prediction of the degree of plasma shielding in plumes induced by nanosecond laser pulses or for processing results of spectroscopic measurements at early stages of expansion of such plumes.

Physics - Plasma physics↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗