Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D.

Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. In the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.

Multiscale Chirping Modes Driven by Thermal Ions in a Plasma with Reactor-Relevant Ion Temperature.

A thermal ion driven bursting instability with rapid frequency chirping, considered as an Alfvénic ion temperature gradient mode, has been observed in plasmas having reactor-relevant temperature in the DIII-D tokamak. The modes are excited over a wide spatial range from macroscopic device size to microturbulence size and the perturbation energy propagates across multiple spatial scales. The radial mode structure is able to expand from local to global in ∼0.1 ms and it causes magnetic topology changes in the plasma edge, which can lead to a minor disruption event. Since the mode is typically observed in the high ion temperature ≳10 keV and high-ß plasma regime, the manifestation of the mode in future reactors should be studied with development of mitigation strategies, if needed. This is the first observation of destabilization of the Alfvén continuum caused by the compressibility of ions with reactor-relevant ion temperature.

Achievement of Reactor-Relevant Performance in Negative Triangularity Shape in the DIII-D Tokamak.

Plasma discharges with a negative triangularity (δ=-0.4) shape have been created in the DIII-D tokamak with a significant normalized beta (ß_{N}=2.7) and confinement characteristic of the high confinement mode (H_{98y2}=1.2) despite the absence of an edge pressure pedestal and no edge localized modes (ELMs). These inner-wall-limited plasmas have a similar global performance as a positive triangularity (δ=+0.4) ELMing H-mode discharge with the same plasma current, elongation and cross sectional area. For cases both of dominant electron cyclotron heating with T_{e}/T_{i}>1 and dominant neutral beam injection heating with T_{e}/T_{i}=1, turbulent fluctuations over radii 0.5<ρ<0.9 were reduced by 10-50% in the negative triangularity shape compared to the matching positive triangularity shape, depending on the radius and conditions.

Quasistationary Plasma Predator-Prey System of Coupled Turbulence, Drive, and Sheared E×B Flow During High Performance DIII-D Tokamak Discharges.

A new, long-lived limit cycle oscillation (LCO) regime has been observed in the edge of near zero torque high performance DIII-D tokamak plasma discharges. These LCOs are localized and composed of density turbulence, gradient drives, and E×B velocity shear damping (E and B are the local radial electric and total magnetic fields). Density turbulence sequentially acts as a predator (via turbulence transport) of profile gradients and a prey (via shear suppression) to the E×B velocity shear. Reported here for the first time is a unique spatiotemporal variation of the local E×B velocity, which is found to be essential for the existence of this system. The LCO system is quasistationary, existing from 3 to 12 plasma energy confinement times (∼30-900 LCO cycles) limited by hardware constraints. This plasma system appears to contribute strongly to the edge transport in these high performance and transient-free plasmas, as evident from oscillations in transport relevant edge parameters at LCO time scale.

Modulation of Core Turbulent Density Fluctuations by Large-Scale Neoclassical Tearing Mode Islands in the DIII-D Tokamak.

We report the first observation of localized modulation of turbulent density fluctuations n[over ˜] (via beam emission spectroscopy) by neoclassical tearing modes (NTMs) in the core of the DIII-D tokamak. NTMs are important as they often lead to severe degradation of plasma confinement and disruptions in high-confinement fusion experiments. Magnetic islands associated with NTMs significantly modify the profiles and turbulence drives. In this experiment n[over ˜] was found to be modulated by 14% across the island. Gyrokinetic simulations suggest that n[over ˜] could be dominantly driven by the ion temperature gradient instability.

Evidence of Toroidally Localized Turbulence with Applied 3D Fields in the DIII-D Tokamak.

New evidence indicates that there is significant 3D variation in density fluctuations near the boundary of weakly 3D tokamak plasmas when resonant magnetic perturbations are applied to suppress transient edge instabilities. The increase in fluctuations is concomitant with an increase in the measured density gradient, suggesting that this toroidally localized gradient increase could be a mechanism for turbulence destabilization in localized flux tubes. Two-fluid magnetohydrodynamic simulations find that, although changes to the magnetic field topology are small, there is a significant 3D variation of the density gradient within the flux surfaces that is extended along field lines. This modeling agrees qualitatively with the measurements. The observed gradient and fluctuation asymmetries are proposed as a mechanism by which global profile gradients in the pedestal could be relaxed due to a local change in the 3D equilibrium. These processes may play an important role in pedestal and scrape-off layer transport in ITER and other future tokamak devices with small applied 3D fields.

