Search NASA⌕ Search

SEARCH · Search NASA

Results for “Collisional Radiative 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

Physics-assisted Latent Space Dynamics Learning for Stiff Collisional-radiative Models

Collisional-radiative (CR) models describe the atomic processes in a plasma by tracking the population density in the ground and excited states for each charge state of the atom/ion. These models predict important plasma properties such as charge state distributions and radiative emissivity and opacity. Accurate descriptions of the CR balance of the plasma are essential in fusion whole device modeling, especially when significant impurities are introduced into the plasmas. In an integrated fusion plasma and CR simulation, a CR model, which is a high-dimensional stiff ODE, needs to be solved on each grid point in the configuration space, and can overwhelm the overall computational cost. In this work, we propose a machine-learning-based method that discovers a latent space and learns its corresponding latent dynamics, which can capture the essential physics to make accurate predictions at much lower online computational cost. The proposed approach is physics-assisted, due to its combination of a physical latent space and a data-driven latent space. It has been demonstrated that the proposed architecture can predict both the full-order CR dynamics and a physical quantity of interest accurately.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters: Experimental Validation and Model Investigations

The recently developed collisional radiative (CR) models for neutral xenon (Xe I) and singly charged xenon ion (Xe II) described by Chaplin et al. in the paper "Insights from Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters" also submitted to the 2021 AIAA Propulsion and Energy Forum are more thoroughly investigated. Comparisons between the existing and recently developed JPL Xe I CR models are done and reasons for discrepancies between the two models are explored. An experiment to validate the Xe II CR model is carried out in the anode region of the UCLA Plasma interactions facility. The Xe II CR model is applied to measured xenon ion emission line intensities to determine electron temperature and density. Spatial resolution of these optical measurements is obtained through an Abel inversion for direct comparison to Langmuir probe results. Although the measured electron densities were on the order of $10^{17}$ m$^{-3}$ and the Xe II CR model is most suitable for use with densities above $10^{18}$ m$^{-3}$ due to the strong dependence of the line ratios on electron temperature in that regime, initial results from the validation experiment show reasonable agreements between the optical and Langmuir probe results. The 441 nm to 605 nm line ratio appears to have the most promise for non-intrusive electron temperature measurements with the Xe II CR model. Shortcomings of the Xe II CR model revealed in the validation experiment are to be explored further.

Lobbia, Robert B.↗

Time-dependent collisional-radiative model for capillary discharge plasmas

A time-dependent, zero-dimensional collisional-radiative model has been developed to study capillary discharge plasmas in view of possibilities for obtaining population inversion leading to laser action in the extreme ultraviolet and soft x-ray spectral region, by amplifying the Balmer-alpha (H-alpha) line of different hydrogenlike ions with nuclear charges from Z = 3 to Z = 6. The model is described in detail, and results for the case of a carbon plasma are presented. Limitations of the model are discussed and comparison is made between our calculations and recent experimental results.

Pockl, M.↗

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗

Investigation of infra-red and nonequilibrium air radiation

This report describes progress on the first year of a research program on the infrared radiation of air plasmas conducted in the High Temperature Gasdynamics Laboratory at Stanford University. This program is intended to investigate the masking of infrared signatures by the air plasma formed behind the bow shock of high velocity missiles. To this date, the radiative emission of air plasmas in the infrared has been the object of few experimental investigations, and although several infrared systems are already modeled in radiation codes such as NEQAIR, measurements are required to validate numerical predictions and indicate whether all transitions of importance are accounted for. The present program is motivated by the fact that 9 excited states (A, B, C, D, B', F, H, and H') of NO radiate in the infrared, especially between 1 and 1.5 microns where at least 9 transitions involving can be observed. Because these IR transitions are relatively well separated from each other, excited NO states concentrations can be easily measured, thus providing essential information on excited-state chemistry for use in optical diagnostics or in electronic excitation model validation. Developing accurate collisional-radiative models for these excited NO states is of importance as the UV-VUV transitions of NO (beta, gamma, epsilon, beta prime, gamma prime) produce a major, if not dominant, fraction of the radiation emitted by air plasmas. During the first year of the program, research has focused on the spectral range 1.0 to 1.5 microns, as detailed in Section 2 of this report. The measurements, conducted in a 50 kW radio-frequency inductively coupled plasma torch operating on air at atmospheric pressure, extend previous shock tube investigations by Wray to a wider spectral range (1.0 to 1.5 microns vs 0.9 to 1.2 microns) and higher temperatures (7600 K in the plasma torch versus 6700 K in the shock-tube). These higher temperatures in the present experiment have made it possible to observe high-lying excited NO states that were previously undetectable. These measurements are currently being extended up to 5 microns, with particular attention paid to the rovibronic bands of ground state NO, molecular continua, CO transitions, and other systems of importance. Publications and presentations resulting from or related to this work are cited in Section 3, and Section 4 lists the personnel who contributed to this report.

