Engineering topics

Marmar, E. S.

Publications and source records attributed to Marmar, E. S..

At least 19 records

DOE OSTI2023

Further rotation reversal studies in C-Mod L-mode plasmas

Studies of core toroidal rotation reversal phenomenology in C-Mod deuterium L-mode plasmas have been expanded to include details of the dependences on plasma current and toroidal magnetic field. Rotation reversal occurs at a critical density, and universal scaling indicates that the product of ncritq95R ∼ BT/2, with ncrit in 1020/m3, R in m, and BT in T. Measurements in H and He plasmas exhibit similar behavior, including a connection with the linear Ohmic confinement/saturated Ohmic confinement transition and the cutoff for non-diffusive heat transport. Electron density and ion cyclotron range of frequencies power modulation experiments suggest that the collisionality ν* is a unifying parameter. Strong impurity puffing causes the critical density to increase, indicating that the situation is more complicated than only collisionality, perhaps involving the details of the effects of dilution on ion temperature gradient mode stability.

Physics↗

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 OSTI2023

Integrated deep learning framework for unstable event identification and disruption prediction of tokamak plasmas

Abstract The ability to identify underlying disruption precursors is key to disruption avoidance. In this paper, we present an integrated deep learning (DL) based model that combines disruption prediction with the identification of several disruption precursors like rotating modes, locked modes, H-to-L back transitions and radiative collapses. The first part of our study demonstrates that the DL-based unstable event identifier trained on 160 manually labeled DIII-D shots can achieve, on average, 84% event identification rate of various frequent unstable events (like H-L back transition, locked mode, radiative collapse, rotating MHD mode, large sawtooth crash), and the trained identifier can be adapted to label unseen discharges, thus expanding the original manually labeled database. Based on these results, the integrated DL-based framework is developed using a combined database of manually labeled and automatically labeled DIII-D data, and it shows state-of-the-art (AUC = 0.940) disruption prediction and event identification abilities on DIII-D. Through cross-machine numerical disruption prediction studies using this new integrated model and leveraging the C-Mod, DIII-D, and EAST disruption warning databases, we demonstrate the improved cross-machine disruption prediction ability and extended warning time of the new model compared with a baseline predictor. In addition, the trained integrated model shows qualitatively good cross-machine event identification ability. Given a labeled dataset, the strategy presented in this paper, i.e. one that combines a disruption predictor with an event identifier module, can be applied to upgrade any neural network based disruption predictor. The results presented here inform possible development strategies of machine learning based disruption avoidance algorithms for future tokamaks and highlight the importance of building comprehensive databases with unstable event information on current machines.

plasma instabilities↗

DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Argon Pumpout by ICRF Waves in C-Mod L- and I-mode Plasmas

Pumpout of argon ions by ICRF waves has been observed in C-Mod deuterium L- and I-mode plasmas that had a substantial hydrogen fraction. The effect is manifested by a reduction of core argon x-ray brightness up to a factor of 90% on time scales of tens of milliseconds following injection of ICRF power. For Ar^16+, the pumpout is strongest for hydrogen minority concentrations between 0.25 and 0.4, when the ICRF waves are not expected to result in minority heating. Modeling with the TORIC code suggests that the pumpout process occurs when the H/D mode conversion layer overlaps with the 2nd harmonic impurity resonance layer. The magnitude of the argon pumpout is independent of ICRF power above an apparent threshold of ~500 kW, independent of electron density and appears to decrease as the plasma current is increased. Potential application as a heavy impurity control tool in reactors is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Argon pumpout by ICRF waves in C-Mod L- and I-mode plasmas

Pumpout of argon ions by ICRF waves has been observed in C-Mod deuterium L- and I-mode plasmas that had a substantial hydrogen fraction. The effect is manifested by a reduction of core argon x-ray brightness up to a factor of 90% on time scales of tens of milliseconds following injection of ICRF power. For Ar 16+ , the pumpout is strongest for hydrogen minority concentrations between 0.25 and 0.4, when the ICRF waves are not expected to result in minority heating. Modeling with the TORIC code suggests that the pumpout process occurs when the H/D mode conversion layer overlaps with the 2nd harmonic impurity resonance layer. Here, the magnitude of the argon pumpout is independent of ICRF power above an apparent threshold of ~500 kW, independent of electron density and appears to decrease as the plasma current is increased. Potential application as a heavy impurity control tool in reactors is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Suppression of first-wall interaction in negative triangularity plasmas on TCV

Magnetically confined fusion plasmas with negative triangularity (d) exhibit greater L-mode confinement than with positive d. Recent experiments in the TCV and DIII-D tokamaks have correlated the confinement improvement to a reduction of fluctuations within the plasma core. We report on fluctuation measurements in the scrape-off layer (SOL) for −0.61 < d < +0.64 in limited and diverted ohmic L-mode plasmas; these reveal a strong reduction in SOL fluctuation amplitudes at d < −0.25, and, surprisingly, an almost full suppression of plasma interaction with the main-chamber first-wall, which could have important implications for the prospects of using negative d plasmas as a reactor solution. An exploration of several physical mechanisms suggests that a reduced connection length—intrinsic to negative d plasmas—plays a critical role in the origin of this phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗

DOE OSTI2021

The Very High n Rydberg Series of Ar^16+ in Alcator C-Mod Tokamak Plasmas

