Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma collisions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Collision operator for electron runaway in cold weakly ionised plasmas

In cold weakly ionised plasmas, the Dreicer generation mechanism can be non-diffusive as demonstrated in a recent study. By expanding the previous letter, we present the detailed description of a proper collision operator to precisely account for the non-diffusive electron kinetics. The operator appropriately combines the Fokker–Planck operator and Boltzmann operator where free-bound collision cross-sections are valid in the low-energy region. The proposed operator is envisaged to predict runaway electron generation in cold weakly ionised plasmas, particularly to design a runaway-free reactor tokamak start-up.

fusion plasma↗
DOE OSTI2025

Search for Quasiparticle Scattering in the Quark-Gluon Plasma with Jet Splittings in 𝑝⁢𝑝 and Pb-Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration reports measurements of the large relative transverse momentum (𝑘 𝑇 ) component of jet substructure in 𝑝⁢𝑝 and Pb-Pb collisions at center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV. Enhancement in the yield of such large-𝑘 T emissions in head-on Pb-Pb collisions is predicted to arise from partonic scattering with quasiparticles of the quark-gluon plasma. The analysis utilizes charged-particle jets reconstructed by the anti-𝑘 T algorithm with resolution parameter 𝑅 = 0.2 in the transverse-momentum interval 60 < 𝑝 T,ch,jet < 80 GeV/𝑐. The soft drop and dynamical grooming algorithms are used to identify high transverse momentum splittings in the jet shower. Comparison of measurements in Pb-Pb and 𝑝⁢𝑝 collisions shows medium-induced narrowing, corresponding to yield suppression of high-𝑘 𝑇 splittings, in contrast to the expectation of yield enhancement due to quasiparticle scattering. The measurements are compared to theoretical model calculations incorporating jet modification due to jet-medium interactions (“jet quenching”), both with and without quasiparticle scattering effects. These measurements provide new insight into the underlying mechanisms and theoretical modeling of jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical simulation of terahertz radiation by laser-driven plasma dipole oscillation

Laser-plasma interactions have many theoretical and technological applications. One is the use of coherently accelerated electrons to provide novel sources of THz radiation. Our research focuses on simulating the cross/self-interactions between two high intensity, ultra-short, counter propagating and detuned laser pulses and an initial neutral target for controlled ionization. Unlike our previous studies of laser-matter interaction over preformed plasma, we explore the injection and collision of laser pulses to induce background plasma driven by the self-guided laser wakefield mechanism, which is then used to perturb the plasma resulting in induced dipole oscillations leading to transverse radiation. Inducing a cylindrical spatial plasma column within the laser beam radius regime provides a stable, spatially localized plasma channel. The emitted radiation from the plasma dipole oscillation (PDO) will not be affected by surrounding plasma absorption, resulting in effective radiation distribution. Finally, results include 3D EM-PIC simulations and a comparison of the self- ionizing plasma against the preformed plasma to assess the efficiency of the mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗
DOE OSTI2025

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗
DOE OSTI2025

Impact of tracking resolution on 𝜙-meson spin alignment measurements

Measurements of global spin alignment of vector mesons in relativistic heavy-ion collisions can provide unique insights into spin-orbit interactions and vector meson dynamics in the quark-gluon plasma produced in those collisions. The global spin alignment is measured by the 00th coefficient of the spin density matrix, 𝜌 00 , via the polar angle (𝜃*) of the decay-kaon momentum in the parent rest frame with respect to the direction of the orbital angular momentum of the collision. Such measurements are affected by the angular and momentum resolutions of the reconstructed tracks in the experiment. Such effects are nontrivial because of kinematic complications caused by the boost to the parent rest frame, and could be important given that the global spin alignment signal is weak. In this paper, we investigate the effects of experimental tracking resolutions on measurements of the 𝜙 (1020) meson 𝜌 00 . We study these effects for two methods of 𝜌 00 measurements, the conventional method analyzing the 𝜙⁢−meson yield versus cos 2 ⁡𝜃* and the invariant mass (𝑚 inv ) method utilizing ⟨cos 2 ⁡𝜃*⟩ versus 𝑚 inv . Using typical resolution values from experiments, we find that the effect of track resolution on 𝜌 00 is small, well within typical measurement uncertainties.

Relativistic heavy-ion collisions↗
DOE OSTI2021

Microtearing Instabilities and Electron Thermal Transport in Low and High Collisionality NSTX Discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. The dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model [T. Rafiq, et al., Phys. Plasmas 23, 062507 (2016)]. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate.In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

fusion plasma technology↗
DOE OSTI2021

Driving rotamak currents with minimal power dissipation

