Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical plasma physics”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Systematic study of flow vector fluctuations in $\sqrt{S_{NN}}$ = 5.02 TeV Pb-Pb collisions

Measurements of the 𝑝 T -dependent flow vector fluctuations in Pb–Pb collisions at $\sqrt{S_{NN}}$ = 5.02 TeV using azimuthal correlations with the ALICE experiment at the Large Hadron Collider are presented. A four-particle correlation approach [ALICE Collaboration, Phys. Rev. C 107, L051901 (2023)] is used to quantify the effects of flow angle and magnitude fluctuations separately. This paper extends previous studies to additional centrality intervals and provides measurements of the 𝑝 T -dependent flow vector fluctuations at $\sqrt{S_{NN}}$ = 5.02 TeV with two-particle correlations. Significant 𝑝 T -dependent fluctuations of the $\vec{𝑉}$ 2 flow vector in Pb–Pb collisions are found across different centrality ranges, with the largest fluctuations of up to ∼15% being present in the 5% most central collisions. In parallel, no evidence of significant 𝑝 T -dependent fluctuations of $\vec{𝑉}$ 3 or $\vec{𝑉}$ 4 is found. Additionally, evidence of flow angle and magnitude fluctuations is observed with more than 5⁢𝜎 significance in central collisions. These observations in Pb–Pb collisions indicate where the classical picture of hydrodynamic modeling with a common symmetry plane breaks down. This has implications for hard probes at high 𝑝 T , which might be biased by 𝑝 T -dependent flow angle fluctuations of at least 23% in central collisions. Given the presented results, existing theoretical models should be reexamined to improve our understanding of initial conditions, quark–gluon plasma properties, and the dynamic evolution of the created system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Theoretical analysis of energetic-ion-driven resistive interchange mode stabilization strategies using a Landau closure model

The aim of the present study is to perform a theoretical analysis of different strategies to stabilize energetic-ion-driven resistive interchange mode (EIC) in Large Helical Device (LHD) plasma. In this work, we use a reduced MHD for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the drift frequency to include their precessional drift. The stabilization trends of the 1/1 EIC observed experimentally with respect to the thermal plasma density and temperature are reproduced by the simulations, showing a reasonable agreement with the data. The LHD operation scenarios with stable 1/1 EIC are identified, leading to the stabilization of the 1/1 EIC if the thermal plasma density and temperature are above a given threshold. Additionally, the 1/1 EIC are also stabilized if the rotational transform is modified in a way that the 1/1 rational surface is located further away than 0.9 times the normalized radius, or the magnetic shear in the plasma periphery is enhanced. Also, LHD discharges with large magnetic fields show a higher EIC destabilization threshold with respect to the thermal plasma density. If the perpendicular NBI deposition region is moved further inward than 0.875 times the normalized radius the 1/1 EIC are also stabilized. Furthermore, increasing the perpendicular NBI voltage such that the EP energy is higher than 30 keV stabilizes the 1/1 EIC. Moreover, deuterium plasmas show a higher stability threshold for the 1/1 EIC than hydrogen plasmas. The experimental data shows a larger time interval between EIC events as the power of the tangential NBI is increased providing that the perpendicular NBI power is at least 13 MW. This implies a stabilizing effect of the tangential NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
DOE OSTI2022

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

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

Estimating cross-field particle transport at the outer midplane of TCV by tracking filaments with machine learning