X-ray transitions of the very high-n Rydberg series in Ar^16+ have been observed from Alcator C-Mod tokamak plasmas. Individual emission lines up to 1s16p - 1s^2 have been resolved and the central chord line brightnesses with principal quantum number n between 7 and 16 are generally found to decay as 1/n^alpha, with alpha slightly larger than 3. In the plasma periphery, emission from 1s9p - 1s^2 and 1s10p - 1s^2 are found to be significantly enhanced relative to this decrease, indicative of selected population of these levels through charge exchange between background neutral deuterium in the ground state and Ar^17+. An unresolved feature between the wavelengths of 1s27p - 1s^2 and 1s30p - 1s^2 is also present, which arises through charge exchange with neutral deuterium in the n^* = 3 excited state. The brightnesses of transitions populated by charge exchange are spatially up/down asymmetric, with an excess on the side of the magnetic surface X-point. The relative brightness of the unresolved very high-n feature compared to 1s7p - 1s^2 is found to increase with electron temperature and decrease with electron density. Simulations of line emission just on the long wavelength side of the Ar^16$ ionization limit indicate that the principal quantum number decay exponent is closer to alpha = 4 at very high n. The brightness dependence on n below 16 is in excellent agreement with calculations from the Flexible Atomic Code package.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Modeling of particle transport, neutrals and radiation in magnetically-confined plasmas with Aurora

In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The very high n Rydberg series of Ar 16+ in Alcator C-Mod tokamak plasmas

We report x-ray transitions of the very high-n Rydberg series in Ar 16+ have been observed from Alcator C-Mod tokamak plasmas. Individual emission lines up to 1s16p-1s 2 have been resolved and the central chord line brightnesses with principal quantum number n between 7 and 16 are generally found to decay as 1/n α , with α slightly larger than 3. In the plasma periphery, emission from 1s9p-1s 2 and 1s10p-1s 2 are found to be significantly enhanced relative to this decrease, indicative of selected population of these levels through charge exchange between background neutral deuterium in the ground state and Ar 17+ . An unresolved feature between the wavelengths of 1s27p-1s 2 and 1s30p-1s 2 is also present, which arises through charge exchange with neutral deuterium in the n* = 3 excited state. The brightnesses of transitions populated by charge exchange are spatially up/down asymmetric, with an excess on the side of the magnetic surface X-point. The relative brightness of the unresolved very high-n feature compared to 1s7p-1s 2 is found to increase with electron temperature and decrease with electron density. Simulations of line emission just on the long wavelength side of the Ar 16+ ionization limit indicate that the principal quantum number decay exponent is closer to α = 4 at very high n. The brightness dependence on n below 16 is in excellent agreement with calculations from the flexible atomic code package.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Scenario adaptive disruption prediction study for next generation burning-plasma tokamaks

Next generation High Performance (HP) tokamaks risk damage from unmitigated disruptions at high current and power. Achieving reliable disruption prediction for a device's HP operation based on its Low Performance (LP) data is a key to its success. In this letter, through explorative data analysis and dedicated numerical experiments on multiple existing tokamaks, we demonstrate how the operational regimes of tokamaks can affect the power of a trained disruption predictor. First, our results suggest data-driven disruption predictors trained on abundant LP discharges work poorly on the HP regime of the same tokamak, which is a consequence of the distinct distributions of the tightly correlated signals related to disruptions in these two regimes. Second, we find that matching operational parameters among tokamaks strongly improves cross-machine accuracy which implies our model learns from the underlying scalings of dimensionless physics parameters like q 95 , β p and confirms the importance of these parameters in disruption physics and cross machine domain matching from the data-driven perspective. Finally, our results show in the absence of HP data from the target devices, the best predictivity of the HP regime for the target machine can be achieved by combining LP data from the target with HP data from other machines. Furthermore, these results provide a possible disruption predictor development strategy for next generation tokamaks, such as ITER and SPARC, and highlight the importance of developing baseline scenario discharges of future tokamaks on existing machines to collect more relevant disruptive data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dependence of alpha-particle-driven Alfven eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conven- tional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfven eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently- developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvenic or nearly sub-Alfvenic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Multi-device Studies of Pedestal Physics and Confinement in the I-mode Regime

This paper describes joint ITPA studies of the I-mode regime, which features an edge thermal barrier but with L-mode-like particle and impurity transport, and without ELMs. The regime has been demonstrated on Alcator C-Mod, ASDEX Upgrade and DIII-D tokamaks, over a wide range of device parameters and pedestal conditions, including low collisionality. Access to I-mode is independent of heating method (neutral beam injection, ICRH and/or ECRH). Normalized energy confinement H98y,2 ≥ 1 has been achieved for a range of 3≤q95≤4.9 and scales favourably with power. Changes in turbulence in the pedestal region accompany the transition from L-mode to I-mode. The L-I threshold increases with plasma density and current, and with device size, but weakly with BT. The upper limit of power for I-modes increases with toroidal field, and is largest on Alcator C-Mod at B > 5 T. Issues for extrapolation to ITER and other future fusion devices are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Inference of Experimental Radial Impurity Transport on Alcator C-Mod: Bayesian Parameter Estimation and Model Selection.

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. We find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dimensionless Parameter Scaling of Intrinsic Torque in C-Mod Enhanced Confinement Plasmas

A dimensionless parameter dependence study of intrinsic torque has been performed on a database of H- and I-mode plasmas from the Alcator C-Mod tokamak. The torque was determined by comparing intrinsic angular momentum density profiles just before and just after L-H and L-I transitions. The intrinsic torque has been found to scale as beta_N^1.5 rho_*^-1.0 nu_$^0.1, with the parameter ranges 0.3 < beta_N <1.5, 0.004 < rho_* < 0.011 and 0.04 < nu_* < 0.9. Comparison with results from other tokamaks suggests that the intrinsic torque should be normalized by some measure of the device size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