In reverse chronological order from the DIII-D weekly: (prep)
A paper titled “Collector Probe Analysis of Tungsten Transport to the Far-SOL from the DIII-D SAS-VW Divertor Experiment” was recently published in the journal Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2023.101566). In this article, experimental results from the SAS-VW campaign coupled with interpretive 3DLIM modelling show opposing trends for core impurity content when compared to impurity deposition on far-SOL Collector Probes (CPs) with increasing main ion density. An experiment was performed involving a series of upper-single-null L-mode discharges in each BT direction with progressively increasing main ion density approaching and slightly exceeding the divertor detachment threshold. The results indicate: a) increased radial W transport with decreasing peak divertor electron temperature; and b) negligible change in W content in the far-SOL at the outer mid-plane with the onset of divertor detachment. Using an interpretive modeling workflow for assessing the transport of W sputtered from the SAS-VW divertor, the analysis suggests that W migration to the main chamber surfaces during the campaign may also contribute to far-SOL deposition.
A paper titled “On the origin of the DIII-D L-H power threshold isotope effect” by K.J. Callahan, L. Schmitz, et al. was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/acf86c). This paper reports on the possible reason for the significantly increased low to high confinement mode (L to H-mode) power threshold P_LH in DIII-D low collisionality ITER similar hydrogen plasmas compared to deuterium. Gyro-fluid (TGLF) and gyro-kinetic (CGYRO) simulations of these experiments indicate the higher impurity (carbon) content in deuterium compared to hydrogen, due to mass-dependent sputtering of graphite, stabilizes ion temperature gradient (ITG) driven turbulence by increased main ion dilution. In the plasma edge, electron non-adiabaticity is also predicted to contribute to the isotope dependence of thermal transport, however its effect is subdominant compared to changes from impurity isotopic behavior. This L-H power threshold reduction with increasing carbon content at low collisionality is found to be in stark contrast to high collisionality results, where additional impurity content appears to increase the power necessary for H-mode access. This trend reversal with respect to impurity content is believed to be due to changing turbulence type, transitioning from an ITG/TEM dominated regime to a regime dominated by Resistive Ballooning Modes (RBM).
A paper titled “Model validation of tungsten erosion and redeposition properties using biased tungsten samples on DiMES” was recently published in the journal Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2023.101551). In this article, five identical samples were inserted into the lower divertor of DIII-D using the DiMES system. Each sample was exposed to a pair of L-mode attached plasma shots and was biased with a DC voltage ranging from -60 V to +25 V. Net erosion was measured with Rutherford Backscattering (RBS), while gross erosion was monitored in situ with the Auburn UV spectrometer (255.14 nm neutral W line). Numerical estimates of net erosion using the ERO2.0 PMI modeling code were in good agreement with measurements, lending confidence to our understanding of how modifying the Debye sheath changes W prompt redeposition. The modelled and measured trends in gross erosion were almost identical, but a scaling factor is required for quantitative agreement, suggesting that the S/XB coefficient for the W-I 255.14 nm line is underestimated.
A paper titled “The onset distribution of rotating m,n = 2,1 tearing modes and its consequences on the stability of high-confinement-mode plasmas in DIII-D” by L. Bardoczi, N. J. Richner and N. C. Logan was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ad0488). This paper reports the statistics of 2,1 NTM onsets of a multi-scenario database of over 13,000 DIII-D H-mode discharges. This analysis shows that the m,n=2,1 magnetic islands are dominantly pressure gradient driven, stochastically triggered non-linear instabilities at all edge safety factor (q95) values, in agreement with the modified Rutherford equation of Neoclassical Tearing Modes (NTMs). The instability onset time closely follows the exponential distribution in intermediate and high q95 scenarios and is characterized by near constant onset rate (λ), in accordance with Poisson-point processes. This implies that the plasmas are operated in marginally stable conditions, characterized by a small threshold for instability growth and variations in the trigger amplitude and/or the stabilizing mechanisms with temporally uniform random distribution in this database. While the majority of the tearing modes occur in the first current-profile relaxation time of the βN flattop, constant λ throughout the βN flattop shows that the tearing onset is insensitive to the evolution of the equilibrium current profile. In low q95 scenarios, where a large fraction of the plasmas are operated at low torque, λ increases over the course of the βN flattop, showing that these plasmas evolve toward more unstable conditions. Overall, the mode onset statistics support that the majority of 2,1 tearing modes are non-linear, neoclassically driven instabilities and classical stability does not play a dominant role in their onset through the current profile relaxation in DIII-D H-mode discharges.
A paper titled “SiC as a core-edge integrated wall solution in DIII-D” by S.A. Zamperini et al. was recently published in Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2023.101535). Silicon carbide (SiC) is a promising material for use in a fusion reactor due to its low hydrogenic diffusivity, high temperature strength and resilience under neutron irradiation. To assess SiC as a main wall material in DIII-D, simulations with TRIM.SP and DIVIMP are performed on a well-diagnosed L-mode discharge. The effective charge, Zeff, across the separatrix is used as a figure of merit in comparing SiC to the current graphite walls. It is found that SiC may reduce Zeff by as much as ∼50%. It is discussed how SiC may be expected to “self-condition” and create wall conditions similar to siliconization, further lowering Zeff due to efficient oxygen gettering (i.e., SiC walls may not require boronizations). This work directly motivates the installation of SiC toroidal limiters within DIII-D, which would be an intermediary step before considering a full-wall changeout.
A paper titled “Tungsten erosion and divertor leakage from the DIII-D SAS-VW tungsten-coated divertor in experiments with neon gas seeding” by M.S. Parsons and many DIII-D contributors was just published in Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2023.101520). In this paper, graphite “collector probes” located in the far Scrape-off-Layer at the outer midplane were used to measure the amount of tungsten escaping the SAS-VW divertor under different conditions with neon seeding. The “divertor leakage” was calculated by taking the ratio between the concentration of tungsten measured on the probe and the spectroscopically measured tungsten gross erosion rate from the divertor. The results show that both the tungsten divertor leakage and tungsten SOL flow pattern are sensitive to the neon seeding conditions.
A paper by Aaron Rosenthal et al., titled “Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D”, was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/acb95a). This paper uses the plasma and neutral density dynamics after an edge localized mode to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization source and neutral density. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport (as shown in figure). Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s−1 and diffusion coefficient of 0.05 m2 s−1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.
A paper titled “Self-consistent simulation of compressional Alfvén eigenmodes excitation and runaway electron diffusion in tokamak disruptions” by C. Liu, A. Lvovskiy, C. Paz-Soldan, S. C. Jardin and A. Bhattacharjee was published in Physical Review Letters (https://doi.org/10.1103/PhysRevLett.131.085102). In this paper the Alfvénic modes driven by high-energy runaway electrons (REs) beams observed in DIII-D disruption experiments were simulated using the kinetic-MHD mode M3D-C1, and both linear and nonlinear results were presented. The study reveals that the precession motion of trapped runaway electrons can exhibit resonance with the n=1 compressional Alfvén eigenmodes (CAEs), which typically have frequencies in the range of a few megahertz. This observation aligns with experimental measurements. The trapped runaway electrons are generated through collisions with high-Z impurities, primarily due to massive gas injection (MGI). In the nonlinear simulations, it is observed that multiple eigenmodes can be excited simultaneously by runaway electrons with different energy levels. The perturbed magnetic fields associated with these modes lead to rapid diffusion of runaway electrons, potentially preventing the formation of a post-disruption runaway electron plateau. These simulation findings validate earlier analyses based on experimental results from DIII-D, offering a novel potential approach to mitigating runaway electron beams during tokamak disruptions.
A paper titled “Empirical probability and machine learning analysis of m,n=2,1 tearing mode onset parameter dependence in DIII-D H-mode scenarios” by L. Bardoczi, N. J. Richner, J. Zhu, C. Rea and N. C. Logan was published in Physics of Plasmas (https://doi.org/10.1063/5.0165859). This paper reports analyses of the 2,1 NTM onsets in a multi-scenario DIII-D database of over 14,000 H-mode discharges. These analyses find that the normalized plasma beta, the rotation profile and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes, while absence of strong correlation with the current profile supports that classical effects are negligible. Interestingly, 2,1 tearing modes are most likely to destabilize when n>1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current & safety factor profile, separatrix shape and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration and reducing the core impurity radiation. In addition, lower triangularity, lower elongation and lower pedestal pressure may also help to improve the stability. The electron and ion temperature, collisionality, resistivity, internal inductance and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database. The formation of n>1 tearing modes and rotation profile flattening appear to be the best predictors of disruptive 2,1 NTM onsets.
The paper titled “Observation of long-radial-range-correlation in turbulence in high-collisionality high-confinement plasmas on DIII-D” by R. Hong (UCLA) et al. has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/acedc2). This research paper delves into the observation of spatially asymmetric turbulent structures in the core of high-collisionality H-mode plasmas on DIII-D tokamak. The key findings include: (1) Streamer-like transport events, with long-radial-range correlation (LRRC), develop in the core region at high collisionality and weak flow shear; (2) the amplitude and radial scale of streamers are regulated by the mean flow shearing rate and are inversely correlated with energy confinement time. These findings shed light on the degrading nature of core confinement at high collisionality on DIII-D tokamak.
The paper titled “Energetic electron transport in magnetic fields with island chains and stochastic regions” by E.G. Kostadinova et al. was published in the Journal of Plasma Physics (https://doi.org/10.1017/S0022377823000879). The paper investigates energetic electron transport in magnetized plasmas with magnetic fields characterized by island chains and regions of stochastic field lines produced by coil perturbations. The paper reports on DIII-D Frontier Science experiments that utilized electron cyclotron heating and current drive pulses to ‘tag’ electron populations within different locations across the discharge. The cross-field transport of these populations was then inferred from electron cyclotron emission measurements and gamma emission signals from scintillator detectors. Two types of energetic particles were distinguished and discussed: non-relativistic suprathermal electrons and relativistic runaway electrons. The magnetic field topology in each discharge was reconstructed with field-line tracing codes, which were also used to determine the location and scale of magnetic islands and stochastic regions. A comparison of simulations and experiments suggested that suprathermal transport was suppressed when the tagging was performed at a smaller radial location than the location of the q = 1 island chain and enhanced otherwise. The paper further demonstrated that increasing the width of the stochastic region within the edge plasma yielded enhancement of the suprathermal electron transport.
A paper by F. Effenberg et al., titled “In-situ coating of silicon-rich films on tokamak plasma-facing components with real-time Si material injection”, was published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/acee98). In this work, silicon-rich films were grown in situ on the lower divertor using real-time small Si pellet injection with the impurity granule injector (IGI) in DIII-D H-mode scenarios. Post-mortem analysis of samples exposed to the plasma through the divertor material evaluation system DiMES revealed Si-rich coatings with thicknesses of at least 0.4–1.2 nm. These layers grew at rates of 0.4–0.7 nm/s and efficiencies of 0.02–0.05 nm/mg on the outer DIII-D divertor. The coatings were found to be composed of silicon oxides, predominantly SiO2, and possibly including SiOC and SiOx, showing that the Si acted as an effective oxygen getter. Si erosion yields determined with the S/XB method showed good agreement with theoretical values for a-SiC and pure Si within error bars, taking surface dilution and preferential sputtering of non-silicon elements such as oxygen and carbon into account. Extrapolations based on the results in this study show that in an FPP, injecting a few hundred kg per year would be necessary to replenish and maintain millimeter-thick low-Z claddings, while less would be needed for occasional surface conditioning. Finer Si or SiC powder injection and heated divertor samples are proposed to optimize the method and promote SiC formation, a promising wall material candidate.
A paper titled “Direct Preemptive Stabilization of m,n=2,1 Neoclassical Tearing Modes by Electron Cyclotron Current Drive in the DIII-D Low-Torque ITER Baseline Scenario” by L. Bardoczi, R. J. La Haye, E. J. Strait, N. C. Logan, S. P. Smith, N. J. Richner and J. D. Callen was published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/acec5c). This paper reports the first preemptive direct stabilization of m,n=2,1 Neoclassical Tearing Modes (NTMs) in DIII-D low-torque ITER Baseline Scenario plasmas by electron cyclotron waves, strongly supporting that the observed stabilization is achieved by replacing the missing bootstrap current in the island by electron cyclotron current (ECCD). These plasmas (with and without ECCD) are stable to 2,1 NTMs when the differential rotation between the q=1 and q=2 surfaces is sustained above a critical value of about 1kHz, but those evolving to low differential rotation routinely develop a disruptive 2,1 island at the time and frequency of a sawtooth precursor. Preemptive, local and direct stabilization by ECCD was tested by scanning the current drive amplitude and location inside and outside of the q=2 rational surface shot by shot. Analysis of only low differential rotation time windows, i.e. when stabilization from the differential rotation is absent, shows focused ECCD at q=2 prevents the onset of 2,1 NTMs. These results give the first demonstration of disruptive NTM control by ECCD in the plasma scenario planned for ITER.
The paper titled “Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor” by R. Perillo was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/acdf02). This paper shows that type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ∼65% in the peak and ∼35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ∼20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ∼0.2 km s−1 in attached to ∼0.8 km s−1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filament fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.
A new high radial resolution 2D multichannel Charge eXchange Imaging (CXI) diagnostic has been developed at the University of Wisconsin-Madison and deployed at DIII-D to measure low-to-intermediate radial wavenumber carbon density fluctuations, as described in M. Major et al., Rev. Sci. Instruments. 93, 113503 (2022), https://doi.org/10.1063/5.0101844. CXI observes the n = 8 − 7 (λ = 529 nm) CVI emission line resulting from charge exchange collisions between heating neutral beam atoms and the intrinsic carbon ion density. The new CXI diagnostic is providing measurements with ΔR ∼ 0.4 cm to access higher kr instabilities (kr < 8 cm^−1) predicted to arise in the steep-gradient region of H-mode pedestals. The CXI system will feature 60 fiber bundles in a 12 × 5 arrangement (current deployment 12 radial x 2 poloidal), with each channel consisting of four 1 mm fibers. A custom optical system was designed to filter and image incoming signals onto an 8 × 8 avalanche photodiode array for high speed, high signal-to-noise measurements. Additionally, a novel electronics suite has been designed and commissioned to amplify and digitize the relatively low-intensity carbon signal at a 2 MHz bandwidth. Initial measurements in DIII-D H-mode pedestals demonstrate broadband fluctuations up to several hundred kHz at sub-cm spatial resolution.
The integration of a high-performance core and a dissipative divertor, or the so-called “core–edge integration,” has been widely identified as a critical gap in the design of future fusion reactors. A recent paper addressing this issue by H. Wang was published in Nuclear Fusion titled “Direct measurement of the electron turbulence-broadening edge transport barrier to facilitate core–edge integration in tokamak fusion plasmas” (https://doi.org/10.1088/1741-4326/acdfe2). This letter reports, for the first time, direct experimental evidence of electron turbulence at the DIII-D H-mode pedestal that correlates with the broadening of the pedestal and thus facilitates core–edge integration. In agreement with gyrokinetic simulations, this electron turbulence is enhanced by high ηe (ηe = Ln/LTe, where Ln is the density scale length and LTe is the electron temperature scale length), which is due to a strong shift between the density and temperature pedestal profiles associated with a closed divertor. The modeled turbulence drives significant heat transport with a lower pressure gradient that may broaden the pedestal to a greater degree than the empirical and theoretically predicted pedestal width scalings. Such a wide pedestal, coupled with a closed divertor, enables us to achieve a good core–edge scenario that integrates a high-temperature low-collisionality pedestal (pedestal top temperature Te,ped > 0.8 keV and a pedestal top collisionality ν*ped < 1) under detached divertor conditions. This paves a new path toward solving the core–edge integration issue in future fusion reactors.
The paper “Validation of density pump-out by pedestal-foot magnetic island formation prior to ELM suppression in KSTAR and DIII-D tokamaks” (Q. Hu et al) was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ace5c0). This work presents experimental validation on density pump-out caused by pedestal-foot magnetic island formation in both KSTAR and DIII-D tokamaks. Pedestal-foot magnetic island formation was proposed to be the responsible mechanism for RMP-caused density pump-out based on TM1 simulations. TM1 predicted several features to validate the mechanism: 1) bifurcated density pump-out by n=1 RMP, 2) sensitivity of density pump-out on q95, 3) the magnitude of density pump-out scale as the square root of RMP coil current. Dedicated experiments are carried out in KSTAR and DIII-D to validate these features and show that: 1) a staircase bifurcation in density pump-out is observed when slowly ramping up the n = 1 RMP current, and the density fluctuations are found to be slightly decreased in the pedestal region; 2) the magnitude of density pump-out becomes weaker when decreasing q95 from 5.5 to 4.9, and a partial recovery of density pump-out is observed when q95 is ramped down to lower than 4.9; 3) analysis of a DIII-D database of n = 3 RMP ELM control experiment finds that the magnitude of density pump-out is proportional to the square root of RMP coil current Δne/ne∝RMP^0.5 . These experimental observations are consistent with TM1 simulations.
It is important to understand how ELM mitigation techniques will affect the erosion rate of the divertor in future devices. A recent paper by A. Cacheris et al., titled “The effectiveness of D2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign” was published in Plasma Physics and Controlled Fusion (https://doi.org/10.1088/1361-6587/acd026). D2 pellets injected at the midplane of the DIII-D tokamak were used to mitigate W divertor erosion caused by type-I ELMs in H-mode plasmas during the Metal Rings Campaign. The overall trend, on average, is a reduction of W erosion per ELM with increasing inter-ELM density and increasing D2 pellet injection rate. Interpretive modeling using the free-streaming recycling model (FSRM) was conducted incorporating mixed-material effects. These results indicate that the interplay between an increasing divertor density and decreasing divertor temperature with increasing D2 injection may explain the similar intra-ELM W erosion at these D2 injection rates and inter-ELM densities.
A paper by Jie Chen et al., titled “Micro-tearing mode dominated electron heat transport in DIII-D H-mode pedestal”, was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/accafb). This paper presents a comprehensive set of evidence which reveals that micro-tearing modes (MTMs) dominate pedestal electron heat transport in an H-mode experiment in the DIII-D tokamak. The experiment investigates the role of MTMs by scanning pedestal collisionality, a main drive of MTM instability, from 0.43 to 0.84 on the pedestal top. Broadband (150-800 kHz) magnetic and density fluctuations originating from the pedestal gradient region and highly consistent with MTMs are observed, with amplitude increasing during the scan. The higher magnetic fluctuation amplitude correlates with a lower pedestal electron temperature gradient, implying MTMs may regulate the pedestal electron heat transport. The collisionality scan results in profile and transport changes consistent with predicted transport capability of MTMs: 1) experimentally-determined electron heat diffusivity increases ~40% at the location where the broadband density fluctuations peak; 2) ion heat diffusivity has less increase (<20%); and 3) a locally flattened region in the electron temperature pedestal is observed at high collisionality. A local, linear gyrokinetic simulation finds MTMs as the most unstable mode in the pedestal gradient region. In addition, local, nonlinear simulations suggest MTMs can dominate and drive experimentally-relevant, megawatt-level electron heat flux. This result establishes MTMs as an effective transport mechanism in the H-mode pedestal, in particular at high collisionality.
The rough surface, caused by sputtering and/or manufacturing process, at the divertor affects the incident ion angle and ion-shadowed area. The incident ion angle is a critical parameter for material erosion, and the ion-shadowed area can be a significant source of redeposited impurities. The paper “Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces”, by S. Abe et al. was recently published in Nuclear Fusion (https://doi.org/10.1063/5.0095155) and featured as a DOE Office of Science highlight (https://www.energy.gov/science/fes/articles/stuck-rough-how-aging-reactor-walls-may-exhibit-lower-erosion). This work presents a computational investigation of the dependence of material erosion on the incident ion angle at rough divertor surfaces. Ion angle distributions (IADs) for D plasmas at DIII-D and NSTX-U divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from DIII-D and NSTX-U experiments. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA). The mean surface inclination angle distribution (SIAD) was found to be a useful parameter to estimate the LIIA and ion-shadowed area from the calculated IADs rather than the root mean square roughness.
The paper “13C surface characterization of midplane and crown collector probes on DIII-D” by J. Duran et al. was recently published in Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2022.101339). Impurity transport outside of the confined plasma region is an important control feature for core impurity mitigation in fusion systems. This paper furthers the physics-based understanding of edge impurity transport through controlled injection of isotopically marked impurities followed by collection with physical probes. These ‘collector probes’ are positioned at multiple locations throughout the tokamak while the isotopic tracer impurities are in transit. The use of an enriched isotopic source of impurities and a multi-collector probe system has enabled the use of an isotopic mixing model that has identified a re-eroded wall source within DIII-D that has contributed to the collected impurities on the collector probe surfaces. These results lay the foundation for modelling efforts that includes accurate representation of the impurity sources which feed the collector probe measurements.
The paper titled “Validation of EDGDE2D-EIRENE predicted 2D distributions of electron temperature and density against divertor Thomson scattering measurements in the low-field side divertor leg in DIII-D“ by M. Groth et al. was published in Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2023.101372). This work shows that the physics models implemented in the edge fluid code EDGE2D-EIRENE predict the measured 2D electron temperature, density and pressure distributions in DIII-D low-confinement mode discharges using Divertor Thomson Scattering (DTS) within their collective uncertainties if the electron temperature in the divertor is 10 eV or higher. The simulations do not predict, however, the peaked electron temperature and density profiles measured adjacent to the target plate when Te,sep, is below 10 eV, i.e., for the plasma downstream from the region of ionization of deuterium atoms. Inclusion of cross-field drifts and a five-fold reduction of radial transport cannot reconcile the discrepancy between the measurements and predictions. Isolating the discrepancy to a region downstream of the ionization front, and identifying the shortcomings in the model likely being a combination of plasma transport and atomic and molecular physics, was made possible through DIII-D's unique, highly reproducible DTS measurements.
The paper “Simultaneous stabilization and control of the n=1 and n=2 resistive wall mode” was published in Nuclear Fusion (Battey et al 2023 Nucl. Fusion https://doi.org/10.1088/1741-4326/accd81) by postdoctoral researcher Alex Battey (Columbia Univ.). A novel controller and observer was demonstrated on DIII-D which was shown to be capable of stabilizing both the n=1 and n=2 resistive wall mode (RWM) up to values of beta_N approaching the n=2 ideal wall limit. This controller was also implemented with a flexible nonzero target which allows the mode to be controlled at a nonzero amplitude and frequency. This allows plasma response data to be obtained in close-loop during active feedback control for comparison with various stability models and to allow for more robust future controllers and observers.
An invited paper titled “First Demonstration of a Fiber Optic Bolometer on a Tokamak Plasma” by S. Lee et al. was published in Review of Scientific Instruments High Temperature Plasma Diagnostics Special Edition (https://doi.org/10.1063/5.0099546). A fiber optic bolometer (FOB) was demonstrated observing a fusion plasma for the first time at the DIII-D tokamak. A FOB uses a fiber optics-based interferometric technique that is designed to have a high sensitivity to temperature changes with a negligible susceptibility to electromagnetic interference (EMI) that can be problematic for resistive bolometers in a tokamak environment. A single-channel FOB was installed on DIII-D and showed a negligible increase in the noise level during typical plasma operations (0.39 mK) compared to the benchtop results (0.38 mK), confirming an insignificant EMI impact to the FOB. Comparisons to DIII-D resistive bolometers showed good agreement with the single-channel FOB, indicating that the FOB is comparable to a resistive bolometer when the impulse calibration is applied. The major potential effect of ionizing radiation on the FOB would be the radiation-induced attenuation, which can be efficiently compensated for by adjusting the probing light power.
The paper, “Low-frequency shear Alfvén waves at DIII-D: theoretical interpretation of experimental observations” by R.R. Ma et al. was recently published in Physics of Plasmas (https://doi.org/10.48550/arXiv.2301.01464). Using a theoretical framework called the “generalized fishbone-like dispersion relation,” the paper calculates the frequencies and growth rates of two DIII-D instabilities, the beta-induced Alfven eigenmode (BAE) and a lower frequency mode called the low-frequency Alfvén mode (LFAM). The calculations successfully reproduce the observed frequencies and the tendency of these two modes to be more unstable at rational values of qmin (see figure).
If the deuterium and tritium nuclei are spin polarized, the D-T fusion cross section is 50% larger. Theoretically, polarized nuclei will remain polarized in a hot plasma but this prediction has never been tested experimentally. Two recently published Nuclear Fusion papers, “Polarized fusion and potential in situ tests of fuel polarization survival in a tokamak plasma” by L. Baylor et al. [DOI 10.1088/1741-4326/acc3ae] and “Conceptual design of DIII-D experiments to diagnose the lifetime of spin polarized fuel” by A.V. Garcia et al. [DOI 10.1088/1741-4326/acaf0d] show that experimental measurements are now feasible in DIII-D. The experiments will use D-3He as a proxy for D-T. The polarized D and 3He fuel is delivered in LiD and 3He shell pellets, respectively. Measurements with an array of charged fusion product detectors that are sensitive to the differential fusion cross section detect the degree of polarization.
The paper titled “Modeling deep slot divertor concepts at DIII-D using SOLPS-ITER with drifts” by R. Maurizio et al. was published in Nuclear Materials and Engineering (https://doi.org/10.1016/j.nme.2022.101356). A staged divertor program is planned to advance DIII-D research on core-edge integration, and the 2nd stage dissipative divertor will address optimization of power and particle exhaust. This paper utilizes SOLPS-ITER simulations with full particle drifts enabled and focuses on the role of slot depth to achieve highly dissipative (detached) divertor conditions, in both Bt directions. With ion B×∇B into the divertor and PSOL = 4 MW, simulations reveal that increasing the depth of a divertor slot from 18 to 50 cm reduces the upstream separatrix electron density needed to detach by 15% primarily due to increased divertor radiation. The deep slot has lower target temperature compared to the shallow slot as a result of lower neutral leakage. Reversing the Bt direction cools and densifies the plasma in the slot, such that both slot options are detached at all simulated densities. Increasing the depth of a slot divertor is beneficial to achieve highly dissipative divertor conditions for both field directions. Additional modeling will build on these results to further inform the design of the Stage 2 divertor.
The paper “Design Study of an Edge Current Density Diagnostic Using New High-Performance Single-Channel Beam Emission Spectrometers at DIII-D” was published in Reviews of Scientific Instruments (https://doi.org/10.1063/5.0101781) by graduate student Ryan Albosta (U. Wisconsin). A new spectral Motional-Stark-Effect manifold diagnostic is being installed at DIII-D to resolve the edge-current density with up to 10 kHz time resolution. The system will simultaneously analyze the Doppler-shifted and Stark split MSE manifold and sigma/pi line ratios of the Balmer-Alpha and Balmer-Beta radiation resulting from collisional fluorescence of the heating deuterium neutral beams. This provides localized information on the magnetic field vector and, using multiple channels, allows inferring current density profiles. Each channel consists of a high-throughput 2-leg Czerny-Turner spectrometer using custom-made aspherical lenses and state-of-the-art EMCCD cameras that are fiber-coupled to the Beam Emission Spectroscopy collection optics viewing the 150-neutral beam line. FIDASIM modeling (see figure) has shown expected B_pol uncertainties of less than 1 mT, yielding a current density accuracy of 80 kA/m^2 at 1 cm radial resolution and spacing of adjacent channels.
The paper “Role of edge-localized neoclassical tearing modes in quiescent H-mode plasmas in the DIII-D tokamak” (Q. Hu et al) has been published in Physics of Plasmas (https://doi.org/10.1063/5.0137123). This work presents an alternative mechanism for the onset of edge-harmonic oscillations (EHO) in quiescent H-mode (QH) plasmas in DIII-D based on nonlinear two-fluid MHD simulations. Using kinetic equilibria constrained by edge profile measurements, the MHD simulations show that the n=1 NTM and its harmonics can be destabilized at the pedestal top of QH-mode plasma by the edge bootstrap current. The simulations further show that the unstable NTMs can saturate either at small or large island width depending on the magnitude of the edge bootstrap current. The onset of the EHO also results in a prompt decrease in the pedestal width and height, consistent with simulation results for the onset of the NTM at the top of the QH-mode pedestal. This suggests that the avoidance of edge-localized modes (ELMs) in QH-mode can be attributed to the enhanced local transport induced by the NTM that is sufficient to prevent the expansion of the pedestal to an unstable width, analogous to the mechanism explored for ELM suppression by resonant magnetic perturbations. Nonlinear MHD simulations scanning the E×B frequency and the ratio of parallel and perpendicular thermal diffusivity at the pedestal top show that edge-localized NTMs are destabilized for conditions of high E×B frequency, high pedestal temperature, and low pedestal density, qualitatively consistent with experimental conditions required for observing the EHO.
The paper “Vertical control of DIII-D discharges with strong negative triangularity” by A.O. Nelson et. al. was published in Plasma Physics and Controlled Fusion (https://doi.org/10.1088/1361-6587/acbe65). This work validates vertical stability calculations included in the TokSys model on Ohmic negative triangularity (NT) configurations on DIII-D, demonstrating pathways towards vertical control on DIII-D that were employed during the NT Armor Campaign. Notably, the paper also identifies a strong stabilizing contribution of a non-conformal wall in NT, suggesting that future NT machines can expand their potential operating space by carefully considering the shape of their stabilizing elements. This finding and the presented low-beta cases will be used to validate flexible MHD codes capable of optimizing the vertical stability problem for next-generation NT reactor designs.
A paper titled “Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal” by M.R. Halfmoon et al. was published in Physics of Plasmas (https://doi.org/10.1063/5.0102152). This paper reports results of gyrokinetic simulations with the GENE code to study pedestal fluctuations for DIII-D discharge 174082. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport (ETG) but at levels insufficient to satisfy power balance. It is shown that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. Several metrics are computed for ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel. This article was featured on the cover of POP's November 2022 issue.
A manuscript titled “Gyrokinetic Simulations Compared with Magnetic Fluctuations Diagnosed with a Faraday-Effect Radial Interferometer-Polarimeter in the DIII-D pedestal” by M. Curie et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac9b76). The work presents the first-time-ever direct comparison between nonlinear gyro kinetic simulations and internal magnetic fluctuation measurements. The shared characteristics of fluctuations such as non-monotonic collision dependence, frequency matching, (dB/B)/(dn_e/n_e) matching, show the microtearing modes (MTM) are present in the DIII-D pedestal. Transport from the simulations shows MTMs contributes a significant amount of electron heat flux. 
A paper entitled “Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma” by A.V. Dudkovskaia et al has been published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/aca48d). The work extends the drift magnetic island theory of [A V Dudkovskaia et al. Plasma Phys. Control. Fusion 63 (2021) 054001] to a realistic tokamak geometry characterized by finite beta values and D-shaped plasmas. The plasma triangularity effects on the neoclassical tearing mode (NTM) threshold are investigated for the first time within the self-consistent kinetic NTM theory of small magnetic islands. In particular, it is found that a higher triangularity plasma is more prone to NTMs, which is in agreement with the recent tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio is extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and successfully benchmarked against the EAST threshold island width measurements.
The paper “Misalignment of magnetic field in DIII-D assessed by post-mortem analysis of divertor targets” by R. Masline et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac9cf1). Measurements of carbon deposition on tungsten-coated lower divertor tile inserts revealed a shift in the toroidal magnetic field of ~5 mm towards the 310 degree toroidal direction. These findings are consistent with previous simulations of the 3D magnetic field distribution performed using the MAFOT field line integration code and Langmuir probe measurements in the small angle slot divertor. Since diagnostics are located at discrete toroidal and poloidal positions, misalignment in the toroidal field may complicate the mapping of these measurements to the same equilibrium. Comparisons with previous measurements of the toroidal field misalignment show that the error fields are gradually changing over time and therefore need to be periodically diagnosed.
A manuscript titled “E×B flow driven electron temperature bifurcation in a closed slot divertor with ion B×∇B away from X-point in the DIII-D tokamak” by X. Ma et al. was published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/ac9743). This work demonstrates a strong interplay between drifts and divertor geometry on divertor dissipation. With the strike point on the inner slanted surface of the SAS divertor, bifurcative transitions were observed near detachment onset resulting in a sharp decrease of Te both experimentally and computationally. SOLPS-ITER modeling with full drifts shows that the reversal of both the poloidal and radial E´B flows near the strike point leads to rapid density accumulation, leading to a step-like transition of divertor conditions resulting in cold plasma across the entire divertor target plate. These results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.
A manuscript titled “Improved Particle Confinement with Resonant Magnetic Perturbations in DIII-D Tokamak H-Mode Plasmas” by N. Logan et al. was published in Physical Review Letters (https://link.aps.org/doi/10.1103/PhysRevLett.129.205001). The work details a newly discovered regime in which the application of nonaxisymmetric (3D) field perturbations improves confinement and leads to a corresponding increase in the plasma density. This happens only in a range of low, counter-current plasma rotation and is opposite the otherwise ubiquitously observed density “pump-out” with the application of 3D fields. The new pump-in phenomenon reduces the frequency of edge localized modes (ELMs) if they are present and is may improve ELMless scenarios like low torque QH-modes that use 3D fields for error field correction and rotation control.
The paper “Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions” by T. Sizyuk and T. Abrams, recently published in the Journal of Nuclear Materials, presents a new modeling approach for the simulation of hydrogen isotope transport and retention in tungsten (https://doi.org/10.1016/j.jnucmat.2022.154095). This model was implemented in the ITMC-DYN simulation package and explains several interesting post-mortem observations from tungsten samples previously irradiated in DIII-D. The variation of temperatures during multi-discharge operation affects the spatial distribution of trapped D atoms. The effect of low material temperature during the cooldown phase, together with high D concentrations, enhances defect formation and growth in the subsurface layers of the D/W material system. Measured thermal desorption spectra of both pristine and pre-damaged tungsten were well-explained by the newly developed model.
The paper “Exploring data-driven models for spatiotemporally local classification of Alfvén eigenmodes” by A. A. Kaptanoglu et al., published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac8a03), uses a manually labeled database of 26 discharges from the DIII-D tokamak in order to train simple neural-network-based models for classifying Alfvén eigenmodes (AEs). The models provide spatiotemporally local identification of four types of AEs by using an array of 40 electron cyclotron emission (ECE) signals as inputs. Despite the minimal dataset, this strategy performs well at spatiotemporally localized classification of AEs, indicating future opportunities for more sophisticated models and incorporation into real-time control strategies. The trained model is then used to generate spatiotemporally-resolved labels for each of the 40 ECE measurements on a much larger database of 1112 DIII-D discharges. This large set of precision labels can be used in future studies for advanced deep predictors and new physical insights.
The paper “Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer” by J. Gonzalez-Martin et al., was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac7406). It discusses the application of different models to reproduce the fast-ion flow induced by Alfvén Eigenmodes that is measured by an imaging neutral particle analyzer (INPA). Ad hoc energetic particle diffusivity modeling by TRANSP significantly deviates from the observations. Comparably, reduced modeling, i.e. a combination of NOVA-K and ASCOT5 codes, reproduce some key features of the observed flow, but largely fail to interpret the observed fast ion depletion near the plasma axis. At last, hybrid simulations predict an RSAE consistent with the experiment that redistributes the injected ions. The resulting synthetic INPA images are in good agreement with the measurement near the injection energy. Nonperturbative effects are required to reproduce that depletion of fast ions near the magnetic axis at the injection energy.
The paper “AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application” by S. Munaretto et al., was recently published in Review of Scientific Instruments (https://doi.org/10.1063/5.0101630). It discusses a time domain algorithm that can be used in real time to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) from the 3D magnetic diagnostic signals. In this algorithm, the signal generated by the coil-induced eddy current is computed as the sum of low-pass filtered time derivatives of the coil current in each coil. While the gain for each low pass filtered signal is fitted, the choice of the low pass filters’ time constant is found to be more complex, since it impacts the uncertainty of the results, i.e. the applicability of the compensation to plasma shots. The code is available both in NSTX-U and DIII-D PCS.
The paper “Modeling the effect of nitrogen recycling on the erosion and leakage of tungsten impurities from the SAS-VW divertor in DIII-D during nitrogen gas injection” by M.S. Parsons et al. has been published in Nuclear Materials and Energy (https://doi.org/10.1016/j.nme.2022.101254). A modeling workflow with the SOLPS-ITER and DIVIMP codes was used to study tungsten behavior in the plasma boundary with a variety of nitrogen injection rates with different assumptions about the nitrogen recycling at the tungsten divertor target. This paper demonstrated that an optimal injection rate may exist in some scenarios which reduces the amount of tungsten leaking out of the divertor and reaching the core (see figure below). However, when the nitrogen recycling rate at the divertor targets is high, the nitrogen redistributes within the SAS-VW slot leading to increased tungsten sputtering, and the range of nitrogen injection rates resulting in tungsten mitigation becomes narrower.
The paper “Computational assessment on impurity sourcing and transport using high-temperature graphite and silicon carbide plasma-facing walls in a tokamak environment” by G. Sinclair et al., published in Fusion Science and Technology (https://doi.org/10.1080/15361055.2022.2099506), utilizes predictive impurity transport simulations to assess the impact of wall temperature on erosion and core contamination in a DIII-D-sized device. DIVIMP impurity tracking simulations conducted on a plasma background generated using SOLPS-ITER showed that the total gross erosion of C from graphite walls could be reduced by a factor of 2 by operating at a wall temperature of 1200 K instead of 400 K and could be further reduced by a factor of 8 by using SiC instead of graphite. The C density at the outer midplane separatrix, which is used as a proxy for core contamination, is 7.5 times lower with a SiC wall than with a graphite wall at a temperature of 800 K. Modeling results motivate near-term upgrades to graphite-clad machines to demonstrate more reactor-relevant wall solutions.
The paper “A method to identify the heat flux from photons and neutrals at the divertor target” by J. Ren et al., was recently published in IEEE Transactions on Plasma Science (https://ieeexplore.ieee.org/document/9831056). For the first time in a tokamak, the fraction of heat flux striking the target surface made up of neutrals and photons relative to that from ions and electrons has been measured. The Surface Eroding Thermocouple (SETC) is an innovative diagnostic designed for fast (~1 ms response time) surface heat flux measurement and has been successfully applied in DIII-D to investigate divertor detachment. Employing the DiMES manipulator, a novel method has been developed of using a combination of surface-facing and recessed SETCs to distinguish between the heat flux from charged particles, and that from neutrals and radiative heating. This technique suggests that at least 20% of the total incident heat flux is attributable to the latter sources that are not trapped in the magnetic field in a fully detached open divertor condition. This not only indicates the non-charged particle heating significantly contributes to the total surface heat flux in detached divertor conditions, but also demonstrates the feasibility of using two operational modes of SETCs to distinguish between the two surface heating mechanisms. Following application in the lower divertor, a combination of both flush and recessed SETCs were installed in the SAS-VW divertor with results expected following the FY22 campaign.
The paper “Local measurements of the pedestal magnetic field profile throughout the ELM cycle on DIII-D” was published in Physics of Plasmas by M. Galen Burke et al. describing the first measurements of the H-mode pedestal magnetic field evolution through the ELM cycle and L-H transition in a high temperature tokamak. The measurement utilizes the separation of the Stark-split neutral beam emission to extract the local magnetic field perpendicular to the beam neutral velocity, B. A new ultra-high optical throughput spectroscopy technique was validated by first measurements that demonstrated excellent agreement with the relative changes in the local magnetic field induced by core current profile changes. Measurements of peeling-limited ELMs show spatially complex and extremely fast (<50 μs) changes in the pedestal magnetic field. These changes imply a fast rearrangement of the edge current peak at the ELM crash, as expected. In some discharges changes in B are observed throughout the ELM cycle. In others, at the edge current peak, B recovers and is relatively stable until 1-2 ms leading up to the next crash. These measurements are made possible by a novel spatial heterodyne spectrometer that was developed under a measurement innovation grant through the University of Wisconsin-Madison.
The paper “Prospects for H-mode inhibition in negative triangularity tokamak reactor plasmas” was published in Nuclear Fusion by Oak Nelson describing access conditions (and the lack thereof) to H-mode in negative triangularity reactor configurations. Using the OMFIT framework, over 9000 equilibria with varying shapes and kinetic profiles were generated and assessed in terms of their stability to infinite-n ballooning modes. Stability calculations indicate that infinite-n ballooning modes are strongly destabilized by changes in the local magnetic shear at strongly negative triangularity, closing access to the second stability region and preventing the growth of an H-mode pedestal. These results indicate that negative triangularity reactors should be able to robustly maintain L-mode operation.
The paper “Faraday-effect polarimetry for current profile measurement in the tokamak plasma edge” by J. Chen et. al. has been published in Review of Scientific Instruments (https://aip.scitation.org/doi/10.1063/5.0101550). The paper shows that a Faraday-effect polarimeter can measure current density profile in DIII-D H-mode pedestal. A low-field-side, vertical-view polarimeter design, covering R=2.15 to 2.27 m on DIII-D with sub-cm chord spacing, is assessed. A synthetic diagnostic simulation, in which current profile is derived from the Faraday rotation angle using Abel inversion, is performed to evaluate impact of various factors, such as diagnostic performance and flux surface quality from kinetic and routine EFIT calculations. It is found that the pedestal current profile can be realistically determined with no greater than 0.12 MA/m^2 uncertainty, or about 10% of the peak pedestal current density, in an investigated DIII-D H-mode plasma.
The paper “Mitigation of plasma-wall interactions with low-Z powders in DIII-D high confinement plasmas” by F. Effenberg et al was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac899d). Experiments with low-Z powder injection in DIII-D high confinement discharges demonstrated increased divertor dissipation and detachment while maintaining good core energy confinement. Lithium (Li), boron (B), and boron nitride (BN) powders were injected in H‑mode plasmas into the upper small-angle slot (SAS) divertor. Multi-species BN powders showed the most substantial increase in divertor neutral compression by more than an order of magnitude and lasting detachment. The application of powders also showed a substantial improvement of wall conditions manifesting in reduced wall fueling source and intrinsic carbon and oxygen content in response to the cumulative injection of non-recycling materials. The results suggest that low-Z powder injection, including mixed element compounds, is a promising new core-edge compatible technique that simultaneously enables divertor detachment and improves wall conditions during high confinement operation.
The paper “Identifying the microtearing modes in the pedestal of DIII-D H-modes using gyrokinetic simulations” by Ehab Hassan et al. has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3be5). This paper reports a study of the instabilities in the pedestal region carried out using gyrokinetic simulations with the GENE code to model an ELMy H-mode DIII-D discharge. The simulations produce MTMs, identified by predominantly by electromagnetic heat flux, small particle flux, and a substantial degree of tearing parity. The magnetic spectrogram from Mirnov coils exhibits three distinct frequency bands—two narrow bands at lower frequency (∼35–55 kHz and ∼70–105 kHz) and a broader band at higher frequency (∼300–500 kHz). Global linear GENE simulations produce MTMs that are centered at the peak of the ω∗ profile and correspond closely with the bands in the spectrogram. The three distinctive frequency bands can be understood from the basic physical mechanisms underlying the instabilities. For example (i) instability of certain toroidal mode numbers (n) is controlled by the alignment of their rational surfaces with the peak in the ω∗ profile, and (ii) MTM instabilities in the lower n bands are the conventional collisional slab MTM, whereas the higher n band depends on curvature drive. While many features of the modes can be captured with the local approximation, a global treatment, as used here, is necessary to quantitatively reproduce the detailed band gaps of the low-n fluctuations.
The paper “Ion thermal transport in the H-mode edge transport barrier on DIII-D” by S. R. Haskey et al. has been published in Physics of Plasmas (https://doi.org/10.1063/5.0072155). The paper describes how the power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ( νi*) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of D+ and the more commonly measured C6+ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low νi* discharges at the radius where the non-linear CGYRO simulations were run. Taken together, the increased fluctuations, power balance calculations, and gyrokinetic simulations show that the ion heat fluxes are in excess of the neoclassical levels, particularly at lower νi*, with the additional fluxes likely caused by weakly suppressed ion scale electrostatic turbulence.
A novel multi-code workflow to interpret collector probe deposition patterns in DIII-D was recently published in a Fusion Science and Technology article entitled “Separatrix-to-Wall Simulations of Impurity Transport with a Fully Three-Dimensional Wall in DIII-D” (https://doi.org/10.1080/15361055.2022.2082791) by S. Zamperini et al. The components of the workflow consist of a detailed computer-aided-design file of the DIII-D vessel wall coupled to the scrape-off-layer (SOL) codes MAFOT, OSM, DIVIMP, and 3DLIM. A special-purpose toolkit enables passing the output of these codes among each other to provide a full-SOL picture of impurity transport. A demonstration of the workflow is described to support evidence of near-SOL tungsten parallel accumulation during trace W impurity experiments on DIII-D. Iteration between simulated deposition patterns in 3DLIM and DIVIMP predicts a region of elevated W density near the separatrix about halfway between the outboard midplane and the top of the plasma. This workflow will be used to better interpret collector probe experiments on DIII-D during the ongoing SAS-VW campaign.
The paper “Influence of triangularity on the plasma response to resonant magnetic perturbations” by S. Gu et al. was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac5f7b). Measurements and modeling from recent RMP-ELM control experiments in DIII-D and other devices find a peak in the resonant plasma response as a function of upper triangularity which corresponds to the most experimentally accessible ELM-suppressed state. As shown in the figure, the plasma response is strongly suppressed at high triangularity compared to that at lower triangularity, while keeping other equilibrium quantities largely unchanged. This result is demonstrated in DIII-D (left), ASDEX Upgrade (middle), and EAST (right) tokamaks, and helps understand the inaccessibility of ELM suppression at high triangularity. Both experiments and modeling find a reduced magnetic plasma response on the high-field side at high triangularity across devices. The amplitude of the dominant mode reveals similar trends with the edge resonance and radial displacement near X-point, which suggests that the multi-mode plasma response provides another way to understand the ELM control physics. These findings indicate that the plasma shape should be taken into consideration when designing an RMP-ELM control strategy in experiment, and that predictive plasma response calculations can be used to maximize access to RMP-ELM control in future devices by maximizing the coupling between coils and the plasma.
The paper “Quantifying heat and particle flux to primary and secondary divertors for various types of edge-localized-modes” by R. Perillo et al. was recently published in Physics of Plasmas and was chosen as Editor’s Pick (https://doi.org/10.1063/5.0092924). This article shows that between 15% and 30% of edge-localized-mode (ELM) heat flux can be deposited to regions outside the main divertor including the far scrape-off-layer (SOL), private flux region (PFR), and secondary divertor inner target. Significantly, the peak heat flux values at the secondary inner divertor can surpass those at the secondary outer target, despite the secondary inner divertor being magnetically isolated from the outer divertor. This behavior is consistently observed for pedestal collisionalities from nu_star,e ∼ 1.5 to ∼ 3.9. Heat flux profiles of the examined ELM types feature rippled structures due to bursts in the outer far-SOL region but not on the secondary inner target, causing long decay lengths in the time-averaged ELM profiles. The contribution of each ELM type to the total time-averaged power deposited to the secondary divertor has been evaluated, showing that grassy ELMs contribute ∼8%, small ELMs ∼67%, and type-I ELMs ∼85%. These findings imply that small ELMs may yet pose a concern for future machines if some regions of the main wall are not designed to withstand significant heat and particle fluxes. Due to the low intra-ELM heat flux contribution, however, the grassy ELM regime is an attractive option for an ELMing scenario in future machines.
The paper “A survey of pedestal magnetic fluctuations using gyrokinetics and a global reduced model for microtearing stability” has recently been published by Max Curie et al in Physics of Plasmas (https://doi.org/10.1063/5.0084842). This article presents a global reduced model for slab-like microtearing modes (MTMs) in the H-mode pedestal, which reproduces distinctive features of experimentally observed magnetic fluctuations, such as discrete frequency bands at noncontiguous mode numbers. The key insight underlying this model is that MTM instability is enabled by the alignment of a rational surface with the peak in the profile of the electron diamagnetic frequency. Conversely, MTMs are strongly stabilized for toroidal mode numbers for which these quantities are misaligned. This property explains the discrete bands in magnetic fluctuation spectra in several DIII-D and JET discharges, which are predicted by the reduced model as well as by global gyrokinetic simulations. As an example, the left figure shows a frequency spectrogram of magnetic fluctuations from a DIII-D discharge with two frequency bands that are candidates for MTMs. The right figure shows output of the SLiM model, which predicts that MTM turbulence at n=3 causes the lower frequency band and turbulence at n=5 and 6 is responsible for the higher frequency band. This fast yet accurate reduced model for MTMs enables rapid interpretation of magnetic fluctuation data from a wide range of experimental conditions to help assess the role of MTM in the pedestal.
The paper “Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D” by Z. Li was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac4acf). This paper reports the recent progress on linear and nonlinear simulations using BOUT++ reduced two fluid MHD model to investigate the broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling-ballooning (PB) mode located near the minimum of Er well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave (DAW) locates in the upper-pedestal propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitatively agreement is likely due to the lack of trapped electron physics in applied fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.
The paper “Deconvolving the Roles of E×B Shear and Pedestal Structure in the Energy Confinement Quality of Super H-Mode Experiments” by A.M. Garofalo et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac4d63). This paper presents results from recent DIII-D experiments investigating the energy confinement quality of super H-mode plasmas and shows that the contribution from high rotation in the E×B shear, not the high pedestal pressure, is responsible for confinement quality in excess of H98y2 ~1.2. The experimental results confirmed quantitatively the predict-first transport modeling results published in a previous paper [S. Ding, et al., Nucl. Fusion 60 (2020) 034001]. By reducing the NBI torque and fixed power, the energy confinement quality during the stationary phase of super H-mode experiments was reduced from H98y2 ~1.6 to H98y2 ~1.2, as predicted. In addition, the experiments showed that the same exceptionally high confinement quality (H98y2 >2) can be achieved at high NBI torque (and nearly identical rotation profiles) with either high or low plasma triangularity, and a 30% difference in the pedestal pressure. The paper also provides new insight into results of previous predictions [W.M. Solomon, et al., Phys. Plasmas 23 (2016) 056105] of ITER Baseline Scenario performance assuming a super H-mode pedestal. Consistent with the experimental results in this paper, no increase of the ITER energy confinement quality is predicted with increasing pedestal pressure, although an improvement in the fusion gain is predicted with the higher core density associated with higher pedestal density.
The magnetic pre-sheath (MPS) width is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In the paper “Determination of the characteristic magnetic pre-sheath length at divertor surfaces using micro-engineered targets on DiMES at DIII-D” recently published by S. Abe et al., the MPS width was determined by measuring the azimuthal and polar D ion incident angle using micro-trench samples on DiMES (https://doi.org/10.1088/1741-4326/ac3cdb). Micro-trenches (30×30×4 µm) fabricated via focused ion beam (FIB) milling on a silicon surface were exposed to steady L-mode plasmas in DIII-D. The D ion incident angles were determined via the areal distribution of C impurity depositions measured by energy-dispersive X-ray spectroscopy (EDS). C deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80° (referenced to the surface normal direction) due to the ion shadowing effect. Those ion incident angles were reproduced by a Monte Carlo equation-of-motion model tracking the ion trajectories through the sheath.
In a recently (Feb 2022) published paper in Physics of Plasmas by E. Hollmann et al (doi https://aip.scitation.org/doi/full/10.1063/5.0080385), the first dynamic (time-dependent) measurements of impurity ion radial (cross-field) diffusion coefficients for post-disruption runaway electron plateaus are described. Small (~ 1 mm diameter) carbon or silicon pellets were fired into the edge of steady-state runaway electron (RE) plateaus created by argon pellet injection. The resulting radial transport of singly-charged impurity ions (C^+or Si^+) was monitored using passive spectroscopy with the DIII-D CER spectrometer system. Radial (cross-field) diffusion coefficients D_⊥≈2-5 m^2/s were obtained, about 10× larger than expected from neo-classical theory. It is speculated at present that these large diffusion coefficients may be due to turbulent transport. Even with these anomalous transport rates, the implications of this work are that impurity ion transport is too slow to mix into a vertically unstable RE plateau in ITER on the vertical drift timescale of ~100 ms, suggesting that RE plateau current dissipation by complete high-Z ion mixing into the RE plateau will be challenging to achieve.
The paper “Disruptive Neoclassical Tearing Mode Seeding in DIII-D with Implications for ITER” was published by R.J. La Haye et. al. in the May 2022 issue of Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac351f). New studies identify critical parameters and physics governing onset of a disruptive neoclassical tearing mode (NTM). A m/n=2/1 mode in DIII-D that grows robustly after a seeding event (edge localized mode ELM or sawtooth crash) is caused by the mode rotation dropping close to the plasma’s Er=0 rest frame; this opens the stabilizing ion-polarization current “gate” and destabilizes an otherwise marginally stable NTM. Our new experimental and theoretical insights and novel toroidal theory-based modelling are benchmarked and scalable to ITER and other future experiments. Mode rotation at q=2 is stabilizing (“gate closed”) except for MHD-induced transients that can “open the gate.” Extrapolating from DIII-D ITER baseline scenario (IBS) discharges, MHD transients are much more likely to destabilize problematic robustly growing 2/1 NTMs in ITER and other burning plasmas.
The paper “Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration” by L. Casali et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3e84). Impurity seeding experiments performed in the SAS graphite divertor at DIII-D using nitrogen and neon as seeded species show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. This work demonstrates the mechanism governing impurity leakage by studying the impurity transport in terms of the parallel momentum balance. Lower divertor enrichment is found for Ne compared to N consistent with its higher ionization potential and longer mean free path. Since carbon is an intrinsic radiator at DIII-D, the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling is also demonstrated. In the cases presented here, plasma drifts, flow reversal caused by high level of closure together with the effect of the thermal force contribute significantly to the shift of the impurity stagnation point leading to impurity leakage. Finally, this work highlights the importance of accompanying experimental studies with numerical modeling with particular attention to plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry.
The Nuclear Fusion Letter “Doubling off-axis electron cyclotron current drive efficiency via velocity space engineering,” by Xi Chen, et al. has been recently published (https://doi.org/10.1088/1741-4326/ac544a). For the first time, experiments on the DIII-D tokamak have demonstrated electron cyclotron current drive with more than double the conventional efficiency by tailoring the wave–particle interactions in velocity space using a novel ‘top launch’ geometry. Steering the EC waves to propagate nearly parallel to the resonance drives current more efficiently by (1) selective damping on electrons with higher parallel velocity v||, and (2) longer absorption path to compensate for inherently weaker absorption at higher v||. Experiments using a fixed-injection top launch system find an optimal velocity space interaction for maximum current drive efficiency at ρ ∼ 0.5 where the ease of drawing out a high v|| electron tail is balanced by sufficient absorption.
The paper “Deconvolving the Roles of E×B Shear and Pedestal Structure in the Energy Confinement Quality of Super H-Mode Experiments” by A.M. Garofalo et al. was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac4d63). This paper presents results from recent DIII-D experiments investigating the energy confinement quality of super H-mode plasmas and shows that the contribution from high rotation in the E×B shear, not the high pedestal pressure, is responsible for confinement quality in excess of H98y2 ~1.2. The experimental results confirmed quantitatively the predict-first transport modeling results published in a previous paper [S. Ding, et al., Nucl. Fusion 60 (2020) 034001]. By reducing the NBI torque and fixed power, the energy confinement quality during the stationary phase of super H-mode experiments was reduced from H98y2 ~1.6 to H98y2 ~1.2, as predicted. In addition, the experiments showed that the same exceptionally high confinement quality (H98y2 >2) can be achieved at high NBI torque (and nearly identical rotation profiles) with either high or low plasma triangularity, and a 30% difference in the pedestal pressure. The paper also provides new insight into results of previous predictions [W.M. Solomon, et al., Phys. Plasmas 23 (2016) 056105] of ITER Baseline Scenario performance assuming a super H-mode pedestal. Consistent with the experimental results in this paper, no increase of the ITER energy confinement quality is predicted with increasing pedestal pressure, although an improvement in the fusion gain is predicted with the higher core density associated with higher pedestal density.
The paper “Design and Physics Basis for the Upcoming DIII-D SAS-VW Campaign to Quantify Tungsten Leakage and Transport in a New Slot Divertor Geometry” by T. Abrams et al. was recently published in Physica Scripta (https://doi.org/10.1088/1402-4896/ac3c5f). This paper describes the physics basis for an upcoming set of DIII-D experiments to exploit the high SOL collisionality enabled by a tightly baffled slot divertor geometry, which is hypothesized to suppress tungsten leakage in DIII-D. A toroidal row of graphite tiles in the Small Angle Slot (SAS) divertor has been coated with 10-15 µm of tungsten. New spectroscopic viewing chords with in-vacuo optics have been developed to measure the W gross erosion source from the divertor surface with high spatial and temporal resolution. SOLPS-DIVIMP simulations conducted with drifts using the planned “V” shape predict a substantial reduction in W sourcing and SOL accumulation in either B×∇B direction relative to either the old SAS divertor shape or the open, lower divertor. DIVIMP-WallDYN modeling indicates that the mixed-material C/W surface layer reaches equilibrium within 10 s of plasma exposure; the surface composition remains >80% W in the common SOL region but is quickly covered by carbon in the private flux region. These physics models will be “stress-tested” in the upcoming FY22 DIII-D campaign to build confidence in advanced tungsten divertor geometries for next-step devices.
The paper “Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D” by R.S. Wilcox et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3b8b). Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. Toroidally asymmetric striations in heat and particle flux are observed in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets, which may be problematic in future devices if heat flux is incident on unintended surfaces.
The paper “Linear Simulation of Magnetohydrodynamic Plasma Response to Three-Dimensional Magnetic Perturbations in High-β_P Plasmas” by Ran Chen et al. was published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/ac39f4). In this paper, linear MHD plasma response to 3D magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β_P plasmas was analyzed numerically using the extended-MHD code M3D-C1. Compared to the vacuum response, single-fluid M3D-C1 simulations predicted a significant decrease of the resonant plasma response to applied n=3 MPs at the maximum-penetration radius as q_95 increases. While single-fluid results showed resonant penetration occurring only near the pedestal top where the E×B toroidal rotation frequency was zero, two-fluid simulations showed comparable resonant penetration at both the top of the pedestal and the foot of the internal transport barrier (ITB), where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot might be responsible for the experimentally observed formation of a two-step “staircase” structure in both electron density and temperature profiles, which leads to a considerable deterioration of global plasma performance when MPs were applied in high-β_P plasmas.
The paper “Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations” by S. Munaretto, recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3b89), analyzes the impact of the poloidal spectrum of the applied RMP on the footprint at the divertor plates that alters the heat and particle flux distributions. Experimentally it is observed that a 3D structure in the IR, Dalpha and langmuir probes data corresponds to a large measured resonant plasma response (m=nq, dominant in the high field side) to applied n=2 RMP, and it is not visible when it is small, while the size of the 3D structure is found to be independent of the amplitude of the non-resonant component of the plasma response (m>nq, dominant in the low field side). Simulations done with M3D-C1 show that an RMP poloidal spectrum that maximizes the resonant plasma response will also maximize the size of the magnetic footprints, and the poloidal spectrum that minimizes the resonant plasma response will also minimize the footprints’ size. One of the consequences of these findings is that it will be challenging to decouple the footprint size from the RMP–ELM control problem, for which the resonant coupling is a requisite.
The paper “Neoclassical toroidal viscosity torque prediction via deep learning” by M. Clement, et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3e83). This paper covers the development of GPECnet, which is a densely connected neural network that has been trained on General Perturbed Equilibrium Code (GPEC) data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs).
The paper “Doppler-Shift Compensated Spatial Heterodyne Spectroscopy (SHS) for Rapidly Moving Sources” by G. Burke, R. Fonck et al., has been published in Applied Optics (https://doi.org/10.1364/AO.426159) describing a novel technique to dramatically improve spectral resolution for the SHS diagnostic while maintaining very high optical throughput required for high-sensitivity studies. SHS measures the n=3-2 Balmer-alpha transition of the full-energy neutral beam emission manifold at high spectral, spatial (dR~1 cm) and temporal (dt~0.1 ms) resolution, allowing for measurements of the local magnetic field (proportional to the spectral separation of the Stark-split pi-components) in the H-Mode pedestal for application to current profile measurements across L-H transitions, ELMs and other fast phenomena. High resolution-luminosity product measurements of neutral beam emission in magnetized plasmas are typically severely limited by the varying Doppler broadening across large diameter collection optics. A broadening compensation method is developed for the spatial heterodyne spectroscopy interferometric technique. The compensation technique greatly reduces this window broadening, thereby enabling high resolution measurements at a significantly higher photon flux than previously available. Compensated and uncompensated measurements of emission generated by impact excitation of 61 keV deuterium neutrals in a tokamak plasma at the DIII-D National Fusion Facility are presented. The spectral width of the compensated measurement is ∼0.13 nm, which is comparable to the instrument resolution. This width is ∼4× smaller than the uncompensated width, which for the 20 cm diameter collection lens system utilized in this study is ∼0.5 nm. Pedestal measurements of the magnetic field evolution across ELMs provided the the basis for G. Burke’s APS-DPP 2021 Invited Talk. Top: Comparison of compensated (SHS) and uncompensated (survey) beam emission spectrum. Note that the SHS integration time is 0.002 s compared to 2 seconds for the survey spectrometer. Bottom: beam emission spectrum model showing SHS eliminates geometric doppler broadening.
A paper “Achievements of Actively Controlled Divertor Detachment Compatible with Sustained High Confinement Core in DIII-D and EAST” by Liang Wang was published in Nuclear Fusion (http://doi.org/10.1088/1741-4326/ac4774). It reports the recent progress made on improving the integration of divertor detachment and core plasma performance jointly in DIII-D and EAST, including the successful development on detachment feedback control techniques for core-edge integration. DIII-D achieved actively controlled full detachment with low plasma temperature (Tet ≤ 5 eV across the entire divertor target plate) and low particle flux simultaneously with excellent core performance (βN ~3, βP >2 and H98~1.5) in the high βP scenario. The high βP scenario facilitates divertor detachment due to its relatively lower ETB, while a very high core energy confinement is enabled by a large-radius ITB. The synergy between ITB and ETB benefits the detachment access without degradation of core performance. The detachment leads to the reduction of ETB, which promotes the development of an even stronger ITB at large radius, and thus improves the plasma confinement eventually. In addition, it is demonstrated nitrogen seeding is more efficient to access full detachment, while neon seeding is more compatible with core plasma and ELM mitigation/suppression. Together with long pulse detachment achievement in EAST, these joint advances show a great potential for achieving a high-performance core plasma suitable for fusion reactors with controllable PWIs.
The paper “The role of BT-dependent flows on W accumulation at the edge of the confined plasma” by S. Zamperini et al. was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac3fe7). Edge codes have long predicted an accumulation of impurities near the separatrix, yet experimental evidence has been lacking. Experiments performed during the 2016 DIII-D Metal Rings Campaign have provided the first experimental evidence of near-SOL accumulation, though only when operating with the ion Bx∇B drift away from the divertor. During the campaign, two toroidally symmetric W rings were installed along the outer divertor to provide a trace source of W impurities into the SOL. Collector probes were inserted into the far-SOL, upon which W impurities were deposited on both sides. Modelling of W transport in the SOL was performed with DIVIMP and the recently upgraded far-SOL impurity transport code 3DLIM. It was found that in order to reproduce the experimentally measured deposition profiles in 3DLIM, the impurities must enter the far-SOL from above the outboard midplane. DIVIMP simulations demonstrate that this only occurs when W accumulates in the near-SOL. Furthermore, it is shown in DIVIMP that often observed fast inner-target parallel flows in the favorable BT direction can “flush out” W that would otherwise accumulate. Ad-hoc addition of these flows was necessary to reproduce the deposition profiles for the probe inserted for the favorable BT direction.
The paper “Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D” by J. Brooks et al., published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac36f4), presents predictive modeling of the erosion properties and surface evolution of the tungsten coating in the upcoming DIII-D SAS-VW divertor. Erosion and redeposition simulations conducted using the fully-kinetic, 3D Monte Carlo code REDEP/WBC, based on a SOLPS-ITER plasma background, estimate a peak net tungsten erosion rate of 0.1 nm/s and a redeposition fraction of ~80% in attached plasma conditions (Te,max = 40 eV), with an 8% leakage flux out of the divertor towards the core. The ITMC-DYN material mixing/response model predicts that within 30 s of cumulative attached plasma exposure (2% carbon flux fraction), migration of eroded carbon to the SAS-VW tungsten-coated surface will yield an equilibrium W surface concentration of ~90%. The simulations presented in the manuscript demonstrate that the SAS-VW divertor can be effectively used to study high-Z impurity sourcing and local transport in a closed slot geometry.
The paper “Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities” by X.D. Du et al., has been recently published in Phys. Rev. Lett. (https://doi.org/10.1103/PhysRevLett.127.235002). For the first time, fast ion phase-space flow, driven by Alfvén eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D experiment. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E−(ω/n)Pζ surfaces, where E, Pζ are the energy and canonical toroidal momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. The measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.
The paper “Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters” by E. U. Zeger, et al. was recently published in IEEE Transactions on Plasma Science (https://doi.org/10.1109/TPS.2021.3114608). An experimental data-driven machine learning approach predicts the pedestal heights and widths of electron density and temperature from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The most important parameters for setting the pedestal structure were plasma current, toroidal magnetic field, neutral beam heating power and shaping quantities. The deep jointly informed neural network (DJINN) algorithm was applied to identify suitable neural network (NN) architectures. The NN outperformed simple linear regression (LR) analysis, indicating nonlinear dependencies in the pedestal structure. [E. U. Zeger (Stanford University) performed her research during the remote 2020 Science Undergraduate Laboratory Internship (SULI) program under the supervision of F. M. Laggner and A. Bortolon and her work received an APS-DPP Outstanding Undergraduate Poster Award.]
The paper “Numerical assessment of the new V-shape small-angle slot divertor on DIII-D” by Roberto Maurizio was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac27c8). This paper reports on the numerical investigations on the effect of V-shaped divertor target in comparison to the small angle slot (SAS) divertor design on DIII-D. Simulations show that the electron temperature near the strike point is reduced in SAS-V compared to SAS, for both magnetic field directions, such that SAS-V achieves divertor detachment at a lower value of the outboard mid-plane separatrix electron density. The detachment threshold is lower because the V-shape focuses recycling neutrals on the V-end, densifying and cooling the plasma in the slot. At sufficiently high density, the V-shape also reduces the radial gradient of the temperature profile at the target, which in turns reduces the radial electric field and the E × B drift velocities, further densifying and cooling the plasma in the slot and leading to detachment.
DIII-D has undergone a major upgrade and successfully injected high power off-axis neutral beams in both co-current and counter-current directions. This capability of high power co/counter steerable off-axis neutral beams on a major tokamak opens a unique parameter space of broad pressure and current profiles for high beta steady-state advanced tokamak (AT) scenarios, while retaining the ability to balance the injected torque for low rotation studies. A recent Nuclear Fusion paper by B.A. Grierson et. al. (https://doi.org/10.1088/1741-4326/ac2872) reports on the critical evaluation of the transmitted power and energetic ion population produced by this unique heating and current drive system, which is assessed through visible imaging, neutron measurements and rotation profile measurements at balanced torque. Minimal losses of neutral beam power have been achieved by optimizing the strongly focused ion sources required to pass through the aperture, and updated calibration coefficients have been implemented to accurately account for the injected power. Good ability to balance the neutral beam torque has been demonstrated by injecting the new off-axis counter injecting beam against the existing off-axis co-injected beam.
The paper “Micro-trench Measurements of the Net Deposition of Carbon Impurity Ions in the DIII-D Divertor and the Resulting Suppression of Surface Erosion” by S. Abe et al has been published in Physica Scripta (https://doi.org/10.1088/1402-4896/ac2af4). This paper reports carbon impurity ion incident angles and deposition rates, along with silicon erosion rates, from measurements of micro-engineered trenches on a silicon surface exposed to L-mode deuterium plasmas at the DIII-D divertor. Post-exposure ex-situ analysis revealed carbon deposition profiles on the trench floor, which showed carbon ion shadowing that was consistent with ERO calculations of average carbon ion angle distributions (IADs) for both polar and azimuthal angles. Measured silicon net erosion rates negatively correlated with the deposited carbon concentration at different locations. Differential erosion of surfaces on two different ion-downstream trench slope structures suggested that carbon deposition rate is affected by the carbon ion incident angle and significantly suppressed the surface erosion. The results suggest the C impurity ion incident angles, determined by the IADs and surface morphology, strongly affect erosion rates as well as the main ion (D, T, He) incident angles.
The paper “Diverted negative triangularity plasmas on DIII-D: the benefit of high confinement without the liability of an edge pedestal” by A. Marinoni et al., published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/ac1f60), describes recent experiments in high confinement L-mode edge plasmas with a negative triangularity shape. The work extends previous results obtained in inner wall limited discharges to a novel LSN configuration and challenges the conventional notion that pedestals are required for operation at high confinement. By reversing triangularity, plasmas do not transition to H-mode despite the auxiliary power far exceeds the L-H power threshold expected from conventional scaling laws. The resulting L-mode edge regime features a power degradation of confinement substantially weaker than the ITER-89P scaling, and sustains normalized confinement and pressure levels typical of standard H-mode scenarios (H_{98,y2}≈1, beta_N≈1). Additional benefits of this regime characterized by relaxed edge pressure profiles are low impurity retention, as displayed by tau_P/tau_E≈1, and a SOL power fall off length exceeding that typically measured in H-mode plasmas.
The paper “Optimizing the differential connection schemes for detecting 3D magnetic perturbations in DIII-D” by S.Munaretto et al., recently published in Review of Scientific Instruments (https://doi.org/10.1063/5.0045453), discusses the methods used at DIII-D to efficiently optimize the number, the location and the connections of a toroidal array of magnetic probes in order to provide the most accurate measurement of stationary or slowly rotating 3D magnetic structures. Acquiring the difference between two sensors highlights these 3D structures, but the desire of simultaneously detecting multiple toroidal harmonics (n) results in a large number of possible connections. Different approaches to optimize the connections by minimizing the condition number (K) of the design matrix -a matrix containing a combination of desired n and sensor locations and connections- are presented, as well as an overview of the updated 3D magnetic diagnostic system of DIII-D.
The paper “Experimental evidence of runaway electron tail generation via localized helical structure in pellet-triggered tokamak disruptions” by X.D. Du et al., has been published in Nucl. Fusion (https://doi.org/10.1088/1741-4326/ac1e5e). The paper describes the development of novel detector, using stacked BGO crystals, for runaway electron studies in the DIII-D tokamak. It is able to resolve fast dynamics of high-energy tail formation of REs with an ultra-high time resolution of ∼ 1 μs. As a cost, the detector estimates the ‘effective’ energy of a given shape of γ-ray spectra and sacrifices the energy resolution. In aid of the new measurement capability, a rapid, inhomogeneous growth of RE tail is observed in detail during a major disruption triggered by an argon pellet. It is found that both the population and energy of a well-confined RE tail significantly oscillate at the early period of the growth. The oscillation phase is locked to a slow rotating magnetohydrodynamic instability, which is briefly destabilized for only ∼ 1 ms at the early period of the current quench. The oscillation ceases promptly, when the mode disappears. The data suggests that the high-energy RE tail is well-confined and accelerated via a localized helical structure in the plasma core.
The paper “Multiscale Chirping Modes Driven by Thermal Ions in a Plasma with Reactor-Relevant Ion Temperature” by X.D. Du et al., has been published in Physical Review Letters (https://doi.org/10.1103/PhysRevLett.127.025001). The paper describes the observation of a thermal ion driven bursting instability with rapid frequency chirping in a plasma with reactor-relevant ion temperature in 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 reconnection in the plasma edge, which can lead to a minor disruption event. The ηi (=∂lnT_i/∂lnn_i) exceeds the theoretically-predicted threshold for the destabilization of Alfvénic continuum modes due to compressibility of core ions. The most unstable modes belong to the strongly coupled kinetic ballooning mode and β-induced Alfvénic eigenmode branch. The key features of the observation are successfully reproduced by linear analysis solving the electromagnetic gyrokinetic equations via CGYRO. Since the mode, assessed to be an Alfvénic ion temperature gradient mode, is typically observed at high ion temperature (>10 keV) and high-β plasma regimes, the manifestation of the mode in future reactors should be studied with the corresponding development of mitigation strategies, if needed.
The Letter “Experimental Signatures of Electron Cyclotron Wave Energy Condensation in Magnetic Islands” by L. Bardoczi and N. C. Logan recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac2413) reports bifurcation of ECCD driven electron temperature perturbation (dTe) in magnetic islands for the first time. In this DIII-D experiment dTe spontaneously increases by 50% in a growing 2,1 island which in turn rapidly shrinks, consistent with improved stabilization by the ECCD due to the temperature bifurcation. This non-linearity of dTe maps out a hysteresis in the island evolution cycle, in qualitative agreement with the theoretically predicted radio frequency current condensation effect by A. H. Reiman and N. J. Fisch [Phys. Rev. Lett. 121 225001 (2018)]. This theory predicts the non-linear amplification of the stabilizing component of the ECCD at the island O-point, following from the facts that (i) the confining magnetic flux surface topology of the island leads to a temperature peak at the O-point when it is subject to a heat source and (ii) the electron cyclotron waves deposit their energy at the electron distribution tail. Therefore, this observed temperature bifurcation can be critical for the active stabilization of tearing modes in tokamaks.
Ion temperature and toroidal rotation velocity fluctuation measurements have been obtained on DIII-D with the Ultra-Fast Charge Exchange Recombination Spectroscopy (UF-CHERS) diagnostic, as described in a recent Review of Scientific Instruments paper, https://doi.org/10.1063/5.0043095 by D. D. Truong et al., which complements an Invited Talk at the 2020 High-Temperature Plasma Diagnostics Conference. The paper describes multiple upgrades to the UF-CHERS diagnostic and optimization of the data analysis procedures that have enabled ion temperature and parallel velocity fluctuation measurements arising from turbulence in DIII-D. Time-averaged UF-CHERS measurements show good agreement with CER. Ensemble averaging during steady plasma conditions (~1-2 sec) provides sufficient signal-to-noise ratio and reveals ion temperature and parallel velocity spectra with broadband fluctuations from 0 to ~300 kHz and the Geodesic Acoustic Mode near 20 kHz. The cross-phase between ion temperature and parallel velocity fluctuations was also measured and provides an important constraint for validating gyrokinetic simulations.
The paper “Physics basis for design of 3D coils in tokamaks” by N.C. Logan et al. recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/abff05) has shown that the use of stellarator tools in designing non-axisymmetric tokamak coils can greatly enhance the efficiency of these coils for rotation and stability control. The paper details how advancements in the physics understanding of exactly what non-axisymmetric or “3D” fields have the largest impact on tokamak plasmas (assisted by DIII-D experimental discovery and validation) has enabled targeting these fields directly with appropriately shaped coil sets. This type of problem is what the stellarator community has been solving for decades, and the tokamak community can now greatly benefit from the advanced optimization schemes originally developed for stellarator coil designs. The result of this new physics basis and cross-cutting optimization work is that these 3D coils can be built further away from the plasma, relaxing the cost and complexity of such systems in reactors. Closer to home, this physics basis is being used to assess future coil upgrades on DIII-D itself.
The paper “ELM and inter-ELM heat and particle flux to a secondary divertor in the DIII-D tokamak” by R. Perillo et. al., has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac0f38). It is shown that the ELM heat flux to the secondary divertor can be up to 1/3 of the total ELM heat flux (primary + secondary divertors), which would result in several tens of MW/m^2 on the upper tiles of the ITER vessel. It is experimentally demonstrated for the first time that, when the distance between the primary and secondary separatrices at the outboard midplane (dRsep) is less than 10mm, significant ELM heat flux is deposited to the secondary inner target, with peak heat flux values that are comparable to those at the secondary outer strike point. This poses a concern for any future tokamak where the secondary inner target is not yet designed to withstand significant heat loads. The ELM heat flux profiles measured at the secondary outer region feature rippled sub-structures and extend into the far-SOL, while at the inner target they are discrete and well defined. It is found that values of dRsep above 25-30mm are needed to reduce the ELM heat flux to the secondary divertor below 10% of the total ELM heat flux deposited to a single-null configuration, implying that the up-down magnetic balance will have to be accurately controlled in order to preserve the plasma facing components at the secondary X-point region in future machines.
Since most reactor plasma facing component (PFC) designs call for mixed-material environments, including ITER's W/Be environment, developing an understanding of how reconstituted surfaces evolve under plasma bombardment is essential. The paper “Modeling of ExB effects on tungsten re-deposition and transport in the DIII-D divertor” by J. Nichols et al. was recently published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac14e6). This work identifies controlling physics that affects material migration patterns in the diverter during an experiment from the DIII-D Metal Rings Campaign in which repeat L-mode discharges were conducted on the W-coated divertor tiles. These simulations indicate that radial and poloidal ExB transport dominates over parallel force balance for high-Z impurities such as W in the divertor region of DIII-D. Time-dependent simulations with scaled ExB impurity drifts and W re-erosion quantitatively reproduced experimentally measured features, including depth-resolved W/C ratios, within a factor of 2 over ~115 seconds of accumulated plasma exposure. The location of co-deposition regions was also shown to be well-represented by an analytic leakage model, driven largely by poloidal ExB drifts.
The paper “Estimate of pre-thermal quench non-thermal electron density profile during Ar pellet shutdowns of low-density target plasmas in DIII-D” by E. Hollmann et al. has been published in Physics of Plasmas (https://doi.org/10.1063/5.0050903). The paper describes first efforts to reconstruct the radial profile of early-time non-thermal electron seeds during disruptions in DIII-D. During disruptions, the thermal temperature collapses due to a combination of radial transport to the wall and impurity radiation. In regions where impurity radiation dominates, non-thermal electron populations can form due to the hot-tail mechanism, where electrons are abandoned at the tail of the distribution function as the bulk of the distribution function cools. These multi-keV non-thermal seeds can subsequently be accelerated to MeV energies during the subsequent toroidal current decay. Observing these early non-thermal seeds is extremely challenging, but has been accomplished in some low-density target plasmas shut down with Ar pellet injection. By combining SXR profiles with ECE profiles, estimates of the non-thermal electron profiles could be made. Interestingly, the non-thermal electrons are observed to appear ahead of the Ar pellet radially; it is not clear yet if this is due to extremely rapid cross-field heat transport or extremely fast cross-field argon ion transport.
A study on the isotope dependence of beta-induced Alfven eigenmode (BAE) and low frequency mode (LFM) stability in DIII-D was published in Nuclear Fusion, https://doi.org/10.1088/1741-4326/ac189d. With deuterium neutral beam injection (NBI), BAEs are at least as unstable in mixed-species plasmas as in deuterium plasmas; however, with hydrogen NBI, the BAEs are stable. In contrast, the LFMs are unaffected by changes in beam species, consistent with the previous observation that LFMs are not driven by high-energy beam ions. As predicted by theory and simulation, the LFMs are more unstable in mixed species plasmas than in pure deuterium plasmas.
The Letter “Neoclassical Tearing Mode Seeding by Nonlinear Three-Wave Interactions in Tokamaks” by Laszlo Bardoczi at al has been published in Physical Review Letters (https://doi.org/10.1103/PhysRevLett.127.055002). The Letter reports the first experimental observation of seed magnetic island formation by nonlinear three-wave coupling of magnetic island triplets. In these DIII-D ITER baseline scenario plasmas, disruptive 2,1 islands are seeded by coupling of 4,3 and 3,2 islands to a central 1,1 sawtooth precursor. These plasmas reach the beta and current flattop, and are robustly stable to classical tearing modes. A number of ELMs and sawtooth crashes undergo without seeding the 2,1 island in the absence of non-linear three-wave coupling. Seeding occurs when high m,n islands couple, at various times in the stationary plasma. Therefore, these observations clarify that magnetic reconnection at q=2 is not caused by a classical current driven instability, but the 2,1 NTM seed island is formed by frequency matching and nonlinearly interacting islands that satisfy the mode number resonance condition. Three-wave interactions are conclusively identified with bi-spectral analysis, indicating fixed phase relationships at the time of 2,1 seeding. These results are general and relevant for future reactors, as 70% of the considered unstable DIII-D ITER baseline scenario discharges of the past decade without ECCD are characterized by frequency matching of resonant tearing modes at the time of 2,1 island seeding. This mechanism is also a candidate to account for TM cascades which has remained unexplained since it's first observation in 1989 by E. J. Strait et al [62(11) 1282 PRL]. As the non-linear three-wave interaction produces seed 2,1 islands in classically stable plasmas, tearing free operation may not be possible by locking in a classically stable current profile. This predicts new challenges for future reactors that must operate in stable plasma equilibria free of disruptive 2,1 islands, calling for high differential rotation at q=2, active control, suppression and avoidance of high m,n modes as much as possible.
A grand challenge for the tokamak approach to fusion energy production is the achievement of a high-performance fusion core integrated with a boundary solution compatible with available first-wall materials. A key issue at this interface is the edge-localized mode (ELM), an instability arising from edge gradients which if uncontrolled delivers an impulsive heat and particle load to the first-wall. A recently published topical review paper “Plasma performance and operational space without ELMs in DIII-D” by C. Paz-Soldan (https://doi.org/10.1088/1361-6587/ac048b) compares the operating space and plasma performance of the stationary no-ELM regimes found in conventional tokamaks through a database study. Normalized performance in DIII-D favores regimes tolerant to power, and peak absolute performance is found to rise with the product of field, current, and radius. The highest absolute performance is found for RMP-ELM suppressed plasmas, but only at high torque and with low pedestal pressure fraction. The highest pedestal pressure fraction is found for QH-mode and EDA-H-mode regimes, though with higher carbon fraction. The normalized performance of negative triangularity (Neg-D) plasmas is comparable to RMP and QH-mode plasmas, though the absolute performance still lags. Considering integration with electron heating, the high pedestal pressure fraction regimes (QH and EDA-H) preserve the highest performance. Considering divertor integration, only the EDA-H, Neg-D, and L-mode scenarios have approached divertor-friendly high separatrix density conditions, with Neg-D preserving the highest performance owing to its compatibility with both high power and density in DIII-D.
The paper “Overview of density pedestal structure: role of fueling versus transport” by S. Mordijck was published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ab8d04). This paper presents a review of physics affecting the pedestal density structure. Evidence from existing machines shows that both fuelling and transport changes can modify the density pedestal. However, due to lack of measurements of the particle source, particle transport is not sufficiently understood to determine if an inward pinch exists. Thus, it is not possible to predict the density pedestal structure in future machines, such as ITER, due to high plasma densities which will screen the pedestal from neutral fuelling. The paper also discusses experiments on DIII-D and C-Mod which were performed to increase pedestal neutral opaqueness na, (product of density and minor radius) toward values expected on ITER. The data sets from the two machines spanned an opaqueness range from 10× to 2× smaller than ITER with DIII-D operating at the smaller values of na and C-Mod at the values closer to ITER. For both devices, the pedestal density increased over the full range of opaqueness whereas the ratio of pedestal to separatrix densities remained constant. These results are consistent with profile measurements showing no significant change of density gradient with opaqueness.
The paper “Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D” by S. R. Haskey et al. has been published in the Review of Scientific Instruments (https://doi.org/10.1063/5.0043836). The paper describes how radially resolved hydrogenic isotope fraction measurement capabilities using the main-ion CER (MICER) spectroscopy system were developed for the recent mixed H/D experiments in the DIII-D isotope thrust. The technique is based on the active (neutral beam on) charge exchange Hα and Dα brightness, which provides a profile of the isotope fraction of the ions from the core to the region approaching the separatrix. Extensive atomic physics calculations were performed using FIDASIM to better understand the relationship between the brightness and the underlying ion densities. 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. These details must be taken into consideration when calculating absolute densities.
The Letter “Experimental inference of flux tunneling between magnetic island chains in tokamaks” by Laszlo Bardoczi and Todd E. Evans has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/ac0411). The Letter reports confinement bifurcations due to flux tunneling in 2,1 islands, triggered by coupling of non-overlapping adjacent 5,2 islands. While subject to ECCD, the temperature at the 2,1 O-point is flat in the coupled phase but a temperature peak forms at the 2,1 island O-point after decoupling. This change in the EC wave energy confinement within the 2,1 island O-point region is shown to be due to a bifurcation from stochastic to nested 2,1 magnetic island topology, caused by flux tunneling through intersecting manifolds of overlapping heteroclinic tangles. The effect of the ECCD on the 2,1 island is shown to correlate with coupling events to a 4,3 island in an ITER baseline scenario plasma: (i) ECCD can not stabilize coupled islands but (ii) it suppresses freely rotating islands to a small size. This shows the critical impact of flux tunneling on disruption avoidance in tokamaks.
A quantitative understanding of RMP ELM suppression is essential for predicting robust ELM suppression in ITER. The recent paper “Predicting operational windows of ELMs suppression by Resonant Magnetic Perturbations in the DIII-D and KSTAR tokamaks” by Q. Hu, et. al. published in Physics of Plasmas as a featured article (https://doi.org/10.1063/5.0043018) reports a newly developed plasma response model quantitatively predicted the narrow isolated q95 windows of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The paper reveals the key physics of close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM permits an applied RMP to produce field penetration due to the lower E×B rotation at the pedestal top rather than being screened and leads to very narrow q95 windows of ELM suppression. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. This can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.
The paper “Evolution of ELMs, pedestal profiles and fluctuations in the inter-ELM period in NBI- and ECH-dominated discharges in DIII-D” by S. Banerjee et. al. has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/abe8b1). In this paper, characterization of how the H-mode pedestal and fluctuations evolve between repeated edge oscillations called edge-localized modes (ELMs) provides insight into the nature of pedestal transport and path towards the ELM trigger in ITER similar shape plasmas. As the fraction of heat transported through the electron is increased relative to the ions, fluctuation measurements reveal distinct modes that appear, which correlate with the evolution of the electron temperature gradient. Transport coefficients obtained from TRANSP show that MTM and/or TEM are plausible candidates for the observed fluctuations. Linear gyrofluid simulations with TGLF corroborate this characterization. When the electron temperature gradient is increased and the modes are present, the frequency of ELMs decreases by 40% with inter-ELM spacing becoming more regular. This reduction in ELM frequency is attributed to the increased fluctuations during electron heating, causing an increase in fluctuation-driven transport in the pedestal and slowing the pedestal recovery between ELMs.
The paper “Noise suppression for MHD characterization with Electron Cyclotron Emission Imaging 1D technique” by G. Yu et al. has been published in Plasma Physics and Controlled Fusion (https://doi.org/10.1088/1361-6587/abe9f2). The paper introduces a data processing technique on ECE-imaging that significantly improves the spectrum quality and avoids biased fitting of the radial structure for coherent MHDs, such as Alfven Eigenmode and Tearing Mode. The technique is particularly useful for transient (~1 ms) coherent modes, as the noise could not be suppressed by a long data length (>10 ms). ECE Radiation modeling shows that the radial resolution of the technique is limited by the finite radiation volume effect and the non-uniform gain of the RF receiver. A narrower bandwidth (<~100 MHz) RF receiver can improve the ECEI capability resolving islands of half-width < 1.5 cm.
The article “Evaluation of silicon carbide as a divertor armor material in DIII-D H-mode discharges” by T. Abrams et al. has been published in Nuclear Fusion (https://doi.org/10.1088/1741-4326/abecee). This paper describes an analytic mixed-material SiC erosion model that uses calculated physical and chemical sputtering yields from SiC, Si, and C. The model was tested using DIII-D DiMES experiments involving graphite-coated SiC samples exposed at the strike-point to repeated H-mode plasma discharges. The qualitative trends from analytic modeling were reproduced by the experimental measurements, obtained via spectroscopic inference using the S/XB method. Post-mortem compositional analysis of the SiC-DiMES samples reveals Si enrichment of about 10%, in line with expectations from the erosion model. Encouragingly, minimal changes to the macroscopic or microscopic surface morphology of the SiC coatings were observed.
In preparation for the upcoming MAST-U campaign, pedestal stability of spherical tokamaks has been revisited by a team of DIII-D and MAST-U researchers by investigating standard H-mode discharges on MAST. The results summarized in the paper “Pedestal stability analysis on MAST in preparation for MAST-U” by M. Knolker et. al. have been published in Nuclear Fusion (https://iopscience.iop.org/article/10.1088/1741-4326/abe804). As a step beyond previous studies, both ion and electron profiles are used for obtaining equilibria and a diverse set of pedestals is evaluated. Stability analysis with the ELITE and CGYRO codes shows that MAST pedestals are constrained by kinetic ballooning modes and medium toroidal mode number peeling-ballooning modes, with most unstable modes ranging from n=25 to n=45. In discharges with a steep q profile at the edge a larger number of poloidal harmonics is excited for each toroidal mode. A comparison with discharges on DIII-D with matched shape and similar non-dimensional parameters indicates that the increased shear at lower aspect ratio stabilizes low n peeling modes.
The paper “Fast modulating electron cyclotron emission (FMECE) diagnostic for tokamaks” by Saeid Houshmandyar et al. (The University of Texas at Austin) has been published in Review of Scientific Instruments (https://doi.org/10.1063/5.0043761). This paper presents the design and development of FMECE, a stand-alone IF section with eight channels, which is a novel application of YIG filters for real-time electron temperature gradient and gradient scale lengths measurements. Implementation and testing of FMECE at DIII-D tokamak, demonstrates its capability in real-time gradient measurements. The data presented in the paper shows that FMECE can identify flattening in the Te-profile; the latter can be used as a sensor for real-time monitoring and control of plasma instabilities. Implementation and application are planned for the EAST tokamak.
The paper “Experimental Observation of Magnetic Island Heteroclinic Bifurcation in Tokamaks” by L. Bardoczi and T. E. Evans was published in Physical Review Letters (https://doi.org/10.1103/PhysRevLett.126.085003). In this new class of magnetic island bifurcations, the second harmonic perturbation grows to a large amplitude causing the number of O-points and X-points of the primary island to double. For example, a 2,1 island's O-points and X-points can double due to a large apparent 4,2 mode. The new O-line is disconnected from the original one, hence called heteroclinic. This bifurcation is observed for the first time in interacting coupled 2,1 tearing modes in DIII-D (see figure). Poincaré maps constrained by measured magnetic amplitudes and phasing show bifurcation from heteroclinic to homoclinic topology in the 2,1 island as the 4,2 relative amplitude decreases. In agreement with external magnetic data, the local electron temperature peak in the 2,1 island splits initially, consistent with two O-points. As the 4,2 relative amplitude decreases a single Te peak forms in the 2,1 island, consistent with one O-point. Heteroclinic island stabilization by ECCD is expected to be significantly harder and strong 4,2 components are common in DIII-D ITER baseline scenario discharges with unstable 2,1 islands. These observations call for developing tearing stability theory and control solutions for heteroclinic islands in tokamaks.EPS Conference Proceedings:
(Earlier papers in separate database)
Link to internal page for preparation of weekly papers is here