Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collisional plasma”

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

Closure and transport theory for plasmas with multiple ion species

This report summarizes a DOE-funded research project on closure and transport theory for plasmas with multiple ion species. The work developed accurate moment-based closure models across the full range of collisionality, incorporating both kinetic and collisional effects. Key contributions include high-fidelity closure relations, multi-temperature models, and transport coefficients applicable to fusion, space, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2021

Fe xvii 2p–3s Line Ratio Diagnostic of Shock Formation Radius in O Stars

The 2p–3s lines of Fe xvii in the X-ray spectrum of the O-type star ζ Puppis exhibit an anomalous (3G + M2)/(3F) line ratio of ∼1.4, in comparison with ∼2.4 for almost all other collisionally excited astrophysical spectra. Based on the work of Mauche et al., we conjectured that the strong UV field of ζ Puppis produces the observed ratio by depopulation of metastable 3s excited states, and that the ratio can potentially be used as an independent diagnostic of plasma formation radius. We used the Flexible Atomic Code collisional-radiative model to model the effect of UV photoexcitation from O stars on the Fe xvii lines. We compared our model calculations to archival spectra of coronal and hot stars from the Chandra HETGS and XMM-Newton RGS to benchmark our calculations for various electron densities and UV field intensities. Our calculations show that UV photoexcitation does not produce a sufficiently large dynamic range in the 3F / (3F + 3G + M2) fraction to explain the difference in the observed ratio between coronal stars and ζ Pup. Thus, this effect likely cannot explain the observed line ratio of ζ Pup, and its origin is still unexplained.

Gabriel J. Grell↗
NASA NTRS2021

Fe XVII 2p–3s Line Ratio Diagnostic of Shock Formation Radius in O Stars

The 2p–3s lines of Fe xvii in the X-ray spectrum of the O-type star ζ Puppis exhibit an anomalous (3G + M2)/(3F) line ratio of ∼1.4, in comparison with ∼2.4 for almost all other collisionally excited astrophysical spectra. Based on the work of Mauche et al., we conjectured that the strong UV field of ζ Puppis produces the observed ratio by depopulation of metastable 3s excited states, and that the ratio can potentially be used as an independent diagnostic of plasma formation radius. We used the Flexible Atomic Code collisional-radiative model to model the effect of UV photoexcitation from O stars on the Fe xvii lines. We compared our model calculations to archival spectra of coronal and hot stars from the Chandra HETGS and XMM-Newton RGS to benchmark our calculations for various electron densities and UV field intensities. Our calculations show that UV photoexcitation does not produce a sufficiently large dynamic range in the 3F / (3F + 3G + M2) fraction to explain the difference in the observed ratio between coronal stars and ζ Pup. Thus, this effect likely cannot explain the observed line ratio of ζ Pup, and its origin is still unexplained.

Gabriel Grell↗
NASA NTRS1977

Oxygen impurity radiation from Tokamak-like plasmas

We have constructed a nonhydrodynamic coronal model for calculating radiation from impurity atoms in a heated plasma. Some recent developments in the calculation of dielectronic recombination rate coefficients and collisional excitation rate coefficients are included. The model is applied to oxygen impurity radiation during the first few milliseconds of a TFR Tokamak plasma discharge, and good agreement with experimental results is obtained. Estimates of total line and continuum radiation from the oxygen impurity are given. It is shown that impurity radiation represents a considerable energy loss.

Rogerson, J. E.↗
DOE OSTI2025

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↗
NASA NTRS1990

Modeling solar flare conduction fronts. I - Homogeneous plasmas and ion-acoustic turbulence. II - Inhomogeneous plasmas and ambipolar electric fields

A one-dimensional, electrostatic, particle-in-cell simulation is used here to model the expansion of a heated electron population in a coronal loop during a solar flare and the characteristics of the associated X-ray emissions. The hot electrons expand outward from the localized region, creating an ambipolar electric field which accelerates a return current of cooler, ambient electrons. Ion-acoustic waves are generated by the return currents as proposed by Brown et al. (1979), but they play little or no role in containing energetic electrons and the conduction front proposed by Brown et al. does not form. The X-ray emission efficiency of the electrons is too low in the corona for them to be the source of hard X-ray bursts. The particle dynamics changes dramatically if the heated plasma is at low altitudes and expands upward into the more tenuous plasma at higher altitudes. Two important applications of this finding are the radio-frequency heating of the corona and the collisional heating of the chromosphere by precipitating energetic electrons. In both cases, the overlying plasma has a density that is too low to supply a balancing return current to the expanding hot electrons. As a result, an ambipolar electric field develops that tends to confine the energetic electrons behind a front that propagate outward at about the speed of sound.

Mckean, M. E.↗
DOE OSTI2022

How sheath properties change with gas pressure: modeling and simulation

Abstract Particle-in-cell simulations are used to study how neutral pressure influences plasma properties at the sheath edge. The high rate of ion–neutral collisions at pressures above several mTorr are found to cause a decrease in the ion velocity at the sheath edge (collisional Bohm criterion), a decrease in the edge-to-center density ratio ( h l factor), and an increase in the sheath width and sheath potential drop. A comparison with existing analytic models generally indicates favorable agreement, but with some distinctions. One is that models for the h l factor need to be made consistent with the collisional Bohm criterion. With this and similar corrections, a comprehensive fluid-based model of the plasma boundary transition is constructed that compares well with the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-dependent collisional radiative modeling of tungsten in the magnetic sheath for erosion diagnosis

