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.

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.

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.

New methodology for measuring electron density perturbations caused by plasma coherent modes using profile reflectometry: Magnitudes and radial profiles in DIII-D.

New capabilities of fast-sweep frequency-modulated profile reflectometry are explored to measure electron density ne perturbation magnitudes and radial profiles due to plasma coherent modes in DIII-D. The first approach is based on the frequency analysis of phase perturbations associated with high frequency (∼MHz) Alfvén eigenmodes (AEs). The measurement of ∼5.5 MHz fast-ion-driven global Alfvén eigenmodes (GAEs) is demonstrated in a neutral beam-heated DIII-D plasma. The GAE induced a broad radial distribution of phase perturbations in the profile reflectometer data. Analysis of these data determined the effective cutoff location displacement and the estimated ne fluctuation profile. In the second approach, high resolution ne profiles are used directly to determine the radial structure of ne perturbations due to a neo-classical tearing mode. These new measurements broaden the application of profile reflectometry and advance the development of AE spectroscopy as a tool for non-invasive diagnosis of fast-ion-driven modes in DIII-D and burning plasmas such as ITER.

Ray-tracing analysis for cross-polarization scattering diagnostic on MAST-upgrade spherical tokamak.

A combined Doppler backscattering/cross-polarization scattering (DBS/CPS) system is being deployed on MAST-U for simultaneous measurements of local density turbulence, turbulence flows, and magnetic turbulence. In this design, CPS shares the probing beam with the DBS and uses a separate parallel-viewing receiver system. In this study, we utilize a modified GENRAY 3D ray-tracing code to simulate the propagation of the probing and scattered beams. The contributions of different scattering locations along the entire beam trajectories are considered, and the corresponding local B̃ wavenumbers are estimated using the wavevector matching criterion. The wavenumber ranges of the local B̃ that are detectable to the CPS system are explored for simulated L- and H-mode plasmas.

Performance demonstration of vacuum microwave components critical for the operation of the ITER low-field side reflectometer.

Final design studies in preparation for manufacturing have been performed for functional components of the vacuum portion of the ITER Low-Field Side Reflectometer (LFSR). These components consist of an antenna array, electron cyclotron heating (ECH) protection mirrors, phase calibration mirrors, and vacuum windows. Evaluation of these components was conducted at the LFSR test facility and DIII-D. The antenna array consists of six corrugated-waveguide antennas for simultaneous profile, fluctuation, and Doppler measurements. A diffraction grating, incorporated into the plasma-facing miter bend, provides protection of sensitive components from stray ECH at 170 GHz. For in situ phase calibration of the LFSR profile reflectometer, an embossed mirror is incorporated into the adjacent miter bend. Measurements of the radiated beam profile indicate that these components have a small, acceptable effect on mode conversion and beam quality. Baseline transmission characteristics of the dual-disk vacuum window are obtained and are used to guide ongoing developments. Preliminary simulations indicate that a surface-relief structure on the window surfaces can greatly improve transmission. The workability of real-time phase measurements was demonstrated on the DIII-D profile reflectometer. The new automated real-time analysis agrees well with the standard post-processing routine.

Turbulence measurements on the high and low magnetic field side of the DIII-D tokamak.

In this paper, we address the challenging question of measuring turbulence levels on the high magnetic field side (HFS) of tokamak plasmas. Although turbulence measurements on the HFS can provide a stringent constraint for the turbulence model validation, to date only low magnetic field side (LFS) measured turbulence has been used in validation studies. To address this issue, an eight channel Correlation Electron Cyclotron Emission (CECE) system at DIII-D was modified to probe both LFS and HFS. In contrast to the second harmonic extraordinary mode electron cyclotron resonance emission that is typically used in CECE, we show that it is possible to probe the HFS using fundamental O-mode electron cyclotron resonance emission. The required hardware modifications for the HFS measurements are presented here, and the potential issues in this measurement are discussed.

Optimized quasi-optical cross-polarization scattering system for the measurement of magnetic turbulence on the DIII-D tokamak.