Pedestal bifurcation and resonant field penetration at the threshold of edge-localized mode suppression in the DIII-D Tokamak.

Rapid bifurcations in the plasma response to slowly varying n=2 magnetic fields are observed as the plasma transitions into and out of edge-localized mode (ELM) suppression. The rapid transition to ELM suppression is characterized by an increase in the toroidal rotation and a reduction in the electron pressure gradient at the top of the pedestal that reduces the perpendicular electron flow there to near zero. These events occur simultaneously with an increase in the inner-wall magnetic response. These observations are consistent with strong resonant field penetration of n=2 fields at the onset of ELM suppression, based on extended MHD simulations using measured plasma profiles. Spontaneous transitions into (and out of) ELM suppression with a static applied n=2 field indicate competing mechanisms of screening and penetration of resonant fields near threshold conditions. Magnetic measurements reveal evidence for the unlocking and rotation of tearinglike structures as the plasma transitions out of ELM suppression.

A Q-band frequency tunable Doppler backscattering (DBS) system for pedestal and scrape-off layer density fluctuation and flow measurements in the DIII-D tokamak.

We present the design and laboratory tests for a new Q-band frequency tunable Doppler backscattering (DBS) system suitable for probing poloidal wavenumber kñ = 6-8 cm-1 density fluctuations and their flow velocities in the pedestal and scape-off layer (SOL) of the DIII-D tokamak. This system will provide new measurements in the increasingly important and under-diagnosed far pedestal and SOL plasma regions. These results are important for experimental transport studies and necessary for the validation of transport models, both of which are important to fusion energy research. The use of a single tunable frequency reduces the complexity and potential failure points as compared to a multichannel system. This new system utilizes a 33-50 GHz tunable source and will be integrated into the current V-band DBS in DIII-D using a broadband Q- and V-band multiplexer. A full-scale mockup of the quasi-optical system was used to test and optimize the performance. These tests include beam profile measurements at different distances (and angles) from a paraboloidal focusing and steering mirror. The measurements cover the full frequency range 33-75 GHz of the integrated/combined Q-V band DBS system and target a large radial coverage of the low-field side of the plasma from ρ = 1.1 to ρ = 0.5, where ρ is the normalized flux surface radial coordinate.

A new synthetic correlation electron cyclotron emission diagnostic for validating nonlinear gyrokinetic simulations of electron temperature turbulence.

To validate nonlinear gyrokinetic simulations of electron temperature turbulence, the experimental correlation electron cyclotron emission (CECE) measurements are to be compared using a synthetic CECE diagnostic, which generates modeled CECE measurement quantities by implementing realistic measurement parameters (e.g., spatial and wavenumber resolutions, radial location, etc.) to nonlinear gyrokinetic simulations. In this work, we calculate the radial and vertical spatial and wavenumber transfer functions, which are defined by the electron cyclotron emission emissivity radial profile and vertical probing antenna pattern, respectively. These transfer functions are applied to nonlinear gyrokinetic simulations of electron temperature turbulence using the continuum gyrokinetic code. A simultaneous comparison of the experimental electron temperature turbulence power spectrum and root-mean-square (RMS) level, as well as the radial correlation length with the new synthetic CECE diagnostic at a core location ρ ∼ 0.75 in an L-mode DIII-D tokamak plasma, is presented. The preliminary result shows that the synthetic CECE output underestimates the RMS level by ∼42% and overestimates the radial correlation length by ∼40%.

Using convolutional neural networks to detect edge localized modes in DIII-D from Doppler backscattering measurements.