Abstract Cross-field transport of particles in the boundary region of magnetically confined fusion plasmas is dominated by turbulence. Blobs, intermittent turbulent structures with large amplitude and a filamentary shape appearing in the scrape-off layer (SOL), are known from theoretical and experimental studies to be the main contributor to the cross-field particle transport. The dynamics of blobs differs depending on various plasma conditions, including triangularity ( δ ). In this work, we analyze triangularity dependence of the cross-field particle transport at the outer midplane of plasmas with δ = + 0.38 , +0.15, −0.14, and −0.26 on the Tokamak à Configuration Variable, using our novel machine learning (ML) blob-tracking approach applied to gas puff imaging data. The cross-field particle flux determined in this way is of the same order as the overall transport inferred from KN1D, GBS, and SOLPS-ITER simulations, suggesting that the blobs identified by the ML blob-tracking account for most of the cross-field particle transport in the SOL. Also, the ML blob-tracking and KN1D show a decrease in the cross-field particle transport as δ becomes more negative. The blob-by-blob analysis of the result from the tracking reveals that the decrease of cross-field particle transport with decreasing δ is accompanied by a decrease in the number of blobs in a fixed time, which tend to have larger area and lower radial speed. Also, the blobs in these plasmas are in the connected sheath regime, and show a velocity scaling consistent with the two-region model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Azimuthal Angle Correlations of Muons Produced via Heavy-Flavor Decays in 5.02 TeV Pb+Pb and pp Collisions with the ATLAS Detector

Angular correlations between heavy quarks provide a unique probe of the quark-gluon plasma created in ultrarelativistic heavy-ion collisions. Results are presented of a measurement of the azimuthal angle correlations between muons originating from semileptonic decays of heavy quarks produced in 5.02 TeV Pb+Pb and pp collisions at the LHC. The muons are measured with transverse momenta and pseudorapidities satisfying p$^μ_T$ > 4 GeV and |η μ | < 2.4, respectively. The distributions of azimuthal angle separation Δ⁢Φ for muon pairs having pseudorapidity separation |Δη| > 0.8, are measured in different Pb +Pb centrality intervals and compared to the same distribution measured in pp collisions at the same center-of-mass energy. Results are presented separately for muon pairs with opposite-sign charges, same-sign charges, and all pairs. A clear peak is observed in all Δ⁢Φ distributions at Δ⁢Φ ~ π, consistent with the parent heavy-quark pairs being produced via hard-scattering processes. The widths of that peak, characterized using Cauchy-Lorentz fits to the Δ⁢Φ distributions, are found to not vary significantly as a function of Pb +Pb collision centrality and are similar for pp and Pb +Pb collisions. This observation will provide important constraints on theoretical descriptions of heavy-quark interactions with the quark-gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Investigating Kinetic Mechanisms of Soot Formation in Plasma Pyrolysis of Methane via Active Learning (Final Technical Report)

Plasma pyrolysis of methane is an effective route for zero-carbon hydrogen production. Yet, soot generated from pyrolysis of hydrocarbons is detrimental to the climate and human health. There is ample experimental and theoretical evidence that suggests polycyclic aromatic hydrocarbons (PAHs) are the molecular precursors to soot particles. The reaction pathways of PAH formation are intricately dependent on a multitude of process parameters, whose kinetic mechanisms are not well-understood in plasma pyrolysis. This project aims to leverage advances in the kinetic modeling of soot formation in combustion, as well as in surrogate modeling and active learning, to systematically investigate the effects of process parameter on the kinetics of PAH formation in plasma pyrolysis of methane. To this end, we propose to use the PAH formation kinetics model developed by the PPPL/PU group based on the well-established ABF and HACA mechanisms, coupled with low-temperature plasma models. We will develop an active learning (AL) framework based on Bayesian optimization to systematically and data-efficiently explore the complex and multivariable parameter space of plasma pyrolysis in order to quantify the effects of plasma and feed parameters on the ABF and HACA kinetic pathways. AL is the branch of machine learning concerned with systematically querying samples from a system (experimental or computational) to train a data-driven model that maps design parameters to a performance criterion. We will use the data generated via AL to perform global sensitivity analysis, combined with uncertainty quantification, to elucidate the impact of different reaction pathways on minimizing formation of soot precursors. This study will result in an improved understanding of kinetics of PAH formation in plasma pyrolysis and can pave the way for more advanced mechanistic studies (e.g., soot nucleation mechanisms). Additionally, the findings will be useful for establishing practical strategies for increasing the pyrolysis efficiency and producing high-grade carbon for synthesis of nanomaterials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma-Assisted Epitaxy of Piezoelectric Sc x Al 1-x N Films on Sapphire for Use in Harsh-Environment Microwave Acoustic Sensors

