J. Walk

Multi-channel transport experiments at Alcator C-Mod and comparison with gyrokinetic simulationsa)

Multi-channel transport experiments have been conducted in auxiliary heated (Ion Cyclotron Range of Frequencies) L-mode plasmas at Alcator C-Mod [Marmar and Alcator C-Mod Group, Fusion Sci. Technol. 51(3), 3261 (2007)]. These plasmas provide good diagnostic coverage for measurements of kinetic profiles, impurity transport, and turbulence (electron temperature and density fluctuations). In the experiments, a steady sawtoothing L-mode plasma with 1.2 MW of on-axis RF heating is established and density is scanned by 20%. Measured rotation profiles change from peaked to hollow in shape as density is increased, but electron density and impurity profiles remain peaked. Ion or electron heat fluxes from the two plasmas are the same. The experimental results are compared directly to nonlinear gyrokinetic theory using synthetic diagnostics and the code GYRO [Candy and Waltz, J. Comput. Phys. 186, 545 (2003)]. We find good agreement with experimental ion heat flux, impurity particle transport, and trends in the fluc...

Improved understanding of physics processes in pedestal structure, leading to improved predictive capability for ITER

Joint experiment/theory/modelling research has led to increased confidence in predictions of the pedestal height in ITER. This work was performed as part of a US Department of Energy Joint Research Target in FY11 to identify physics processes that control the H-mode pedestal structure. The study included experiments on C-Mod, DIII-D and NSTX as well as interpretation of experimental data with theory-based modelling codes. This work provides increased confidence in the ability of models for peeling–ballooning stability, bootstrap current, pedestal width and pedestal height scaling to make correct predictions, with some areas needing further work also being identified. A model for pedestal pressure height has made good predictions in existing machines for a range in pressure of a factor of 20. This provides a solid basis for predicting the maximum pedestal pressure height in ITER, which is found to be an extrapolation of a factor of 3 beyond the existing data set. Models were studied for a number of processes that are proposed to play a role in the pedestal ne and Te profiles. These processes include neoclassical transport, paleoclassical transport, electron temperature gradient turbulence and neutral fuelling. All of these processes may be important, with the importance being dependent on the plasma regime. Studies with several electromagnetic gyrokinetic codes show that the gradients in and on top of the pedestal can drive a number of instabilities.

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-modes 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.

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.

Changes in core electron temperature fluctuations across the ohmic energy confinement transition in Alcator C-Mod plasmas

The first measurements of long wavelength (ky?s?<?0.3) electron temperature fluctuations in Alcator C-Mod made with a new correlation electron cyclotron emission diagnostic support a long-standing hypothesis regarding the confinement transition from linear ohmic confinement (LOC) to saturated ohmic confinement (SOC). Electron temperature fluctuations decrease significantly (?40%) crossing from LOC to SOC, consistent with a change from trapped electron mode (TEM) turbulence domination to ion temperature gradient (ITG) turbulence as the density is increased. Linear stability analysis performed with the GYRO code (Candy and Waltz 2003 J. Comput. Phys. 186 545) shows that TEMs are dominant for long wavelength turbulence in the LOC regime and ITG modes are dominant in the SOC regime at the radial location (????0.8) where the changes in electron temperature fluctuations are measured. In contrast, deeper in the core (??<?0.8), linear stability analysis indicates that ITG modes remain dominant across the LOC/SOC transition. This radial variation suggests that the robust global changes in confinement of energy and momentum occurring across the LOC/SOC transition are correlated to local changes in the dominant turbulent mode near the edge.

Pedestal structure and stability in H-mode and I-mode: a comparative study on Alcator C-Mod

New experimental data from the Alcator C-Mod tokamak are used to benchmark predictive modelling of the edge pedestal in various high-confinement regimes, contributing to greater confidence in projection of pedestal height and width in ITER and reactors. ELMy H-modes operate near stability limits for ideal peeling–ballooning modes, as shown by calculations with the ELITE code. Experimental pedestal width in ELMy H-mode scales as the square root of βpol at the pedestal top, i.e. the dependence expected from theory if kinetic ballooning modes (KBMs) were responsible for limiting the pedestal width. A search for KBMs in experiment has revealed a short-wavelength electromagnetic fluctuation in the pedestal that is a candidate driver for inter-edge localized mode (ELM) pedestal regulation. A predictive pedestal model (EPED) has been tested on an extended set of ELMy H-modes from C-Mod, reproducing pedestal height and width reasonably well across the data set, and extending the tested range of EPED to the highest absolute pressures available on any existing tokamak and to within a factor of three of the pedestal pressure targeted for ITER. In addition, C-Mod offers access to two regimes, enhanced D-alpha (EDA) H-mode and I-mode, that have high pedestals, but in which large ELM activity is naturally suppressed and, instead, particle and impurity transport are regulated continuously. Pedestals of EDA H-mode and I-mode discharges are found to be ideal magnetohydrodynamic (MHD) stable with ELITE, consistent with the general absence of ELM activity. Invocation of alternative physics mechanisms may be required to make EPED-like predictions of pedestals in these kinds of intrinsically ELM-suppressed regimes, which would be very beneficial to operation in burning plasma devices.