In H-mode tokamak plasmas, the plasma is sometimes ejected beyond the edge transport barrier. These events are known as edge localized modes (ELMs). ELMs cause a loss of energy and damage the vessel walls. Understanding the physics of ELMs, and by extension, how to detect and mitigate them, is an important challenge. In this paper, we focus on two diagnostic methods-deuterium-alpha (Dα) spectroscopy and Doppler backscattering (DBS). The former detects ELMs by measuring Balmer alpha emission, while the latter uses microwave radiation to probe the plasma. DBS has the advantages of having a higher temporal resolution and robustness to damage. These advantages of DBS diagnostic may be beneficial for future operational tokamaks, and thus, data processing techniques for DBS should be developed in preparation. In sight of this, we explore the training of neural networks to detect ELMs from DBS data, using Dα data as the ground truth. With shots found in the DIII-D database, the model is trained to classify each time step based on the occurrence of an ELM event. The results are promising. When tested on shots similar to those used for training, the model is capable of consistently achieving a high f1-score of 0.93. This score is a performance metric for imbalanced datasets that ranges between 0 and 1. We evaluate the performance of our neural network on a variety of ELMs in different high confinement regimes (grassy ELM, RMP mitigated, and wide-pedestal), finding broad applicability. Beyond ELMs, our work demonstrates the wider feasibility of applying neural networks to data from DBS diagnostic.

Design of a low frequency, density profile reflectometer system for the MAST-U spherical tokamak.

Validated and accurate edge profiles (temperature, density, etc.) are vitally important to the Mega Ampere Spherical Tokamak Upgrade (MAST-U) divertor and confinement effort. Density profile reflectometry has the potential to significantly add to the measurement capabilities currently available on MAST-U (e.g., Thomson scattering and Langmuir probes). This work presents the diagnostic requirements, problems, and solutions facing profile reflectometry in spherical tokamaks and MAST-U in particular. Requirements include density measurements near zero electron density in the scrape off layer region, coverage for a broad range of MAST-U plasma parameters, high time (≤10 microseconds) and spatial resolutions (≤1 cm), reliability, and identification of the plasma start frequency.

New millimeter-wave diagnostics to locally probe internal density and magnetic field fluctuations in National Spherical Torus Experiment-Upgrade (invited).

A set of new millimeter-wave diagnostics will deliver unique measurement capabilities for National Spherical Torus Experiment-Upgrade to address a variety of plasma instabilities believed to be important in determining thermal and particle transport, such as micro-tearing, global Alfvén eigenmodes, kinetic ballooning, trapped electron, and electron temperature gradient modes. These diagnostics include a new integrated intermediate-k Doppler backscattering (DBS) and cross-polarization scattering (CPS) system (four channels, 82.5-87 GHz) to measure density and magnetic fluctuations, respectively. The system can access reasonably large normalized wavenumbers kθρs ranging from ≤0.5 to 15 (where ion sound gyroradius ρs = 1 cm and kθ is the binormal density turbulence wavenumber). The system addresses the challenges for making useful DBS/CPS measurements with a remote control of launch polarization (X- or O-mode), probed wavenumber, polarization match of the launch beam with the edge magnetic field pitch angle, and beam steering of the launched beam for wave-vector alignment. In addition, a low-k DBS system consisting of eight fixed frequencies (34-52 GHz) and four tunable frequencies (55-75 GHz) for low-k density turbulence and fast ion physics will be located at a nearby port location. The combined systems cover the near LCFS and pedestal regions (34-52 GHz), the pedestal or mid-radius (50-75 GHz), and core plasmas (82.5-87 GHz).

Observation of a critical gradient threshold for electron temperature fluctuations in the DIII-D Tokamak.

A critical gradient threshold has been observed for the first time in a systematic, controlled experiment for a locally measured turbulent quantity in the core of a confined high-temperature plasma. In an experiment in the DIII-D tokamak where L(T(e))(-1) = |∇T(e)|/T(e) and toroidal rotation were varied, long wavelength (k(θ)ρ(s) ≲ 0.4) electron temperature fluctuations exhibit a threshold in L(T(e))(-1): below, they change little; above, they steadily increase. The increase in δT(e)/T(e) is concurrent with increased electron heat flux and transport stiffness. Observations were insensitive to rotation. Accumulated evidence strongly enforces the identification of the experimentally observed threshold with ∇T(e)-driven trapped electron mode turbulence.

A novel Doppler backscattering (DBS) system to simultaneously measure radio frequency plasma fluctuations and low frequency turbulence.