The Sc x Al 1-x N wurtzite structure has been shown theoretically and experimentally to exhibit significantly higher piezoelectric coupling compared to pure AlN. In this work, a plasma-assisted epitaxial growth method has been used to synthesize epitaxial (0002) Sc x Al 1-x N films on c-sapphire substrates from x = 0.07 to 0.30 by co-evaporating high-purity Sc and Al sources in the presence of a nitrogen plasma generated by an RF plasma source. Epitaxial Sc x Al 1-x N films with highly oriented (0002) grains and in-plane registry were produced on c-sapphire substrates that were pre-exposed to the nitrogen plasma to form an oxynitride seed layer. Growth of Sc x Al 1-x N films was carried out at 930°C under both metal-rich and N-rich conditions using precisely controlled Sc, Al, and N-plasma fluxes. Metal-rich depositions yielded non-(0002)-oriented Sc x Al 1-x N grains and intermetallic ScAl grains. Nitrogen-rich growth with a Sc/Al flux ratio of 1/3 produced the best (0002) epitaxy as determined by x-ray diffraction analysis. Surface acoustic wave resonator (SAWR) devices were fabricated from 500-nm-thick Sc x Al 1-x N and AlN films to extract their electromechanical coupling coefficients, k 2 . As the Sc concentration in the films increases, the degree of (0002) epitaxy is reduced, yet the value of k 2 increases becasue there is more Sc in the wurzite lattice despite the decreased level of (0002) grain alignment. As a result, the use of a 10-nm-thick Si x N y capping layer on top of the Sc x Al 1-x N films aids in preventing etching during SAWR device photolithography and also helps hinder film oxidation up to 800°C.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of the electron density and temperature gradient on drift-wave turbulence in the Large Plasma Device

In this paper we present an experimental study of edge turbulence in the Large Plasma Device at UCLA. We utilize a scan of discharge power and prefill pressure (neutral density) to show experimentally that turbulent density fluctuations decrease with decreasing density gradient, as predicted for resistive drift-wave turbulence (RDWT). As expected for RDWT, we observe that the cross-phase between the density and potential fluctuations is close to 0. Moreover, the addition of an electron temperature gradient leads to a reduction in the amplitude of the density fluctuations, as expected for RDWT. However, counter to theoretical expectations, we find that the potential fluctuations do not follow the same trends as the density fluctuations for changes either in density gradients or the addition of a temperature gradient. The disconnect between the density and potential fluctuations is connected to changes in the parallel flows as a result of differences in the prefill pressure, i.e. neutral density. Further analysis of the density and potential fluctuation spectra show that the electron temperature gradient reduces the low frequency fluctuations up to $10 \,{\rm kHz}$ and the introduction of a temperature gradient leads to an unexpected ${\sim }{\rm \pi}$ shift of the density–potential cross-phase at ${\sim }10\,{\rm kHz}$ , while maintaining the typical resistive drift-wave cross-phase at lower frequencies. These experiments partly confirm existing knowledge on resistive drift-wave turbulence, but also introduce new observations that indicate a need for dedicated nonlinear three-dimensional turbulence simulations that include neutrals.

Physics↗
DOE OSTI2022

Measurement of prompt D s + -meson production and azimuthal anisotropy in Pb–Pb collisions at s NN = 5.02 TeV

