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 73 records · Page 4

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

Physics↗

Study of quasi-collisional effects in laboratory and astrophysical plasmas

High-amplitude turbulence excited in plasmas at small-scales can leave imprint on the radiation produced by (accelerated) plasma particles -- mostly electrons in electron-ion plasmas and both electrons and positrons in lepton pair plasmas. Furthermore, turbulence is known to introduce "effective collisional" effects -- anomalous resistivity, dissipation, diffusion, etc. -- in the otherwise collisionless plasmas. These effective collisions affect radiation transfer and modify optical and magneto-optic effects -- the transmittance and reflectivity coefficients, Faraday rotation, etc. Note that the effective collisions considered in this project are not resonant wave-particle interactions, but much less conventional randomization of particles’ paths. Therefore, we colloquially refer them to as "quasi-collisions". The study of quasi-collisions on transport, radiative and optical properties of plasmas is the main focus of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Noninvasive THz-TDS measurements of plasma bounded and optically shielded by Hall thruster wall material

We experimentally demonstrate the capability of terahertz time-domain spectroscopy (THz-TDS) to noninvasively measure the electron density and collision frequency of plasma bounded and optically shielded by Hall thruster wall material. This paper augments the standard THz-TDS plasma diagnostic theory to account for plasma boundaries, presents THz optical property measurements of three different wall materials (grades M, M26, and HP boron nitride composite), and provides electron density and collision frequency measurements of an inductively coupled plasma bounded and optically shielded by each wall material. We find that the electron density measurement capability is weakly impacted by the boundaries, whereas the electron collision frequency measurement capability is strongly reduced by the boundaries. Here, the bounded plasma electron density trends deviate substantially from those of the unbounded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anomalous cross-field transport in a Hall thruster inferred from direct measurement of instability growth rates

The contribution of the electron drift instability to anomalous electron transport is experimentally assessed in a Hall effect discharge. The transport is represented by an anomalous collision frequency, which is related through quasilinear theory to the energy and growth rate of the instability. The wave energy is measured directly with ion saturation probes, while estimates of the growth rate are employed based on both linearized theory and direct measurement. The latter measurement is performed with a bispectral analysis method. The wave-driven collision frequency is compared to measurements of the actual collision frequency inferred from a method based on laser- induced fluorescence. It is found that estimates for transport using linearized theory for the growth differ by over an order of magnitude from the actual anomalous collision frequency in the plasma. The wave-driven anomalous collision frequency with measured growth, however, is shown to agree with the electron collision frequency in magnitude and capture aspects of the trends in spatial variation. This result demonstrates experimentally that wave-driven effects ultimately can explain the observed cross-field transport in these devices. As a result, the implications of this finding are discussed in the context of the key lengthscales that drive the transport as well as the implications identifying reduced fidelity models that could be used to predict anomalous collision frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anisotropic Flow in Fixed-Target 208 Pb + 20 Ne Collisions as a Probe of Quark-Gluon Plasma

The System for Measuring Overlap with Gas (SMOG2) at the LHCb detector enables the study of fixed-target ion-ion collisions at relativistic energies ($\sqrt{𝑠_{NN}}$ ∼ 100 GeV in the center of mass). Here, with input from ab initio calculations of the structure of 16 O and 20 Ne , we compute 3+1⁢D hydrodynamic predictions for the anisotropic flow of Pb+Ne and Pb+O collisions to be tested with upcoming LHCb data. This will allow the detailed study of quark-gluon plasma formation as well as experimental tests of the predicted nuclear shapes. Elliptic flow (𝑣 2 ) in Pb + Ne collisions is greatly enhanced compared to the Pb + O baseline due to the shape of 20 Ne , which is deformed in a bowling-pin geometry. Owing to the large 208 Pb radius, this effect is seen in a broad centrality range, a unique feature of this collision configuration. Larger elliptic flow further enhances the quadrangular flow (𝑣 4 ) of Pb + Ne collisions via nonlinear coupling, and impacts the sign of the kurtosis of the elliptic flow vector distribution (𝑐 2 ⁡{4}). Exploiting the shape of 20 Ne proves thus an ideal method to investigate the formation of quark-gluon plasma in fixed-target experiments at LHCb, and demonstrates the power of System for Measuring Overlap with Gas as a tool to image nuclear ground states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

$$\Lambda $$ polarization in very high energy heavy ion collisions as a probe of the quark–gluon plasma formation and properties

Abstract We have studied the spin polarization of $$\Lambda $$ Λ hyperons in heavy ion collisions at center-of-mass energies $$\sqrt{s_\textrm{NN}}= 200$$ s NN = 200 GeV and $$\sqrt{s_\textrm{NN}}= 5.02$$ s NN = 5.02 TeV carried out at RHIC and LHC colliders. We have calculated the mean spin vector at local thermodynamic equilibrium, including all known first-order terms in the gradients of the thermo-hydrodynamic fields, assuming that the hadronization hypersurface has a uniform temperature. We have also included the feed-down contributions to the polarization of $$\Lambda $$ Λ stemming from the decays of polarized $$\Sigma ^*$$ Σ ∗ and $$\Sigma ^0$$ Σ 0 hyperons. The obtained results are in good agreement with the data. In general, the component of the spin vector along the global angular momentum, orthogonal to the reaction plane, shows strong sensitivity to the initial longitudinal flow velocity. Furthermore, the longitudinal component of the spin vector turns out to be very sensitive to the bulk viscosity of the plasma at the highest LHC energy. Therefore, the azimuthal dependence of spin polarization can effectively constrain the initial hydrodynamic conditions and the transport coefficients of the quark gluon plasma.

