Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast ions”

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 271 records · Page 15

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Interpretation of electromagnetic modes in the sub-TAE frequency range in JET plasmas with elevated monotonic q-profiles

Recent JET deuterium experiments with an advanced tokamak scenario using an internal transport barrier (ITB) have been performed to clearly observe destabilized toroidicity-induced Alfvén eigenmodes (TAEs) by fast ions; interestingly, these also exhibit unstable electromagnetic (EM) perturbations in the sub-TAE frequency range. We identify such EM perturbations to be beta-induced ion temperature gradient (BTG) eigenmodes and not beta-induced Alfvén eigenmodes nor beta-induced Alfvén acoustic eigenmodes, which are often unstable in such high-beta plasmas with high-power neutral beam injection (NBI). The BTG modes are the most unstable modes due to the high thermal ion temperature gradient related to the ITB, high thermal ion temperature compared to thermal electron temperature (high T i /T⁠ e ), and a high ion beta regime. BTG mode experimental characteristics match the analytical theory, i.e., location in the vicinity of a rational magnetic surface with a low magnetic shear, mode frequency scaling with the ion diamagnetic frequency (⁠ω$^*_i$⁠), and a coupling among Alfvén and drift waves. We also perform linear gyrokinetic simulations with validated plasma profiles and equilibrium, and find a mode kinetically driven by thermal ions with similar characteristics as the experimental BTG modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Installation of a solid state neutral particle analyzer array on mega ampere spherical tokamak upgrade

A compact solid state neutral particle analyzer (SSNPA) diagnostic, previously installed at NSTX-U, has been moved to MAST-U and successfully operated in the first physics campaign (MU01). The SSNPA operates by detecting the flux of fast neutral particles produced by charge exchange (CX) reactions to diagnose the fast ion distribution. The diagnostic consists of three 16-channel sensors, which provide a radial view of the plasma and have a sightline intersection with the South–South neutral beam line. From this radial geometry, active CX signals from mostly trapped particles are observed. Each channel has a spatial resolution of 3–4 cm, a temporal resolution of 200 kHz, and an average pitch angle resolution of a few degrees. The three-sensor configuration allows for coarse energy resolution of the CX signals; each sensor sees similar sightlines but different filter thicknesses alter the energy cutoffs by known amounts. Experimental data show that all channels are collecting data as intended. The signal to noise ratio is typically around 15. Preliminary data analysis shows a correlation between the SSNPA signal and magnetohydrodynamic activity in the plasma as expected.

Instruments & Instrumentation↗

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of neutron emission and neutron energy spectrum calculations on a Mega Ampere Spherical Tokamak with directional relativistic spectrum simulator

The recently developed directional relativistic spectrum simulator (DRESS) code has been validated for the first time against numerical calculations and experimental measurements performed on the Mega Ampere Spherical Tokamak . In this validation, the neutron emissivities and rates computed by DRESS are benchmarked against TRANSP/NUBEAM predictions while the neutron energy spectra provided by DRESS taking as input TRANSP/NUBEAM and ASCOT/BBNBI in Gyro-Orbit mode fast ion distributions are validated against proton pulse height spectra measured by the neutron flux monitor. Finally, excellent agreement was found between DRESS and TRANSP/NUBEAM predictions of local and total neutron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predicting the rotation profile in ITER

Determining the toroidal rotation for future tokamaks like ITER is a challenging and important problem. By combining empirical scalings for the intrinsic rotation at the top of the pedestal with the expected neutral beam torque and modeling of momentum transport, the toroidal rotation profile for ITER is predicted with TGYRO using TGLF (SAT0 and SAT1). On axis rotation exceeds 20 krad s –1 and the E x B shear is significant enough to reduce turbulent transport and significantly increase confinement and fusion power when comparing to cases that ignore the effect of rotation. Here, the prediction of the rotation at the top of the pedestal is made with increased confidence due to experiments and modeling in DIII-D that have determined the importance of fast-ion and neutral particle transport effects on intrinsic rotation at this location. In particular, the effect of neutral particles on momentum transport in the pedestal region is found to be insignificant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic simulation of low-frequency Alfvénic modes in DIII-D tokamak

Global gyrokinetic simulations find that a beta-induced Alfvén eigenmode (BAE) and a Low frequency mode (LFM) coexist in the DIII-D tokamak experiments. The LFM can be excited without fast ions and has a frequency inside the gap of the beta-induced Alfvén-acoustic eigenmode (BAAE). The mode structure of the LFM from simulations and many of its parametric dependencies are consistent with experimental observations. For the BAE, the simulated mode structure is consistent with experiment but the frequency is lower than experiment. An antenna scan shows that the LFM is NOT the conventioanl BAAE; rather, it is an interchange-like electromagnetic mode excited by non-resonant drive of pressure gradients. The compressible magnetic perturbations significantly increase the growth rates of the BAE and LFM. On the other hand, effects of trapped electrons and equilibrium current have modest effects on the BAE and LFM. Finally, radial transfer of wave energy is different between the BAE and LFM, resulting in their respective triangular and elliptic mode shapes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stability of beta-induced Alfvén eigenmodes (BAE) in DIII-D