Kruger, Charles H.↗

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

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↗

Laboratory astrophysics on ASDEX Upgrade: Measurements and analysis of K-shell O, F, and Ne spectra in the 9 - 20 A region

High-resolution measurements of K-shell emission from O, F, and Ne have been performed at the ASDEX Upgrade tokamak in Garching, Germany. Independently measured temperature and density profiles of the plasma provide a unique test bed for model validation. We present comparisons of measured spectra with calculations based on transport and collisional-radiative models and discuss the reliability of commonly used diagnostic line ratios.

Hansen, S. B.↗

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomic Data for Plasma Spectroscopy: The CHIANTI Database, Improvements and Challenges

CHIANTI is an atomic database and software package for modeling emission lines and continua from hot astrophysical plasmas. It is freely available to all researchers and has been widely used in the Heliophysics and Astrophysics communities for almost 25 years. In this review, we summarize the properties of the current version of the database and give an overview of the relevant atomic processes. We also discuss progress towards a complete implementation of collisional-radiative modeling, simultaneously solving for atomic level and ion populations for individual elements.

A&M databases↗

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward electron temperature profiles in hot-dense plasmas from x-ray spectral ensembles

High repetition rate laser systems enable new strategies for diagnosing plasma behavior with large datasets. Here, we define an ensemble technique that relies on randomized targeting of x-ray tracer micro-stripes. On each shot, a high-intensity laser pulse is focused on a solid target with Ti tracer stripes embedded in an Al foil, randomly targeting a micro-stripe, a portion of a stripe, or a gap between stripes. High-resolution, time-integrated x-ray spectrometers capture line emission from the portion of the micro-stripe that is heated to sufficiently high electron temperatures. Accumulation of many such cases is used to construct ensemble distributions of x-ray line intensities that encompass all relative offsets of the laser focus to the micro-stripe centers. Synthetic intensity distributions are likewise generated using collisional-radiative modeling. Bayesian fitting of modeled to measured intensity distributions establishes the most likely radial temperature profiles, enabling comparison to hydrodynamic models and calling into question the cylindrical symmetry of these micro-stripe-embedded systems. Ensemble techniques have significant potential for high-energy-density plasma diagnostics, especially with the advent of high repetition rate experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Helium metastable density determination in the COST reference source by absolutely calibrated optical emission spectroscopy

Helium metastable densities in the COST Reference Microplasma Jet are estimated for a variety of He/N2 admixtures and dissipated powers by applying a collisional-radiative model to absolutely calibrated optical emission spectroscopy measurements. This is accomplished by delineating the excitation mechanisms that result in the N2(C–B) and N2+(B–X) emission bands, the latter of which is strongly coupled to the presence of helium metastables. A number of other plasma parameters are established and discussed for each operating condition including the electron energy distribution function, reduced electric field, rate constants, and electron density. With these parameters, the reaction rates for the primary ionization pathways are also calculated, emphasizing the importance of helium metastables for discharge sustainment. Good agreement with the existing literature is found for most plasma parameters and for helium metastable densities, in particular. A clear [N2]−1 relationship between the nitrogen concentration and density of helium metastables is demonstrated, as has been identified in previous studies in analogous atmospheric pressure plasma jets. This validates the efficacy of this optical technique for determining helium metastable densities and establishes it as a viable, and in many cases, more accessible alternative to other means of quantifying helium metastables in low-temperature plasmas.

Physics↗

Experimental validation of a collision-radiation dataset for molecular hydrogen in plasmas