The Rotamak is a proposed thermonuclear fusion device which employs rotating magnetic fields (RMF) to generate an azimuthal current to produce a field-reversed configuration. The efficiency of the currents that produce the field reversal by RMFs was debated some 40 years ago. The debate revolved around whether the currents would incur dissipation by the conventional Spitzer perpendicular resistivity, or whether some other relation between current and dissipation would be more appropriate. By employing an electron–ion pitch-angle scattering model, we find that the dissipation is non-Spitzer in nature. However, curiously, there appears to exist a regime where the power dissipated to maintain the current becomes vanishingly small.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A scaling law of the neutral penetration length and Balmer- α wing shape in high-temperature plasmas

Hydrogen atoms penetrating deep inside high-temperature magnetically confined plasmas by repetitive charge-exchange collisions result in a particle source of the plasma, which affects the plasma transport significantly. In this paper, we present an approximate solution of the fluid equations for neutral transport and an analytical representation of the neutral penetration length, in a simplified plasma geometry. This analysis predicts a power-law decay in the Balmer-α line wings which reflects the velocity distribution of the neutral atoms, with the power-law exponent analytically represented as well. These scaling laws are compared with a simple Monte–Carlo simulation and spectroscopic observations of Large Helical Device plasmas. Since the Balmer-α line wings are experimentally accessible, our formulation opens the possibility to quickly estimate the neutral penetration length from spectroscopic observations.

neutral opacity↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2021

The smallest fluid on earth

High energy heavy ion collisions create quark gluon plasmas that behave like almost perfect fluids. Very similar features to those that led to this insight have also been observed in experimental data from collisions of small systems, involving protons or other light nuclei. Here we describe recent developments aimed at understanding whether, and if so how, systems that produce relatively few particles (orders of magnitude less than in typical heavy ion collisions) and are only one to a few times the size of a proton, can behave like fluids. This involves a deeper understanding of fluid dynamics and its applicability, improvements of our understanding of the initial geometry of the collisions by considering fluctuations of the proton shape, as well as advancements in the calculation of initial state effects within an effective theory of quantum chromodynamics, which can affect the observables that are used to study fluid behavior. We further address open questions and discuss future directions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fall SULI Report

This fall, I interned under Fan Guo and Xiaocan Li in T-2 at Los Alamos National Laboratory through the SULI program. I investigated how shock thickness is impacted by collisions between particles in plasmas and how particle acceleration in shocks is impacted by the mass ratio for protons and electrons used in simulations. Both investigations were carried out using the VPIC simulation code. I found that collisions do not meaningfully impact shock width and that higher mass ratios lead to somewhat more electron acceleration in shocks. This write-up is in three parts, the first discussing a description of what I was up to personally throughout the fall, the second discussing the results of the collisional shock study, and the third discussing the results of the high mass ratio shock study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Viscosities of the Baryon-Rich Quark-Gluon Plasma from Beam Energy Scan Data

Here, this work presents the first Bayesian inference study of the (3 + 1)D dynamics of relativistic heavy-ion collisions and quark-gluon plasma viscosities using an event-by-event (3 + 1)D hydrodynamics + hadronic transport theoretical framework and data from the Relativistic Heavy Ion Collider Beam energy scan program. Robust constraints on initial state nuclear stopping and the baryon chemical potential-dependent shear viscosity of the produced quantum chromodynamic (QCD) matter are obtained. The specific bulk viscosity of the QCD matter is found to exhibit a preferred maximum around $\sqrt{s_{NN}}$ = 19.6 GeV . This result allows for the alternative interpretation of a reduction (and/or increase) of the speed of sound relative to that of the employed lattice-QCD based equation of state for net baryon chemical potential μ B ~0.2(0.4) GeV .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of charged hadron multiplicity in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV with the sPHENIX detector

The pseudorapidity distribution of charged hadrons produced in Au + Au collisions at a center-of-mass energy of $\sqrt{s_{NN}}$ = 200 GeV is measured using data collected by the sPHENIX detector. Charged hadron yields are extracted by counting cluster pairs in the inner and outer layers of the Intermediate Silicon Tracker, with corrections applied for detector acceptance, reconstruction efficiency, combinatorial pairs, and contributions from secondary decays. The measured distributions cover |η| < 1.1 across various centralities, and the average pseudorapidity density of charged hadrons at mid-rapidity is compared to predictions from Monte Carlo heavy-ion event generators. This result, featuring full azimuthal coverage at mid-rapidity, is consistent with previous experimental measurements at the Relativistic Heavy Ion Collider, thereby supporting the broader sPHENIX physics program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