Although the stability of ellipticity, toroidal and reversed-shear Alfvén eigenmodes (EAE, TAE, RSAE) are relatively well understood, less is known about the stability of lower-frequency modes such as the beta-induced Alfvén eigenmode (BAE) but, because they are often unstable in present devices and are implicated in fast-ion transport, understanding their stability is vital. BAE stability is studied in primarily weak or reversed shear DIII-D plasmas with sub-Alfvénic deuterium beams. Modes are classified based on electron cyclotron emission, beam emission spectroscopy, magnetics, and interferometer data. The study is limited to the initial two seconds of the discharge, where the evolving q profile provides an effective scan of the dependence of stability upon q. In a dedicated experiment, BAEs are unstable at times in the discharge when the minimum of the safety factor q min is close to a rational number. The observed mode frequencies are usually close to analytic estimates of the BAE accumulation point and the eigenfunction peaks in the vicinity of q min . Unstable BAEs usually occur in bursts that chirp rapidly in frequency. To isolate the importance of thermal and beam gradients in driving the modes, the beam and electron cyclotron heating power is altered for 50–100 ms durations in reproducible discharges. As expected from the resonance condition, BAEs depend sensitively on the beam power and injection geometry. Modes only persist for ~25 ms because the anisotropic beam population only interacts strongly with the modes over a relatively narrow range of q. Lastly, a database of over 1000 beam-heated discharges shows that BAEs are more likely to be unstable when the poloidal beta exceeds 0.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic particle-induced geodesic acoustic modes on DIII-D

Various properties of the energetic particle-induced geodesic acoustic mode (EGAM) are explored in this large database analysis of DIII-D experimental data. EGAMs are n = 0 modes with m = 0 electrostatic potential fluctuations (where n/m = toroidal/poloidal mode number), m = 1 density fluctuations, and m = 2 magnetic fluctuations. The fundamental frequency (~20–40 kHz) of the mode is typically below that of the traditional geodesic acoustic mode frequency. EGAMs are most easily destabilized by beams in the counter plasma current (counter-I p ) direction as compared to co-Ip and off-axis beams. During counter beam injection, the mode frequency is found to have the strongest linear dependence (correlation coefficient r = –0.71) with the safety factor (q). Here, the stability of the mode in the space of q and poloidal beta (β p ) shows a clear boundary for the mode stability. The stability of the mode depends more strongly on damping rate than on fast-ion drive for a given injection geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perturbative analysis of low-frequency instabilities in high-field ST40 experiments

Abstract Linear and nonlinear perturbative analyses of low-frequency instabilities are presented for ST40 spherical torus high field plasmas, with simultaneous neutral beam heating due to two tangential sources at 25 keV and 55 keV. Initial plasma profiles on which the TRANSP and NUBEAM codes were based were obtained using an integrated analysis of several diagnostics including line-of-sight integrated and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. Magnetohydrodynamics (MHD) analysis using the NOVA/NOVA-K codes suggests that the modes in the experiments are core-localized n = 1 eigenmodes with mixed acoustic and electromagnetic Alfvénic polarization. The global transport combined with the perturbative mode analysis indicates that there is greater tendency for bursty chirping response as the level of background micro-turbulence decreases in time and those predictions are shown to be consistent with the experimental data. This finding is interpreted in terms of the suppression of coherent phase-space structures that support chirping when the resonant fast ion dynamics is dominated by frequent stochastic, orbit-decorrelating events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of performance enhancement by balanced double-null shaping in KSTAR

We report experimental observations on the effect of plasma boundary shaping towards balanced double-null (DN) configuration on the plasma performance in KSTAR. The transition from a single-null to a DN configuration resulted in improved plasma performance, manifested through changes in the pedestal region, decreased density, and core MHD activity variation. Specifically, the DN transition led to a wider and higher pedestal structure, accompanied by grassy edge-localized modes (ELMs) characteristics. The density decrease was a prerequisite for performance enhancement during DN shaping, increasing fast ion confinement. Optimizing the plasma near the core region was associated with the suppression of sawtooth instabilities and the occurrence of fishbone modes during the DN transition. Integrated modeling demonstrated that secondary effects of the DN shaping could increase core thermal energy confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