Details of the neutral energy distribution and ionization source using spectrally resolved Balmer-alpha measurements on DIII-D

Spectrally resolved passive Balmer-α (D-α, H-α) measurements from the DIII-D 16 channel edge main-ion charge exchange recombination system confirm the presence of higher energy neutrals ("thermal" neutrals) in addition to the cold neutrals that recycle off the walls in the edge region of DIII-D plasmas. Charge exchange between thermal ions and edge neutrals transfers energy and momentum between the populations giving rise to thermal neutrals with energies approximating the ions in the pedestal region. Multiple charge exchange events in succession allow an electron to effectively take a random walk, transferring from ion to ion, providing a pathway of increasing energy and velocity, permitting a neutral to get deeper into the plasma before a final ionization event that contributes to the ion and electron particle fueling. Spectrally resolved measurements provide information about the density and velocity distribution of these neutrals, which has been historically valuable for validating Monte Carlo neutral models, which include the multi stage charge exchange dynamics. Here, a multi-channel set of such measurements is used to specifically isolate the details of the thermal neutrals that are responsible for fueling inside the pedestal top. Being able to separate the thermal from the cold emission overcomes several challenges associated with optical filter-based neutral density measurements. The neutral dynamics, deeper fueling by the thermal neutrals, and spectral measurement are modeled with the FIDASIM Monte Carlo collisional radiative code, which also produces synthetic spectra with a shape that is in close agreement with the measurements. By scaling the number of neutrals in the simulation to match the intensity of the thermal emission, we show it is possible to obtain local neutral densities and ionization source rates.

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. These results demonstrate that for the same plasma parameters, the Dα emission is 20%-30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n_H/(n_H + n_D)] by using the Hα radiance fraction [L_Hα/(L_Hα + L_Dα)] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

Charge exchange recombination spectroscopy measurements of DIII-D poloidal rotation with poloidal asymmetry in angular rotation

Sixteen new tangential views for the charge exchange recombination (CER) spectroscopy diagnostic at DIII-D were installed in 2019 on the high-field side (HFS) of the tokamak with the main goal being the measurement of main-ion (deuterium) poloidal rotation. Eight of the new views are connected to spectrometers, which view the main-ion spectrum, adding main-ion measurements where there were previously none, and another eight new views increased the spatial resolution of existing impurity (carbon) measurements on the HFS. When combined with the existing low-field side measurements, measurements at two locations on flux surfaces out to a normalized minor radius of ≈0.6 are possible. The new tangential views have been used to measure the deuterium poloidal rotation directly for the first time using the Poloidal Asymmetry in Angular Rotation (PAAR) method. These new measurements enable further testing of the validity of neoclassical poloidal rotation predictions. Separate measurements of the radial electric field can be made for an impurity ion and the main-ion by combining the PAAR measurements with additional CER measurements of toroidal rotation, temperature, and density. These independent measurements of the radial electric field agree reasonably well.

Formation of a High Pressure Staircase Pedestal with Suppressed Edge Localized Modes in the DIII-D Tokamak

We observe the formation of a high-pressure staircase pedestal (≈16-20  kPa) in the DIII-D tokamak when large amplitude edge localized modes are suppressed using resonant magnetic perturbations. The staircase pedestal is characterized by a flattening of the density and temperature profiles in midpedestal creating a two-step staircase pedestal structure correlated with the appearance of midpedestal broadband fluctuations. The pedestal oscillates between the staircase and single-step structure every 40-60 ms, correlated with oscillations in the heat and particle flux to the divertor. Gyrokinetic analysis using the cgyro code shows that when the heat and particle flux to the divertor decreases, the pedestal broadens and the E×B shear at the midpedestal decreases, triggering a transport bifurcation from the kinetic ballooning mode (KBM) to trapped electron mode (TEM) limited transport that flattens the density and temperature profiles at midpedestal and results in the formation of the staircase pedestal. As the heat flux to the divertor increases, the pedestal narrows and the E×B shear at the midpedestal increases, triggering a back transition from TEM to KBM limited transport. The pedestal pressure increases during the staircase phase, indicating that enhanced midpedestal turbulence can be beneficial for confinement.

Relative intensity calibration of the DIII-D charge-exchange recombination spectroscopy system using neutral beam injection into gas

A new calibration method for the DIII-D charge-exchange spectroscopy system produces a smoother impurity density profile compared to previous techniques, improving the accuracy of the impurity density profile reconstruction. The relative intensity calibration between the chords of the DIII-D charge-exchange recombination spectroscopy system is performed by firing neutral beams into the evacuated vacuum vessel pre-filled with neutral gas. Relative calibration is required in order to account for uncertainty in the 3D geometry of the neutral beam. Previous methods using helium gas have been improved by using xenon, which emits an emission line close to the commonly used carbon wavelength 5290.5 Å, as well as improved timing of the gas injection, inclusion of variations in the vessel pressure, and timing of neutral beam injection. Photoemission spectra recorded by 112 sightlines viewing 6 neutral beams are compared and used to form a relative calibration factor for each sightline. This relative calibration is shown to improve the quality of the measured ion density profile.