The production yield and angular anisotropy of prompt $D^{+}_{s}$ mesons were measured as a function of transverse momentum (p T ) in Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV collected with the ALICE detector at the LHC. $D^{+}_{s}$ mesons and their charge conjugates were reconstructed at midrapidity (|y| < 0.5) from their hadronic decay channel $D^{+}_{s}$ → φπ + , with φ → K - K + , in the p T intervals 2 < p T < 50 GeV/c and 2 < p T < 36 GeV/c for the 0–10% and 30–50% centrality intervals. For p T > 10 GeV/c, the measured $D^{+}_{s}$ -meson nuclear modification factor R AA is consistent with the one of non-strange D mesons within uncertainties, while at lower pT a hint for a $D^{+}_{s}$ -meson R AA larger than that of non-strange D mesons is seen. The enhanced production of $D^{+}_{s}$ relative to non- strange D mesons is also studied by comparing the pT-dependent $D^{+}_{s}$ /D 0 production yield ratios in Pb–Pb and in pp collisions. The ratio measured in Pb–Pb collisions is found to be on average higher than that in pp collisions in the interval 2 < p T < 8 GeV/c with a significance of 2.3σ and 2.4σ for the 0–10% and 30–50% centrality intervals. The azimuthal anisotropy coefficient v2 of prompt $D^{+}_{s}$ mesons was measured in Pb–Pb collisions in the 30–50% centrality interval and is found to be compatible with that of non-strange D mesons. The main features of the measured RAA, $D^{+}_{s}$ /D 0 ratio, and v 2 as a function of p T are described by theoretical calculations of charm-quark transport in a hydrodynamically expanding quark–gluon plasma including hadronisation via charm-quark recombination with light quarks from the medium. The p T -integrated production yield of $D^{+}_{s}$ mesons is compatible with the prediction of the statistical hadronisation model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Anisotropic jet broadening and jet shape

In this paper, we explore the use of jet substructure as a way of probing phenomena which break the isotropic behavior of jets, such as jet propagation through an anisotropically flowing quark-gluon plasma or spin correlations. We introduce two novel observables for this purpose: the azimuthal-dependent jet broadening and the azimuthal-dependent jet shape, which generalize the traditional isotropic substructure studies. Using Soft-Collinear Effective Theory, we explicitly calculate the jet functions associated with these observables with a standard jet axis and with a Winner-Take-All jet axis in both the resummed and fixed order limits. While our analysis first and foremost establishes the formalism for the azimuthal-dependent jet substructure, it also brings to light new results for jet substructure in the azimuthally integrated case, such as the semi-inclusive jet function and the exclusive jet shape for the Winner-Take-All axis, and the jet broadening in the fixed order region. As an illustrative example for the new formalism we demonstrate that the azimuthal-dependent jet broadening can be used as a direct probe of the transversity parton distribution function in deep inelastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Study of the Effect of Reflections on High-Power, 110-GHz Pulsed Gyrotron Operation

The effect of reflection is studied experimentally and theoretically on a high-power 110-GHz gyrotron operating in the TE 22,6 mode in 3 μs pulses at 96 kV, 40 A. The experimental setup allows variation of the reflected power from 0 to 33% over a range of gyrotron operating conditions. The phase of the reflection is varied by translating the reflector along the axis. Operating at a higher efficiency point, at 4.40 T with 940kW of output power, reflected power exceeding 11% causes a switch from operation in the TE 22,6 to simultaneous operation in the TE 22,6 and TE 21,6 modes with a large decrease of the total gyrotron output power. This switching effect is in good agreement with simulations using the MAGY code. Operating at a more stable point, 4.44 T with 580kW of output power, when the reflection is increased, the output power remains in the TE 22,6 mode but it decreases monotonically with increasing reflection, dropping to 200kW at 33% reflection. Furthermore, at a reflection above 22%, a power modulation at 25 to 30MHz is observed, independent of the phase of the reflected wave. Such a modulated signal may be useful in spectroscopic and other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurement of azimuthal anisotropy of muons from charm and bottom hadrons in Pb+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