Physics↗

Phenomenological study on correlation between flow harmonics and mean transverse momentum in nuclear collisions

To assess the properties of the quark-gluon plasma formed in nuclear collisions, the Pearson correlation coefficient between flow harmonics and mean transverse momentum, \rho\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v n 2 , [ p T ] ) , reflecting the overlapped geometry of colliding atomic nuclei, is measured. \rho\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v 2 2 , [ p T ] ) was found to be particularly sensitive to the quadrupole deformation of the nuclei. We study the influence of the nuclear quadrupole deformation on \rho\left(v_{n}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v n 2 , [ p T ] ) in Au+Au and U+U collisions at RHIC energy using AMPT transport model, and show that the \rho\left(v_{2}^{2},\left[p_{\mathrm{T}}\right]\right) ρ ( v 2 2 , [ p T ] ) is reduced by the quadrupole deformation \beta_2 β 2 and turns to change sign in ultra-central collisions (UCC).

Zhang, Chunjian↗

Neural network model of neutral beam injection in the EAST tokamak to enable fast transport simulations

The neutral beam injection (NBI) system in EAST produces energetic neutral particles, which collide with electrons and ions in tokamak plasmas and heat the plasmas through Coulomb collisions. Moreover, it drives a non-inductive source of current, due to the charge-exchange collision between neutral particles and ions, and injects toroidal torque, which generates a toroidal rotation of the plasma. The effect caused by the NBI system, such as plasma heating, current drive, total neutron rate, momentum transfer, and shine-through, are modeled by a comprehensive module called NUBEAM. However, NUBEAM is computationally intensive since it relies on Monte Carlo methods. In this work, a neural network model has been developed as a surrogate model for NUBEAM in EAST. The database for neural-network model training, validation and testing is generated by running TRANSP for experimental discharges from recent EAST campaigns (after the latest NBI upgrade) while using the NUBEAM module. Simulation results illustrate that the trained neural network has the capability of replicating the predictions made by NUBEAM while demanding a significantly shorter execution time. Finally, these results indicate that surrogate models like the one proposed in this work could enable fast transport simulations for EAST after integrating them into a control-oriented predictive code such as COTSIM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗

Measuring dynamics of differentially rotating, unmagnetized, free-boundary plasma produced by soft x-ray irradiation

We present results from a new experiment, fielded on the MAGPIE pulsed power generator, producing a rotating, boundary-free plasma with no pre-imposed magnetic field. Angular momentum is introduced to the rotating system by the oblique collision of multiple plasma flows, which are driven by x-ray ablation of solid targets using the soft x-rays emitted from the implosion of wire array z-pinches. This produces a hydrodynamically stable plasma which undergoes ~2-3 rotations over the duration of the experiment, significantly more than previous pulsed-power platforms. Estimating the angular frequency from the electron density profile in the inner part of the rotating plasma shows that it has a quasi-Keplerian rotation profile. The system also allows for the addition of a controllable magnetic field. This, combined with the sufficiently large Reynolds (~10 5 ) and magnetic Reynolds (~10) numbers, will enable investigation of the effect of magnetic field on the structure and stability of the rotating plasma.

Physics - Plasma physics↗

A theoretical model for sheath dynamics in the two-plasma mode of a cylindrical filament discharge

We present a refined understanding of cylindrical plasma-facing emitter cathodes by examining the limitations of the conventional space charge limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge limited regimes based on electron current constraints, our findings reveal that trapped ions can form a quasi-neutral “second plasma,” expanding upstream and interacting with the primary plasma. In cylindrical and other non-planar geometries, an expanding trapped-ions layer produces current enhancement by converting more of the emitted electron flux to passing flux. The current enhancement leads to a complex coupling between the trapped-ions plasma and upstream plasma through ionizing collisions in the upstream region and charge exchange collisions in the virtual cathode. The coupling results in mutual density enhancement and at the same time, antagonistic expansion dynamics between the two plasma layers. The process is understood using a new “aid-compete” model. Numerical solutions, validated against simulations, reveal the model's capacity to predict system evolution and complex discharge behaviors, albeit with qualitative approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of Suppressed Charged-Particle Production in Ultrarelativistic Oxygen-Oxygen Collisions

A hot and dense state of nuclear matter, known as the quark-gluon plasma, is created in collisions of ultrarelativistic heavy nuclei. Highly energetic quarks and gluons, collectively referred to as partons, lose energy as they travel through this matter, leading to suppressed production of particles with large transverse momenta (𝑝 T ). Conversely, high-𝑝 T particle suppression has not been seen in proton-lead collisions, raising questions regarding the minimum system size required to observe parton energy loss. Oxygen-oxygen (OO) collisions examine a region of effective system size that lies between these two extreme cases. The CMS detector at the CERN LHC has been used to quantify charged-particle production in inclusive OO collisions for the first time via measurements of the nuclear modification factor (𝑅 AA ). The 𝑅 AA is derived by comparing particle production to expectations based on proton-proton (𝑝⁢𝑝) data and has a value of unity in the absence of nuclear effects. The data for OO and 𝑝⁢𝑝 collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{NN}}$ = 5.36 TeV correspond to integrated luminosities of 6.1 nb −1 and 1.02 pb −1 , respectively. The 𝑅 AA is below unity with a minimum of 0.69 ± 0.04 around 𝑝 T = 6 GeV. The data exhibit better agreement with theoretical models incorporating parton energy loss as compared to baseline models without energy loss.

FOS: Physical sciences↗

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