A novel quadrature Doppler Backscattering (DBS) system has been developed and optimized for the E-band (60-90 GHz) frequency range using either O-mode or X-mode polarization in DIII-D plasmas. In general, DBS measures the amplitude of density fluctuations and their velocity in the lab frame. The system can simultaneously monitor both low-frequency turbulence (f < 10 MHz) and radiofrequency plasma density fluctuations over a selectable frequency range (20-500 MHz). Detection of high-frequency fluctuations has been demonstrated for low harmonics of the ion cyclotron frequency (e.g., 2fci ∼ 23 MHz) and externally driven high-frequency helicon waves (f = 476 MHz) using an adjustable frequency down conversion system. Importantly, this extends the application of DBS to a high-frequency spectral domain while maintaining important turbulence and flow measurement capabilities. This unique system has low phase noise, good temporal resolution (sub-millisecond), and excellent wavenumber coverage (kθ ∼ 1-20 cm-1 and kr ≲ 30 cm-1). As a demonstration, localized internal DIII-D plasma measurements are presented from turbulence (f ≤ 5 MHz), Alfvenic waves (f ∼ 6.5 MHz), ion cyclotron waves (f ≥ 20 MHz), as well as fluctuations around 476 MHz driven by an external high-power 476 MHz helicon wave antenna. In the future, helicon measurements will be used to validate GENRAY and AORSA modeling tools for prediction of helicon wave propagation, absorption, and current drive location for the newly installed helicon current drive system on DIII-D.

Role of zonal flow predator-prey oscillations in triggering the transition to H-mode confinement.

Direct evidence of zonal flow (ZF) predator-prey oscillations and the synergistic roles of ZF- and equilibrium E×B flow shear in triggering the low- to high-confinement (L- to H-mode) transition in the DIII-D tokamak is presented. Periodic turbulence suppression is first observed in a narrow layer at and just inside the separatrix when the shearing rate transiently exceeds the turbulence decorrelation rate. The final transition to H mode with sustained turbulence and transport reduction is controlled by equilibrium E×B shear due to the increasing ion pressure gradient.

Analysis method for calculating radial correlation length of electron temperature turbulence from correlation electron cyclotron emission radiometer.

The radial correlation length (Lr) is one of the essential quantities to measure in order to more fully characterize and understand turbulence and anomalous transport in magnetic fusion plasmas. The analysis method for calculating Lr of electron temperature (Te) turbulence from correlation electron cyclotron emission (correlation ECE or CECE) radiometer measurements has not been fully developed partly due to the fact that the turbulent electron temperature fluctuations are generally imbedded in much larger amplitude thermal noise, which leads to a greatly reduced cross correlation coefficient (ϱ) between two spatially separated ECE signals. This work finds that this ϱ reduction factor due to thermal noise is a function of the local relative temperature fluctuation power and CECE system bandwidths of intermediate and video frequencies, independent of radial separations. This indicates that under the approximation of constant relative temperature fluctuation power for a small radial range of local CECE measurements, the original shape of ϱ as a function of radial separation without thermal noise is preserved in the CECE data with thermal noise present. For Te turbulence with a Gaussian radial wavenumber spectrum, a fit function using the product of Gaussian and sinusoidal functions is derived for calculating Lr. This analysis method has been numerically tested using simulated ECE radiometer data over a range of parameters. Using this method, the experimental temperature turbulence correlation length Lr in a DIII-D L-mode plasma is found to be ∼10 times the local ion gyroradius.

Evaluation of a new DIII-D Doppler backscattering system for higher wavenumber measurement and signal enhancement.

The high density fluctuation poloidal wavenumber, kθ (kθ > 8 cm-1, kθρs > 5, ρs is the ion gyro radius using the ion sound velocity), measurement capability of a new Doppler backscattering (DBS) system at the DIII-D tokamak has been experimentally evaluated. In DBS, wavenumber (k) matching becomes more important at higher wavenumbers, owing to the exponential dependence of the measured signal loss factor on wave vector mismatch. Wave vector matching allows for the Bragg scattering condition to be satisfied, which minimizes the signal loss at higher k's. In the previous DBS system, without toroidal wave vector matching, the measured DBS signal-to-noise ratio at higher kθ (>8 cm-1) is substantially reduced, making it difficult to measure higher kθ turbulence. The new DBS system has been optimized to access higher wavenumber, kθ ≤ 20 cm-1, density turbulence measurement. The optimization hardware addresses fluctuation wave vector matching using toroidal steering of the launch mirror to produce a backscattered signal with improved intensity. The probe's sensitivity to high-k density fluctuations has been increased by approximately an order of magnitude compared to the old system that has been in use at DIII-D. Note that typical measurement locations are above or below the tokamak midplane on the low field side with normalized radial ranges of 0.5-1.0. The new DBS probe system with the toroidal matching of fluctuation wave vectors is thought to be critical to understanding high-k turbulent transport in fusion-relevant research at DIII-D.