Simulations and laboratory tests are used to design and optimize a quasi-optical system for cross-polarization scattering (CPS) measurements of magnetic turbulence on the DIII-D tokamak. The CPS technique uses a process where magnetic turbulence scatters electromagnetic radiation into the perpendicular polarization enabling a local measurement of the perturbing magnetic fluctuations. This is a challenging measurement that addresses the contribution of magnetic turbulence to anomalous thermal transport in fusion research relevant plasmas. The goal of the new quasi-optical design is to demonstrate the full spatial and wavenumber capabilities of the CPS diagnostic. The approach used consists of independently controlled and in vacuo aiming systems for the probe and scattered beams (55-75 GHz).

First step toward a synthetic diagnostic for magnetic fluctuation measurements using cross-polarization scattering on DIII-D.

Cross-polarization scattering (CPS) provides localized magnetic fluctuation ( B ̃ ) measurements in fusion plasmas based on the process where B ̃ scatters electromagnetic radiation into the orthogonal polarization. The CPS system on DIII-D utilizes the probe beam of a Doppler backscattering (DBS) diagnostic combined with a cross-view CPS receiver system, which allows simultaneous density and B ̃ measurements with good spatial and wavenumber coverage. The interpretation of the signals is challenging due to the complex plasma propagation of the DBS probe beam and CPS receive beams. A synthetic diagnostic for CPS is therefore essential to interpret data and perform detailed validation tests of non-linear turbulence simulations. This work reports a first step toward a synthetic diagnostic for CPS utilizing GENRAY, a 3-D ray tracing code, to simulate the propagation of the probe and scattered rays. The local B ̃ wavenumber is calculated from the local O- and X-mode wavenumbers using the wave vector matching scattering condition. The CPS wavenumber values and spatial locations are determined by a complex consideration that includes the local density and B ̃ level, receive antenna pattern and orientation, scattering volume, wavenumber values detected at the various scattering centers, and alignment of the magnetic wave vector with the plane perpendicular to the magnetic field. The issue of a spurious CPS signal due to polarization mismatches for launch and receive is also discussed. It is suggested that simultaneous O- and X-mode DBS measurements should be utilized for better understanding of the CPS signal contamination when the cutoff locations for both polarizations are close.

A free-standing wire scattering technique to monitor calibration variations of the DIII-D density profile reflectometer.

Real-time phase calibration of the ITER profile reflectometer is essential due to the long plasma duration and expected waveguide path length changes during a discharge. Progress has been recently made in addressing this issue by employing a phase calibration technique on DIII-D that monitors calibration variations that occur during each plasma discharge. By installing a thin free-standing metallic wire (1 mm diameter) near the end of the overmoded waveguide transmission system (oriented perpendicular to the waveguide axis), the round-trip phase shift from the wire is detected simultaneously with the plasma phase shifts. Variations in the reflectometer round trip path length (∼26 m) are then calculated after each DIII-D plasma discharge, allowing the calibration phase to be accurately monitored and updated. The round-trip reflectometer path length is observed to vary by ∼3 mm (root mean square value) during a typical DIII-D discharge. Using the variations in calibration phase, the density profile measurement accuracy can be improved. Since the wire retro-reflected power is ∼0.01 of the plasma signal, minimal effect is observed on the reflected signal from the plasma. Importantly, through a suitable choice in wire diameter, the calibration signal can be made approximately independent of the V-band reflectometer launch polarization. This is particularly important on DIII-D since orthogonal X- and O-mode polarized beams are coupled into the same transmission waveguide and launch antenna.

Millimeter-wave interferometry and far-forward scattering for density fluctuation measurements on LTX-ß.

The λ ≈ 1 mm (f = 288 GHz) interferometer for the Lithium Tokamak Experiment-ß (LTX-ß) will use a chirped-frequency source and a centerstack-mounted retro-reflector mirror to provide electron line density measurements along a single radial chord at the midplane. The interferometer is unique in the use of a single source (narrow-band chirped-frequency interferometry) and a single beam splitter for separating and recombining the probe and reference beams. The current work provides a documentation of the interferometry hardware and evaluates the capabilities of the system as a far-forward collective scattering diagnostic. As such, the current optical setup is estimated to have a detection range of 0.4 ≲ k ⊥ ≲ 1.7 cm-1, while an improved layout will extend the upper k ⊥ limit to ∼3 cm-1. Measurements with the diagnostic on LTX are presented, showing interferometry results and scattered signal data. These diagnostics are expected to provide routine measurements on LTX-ß for high frequency coherent density oscillations (e.g., Alfvénic modes during neutral beam injection) as well as for broadband turbulence.