Active spectroscopy measurements of the deuterium temperature, rotation, and density from the core to scrape off layer on the DIII-D tokamak (invited)

Main-ion charge exchange recombination spectroscopy (MICER) uses the neutral beam induced D _α spectrum to measure the local deuterium ion (D⁺) temperature, rotation, and density, as well as parameters related to the neutral beams, fast ions, and magnetic field. An edge MICER system consisting of 16 densely packed chords was recently installed on DIII-D, extending the MICER technique from the core to the pedestal and steep gradient region of H-mode plasmas where the D⁺ and commonly measured impurity ion properties can differ significantly. A combination of iterative collisional radiative modeling techniques and greatly accelerated spectral fitting allowed the extension of this diagnostic technique to the plasma edge where the steep gradients introduce significant diagnostic challenges. The importance of including the fast ion D _α emission in the fit to the spectrum for the edge system is investigated showing that it is typically not important except for cases which can have significant fast ion fractions near the plasma edge such as QH-mode. Example profiles from an Ohmic L-mode and a high power ITER baseline case show large differences in the toroidal rotation of the two species near the separatrix including a strong co-current D⁺ edge rotation. The measurements and analysis demonstrate the state of the art in active spectroscopy and integrated modeling for diagnosing fusion plasmas and the importance of direct main ion measurements.

Synthetic diagnostic for assessing spatial averaging of charge exchange recombination spectroscopy measurements

A synthetic charge exchange recombination spectroscopy diagnostic based on the FIDASIM modeling suite has been created for the DIII-D tokamak. This synthetic diagnostic assumes that the ions have Maxwellian distribution functions on each flux surface and models emission from charge exchange events between the beam neutrals and a fully ionized impurity. This work was motivated by the observation of non-Gaussian spectra that may be caused by spatial averaging, atomic physics, or non-Maxwellian distribution functions. Measurements of non-Gaussian spectra commonly observed in the high confinement mode pedestal and in plasmas with large core gradients are compared to the synthetic diagnostic. Spatial averaging alone cannot account for the observations in these two cases, opening up the possibility of there being other causes such as non-Maxwellian distribution functions. The synthetic diagnostic has also been used to resolve a long-standing issue: it is shown that the lower temperatures measured by using vertical view chords relative to tangential view chords are due to increased spatial averaging for vertical views due to the DIII-D neutral beams being approximately twice as tall as they are wide.

Measurement of deuterium density profiles in the H-mode steep gradient region using charge exchange recombination spectroscopy on DIII-D

Recent completion of a thirty two channel main-ion (deuterium) charge exchange recombination spectroscopy (CER) diagnostic on the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] enables detailed comparisons between impurity and main-ion temperature, density, and toroidal rotation. In a H-mode DIII-D discharge, these new measurement capabilities are used to provide the deuterium density profile, demonstrate the importance of profile alignment between Thomson scattering and CER diagnostics, and aid in determining the electron temperature at the separatrix. Sixteen sightlines cover the core of the plasma and another sixteen are densely packed towards the plasma edge, providing high resolution measurements across the pedestal and steep gradient region in H-mode plasmas. Extracting useful physical quantities such as deuterium density is challenging due to multiple photoemission processes. These challenges are overcome using a detailed fitting model and by forward modeling the photoemission using the FIDASIM code, which implements a comprehensive collisional radiative model.

High resolution main-ion charge exchange spectroscopy in the DIII-D H-mode pedestal

A new high spatial resolution main-ion (deuterium) charge-exchange spectroscopy system covering the tokamak boundary region has been installed on the DIII-D tokamak. Sixteen new edge main-ion charge-exchange recombination sightlines have been combined with nineteen impurity sightlines in a tangentially viewing geometry on the DIII-D midplane with an interleaving design that achieves 8 mm inter-channel radial resolution for detailed profiles of main-ion temperature, velocity, charge-exchange emission, and neutral beam emission. At the plasma boundary, we find a strong enhancement of the main-ion toroidal velocity that exceeds the impurity velocity by a factor of two. The unique combination of experimentally measured main-ion and impurity profiles provides a powerful quasi-neutrality constraint for reconstruction of tokamak H-mode pedestals.

Improved edge charge exchange recombination spectroscopy in DIII-D