Design elements and first data from a new Doppler backscattering system on the MAST-U spherical tokamak.

A new Doppler backscattering (DBS) system has been installed and tested on the MAST-U spherical tokamak. It utilizes eight simultaneous fixed frequency probe beams (32.5, 35, 37.5, 40, 42.5, 45, 47.5, and 50 GHz). These frequencies provide a range of radial positions from the edge plasma to the core depending on plasma conditions. The system utilizes a combination of novel features to provide remote control of the probed density wavenumber, the launched polarization (X vs O-mode), and the angle of the launched DBS to match the magnetic field pitch angle. The range of accessible density turbulence wavenumbers (kθ) is reasonably large with normalized wavenumbers kθρs ranging from ≤0.5 to 9 (ion sound gyroradius ρs = 1 cm). This wavenumber range is relevant to a variety of instabilities believed to be important in establishing plasma transport (e.g., ion temperature gradient, trapped electron, electron temperature gradient, micro-tearing, kinetic ballooning modes). The system is specifically designed to address the requirement of density fluctuation wavevector alignment which can significantly reduce the SNR if not accounted for.

Validating and optimizing mismatch tolerance of Doppler backscattering measurements with the beam model (invited).

We use the beam model of Doppler backscattering (DBS), which was previously derived from beam tracing and the reciprocity theorem, to shed light on mismatch attenuation. This attenuation of the backscattered signal occurs when the wavevector of the probe beam's electric field is not in the plane perpendicular to the magnetic field. Correcting for this effect is important for determining the amplitude of the actual density fluctuations. Previous preliminary comparisons between the model and Mega-Ampere Spherical Tokamak (MAST) plasmas were promising. In this work, we quantitatively account for this effect on DIII-D, a conventional tokamak. We compare the predicted and measured mismatch attenuation in various DIII-D, MAST, and MAST-U plasmas, showing that the beam model is applicable in a wide variety of situations. Finally, we performed a preliminary parameter sweep and found that the mismatch tolerance can be improved by optimizing the probe beam's width and curvature at launch. This is potentially a design consideration for new DBS systems.

New, improved analysis of correlation ECE data to accurately determine turbulent electron temperature spectra and magnitudes (invited).

Turbulent electron temperature fluctuation measurement using a correlation electron cyclotron emission (CECE) radiometer has become an important diagnostic for studying energy transport in fusion plasmas, and its use is widespread in tokamaks (DIII-D, ASDEX Upgrade, Alcator C-Mod, Tore Supra, EAST, TCV, HL-2A, etc.). The CECE diagnostic typically performs correlation analysis between two closely spaced (within the turbulent correlation length) ECE channels that are dominated by uncorrelated thermal noise emission. This allows electron temperature fluctuations embedded in the thermal noise to be revealed and fluctuation level and spectra determined. We have demonstrated a new, improved CECE coherency-based analysis for calculating the temperature fluctuation frequency spectrum and level, which has been verified both numerically through the simulation of synthetic ECE radiometer data and through analysis of experimental data from the CECE system on DIII-D. The new formulation places coherency-based analysis on a firm foundational footing and corrects some currently published methodologies. This new method accurately accounts for bias error in the coherence function and correctly calculates noise levels for a fixed data record length. It provides excellent accuracy in determining temperature fluctuation level (e.g., <10% error) even for a small realization number in the ensemble average. The method also has a smaller uncertainty (i.e., error bar) in the power spectrum when compared to the more standard cross-power method when evaluated at low coherency. Direct calculation of system noise level using correlation between randomized intermediate frequency signals is recommended.