Evaluation of low-frequency operational limit of proposed ITER low-field-side reflectometer waveguide run including miter bends.

The present design concept for the ITER low-field-side reflectometer transmission line (TL) consists of an ∼40 m long, 6.35 cm diameter helically corrugated waveguide (WG) together with ten 90° miter bends. This paper presents an evaluation of the TL performance at low frequencies (33-50 GHz) where the predicted HE11 mode ohmic and mode conversion losses start to increase significantly. Quasi-optical techniques were used to form a near Gaussian beam to efficiently couple radiation in this frequency range into the WG. It was observed that the output beams from the guide remained linearly polarized with cross-polarization power levels of ∼1.5%-3%. The polarization rotation due to the helical corrugations was in the range ∼1°-3°. The radiated beam power profiles typically show excellent Gaussian propagation characteristics at distances >20 cm from the final exit aperture. The round trip propagation loss was found to be ∼2.5 dB at 50 GHz and ∼6.5 dB at 35 GHz, showing an inverse increase with frequency. This was consistent with updated calculations of miter bend and ohmic losses. At low frequencies (33-50 GHz), the mode purity remained very good at the exit of the waveguide, and the losses are perfectly acceptable for operation in ITER. The primary challenge may come from the future addition of a Gaussian telescope and other filter components within the corrugated guide, which will likely introduce additional perturbations to the beam profile and an increase in mode-conversion loss.

A frequency-modulated continuous-wave reflectometer for the Lithium Tokamak Experiment.

The frequency-modulated continuous-wave reflectometer on LTX (Lithium Tokamak Experiment) and the data analysis methods used for determining electron density profiles are described. The diagnostic uses a frequency range of 13.1-33.5 GHz, for covering a density range of 0.21-1.4×1013 cm-3 (in O-mode polarization) with a time resolution down to 8 µs. The design of the diagnostic incorporates the concept of an "optimized" source frequency sweep, which minimizes the large variation in the intermediate frequency signal due to a long dispersive transmission line. The quality of the raw data is dictated by the tuning characteristics of the microwave sources, as well as the group delay ripple in the transmission lines, which can generate higher-order nonlinearities in the frequency sweep. Both effects are evaluated for our diagnostic and best practices are presented for minimizing "artifacts" generated in the signals. The quality of the reconstructed profiles is also improved using two additional data analysis methods. First, the reflectometer data are processed as a radar image, where clutter due to echoes from the wall and backscattering from density fluctuations can be easily identified and removed. Second, a weighed least-squares lamination algorithm POLAN (POLynomial ANalysis) is used to reconstruct the electron density profile. Examples of density profiles in LTX are presented, along with comparisons to measurements from the Thomson scattering and the λ = 1 mm interferometer diagnostics.

An eight-channel Doppler backscattering system in the experimental advanced superconducting tokamak.

Doppler backscattering system can measure the perpendicular velocity and fluctuation amplitude of the density turbulence with intermediate wavenumber. An eight-channel Doppler backscattering system has been installed in the Experimental Advanced Superconducting Tokamak (EAST), which can probe eight different radial locations simultaneously by launching eight fixed frequencies (55, 57.5, 60, 62.5, 67.5, 70, 72.5, 75 GHz) into plasma. The quasi-optical system consists of circular corrugated waveguide transmission, a fixed parabolic mirror, and a rotatable parabolic mirror which are integrated with quasi-optics front-end of the profile reflectometer inside the vacuum vessel. The incidence angle can be chosen from 5° to 12°, and the wavenumber range is 2-15/cm with the wavenumber resolution Δk/k≤0.21. Ray tracing simulations are used to calculate the scattering locations and the perpendicular wavenumber. The dynamic range of this new eight-channel Doppler backscattering system can be as large as 40 dB in the EAST. In this article, the hardware design, the ray tracing, and the preliminary experimental results in the EAST will be presented.