Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasma modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Experimental methods to test photoionized plasma models in emission

The expanding foil photoionized plasma platform was developed on the Z-machine at Sandia National Laboratories to produce terrestrial photoionized plasmas at some of the same conditions as those observed in astrophysical black hole accretion disks. The specific scientific goals of the platform are motivated by long-standing astrophysical puzzles related to accreting black hole systems. Ultimately, the experimental goal is to inform the supersolar Fe abundance problem by testing photoionized plasma models in emission with at-parameter laboratory data. High quality high resolution emission data with absolute intensity calibration from a laboratory photoionized plasma have never been collected prior to the experiments described here, providing broader motivation. A robust test of the model predictions requires the successful production of a laboratory photoionized plasma along with measurements of the supporting model inputs including the radiation drive spectrum used to create the photoionized plasma, the electron density, and the temperature. The measured absolute spectral radiance in emission collected from the independently diagnosed plasma can be compared against model calculations performed using the measured inputs. In this paper, we focus on the experimental platform and the model inputs, specifically detailing the methodology used to measure the plasma conditions. The moderate electron temperature of ∼41 eV ± 15 eV and electron density of ∼1e19 ± 1.6e18 e−/cm3 and the ∼3e12 W/cm2 x-ray irradiance at the sample confirm that photoionization dominates the plasma ionization and spectra.

Cho, P. B. (ORCID:0000000291632905)↗
DOE OSTI2023

Prediction of Breakdown in Air and Solid Dielectrics: A Complete Plasma Model from Discharge Initiation to Flashover (Final Technical Report)

A major challenge for most of the technologies of today, and the technologies of the future, is the compact and reliable design of electrical insulation systems for their electrical components. Although controlled low-temperature plasmas (LTPs) or non-thermal plasmas in air at atmospheric pressure are of interest for a wide range of technologies, the appearance of such discharges is considered highly undesirable in insulation systems because they are precursors of a breakdown or flashover, where complete insulation failure takes place. For non-thermal plasmas, often only one mechanism, or the transition from one discharge mechanism to another, has typically been studied to date. On the other hand, in the case of thermal plasma applications, the conditions following the occurrence of an electric arc or flashover have always been considered. Here, there is a big technical gap: there is no paper describing a complete plasma model simulating all discharge mechanisms and their transitions from discharge initiation to electrical arc and flashover. To address this gap, a complete plasma model for both thermal and non-thermal plasmas should be developed using a hydrodynamic approach. The project aimed to build the model. The work has the potential to address industrial applications and is of great importance to higher voltage systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Formulation of a one-dimensional electrostatic plasma model for testing the validity of kinetic theory

Here, we present a one-dimensional (1-D) model composed of aligned, electrostatically interacting charged disks, conceived to address in a computable model the validity of the Bogoliubov assumption on the decay of particle correlations in the Born–Bogoliubov–Green–Kirkwood–Yvon hierarchy. This assumption is a basic premise of plasma kinetic theory. The disk model exhibits spatially 1-D features at short distances, but retains 3-D features at large distances. Here the collective dynamics of this model plasma is investigated by solving the corresponding Vlasov equation. In addition, the implementation of the model for the numerical validation of the Bogoliubov assumption is formulated.

1-D plasma model↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗
DOE OSTI2022

Modeling plasma response to non-axisymmetric magnetic field perturbations in tokamak boundaries

The primary goal of this project is to provide physics insight into models that elucidate the mechanisms responsible for controlling the performance of H-mode plasmas and edge localized mode (ELM) stability when subjected to intrinsic and applied non-axisymmetric magnetic perturbations. The overall goal is to facilitate the development of physics-based criteria for the design of magnetic coils to be used for ELM suppression and H-mode pedestal control in future burning plasma devices such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of MARS-F modeling of plasma response to RMPs using internal measurements on DIII-D

The linear resistive plasma response model is validated against the plasma internal measurement data from DIII-D edge-localized mode (ELM) control experiments with applied resonant magnetic perturbation (RMP). Considered are DIII-D discharges where the n = 1, 2, and 3 (n is the toroidal mode number) RMP field was applied. Experimental data for the plasma boundary displacement, as well as the three-dimensional (3D) pressure perturbation in the edge pedestal region, are deduced from the vertical Thomson scattering (TS) system and the horizontal charge exchange recombination (CER) system on DIII-D. The linear response model produces results that are in reasonable quantitative agreement with the DIII-D internal measurements. The plasma boundary displacement of up to 15 mm is modeled, with the pedestal pressure perturbation reaching 3 kPa. As an important insight, the larger plasma displacement measured by the vertical TS system, as compared to that measured by the horizontal CER system, is due to the contribution from the tangential component of the plasma displacement to the former. This mixing of displacement components is also found to influence the sensitivity of the CER measurement comparisons. The results of this study provide further confidence in the linear resistive plasma response model for analyzing ELM control experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗
DOE OSTI2026

Multi-machine validation of plasma initiation modelling and prospects for future devices: Predicting plasma initiation using only hardware design and control room input data

This paper reports on the generic prediction capability of full electromagnetic plasma initiation modelling with DYON, which was carried out for the first time in fusion research by the joint modelling of the International Tokamak Physics Activity—Integrating Operation Scenario group. The following devices were included in the experiment database: VEST (spherical torus, copper coils, Stainless steel wall, R/a = 0.3 m/0.2 m, V v = 3.7 m 3 ), MAST-U (spherical torus, copper coils, C wall, R/a = 0.7 m/0.5 m, V v = 55 m 3 ), EAST (conventional tokamak, superconducting coils, metallic wall, R/a = 1.85 m/0.5 m, V v = 38 m 3 ), DIII-D (conventional tokamak, copper coils, C wall, R/a = 1.67 m/0.65 m, V v = 35 m 3 ), and KSTAR (conventional tokamak, superconducting coils, C wall, R/a = 1.8 m/0.5 m, V v = 55 m 3 ). Despite the different hardware features of the devices, the required operating spaces of the loop voltage induction and prefill gas pressure for inductive plasma initiation in each device were successfully reproduced by the predictive simulations with DYON using only the individual hardware design and the control room input data for each discharge. This successful validation across multiple machines demonstrates that the full electromagnetic DYON modelling can capture the essential physics of inductive plasma initiation. The simulation settings commonly employed for all modelling and the modifications necessary to account for the discrepancies between individual devices are reported. Predictions for ITER based on the multi-machine validation indicate that a wide range of prefill gas pressures exists for the Townsend breakdown and the plasma burn-through (0.01–1.5 mPa).

DYON↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