The charge exchange recombination spectroscopy diagnostic on the DIII-D tokamak has been upgraded with the addition of more high radial resolution view chords near the edge of the plasma (r/a > 0.8). The additional views are diagnosed with the same number of spectrometers by placing fiber optics side-by-side at the spectrometer entrance with a precise separation that avoids wavelength shifted crosstalk without the use of bandpass filters. The new views improve measurement of edge impurity parameters in steep gradient, H-mode plasmas with many different shapes. The number of edge view chords with 8 mm radial separation has increased from 16 to 38. New fused silica fibers have improved light throughput and clarify the observation of non-Gaussian spectra that suggest the ion distribution function can be non-Maxwellian in low collisionality plasmas.

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.

Observation of a multimode plasma response and its relationship to density pumpout and edge-localized mode suppression

Density pumpout and edge-localized mode (ELM) suppression by applied n=2 magnetic fields in low-collisionality DIII-D plasmas are shown to be correlated with the magnitude of the plasma response driven on the high-field side (HFS) of the magnetic axis but not the low-field side (LFS) midplane. These distinct responses are a direct measurement of a multimodal magnetic plasma response, with each structure preferentially excited by a different n=2 applied spectrum and preferentially detected on the LFS or HFS. Ideal and resistive magneto-hydrodynamic (MHD) calculations find that the LFS measurement is primarily sensitive to the excitation of stable kink modes, while the HFS measurement is primarily sensitive to resonant currents (whether fully shielding or partially penetrated). The resonant currents are themselves strongly modified by kink excitation, with the optimal applied field pitch for pumpout and ELM suppression significantly differing from equilibrium field alignment.

An upgrade of the magnetic diagnostic system of the DIII-D tokamak for non-axisymmetric measurements

The DIII-D tokamak magnetic diagnostic system [E. J. Strait, Rev. Sci. Instrum. 77, 023502 (2006)] has been upgraded to significantly expand the measurement of the plasma response to intrinsic and applied non-axisymmetric "3D" fields. The placement and design of 101 additional sensors allow resolution of toroidal mode numbers 1 ≤ n ≤ 3, and poloidal wavelengths smaller than MARS-F, IPEC, and VMEC magnetohydrodynamic model predictions. Small 3D perturbations, relative to the equilibrium field (10(-5) < δB/B0 < 10(-4)), require sub-millimeter fabrication and installation tolerances. This high precision is achieved using electrical discharge machined components, and alignment techniques employing rotary laser levels and a coordinate measurement machine. A 16-bit data acquisition system is used in conjunction with analog signal-processing to recover non-axisymmetric perturbations. Co-located radial and poloidal field measurements allow up to 14.2 cm spatial resolution of poloidal structures (plasma poloidal circumference is ~500 cm). The function of the new system is verified by comparing the rotating tearing mode structure, measured by 14 BP fluctuation sensors, with that measured by the upgraded B(R) saddle loop sensors after the mode locks to the vessel wall. The result is a nearly identical 2/1 helical eigenstructure in both cases.

Synchronous imaging of coherent plasma fluctuations

A new method for imaging high frequency plasma fluctuations is described. A phase locked loop and field programmable gate array are used to generate gating triggers for an intensified CCD camera. A reference signal from another diagnostic such as a magnetic probe ensures that the triggers are synchronous with the fluctuation being imaged. The synchronous imaging technique allows effective frame rates exceeding millions per second, good signal to noise through the accumulation of multiple exposures per frame, and produces high resolution images without generating excessive quantities of data. The technique can be used to image modes in the MHz range opening up the possibility of spectrally filtered high resolution imaging of MHD instabilities that produce sufficient light fluctuations. Some examples of projection images of plasma fluctuations on the H-1NF heliac obtained using this approach are presented here.

A multichannel magnetic probe system for analysing magnetic fluctuations in helical axis plasmas

The need to understand the structure of magnetic fluctuations in H-1NF heliac [S. Hamberger et al., Fusion Technol. 17, 123 (1990)] plasmas has motivated the installation of a sixteen former, tri-axis helical magnetic probe Mirnov array (HMA). The new array complements two existing poloidal Mirnov arrays by providing polarisation information, higher frequency response, and improved toroidal resolution. The helical placement is ideal for helical axis plasmas because it positions the array as close as possible to the plasma in regions of varying degrees of favourable curvature in the magnetohydrodynamic sense, but almost constant magnetic angle. This makes phase variation with probe position near linear, greatly simplifying the analysis of the data. Several of the issues involved in the design, installation, data analysis, and calibration of this unique array are presented including probe coil design, frequency response measurements, mode number identification, orientation calculations, and mapping probe coil positions to magnetic coordinates. Details of specially designed digitally programmable pre-amplifiers, which allow gains and filters to be changed as part of the data acquisition initialisation sequence and stored with the probe signals, are also presented. The low shear heliac geometry [R. Jiménez-Gómez et al., Nucl. Fusion 51, 033001 (2011)], flexibility of the H-1NF heliac, and wealth of information provided by the HMA create a unique opportunity for detailed study of Alfvén eigenmodes, which could be a serious issue for future fusion reactors.