Quantitative spectroscopy of molecular hydrogen has generated substantial demand, leading to the accumulation of diverse elementary process data encompassing radiative transitions, electron-impact transitions, predissociations, and quenching. However, their rates currently available are still sparse, and there are inconsistencies among those proposed by different authors. In this study, we demonstrate an experimental validation of such a molecular dataset by composing a collisional-radiative model (CRM) for molecular hydrogen and comparing experimentally obtained vibronic populations across multiple levels. From the population kinetics of molecular hydrogen, the importance of each elementary process in various parameter space is studied. In low-density plasmas (electron density ne≲1017 m−3) the excitation rates from the ground states and radiative decay rates, both of which have been reported previously, determine the excited state population. The inconsistency in the excitation rates affects the population distribution the most significantly in this parameter space. However, in higher density plasmas (ne≳1018 m−3), the excitation rates from excited states become important, which have never been reported in the literature, and may need to be approximated in some way. In order to validate these molecular datasets and approximated rates, we carried out experimental observations for two different hydrogen plasmas; a low-density radio frequency heated plasma (ne≈1016 m−3) and the Large Helical Device (LHD) divertor plasma (ne≳1018 m−3). The visible emission lines from EF1Σg+, HH¯1Σg+, D1Πu±, GK1Σg+, I1Πg±, J1Δg±, h3Σg+, e3Σu+, d3Πu±,g3Σg+, i3Πg±, and j3Δg± states were observed simultaneously and their population distributions were obtained from their intensities. We compared the observed population distributions with the CRM prediction, in particular the CRM with the rates compiled by Janev et al., Miles et al., and those calculated with the molecular convergent close-coupling (MCCC) method. The MCCC prediction gives the best agreement with the experiment, particularly for the emission from the low-density plasma. However, the population distribution in the LHD divertor shows a worse agreement with the CRM than those from low-density plasma, indicating the necessity of the precise excitation rates from excited states. We also found that the rates for the electron attachment is inconsistent with experimental results. This requires further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron-impact excitation data for W 2+ in support of tungsten spectroscopy and re-deposition measurements for magnetically-confined plasmas

Abstract To better understand plasma wall interactions involving tungsten, accurate atomic structure and electron-impact driven collisional processes for near-neutral ion stages of tungsten are required. Complementing existing work on neutral and singly ionised tungsten, atomic structure and collisional calculations for W 2+ electron-impact excitation have been completed. These excitation calculations are an important component of S/XB coefficients for near-neutral charge states, which may be used to spectroscopically infer re-deposition of tungsten at the plasma-solid boundary of fusion relevant devices. With W 2+ in particular having emission lines that can be observed at ultraviolet (UV) wavelengths, while higher charge states of tungsten are unlikely to have lines possible to observe outside of the vacuum UV range. The atomic structure was generated using the General-purpose Relativistic Atomic Structure Package (GRASP 0 ), implementing the Multi-configuration Dirac Fock approach. This structure was the basis for a subsequent Dirac R -matrix electron-impact excitation calculation to provide Maxwellian averaged rate coefficients. A synthetic spectrum was generated from this data using a collisional-radiative model to predict the strongest W III spectral lines and these lines were compared to emission from the Compact Toroidal Hybrid (CTH) plasma device. Several of the strongest W III lines are observed in CTH and agree well with the modelled line wavelengths and intensities, a table of these lines is provided that could be observed in other devices.

McCann, M. (ORCID:0000000215321240)↗

Estimation of mean electron energy in helium surface ionization waves on dielectric substrates

Abstract The determination of basic plasma parameters in atmospheric pressure discharges is critical to advancing their use in applications. Atmospheric pressure plasma jets have found use in the fields of medicine, agriculture, material modification and others. Atmospheric pressure plasma jets often generate plasma surface ionization waves (SIW) which interact with and propagate over surfaces. Electrical diagnostics are challenging in SIW due to high collision frequencies and small scale of the plasma discharge. This work employs a passive optical emission line ratio technique to estimate the mean electron energy in SIW over planar dielectric substrates. The method uses an intensity ratio of two helium triplet lines: He( 3 3 S ) at 706.5 nm and He( 3 3 D ) at 587.56 nm. A collisional-radiative model is used to correlate line ratio to mean electron energy and determine dependencies on electron density and He/air gas mixture. Mean electron energies ranging from 3–8 eV are determined in He/air mixtures and are found to remain constant as the surface wave propagates radially. This work provides a 2D, time-resolved, mean electron energy diagnostic for surface ionization wave propagation and validation of numerical modeling in atmospheric pressure systems with spatially varying He/air gas mixtures. The model in question is designed for use with any He line ratio in the n = 3 excitation level.

Morsell, Joshua (ORCID:0000000331073002)↗