Tungsten impurity transport experiments in Alcator C-Mod to address high priority research and development for ITER

Experiments in Alcator C-Mod tokamak plasmas in the Enhanced D-alpha H-mode regime with ITER-like mid-radius plasma density peaking and Ion Cyclotron Resonant heating, in which tungsten is introduced by the laser blow-off technique, have demonstrated that accumulation of tungsten in the central region of the plasma does not take place in these conditions. The measurements obtained are consistent with anomalous transport dominating tungsten transport except in the central region of the plasma where tungsten transport is neoclassical, as previously observed in other devices with dominant neutral beam injection heating, such as JET and ASDEX Upgrade. In contrast to such results, however, the measured scale lengths for plasma temperature and density in the central region of these Alcator C-Mod plasmas, with density profiles relatively flat in the core region due to the lack of core fuelling, are favourable to prevent inter and intra sawtooth tungsten accumulation in this region under dominance of neoclassical...

Inboard and outboard radial electric field wells in the H- and I-mode pedestal of Alcator C-Mod and poloidal variations of impurity temperature

We present inboard (HFS) and outboard (LFS) radial electric field (Er) and impurity temperature (Tz) measurements in the I-mode and H-mode pedestal of Alcator C-Mod. These measurements reveal strong Er wells at the HFS and the LFS midplane in both regimes and clear pedestals in Tz, which are of similar shape and height for the HFS and LFS. While the H-mode Er well has a radially symmetric structure, the Er well in I-mode is asymmetric, with a stronger ExB shear layer at the outer edge of the Er well, near the separatrix. Comparison of HFS and LFS profiles indicates that impurity temperature and plasma potential are not simultaneously flux functions. Uncertainties in radial alignment after mapping HFS measurements along flux surfaces to the LFS do not, however, allow direct determination as to which quantity varies poloidally and to what extent. Radially aligning HFS and LFS measurements based on the Tz profiles would result in substantial inboard-outboard variations of plasma potential and electron density. Aligning HFS and LFS Er wells instead also approximately aligns the impurity poloidal flow profiles, while resulting in a LFS impurity temperature exceeding the HFS values in the region of steepest gradients by up to 70%. Considerations based on a simplified form of total parallel momentum balance and estimates of parallel and perpendicular heat transport time scales seem to favor an approximate alignment of the Er wells and a substantial poloidal asymmetry in impurity temperature.

Effects of LHRF on toroidal rotation in Alcator C-Mod plasmas

Application of lower hybrid range of frequencies (LHRF) waves can induce both co- and counter-current directed changes in toroidal rotation in Alcator C-Mod plasmas, depending on the target plasma current, electron density, confinement regime and magnetic shear. For ohmic L-mode discharges with good core LH wave absorption, and significant current drive at a fixed LH power near 0.8 MW, the interior (r/a q95/11.5, and in the co-current direction if ne(1020 m−3) 1, indicating a good correlation with driven current fraction, unifying the results observed on various tokamaks. For high density (ne ≥ 1.2 × 1020 m−3) L-mode target discharges, where core LH wave absorption is low, the rotation change is in the co-current direction, but evolves on a shorter momentum transport time scale, and is seen across the entire spatial profile. For H-mode target plasmas, both co- and counter-current direction increments have been observed with LHRF. The H-mode co-rotation is correlated with the pedestal temperature gradient, which itself is enhanced by the LH waves absorbed in the plasma periphery. The H-mode counter-rotation increment, a flattening of the peaked velocity profile in the core, is consistent with a reduction in the momentum pinch correlated with a steepening of the core density profile. Most of these rotation changes must be due to indirect transport effects of LH waves on various parameters, which modify the momentum flux.

Characterization of the Pedestal in Alcator C-Mod ELMing H-Modes and Comparison to the EPED Model

A dedicated series of ELMing H-mode discharges on Alcator C-Mod spanning a broad range of plasma parameters, including plasma current (400–1000 kA), magnetic field (3.5–8 T), and plasma shaping, are presented with experimental scalings of the plasma pedestal with bulk plasma and engineering parameters. The H-modes presented achieve pedestals with densities spanning 5 × 1019–2.5 × 1020 m−3 and temperatures of 150–1000 eV (corresponding to 5–40 kPa in the pressure pedestal), over a width of 3–5% of poloidal flux. The observed pedestal structure is compared with the most recent iteration of the EPED class of models, which uniquely predict the pedestal width and height for a set of scalar input parameters via a combination of stability calculations for peeling–ballooning MHD modes and kinetic-ballooning modes.