Azimuthal anisotropies of muons from charm and bottom hadron decays are measured in Pb+Pb collisions at s NN = 5.02 TeV . The data were collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2018 with integrated luminosities of 0.5 nb - 1 and 1.4 n b - 1 , respectively. The kinematic selection for heavy-flavor muons requires transverse momentum 4 <lt; p T < 30 GeV and pseudorapidity | η | < 2.0 . The dominant sources of muons in this p T range are semi-leptonic decays of charm and bottom hadrons. These heavy-flavor muons are separated from light-hadron decay muons and punch-through hadrons using the momentum imbalance between the measurements in the tracking detector and in the muon spectrometers. Azimuthal anisotropies, quantified by flow coefficients, are measured via the event-plane method for inclusive heavy-flavor muons as a function of the muon p T and in intervals of Pb+Pb collision centrality. Heavy-flavor muons are separated into contributions from charm and bottom hadron decays using the muon transverse impact parameter with respect to the event primary vertex. Non-zero elliptic ( v 2 ) and triangular ( v 3 ) flow coefficients are extracted for charm and bottom muons, with the charm muon coefficients larger than those for bottom muons for all Pb+Pb collision centralities. The results indicate substantial modification to the charm and bottom quark angular distributions through interactions in the quark-gluon plasma produced in these Pb+Pb collisions, with smaller modifications for the bottom quarks as expected theoretically due to their larger mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theoretical study of the Alfven eigenmode stability in CFETR steady state discharges

The aim of this study is to analyze the stability of Alfven eigenmodes (AE) in the China Fusion Engineering Test Reactor (CFETR) plasma for steady state operations. The analysis is done using the gyro-fluid code FAR3d including the effect of the acoustic modes, EP finite Larmor radius damping effects and multiple energetic particle populations. Two high poloidal β scenarios are studied with respect to the location of the internal transport barrier (ITB) at r/a ≈ 0.45 (case A) and r/a ≈ 0.6 (case B). Both operation scenarios show a narrow TAE gap between the inner-middle plasma region and a wide EAE gap all along the plasma radius. In this work, the AE stability of CFETR plasmas improves if the ITB is located inwards, case A, showing AEs with lower growth rates with respect to the case B. The AEs growth rate is smaller in the case A because the modes are located in the inner-middle plasma region where the stabilizing effect of the magnetic shear is stronger with respect to the case B. Multiple EP populations effects (NBI driven EP + alpha articles) are negligible for the case A, although the simulations for the case B show a stabilizing effect of the NBI EP on the n = 1 BAE caused by α particles during the thermalization process. If the FLR damping effects are included in the simulations, the growth rate of the EAE/NAE decreases up to 70%, particularly for n > 3 toroidal families. Low n AEs (n < 6) show the largest growth rates. On the other hand, high n modes (n = 6 to 15) are triggered in the frequency range of the NAE, strongly damped by the FLR effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Onset of Constituent Quark Number Scaling in Heavy-Ion Collisions at RHIC

Partonic collectivity is one of the necessary signatures for the formation of quark-gluon plasma in high-energy nuclear collisions. Number of constituent quarks (NCQ) scaling has been observed for hadron elliptic flow 𝑣2 in top energy nuclear collisions at the Relativistic Heavy Ion Collider and the LHC, and this has been theoretically suggested as strong evidence for partonic collectivity. In this Letter, a systematic analysis of 𝑣 2 of 𝜋 ± , 𝐾 ± , 𝐾$^{0}_{𝑆}$, 𝑝, and Λ in Au+Au collisions at $\sqrt{s_{NN}}$ = 3.2, 3.5, 3.9, and 4.5 GeV, with the STAR experiment at the Relativistic Heavy Ion Collider, is presented. NCQ scaling is markedly violated at 3.2 GeV, consistent with a hadronic-interaction dominated equation of state. However, as the collision energy increases, a gradual evolution to NCQ scaling is observed. In conclusion, this beam-energy dependence of 𝑣 2 for all hadrons studied provides evidence for the onset of dominant partonic interactions by $\sqrt{s_{NN}}$ = 4.5 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