Abstract Tungsten is the material of choice for the divertors in ITER, SPARC and future fusion reactors. Accurate diagnosis of tungsten erosion and migration is important for first wall life time, slag production and core performance. The addition of a magnetic presheath requires time-dependent collisional radiative effects to be included for accurate neutral tungsten collisional radiative modeling. Gross erosion measurements could be modified by a factor of 10 due to the inclusion of time-dependent effects for ITER relevant divertor conditions. A simple sputtering model and sheath density model are developed to investigate time-dependent collisional radiative effects. Neutral tungsten spectral lines populated from different metastable levels depend on model parameters leading to potential spectroscopic diagnostics of plasma parameters. Electron temperatures inferred from spectroscopic line ratios are in agreement with Langmuir probe measurements in the Compact Toroidal Hybrid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fast correlation heating in moderately coupled electron–ion plasmas

If the electrons in a plasma are suddenly heated, the resulting change in Debye shielding causes the ion kinetic energy to quickly increase. For the first time, this correlation heating , which is much faster than collisional energy exchange, is rigorously derived for a moderately coupled, electron–ion plasma. The electron–ion mass ratio is taken to be the smallest parameter in the Bogoliubov–Born–Green–Kirkwood–Yvon hierarchy, smaller even than the reciprocal of the plasma parameter. This ordering differs from conventional kinetic theory by making the electron collision rates faster than the ion plasma frequency, which allows stronger coupling and makes the ion heating a function only of the total energy supplied to the electrons. The calculation uses known formulae for correlations in a two-temperature plasma, for which a new, elementary derivation is presented. Suprathermal ions may be created more rapidly by this mechanism than by ion–electron Coulomb collisions. This means that the use of a femtosecond laser pulse could potentially help to achieve ignition in certain fast ignition approaches to inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Implementation and assessment of an extended hydrogenic molecular model in UEDGE

Abstract Two sets of pure deuterium plasmas are simulated using UEDGE, one including atoms and molecules and the other including atoms only as separate fluid species. Simulation results in one and two dimensions are reported to assess the role of molecules in tokamak plasmas. It is shown that thermal coupling of the molecules to the plasma can be relevant for the onset of detachment under highly collisional conditions in simple geometries, but play a modest role at low collisionality. Ion‐molecule equipartition presents an additional ion energy loss channel, dissipating ion energy that would otherwise heat the electrons by equipartition. The resulting reduced electron temperature increases the ion‐electron recombination by an order of magnitude as the plasma density is increased, yielding earlier and deeper detachment.

Holm, A.↗
DOE OSTI2022

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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)↗
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

Ion orbit loss is considered important for generating the radially inward electric field Er in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of Er in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from Er approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on Er is found to be only mild.

gyrokinetic simulations↗
DOE OSTI2026

Upgrades to the Thomson Scattering Diagnostic on the Lithium Tokamak eXperiment-Beta (LTX-$β$)

We describe the upgrades to the Thomson scattering (TS) diagnostic on the Lithium Tokamak eXperiment-beta (LTX- β ) to provide complete radial coverage of the plasma. LTX- β is a spherical tokamak with plasma limited at a major radius of R = 14 cm and R = 66 cm. Flat temperature profiles with peaked density were observed in LTX- β using an existing multipoint TS diagnostic with radial coverage of the plasma core. There is a need for more comprehensive core and edge TS measurements to thoroughly examine beam fueling, beam heating, and transport under low recycling conditions with liquid lithium walls and a lower collisionality scrape-off layer (SOL). The core TS set up provides plasma profiles from R = 40–60 cm. To extend the TS diagnostic coverage: 1) we added a TS viewing optics to cover the SOL region at five radial points from R = 57–64 cm, and 2) another TS viewing optics in the HFS region from R = 15–45 cm to obtain TS data from ten radial positions while making use of the original laser beam. The main focus of this article is on the new HFS TS setup that is equipped with better-matched collection optics located in the horizontal midplane where the laser polarization lies and views from a scattering angle of 130°–150° with respect to the laser path. A set of unequal-length Y-branched fiber bundles transports the TS light from two radial points to each polychromator. The new HFS view is better matched in spot size, scattering angle, and throughput to the laser and polychromators and has new V-grooved style viewing dumps. The new HFS TS view, along with the core and SOL TS views, will provide a full radial measurement coverage from R = 15–64 cm to explore the flat temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparing particle transport in JET and DIII-D plasmas: gyrokinetic and gyrofluid modelling

The paper models the impact of fueling versus particle transport in the DIII-D and JET tokamaks. Experimentally, it was found that in JET an increase in the density peaking at lower collisionality could be attributed to an increase in core fueling, however on DIII-D, the opposite was observed. In DIII-D, the increase in density peaking was the result of an increase in inward transport, not an increase in core fueling. The paper shows that the difference in temperature gradients and thus turbulence dynamics between the two experiments can explain the experimentally observed discrepancy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