Improved profile fitting and quantification of uncertainty in experimental measurements of impurity transport coefficients using Gaussian process regression

The need to fit smooth temperature and density profiles to discrete observations is ubiquitous in plasma physics, but the prevailing techniques for this have many shortcomings that cast doubt on the statistical validity of the results. This issue is amplified in the context of validation of gyrokinetic transport models (Holland et al 2009 Phys. Plasmas 16 052301), where the strong sensitivity of the code outputs to input gradients means that inadequacies in the profile fitting technique can easily lead to an incorrect assessment of the degree of agreement with experimental measurements. In order to rectify the shortcomings of standard approaches to profile fitting, we have applied Gaussian process regression (GPR), a powerful non-parametric regression technique, to analyse an Alcator C-Mod L-mode discharge used for past gyrokinetic validation work (Howard et al 2012 Nucl. Fusion 52 063002). We show that the GPR techniques can reproduce the previous results while delivering more statistically rigorous fits and uncertainty estimates for both the value and the gradient of plasma profiles with an improved level of automation. We also discuss how the use of GPR can allow for dramatic increases in the rate of convergence of uncertainty propagation for any code that takes experimental profiles as inputs. The new GPR techniques for profile fitting and uncertainty propagation are quite useful and general, and we describe the steps to implementation in detail in this paper. These techniques have the potential to substantially improve the quality of uncertainty estimates on profile fits and the rate of convergence of uncertainty propagation, making them of great interest for wider use in fusion experiments and modelling efforts.United States. Dept. of Energy. Office of Fusion Energy Sciences (Award DE-FC02-99ER54512)United States. Dept. of Energy. Office of Science (Contract DE-AC05-06OR23177)United States. Dept. of Energy. Office of Advanced Scientific Computing Research (Award DE-SC0007099

Measurement of LHCD edge power deposition through modulation techniques on Alcator C-Mod

The efficiency of LHCD on Alcator C-Mod drops exponentially with line average density. At reactor relevant densities (> 1 · 1020 [m[-3 superscript]]) no measurable current is driven. While a number of causes have been suggested, no specific mechanism has been shown to be responsible for the loss of current drive at high density. Fast modulation of the LH power was used to isolate and quantify the LHCD deposition within the plasma. Measurements from these plasmas provide unique evidence for determining a root cause. Modulation of LH power in steady plasmas exhibited no correlated change in the core temperature. A correlated, prompt response in the edge suggests that the loss in efficiency is related to a edge absorption mechanism. This follows previous results which found the generation of n||-independent SOL currents. Multiple Langmuir probe array measurements of the conducted heat conclude that the lost power is deposited near the last closed flux surface. The heat flux induced by LH waves onto the outer divertor is calculated. Changes in the neutral pressure, ionization and hard X-ray emission at high density highlight the importance of the active divertor in the loss of efficiency. Results of this study implicate a mechanism which may occur over multiple passes, leading to power absorption near the LCFS

Poloidal asymmetries in edge transport barriers

Measurements of impurities in Alcator C-Mod indicate that in the pedestal region, significant poloidal asymmetries can exist in the impurity density, ion temperature, and main ion density. In light of the observation that ion temperature and electrostatic potential are not constant on a flux surface [Theiler et al., Nucl. Fusion 54, 083017 (2014)], a technique based on total pressure conservation to align profiles measured at separate poloidal locations is presented and applied. Gyrokinetic neoclassical simulations with XGCa support the observed large poloidal variations in ion temperature and density, and that the total pressure is approximately constant on a flux surface. With the updated alignment technique, the observed in-out asymmetry in impurity density is reduced from previous publishing [Churchill et al., Nucl. Fusion 53, 122002 (2013)], but remains substantial (nz,H/nz,L∼6). Candidate asymmetry drivers are explored, showing that neither non-uniform impurity sources nor localized fluctuation-driven transport are able to explain satisfactorily the impurity density asymmetry. Since impurity density asymmetries are only present in plasmas with strong electron density gradients, and radial transport timescales become comparable to parallel transport timescales in the pedestal region, it is suggested that global transport effects relating to the strong electron density gradients in the pedestal are the main driver for the pedestal in-out impurity density asymmetry.United States. Department of Energy (DE-FC02-99ER54512)United States. Department of Energy (DE-FG02-06ER54845)United States. Department of Energy (DE-FG02-86ER53223)United States. Department of Energy (DE-AC02-09CH11466

The physics mechanisms of the weakly coherent mode in the Alcator C-Mod Tokamak

The weakly coherent mode (WCM) in I-mode has been studied by a six-field two-fluid model based on the Braginskii equations under the BOUT++ framework for the first time. The calculations indicate that a tokamak pedestal exhibiting a WCM is linearly unstable to drift Alfven wave (DAW) instabilities and the resistive ballooning mode. The nonlinear simulation shows promising agreement with the experimental measurements of the WCM. The shape of the density spectral and location of the spectral peak of the dominant toroidal number mode n = 20 agrees with the experimental data from reflectometry. The simulated mode propagates in electron diamagnetic direction is consistent with the results from the magnetic probes in the laboratory frame, a large ratio of particle to heat diffusivity is consistent with the distinctive experimental feature of I-mode, and the value of the simulated χeat the edge is in the range of experimental errors of χefffrom the experiment. The prediction of the WCM shows that free energy is mainly provided by the electron pressure gradient, which gives guidance for pursuing future I-mode studies

Correlation ECE diagnostic in Alcator C-Mod

Correlation ECE (CECE) is a diagnostic technique that allows measurement of small amplitude electron temperature, T[subscript e], fluctuations through standard cross-correlation analysis methods. In Alcator C-Mod, a new CECE diagnostic has been installed[Sung RSI 2012], and interesting phenomena have been observed in various plasma conditions. We find that local T[subscript e] fluctuations near the edge (ρ ~ 0:8) decrease across the linearto- saturated ohmic confinement transition, with fluctuations decreasing with increasing plasma density[Sung NF 2013], which occurs simultaneously with rotation reversals[Rice NF 2011]. T[subscript e] fluctuations are also reduced across core rotation reversals with an increase of plasma density in RF heated L-mode plasmas, which implies that the same physics related to the reduction of T[subscript e] fluctuations may be applied to both ohmic and RF heated L-mode plasmas. In I-mode plasmas, we observe the reduction of core T[subscript e] fluctuations, which indicates changes of turbulence occur not only in the pedestal region but also in the core across the L/I transition[White NF 2014]. The present CECE diagnostic system in C-Mod and these experimental results are described in this paper

High-field side scrape-off layer investigation: Plasma profiles and impurity screening behavior in near-double-null configurations

New experiments on Alcator C-Mod reveal that the favorable impurity screening characteristics of the high-field side (HFS) scrape-off layer (SOL), previously reported for single null geometries, is retained in double null configurations, despite the formation of an extremely thin SOL. In balanced double-null, nitrogen injected locally into the HFS SOL is better screened by a factor of 2.5 compared to the same injection into the low field side (LFS) SOL. This result is insensitive to plasma current and Greenwald fraction. Nitrogen injected into the HFS SOL is not as well screened (only a factor of 1.5 improvement over LFS) in unbalanced double-null discharges, when the primary divertor is in the direction of B×∇B. In this configuration, impurity ‘plume’ emission patterns indicate that an opposing E × B drift competes with the parallel impurity flow to the divertor. In balanced double-null plasmas, the dispersal pattern exhibits a dominant E × B motion. Unbalanced discharges with the primary divertor opposite the direction of B×∇B exhibit excellent HFS screening characteristics – a factor of 5 enhancement compared to LFS. These data support the idea that future tokamaks should locate all RF actuators and close-fitting wall structures on the HFS and employ near-double-null magnetic topologies, both to precisely control plasma conditions at the antenna/plasma interface and to maximally mitigate the impact of local impurity sources arising from plasma-material interactions. Keywords: Alcator C-Mod; Impurity screening; Double null; High field side scrape-off layerUnited States. Department of Energy (Contract DE-FC02-99ER54512

Quasi-coherent fluctuations limiting the pedestal growth on Alcator C-Mod: experiment and modelling

Performance predictions for future fusion devices rely on an accurate model of the pedestal structure. The candidate for predictive pedestal structure is EPED, and it is imperative to test the underlying hypotheses to further gain confidence for ITER projections. Here, we present experimental work testing one of the EPED hypotheses, namely the existence of a soft limit set by microinstabilities such as the kinetic ballooning mode. This work extends recent work on Alactor C-Mod (Diallo et al 2014 Phys. Rev. Lett. 112 115001), to include detailed measurements of the edge fluctuations and comparisons of edge simulation codes and experimental observations

Nonlinear gyrokinetic simulations of the I-mode high confinement regime and comparisons with experimenta)

For the first time, nonlinear gyrokinetic simulations of I-mode plasmas are performed and compared with experiment. I-mode is a high confinement regime, featuring energy confinement similar to H-mode, but without enhanced particle and impurity particle confinement [D. G. Whyte et al., Nucl. Fusion 50, 105005 (2010)]. As a consequence of the separation between heat and particle transport, I-mode exhibits several favorable characteristics compared to H-mode. The nonlinear gyrokinetic code GYRO [J. Candy and R. E. Waltz, J Comput. Phys. 186, 545 (2003)] is used to explore the effects of E × B shear and profile stiffness in I-mode and compare with L-mode. The nonlinear GYRO simulations show that I-mode core ion temperature and electron temperature profiles are more stiff than L-mode core plasmas. Scans of the input E × B shear in GYRO simulations show that E × B shearing of turbulence is a stronger effect in the core of I-mode than L-mode. The nonlinear simulations match the observed reductions in long wavelength density fluctuation levels across the L-I transition but underestimate the reduction of long wavelength electron temperature fluctuation levels. The comparisons between experiment and gyrokinetic simulations for I-mode suggest that increased E × B shearing of turbulence combined with increased profile stiffness are responsible for the reductions in core turbulence observed in the experiment, and that I-mode resembles H-mode plasmas more than L-mode plasmas with regards to marginal stability and temperature profile stiffness.United States. Department of Energy (Contract No. DE-FC02-99ER54512-CMOD)United States. Department of Energy. Office of Science (Contract No. DE-AC02- 05CH11231

Non-local heat transport in Alcator C-Mod ohmic L-mode plasmas

Non-local heat transport experiments were performed in Alcator C-Mod ohmic L-mode plasmas by inducing edge cooling with laser blow-off impurity (CaF2) injection. The non-local effect, a cooling of the edge electron temperature with a rapid rise of the central electron temperature, which contradicts the assumption of 'local' transport, was observed in low collisionality linear ohmic confinement (LOC) regime plasmas. Transport analysis shows this phenomenon can be explained either by a fast drop of the core diffusivity, or the sudden appearance of a heat pinch. In high collisionality saturated ohmic confinement (SOC) regime plasmas, the thermal transport becomes 'local': the central electron temperature drops on the energy confinement time scale in response to the edge cooling. Measurements from a high resolution imaging x-ray spectrometer show that the ion temperature has a similar behaviour as the electron temperature in response to edge cooling, and that the transition density of non-locality correlates with the rotation reversal critical density. This connection may indicate the possible connection between thermal and momentum transport, which is also linked to a transition in turbulence dominance between trapped electron modes (TEMs) and ion temperature gradient (ITG) modes. Experiments with repetitive cold pulses in one discharge were also performed to allow Fourier analysis and to provide details of cold front propagation. These modulation experiments showed in LOC plasmas that the electron thermal transport is not purely diffusive, while in SOC the electron thermal transport is more diffusive like. Linear gyrokinetic simulations suggest the turbulence outside r/a = 0.75 changes from TEM dominance in LOC plasmas to ITG mode dominance in SOC plasmas.United States. Dept. of Energy (DoE Contract No DE-FC02-99ER54512)Oak Ridge Institute for Science and Education (DOE Fusion Energy Postdoctoral Research Program

20 years of research on the Alcator C-Mod tokamak

The object of this review is to summarize the achievements of research on the Alcator C-Mod tokamak [Hutchinson et al., Phys. Plasmas 1, 1511 (1994) and Marmar, Fusion Sci. Technol. 51, 261 (2007)] and to place that research in the context of the quest for practical fusion energy. C-Mod is a compact, high-field tokamak, whose unique design and operating parameters have produced a wealth of new and important results since it began operation in 1993, contributing data that extends tests of critical physical models into new parameter ranges and into new regimes. Using only high-power radio frequency (RF) waves for heating and current drive with innovative launching structures, C-Mod operates routinely at reactor level power densities and achieves plasma pressures higher than any other toroidal confinement device. C-Mod spearheaded the development of the vertical-target divertor and has always operated with high-Z metal plasma facing components—approaches subsequently adopted for ITER. C-Mod has made ground-breaking discoveries in divertor physics and plasma-material interactions at reactor-like power and particle fluxes and elucidated the critical role of cross-field transport in divertor operation, edge flows and the tokamak density limit. C-Mod developed the I-mode and the Enhanced Dα H-mode regimes, which have high performance without large edge localized modes and with pedestal transport self-regulated by short-wavelength electromagnetic waves. C-Mod has carried out pioneering studies of intrinsic rotation and demonstrated that self-generated flow shear can be strong enough in some cases to significantly modify transport. C-Mod made the first quantitative link between the pedestal temperature and the H-mode's performance, showing that the observed self-similar temperature profiles were consistent with critical-gradient-length theories and followed up with quantitative tests of nonlinear gyrokinetic models. RF research highlights include direct experimental observation of ion cyclotron range of frequency (ICRF) mode-conversion, ICRF flow drive, demonstration of lower-hybrid current drive at ITER-like densities and fields and, using a set of novel diagnostics, extensive validation of advanced RF codes. Disruption studies on C-Mod provided the first observation of non-axisymmetric halo currents and non-axisymmetric radiation in mitigated disruptions. A summary of important achievements and discoveries are included.United States. Dept. of Energy (Cooperative Agreement DE-FC02-99ER54512)United States. Dept. of Energy (Cooperative Agreement DE-FG03-94ER-54241)United States. Dept. of Energy (Cooperative Agreement DE-AC02-78ET- 51013)United States. Dept. of Energy (Cooperative Agreement DE-AC02-09CH11466)United States. Dept. of Energy (Cooperative Agreement DE-FG02-95ER54309)United States. Dept. of Energy (Cooperative Agreement DE-AC02-05CH11231)United States. Dept. of Energy (Cooperative Agreement DE-AC52-07NA27344)United States. Dept. of Energy (Cooperative Agreement DE-FG02- 97ER54392)United States. Dept. of Energy (Cooperative Agreement DE-SC00-02060