Search

The LLAMA (Lyman Alpha Measurement Apparatus) pinhole camera diagnostic had previously been deployed on DIII-D to measure radial profiles of the Lyman-α (Ly-α) deuterium neutral line brightness across the plasma boundary in the lower chamber to infer neutral deuterium density and ionization rate profiles. This system has recently been upgraded with a new diagnostic head, named ALPACA, that also encloses two pinhole cameras and duplicates the LLAMA views in the upper chamber. Similar to LLAMA, ALPACA provides two times 20 lines of sight, viewing the plasma edge on the inboard and outboard sides with a radial resolution of ∼ 2.5 cm (FWHM) and an effective time resolution of ∼ 1 ms that allows for the investigation of inter-ELM dynamics. The extended Ly-α system provides better coverage to study neutrals in experiments with various plasma shapes utilizing both the upper and lower divertors. Furthermore, post-campaign calibration of the LLAMA diagnostic has successfully been demonstrated for the first time. This was facilitated by various upgrades to the calibration set-up and detailed measurements of the emissivity distribution of the Ly-α calibration source using a pinhole collimator. It was found that the sensitivity of the inboard LLAMA pinhole camera was reduced by a factor of 2.0 ± 0.2 over the course of six months of plasma operation in 2021. The upgraded Ly-α system, equipped with improved absolute calibration, will provide key input for neutral fueling and pedestal particle transport studies and for 2D edge transport code validation on the DIII-D tokamak. [ABSTRACT FROM AUTHOR]

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions. [ABSTRACT FROM AUTHOR]

We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments. [ABSTRACT FROM AUTHOR]

The recently developed integrated model based on engineering parameters (IMEP) (Luda et al 2020 Nucl. Fusion 61 126048; Luda et al 2021 Nucl. Fusion 60 036023), so far validated on ASDEX Upgrade, has been tested on a database of 3 Alcator C-Mod and 55 JET-ILW ELMy (type I) H-mode stationary phases. The empirical pedestal transport model included in IMEP, consisting now of imposing a fixed value of R < ∇ T e > / T e , t o p = − 82.5 , allows an accurate prediction of the pedestal top temperature (when the pedestal top density is fixed to the experimental measurements) across these three machines with different sizes, when the pedestal is peeling–ballooning (PB) limited. Cases far from the ideal PB boundary, corresponding to high edge Spitzer resistivity, are instead strongly overpredicted by IMEP. A comparison between the predictions of Europed and IMEP for a subset of JET-ILW cases shows that IMEP can more accurately reproduce the experimental pedestal width. This allows IMEP to better capture profile effects on the pedestal stability, and therefore to correctly describe the negative effect of fueling on the pedestal pressure for PB limited cases. A strong correlation between the separatrix density and the fueling rate has been identified for a subset of JET-ILW cases, when taking into account different divertor configurations. Overall, these promising results encourage further developments of integrated models to obtain reliable predictions of pedestal and global confinement using only engineering parameters for present and future machines. [ABSTRACT FROM AUTHOR]

Electromagnetic pickup noise in the tokamak environment imposes an imminent challenge for measuring weak diagnostic photocurrents in the nA range. The diagnostic signal can be contaminated by an unknown mixture of crosstalk signals from coils powered by currents in the kA range. To address this issue, an algorithm for robust identification of linear multi-input single-output (MISO) systems has been developed. The MISO model describes the dynamic relationship between measured signals from power sources and observed signals in the diagnostic and allows for a precise subtraction of the noise component. The proposed method was tested on experimental diagnostic data from the DIII-D tokamak, and it has reduced noise by up to 20 dB in the 1–20 kHz range. [ABSTRACT FROM AUTHOR]

The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the 33D − 23P transition in atomic helium with unbounded correlation constraints between the electron density and temperature. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based on fast camera data on the Alcator C-Mod tokamak, we present the first two-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density, revealing shadowing effects in a fusion plasma using a single spectral line. [ABSTRACT FROM AUTHOR]

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modeling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models. [ABSTRACT FROM AUTHOR]

Emission lines which appear in the spectral ranges of ground state transitions from n = 2 levels in He- and H-like argon ions are discussed. X-ray transitions from elements commonly found in tokamaks (tungsten, molybdenum, iron and sulphur) which radiate in the wavelength range from 3700–4000 mĹ are identified by comparison with atomic structure calculations. Individual lines from tungsten charge states in the vicinity of Zn-like W44+ are documented, along with B-like Mo37+. The behaviour of line ratios as a function of electron temperature is examined, in support of the identifications. [ABSTRACT FROM AUTHOR]

A new pseudolocal tomography algorithm is developed for soft X-ray(SXR) imaging measurements of the turbulent electron temperature fluctuations (δ Te) in tokamaks and stellarators. The algorithm overcomes the constraints of limited viewing ports on the vessel wall (viewing angle) and limited number of lines of sight (LOS). This is accomplished by increasing the number of LOS locally in a region of interest. Numerical modeling demonstrates that the wavenumber spectrum of the turbulence can be reliably reconstructed, with an acceptable number of viewing angles and LOS and suitable low SNR detectors. We conclude that a SXR imaging diagnostic for measurements of turbulent δ Te using a pseudolocal reconstruction algorithm is feasible. [ABSTRACT FROM AUTHOR]

A one dimensional, absolutely calibrated pinhole camera system was installed on the DIII-D tokamak to measure edge Lyman-alpha (Ly-α) emission from hydrogen isotopes, which can be used to infer neutral density and ionization rate profiles. The system is composed of two cameras, each providing a toroidal fan of 20 lines of sight, viewing the plasma edge on the inboard and outboard side of DIII-D. The cameras' views lie in a horizontal plane 77 cm below the midplane. At its tangency radius, each channel provides a radial resolution of ∼2 cm full width at half maximum (FWHM) with a total coverage of 22 cm. Each camera consists of a rectangular pinhole, Ly-α reflective mirror, narrow-band Ly-α transmission filter, and a 20 channel AXUV photodetector. The combined mirror and transmission filter have a FWHM of 5 nm, centered near the Ly-α wavelength of 121.6 nm and is capable of rejecting significant, parasitic carbon-III (C-III) emission from intrinsic plasma impurities. To provide a high spatial resolution measurement in a compact footprint, the camera utilizes advanced engineering and manufacturing techniques including 3D printing, high stability mirror mounts, and a novel alignment procedure. Absolutely calibrated, spatially resolved Ly-α brightness measurements utilize a bright, isolated line with low parasitic surface reflections and enable quantitative comparison to modeling to study divertor neutral leakage, main chamber fueling, and radial particle transport. [ABSTRACT FROM AUTHOR]

The LLAMA (Lyman-Alpha Measurement Apparatus) diagnostic was recently installed on the DIII-D tokamak [Rosenthal et al., Rev. Sci. Instrum. (submitted) (2020)]. LLAMA is a pinhole camera system with a narrow band Bragg mirror, a bandpass interference filter, and an absolute extreme ultraviolet photodiode detector array, which measures the Ly-α brightness in the toroidal direction on the inboard, high field side (HFS) and outboard, low field side (LFS). This contribution presents a setup and a procedure for an absolute calibration near the Ly-α line at 121.6 nm. The LLAMA in-vacuum components are designed as a compact, transferable setup that can be mounted in an ex situ vacuum enclosure that is equipped with an absolutely calibrated Ly-α source. The spectral purity and stability of the Ly-α source are characterized using a vacuum ultraviolet spectrometer, while the Ly-α source brightness is measured by a NIST-calibrated photodiode. The non-uniform nature of the Ly-α source emission was overcome by performing a calibration procedure that scans the Ly-α source position and employs a numerical optimization to determine the emission pattern. Nominal and measured calibration factors are determined and compared, showing agreement within their uncertainties. A first conversion of the measured signal obtained from DIII-D indicates that the Ly-α brightness on the HFS and LFS is on the order of 1020 Ph sr−1 m−2 s−1. The established calibration setup and procedure will be regularly used to re-calibrate the LLAMA during DIII-D vents to monitor possible degradation of optical components and detectors. [ABSTRACT FROM AUTHOR]

Analysis and modeling of rotation reversal hysteresis experiments show that a single turbulent bifurcation is responsible for the Linear to Saturated Ohmic Confinement (LOC/SOC) transition and concomitant intrinsic rotation reversal on Alcator C-Mod. Plasmas on either side of the reversal exhibit different toroidal rotation profiles and therefore different turbulence characteristics despite the profiles of density and temperature, which are indistinguishable within measurement uncertainty. Elements of this bifurcation are also shown to persist for auxiliary heated L-modes. The deactivation of subdominant (in the linear growth rate and contribution to heat transport) ion temperature gradient and trapped electron mode instabilities is identified as the only possible change in turbulence within a reduced quasilinear transport model across the reversal, which is consistent with the measured profiles and inferred heat and particle fluxes. Experimental constraints on a possible change from strong to weak turbulence, outside the description of the quasilinear model, are also discussed. These results indicate an explanation for the LOC/SOC transition that provides a mechanism for the hysteresis through the dynamics of subdominant modes and changes in their relative populations and does not involve a change in the most linearly unstable ion-scale drift-wave instability. [ABSTRACT FROM AUTHOR]

The fluid transport code [trans-electric field (Er) module] under the BOUT++ framework has been used to simulate divertor heat flux width and boundary Er with all drifts and the sheath potential in the scrape-off layer. The calculated steady state radial Er in the pedestal region has been compared with that of experimental measurements from the Alcator C-Mod tokamak. The magnitude and shape of Er are similar to those of the experimental data. In order to understand the relative role of cross-field drifts vs turbulent transport in setting the heat flux width, four C-Mod enhanced Dα H-mode discharges with a lower single null divertor configuration should be simulated. BOUT++ transport simulations with cross-field drifts included yield similar heat flux width λq to that of experimental measurements (within a factor of 2) from both the probe and the surface thermocouple diagnostics and show a similar trend with plasma current to that of the Eich experimental scaling. The simulations show that both drifts and turbulent transport compete to determine the heat flux width. The magnetic drifts play a dominant role in setting the divertor heat-flux width, while the E × B drift decreases the heat flux width by 10%–25%, leading to improved agreement with the experiment relative to Goldston's model. A turbulence diffusivity scan (χ⊥) identifies two distinct regimes: a drift dominant regime when χ⊥ is small and a turbulence dominant regime when χ⊥ is large. The Goldston heuristic drift model yields a lower limit of the width λq. [ABSTRACT FROM AUTHOR]

The heat flux distributions on divertor targets in H-mode plasmas are serious concerns for future devices. We seek to simulate the tokamak boundary plasma turbulence and heat transport in the edge localized mode-suppressed regimes. The improved BOUT++ model shows that not only Ip but also the radial electric field Er plays an important role on the turbulence behavior and sets the heat flux width. Instead of calculating Er from the pressure gradient term (diamagnetic Er), it is calculated from the plasma transport equations with the sheath potential in the scrape-off layer and the plasma density and temperature profiles inside the separatrix from the experiment. The simulation results with the new Er model have better agreement with the experiment than using the diamagnetic Er model: (1) The electromagnetic turbulence in enhanced Dα H-mode shows the characteristics of quasi-coherent modes (QCMs) and broadband turbulence. The mode spectra are in agreement with the phase contrast imaging data and almost has no change in comparison to the cases which use the diamagnetic Er model; (2) the self-consistent boundary Er is needed for the turbulence simulations to get the consistent heat flux width with the experiment; (3) the frequencies of the QCMs are proportional to Er, while the divertor heat flux widths are inversely proportional to Er; and (4) the BOUT++ turbulence simulations yield a similar heat flux width to the experimental Eich scaling law and the prediction from the Goldston heuristic drift model. [ABSTRACT FROM AUTHOR]

A profile for the critical gradient scale length (Lc) has been measured in L-mode discharges at the Alcator C-Mod tokamak, where electrons were heated by an ion cyclotron range of frequency through minority heating with the intention of simultaneously varying the heat flux and changing the local gradient. The electron temperature gradient scale length (LTe−1 = |∇Te|/Te) profile was measured via the BT-jog technique [Houshmandyar et al., Rev. Sci. Instrum. 87, 11E101 (2016)] and it was compared with electron heat flux from power balance (TRANSP) analysis. The Te profiles were found to be very stiff and already above the critical values, however, the stiffness was found to be reduced near the q = 3/2 surface. The measured Lc profile is in agreement with electron temperature gradient (ETG) models which predict the dependence of Lc−1 on local Zeff, Te/Ti, and the ratio of the magnetic shear to the safety factor. The results from linear Gene gyrokinetic simulations suggest ETG to be the dominant mode of turbulence in the electron scale (k⊥ρs > 1), and ion temperature gradient/trapped electron mode modes in the ion scale (k⊥ρs < 1). The measured Lc profile is in agreement with the profile of ETG critical gradients deduced from Gene simulations. [ABSTRACT FROM AUTHOR]

New validation of global, nonlinear, ion-scale gyrokinetic simulations (GYRO) is carried out for L- and I-mode plasmas on Alcator C-Mod, utilizing heat fluxes, profile stiffness, and temperature fluctuations. Previous work at C-Mod found that ITG/TEM-scale GYRO simulations can match both electron and ion heat fluxes within error bars in I-mode [White PoP 2015], suggesting that multi-scale (cross-scale coupling) effects [Howard PoP 2016] may be less important in I-mode than in L-mode. New results presented here, however, show that global, nonlinear, ion-scale GYRO simulations are able to match the experimental ion heat flux, but underpredict electron heat flux (at most radii), electron temperature fluctuations, and perturbative thermal diffusivity in both L- and I-mode. Linear addition of electron heat flux from electron scale runs does not resolve this discrepancy. These results indicate that single-scale simulations do not sufficiently describe the I-mode core transport, and that multi-scale (coupled electron- and ion-scale) transport models are needed. A preliminary investigation with multi-scale TGLF, however, was unable to resolve the discrepancy between ion-scale GYRO and experimental electron heat fluxes and perturbative diffusivity, motivating further work with multi-scale GYRO simulations and a more comprehensive study with multi-scale TGLF. [ABSTRACT FROM AUTHOR]

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 ve at the edge is in the range of experimental errors of veff from 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. [ABSTRACT FROM AUTHOR]

Two-dimensional scrape-off layer turbulence (SOLT) code simulations are compared with an L-mode discharge on the Alcator C-Mod tokamak [Greenwald et al., Phys. Plasmas 21, 110501 (2014)]. Density and temperature profiles for the simulations were obtained by smoothly fitting Thomson scattering and mirror Langmuir probe (MLP) data from the shot. Simulations differing in turbulence intensity were obtained by varying a dissipation parameter. Mean flow profiles and density fluctuation amplitudes are consistent with those measured by MLP in the experiment and with a Fourier space diagnostic designed to measure poloidal phase velocity. Blob velocities in the simulations were determined from the correlation function for density fluctuations, as in the analysis of gas-puff-imaging (GPI) blobs in the experiment. In the simulations, it was found that larger blobs moved poloidally with the E × B flow velocity, vE, in the near-SOL, while smaller fluctuations moved with the group velocity of the dominant linear (interchange) mode, vE+1/2 vdi, where vdi is the ion diamagnetic drift velocity. Comparisons are made with the measured GPI correlation velocity for the discharge. The saturation mechanisms operative in the simulation of the discharge are also discussed. It is found that neither sheared flow nor pressure gradient modification can be excluded as saturation mechanisms. [ABSTRACT FROM AUTHOR]

For the first time, the power deposition of lower hybrid RF waves into the edge plasma of a diverted tokamak has been systematically quantified. Edge deposition represents a parasitic loss of power that can greatly impact the use and efficiency of Lower Hybrid Current Drive (LHCD) at reactor-relevant densities. Through the use of a unique set of fast time resolution edge diagnostics, including innovative fast-thermocouples, an extensive set of Langmuir probes, and a Lyα ionization camera, the toroidal, poloidal, and radial structure of the power deposition has been simultaneously determined. Power modulation was used to directly isolate the RF effects due to the prompt (t < τE) response of the scrape-off-layer (SOL) plasma to Lower Hybrid Radiofrequency (LHRF) power. LHRF power was found to absorb more strongly in the edge at higher densities. It is found that a majority of this edge-deposited power is promptly conducted to the divertor. This correlates with the loss of current drive efficiency at high density previously observed on Alcator C-Mod, and displaying characteristics that contrast with the local RF edge absorption seen on other tokamaks. Measurements of ionization in the active divertor show dramatic changes due to LHRF power, implying that divertor region can be a key for the LHRF edge power deposition physics. These observations support the existence of a loss mechanism near the edge for LHRF at high density (ne > 1:0 × 1020 (m -3)). Results will be shown addressing the distribution of power within the SOL, including the toroidal symmetry and radial distribution. These characteristics are important for deducing the cause of the reduced LHCD efficiency at high density and motivate the tailoring of wave propagation to minimize SOL interaction, for example, through the use of high-field-side launch. [ABSTRACT FROM AUTHOR]

The high-fidelity BOUT++ two-fluid code suite has demonstrated significant recent progress toward integrated multi-scale simulations of tokamak pedestal, including Edge-Localized-Mode (ELM) dynamics, evolution of ELM cycles, and continuous fluctuations, as observed in experiments. Nonlinear ELM simulations show three stages of an ELM event: (1) a linear growing phase; (2) a fast crash phase; and (3) a slow inward turbulence spreading phase lasting until the core heating flux balances the ELM energy loss and the ELM is terminated. A new coupling/splitting model has been developed to perform simulations of multi-scale ELM dynamics. Simulation tracks five ELM cycles for 10 000 Alfvén times for small ELMs. The temporal evolution of the pedestal pressure is similar to that of experimental measurements for the pedestal pressure profile collapses and recovers to a steep gradient during ELM cycles. To validate BOUT++ simulations against experimental data and develop physics understanding of the fluctuation characteristics for different tokamak operation regimes, both quasi-coherent fluctuations (QCFs) in ELMy H-modes and Weakly Coherent Modes in I-modes have been simulated using three dimensional 6-field 2-fluid electromagnetic model. The H-mode simulation results show that (1) QCFs are localized in the pedestal region having a predominant frequency at f ≃ 300 - 400 kHz and poloidal wavenumber at kθ 0.7 cm-1, and propagate in the electron diamagnetic direction in the laboratory frame. The overall signatures of simulation results for QCFs show good agreement with C-Mod and DIII-D measurements. (2) The pedestal profiles giving rise to QCFs are near the marginal instability threshold for ideal peeling-ballooning modes for both C-Mod and DIII-D, while the collisional electromagnetic drift-Alfvén wave appears to be dominant for DIII-D. (3) Particle diffusivity is either smaller than the heat diffusivity for DIIID or similar to the heat diffusivity for C-Mod. Key I-mode simulation results are that (1) a strong instability exists at n ⩾ 20 for resistive ballooning mode and drift-Alfvén wave; (2) the frequency spectrum of nonlinear BOUT++ simulation features a peak around 300 kHz for the mode number n = 20, consistent with a reflectometer measurement at nearby position; (3) the calculated particle diffusivity is larger than the calculated heat diffusivity, which is consistent with a key feature of the I-mode pedestal with no particle barrier. [ABSTRACT FROM AUTHOR]

Peaked density profiles in low-collisionality AUG and JET H-mode plasmas are probably caused by a turbulently driven particle pinch, and Alcator C-Mod experiments confirmed that collisionality is a critical parameter. Density peaking in reactors could produce a number of important effects, some beneficial, such as enhanced fusion power and transport of fuel ions from the edge to the core, while others are undesirable, such as lower beta limits, reduced radiation from the plasma edge, and consequently higher divertor heat loads. Fundamental understanding of the pinch will enable planning to optimize these impacts. We show that density peaking is predicted by nonlinear gyrokinetic turbulence simulations based on measured profile data from low collisionality H-mode plasma in Alcator C-Mod. Multiple ion species are included to determine whether hydrogenic density peaking has an isotope dependence or is influenced by typical levels of low-Z impurities, and whether impurity density peaking depends on the species. We find that the deuterium density profile is slightly more peaked than that of hydrogen, and that experimentally relevant levels of boron have no appreciable effect on hydrogenic density peaking. The ratio of density at r/a0.44 to that at r/a0.74 is 1.2 for the majority D and minority H ions (and for electrons), and increases with impurity Z: 1.1 for helium, 1.15 for boron, 1.3 for neon, 1.4 for argon, and 1.5 for molybdenum. The ion temperature profile is varied to match better the predicted heat flux with the experimental transport analysis, but the resulting factor of two change in heat transport has only a weak effect on the predicted density peaking. [ABSTRACT FROM AUTHOR]

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 transport. Simulations of ITER H-mode plasmas, including both anomalous (modelled by the Gyro-Landau-Fluid code GLF23) and neoclassical transport for main ions and tungsten and with density profiles of similar peaking to those obtained in Alcator C-Mod show that accumulation of tungsten in the central plasma region is also unlikely to occur in stationary ITER H-mode plasmas due to the low fuelling source by the neutral beam injection (injection energy ~ 1 MeV), which is in good agreement with findings in the Alcator C-Mod experiments. [ABSTRACT FROM AUTHOR]

Injecting Lower Hybrid Range of Frequency (LHRF) waves into Alcator C-Mod's high-density H-mode plasmas has led to enhanced global energy confinement by increasing pedestal temperature and pressure gradients, decreasing the separatrix density, modifying the pedestal radial electric field and rotation, and decreasing edge turbulence. These experiments indicate that edge LHRF can be used as an actuator to increase energy confinement via modification of boundary quantities. H98-factor increases of up to ~35% (e.g., H98 from 0.75 to 1.0) are seen when moderate amounts of LH power (PLH/Ptot~0.15) are applied to H-modes of densities e ~3 × 1020 m−3, corresponding to values ~0.5 of the Greenwald density. However, the magnitude of the improvement is reduced if the confinement quality of the target H-mode plasma is already good (i.e., H98 target ~ 1). Ray-tracing modeling and accessibility calculations for the LH waves indicate that they do not penetrate to the core. The LHRF power appears to be deposited in plasma boundary region, with a large fraction of the injected power increment appearing promptly on the outer divertor target. There is no evidence that the LH waves are driving current in these plasmas. The LHRF-actuated improvements are well correlated with suppressed pedestal density fluctuations in the 100–300kHz range. There is also a correlation between the improved confinement and a drop in separatrix density, a correlation that is consistent with previous H-mode results with no LHRF. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

A new "Mirror Langmuir Probe" diagnostic, combined with a double-coil scanning magnetic probe, is used to interrogate Alcator C-Mod's quasi-coherent mode (QCM) with unprecedented detail. In ohmic EDA H-modes, the QCM is found to reside in a region of positive radial electric field, with a radial width (~3 mm) that spans open and closed field line regions. Large amplitude, in-phase sinusoidal bursts (~100 kHz) in density, electron temperature, and plasma potential are observed, with potential lagging density by ~16°, producing an outward radial transport velocity of ~10 m/s. Mode propagation corresponds to the sum of local E × B and electron diamagnetic drift velocities. Poloidal magnetic field fluctuations project to current filaments carrying peak current densities of ~25 A/cm2. An evaluation of parallel electron force balance (Ohm's law) over a fluctuation cycle indicates a significant electromotive component. Interchange drive is also a contributor in the current continuity (vorticity) equation. Thus, the QCM is primarily a separatrix-spanning electron drift-wave with interchange and electromagnetic contributions. [ABSTRACT FROM AUTHOR]

I-mode is a high-performance tokamak regime characterized by the formation of a temperature pedestal and enhanced energy confinement, without an accompanying density pedestal or drop in particle and impurity transport. I-mode operation appears to have naturally occurring suppression of large Edge-Localized Modes (ELMs) in addition to its highly favorable scalings of pedestal structure and overall performance. Extensive study of the ELMy H-mode has led to the development of the EPED model, which utilizes calculations of coupled peeling-ballooning MHD modes and kinetic-ballooning mode (KBM) stability limits to predict the pedestal structure preceding an ELM crash. We apply similar tools to the structure and ELM stability of I-mode pedestals. Analysis of I-mode discharges prepared with high-resolution pedestal data from the most recent C-Mod campaign reveals favorable pedestal scalings for extrapolation to large machines-pedestal temperature scales strongly with power per particle Pnet/e, and likewise pedestal pressure scales as the net heating power (consistent with weak degradation of confinement with heating power). Matched discharges in current, field, and shaping demonstrate the decoupling of energy and particle transport in I-mode, increasing fueling to span nearly a factor of two in density while maintaining matched temperature pedestals with consistent levels of Pnet/-e. This is consistent with targets for increased performance in I-mode, elevating pedestal βp and global performance with matched increases in density and heating power. MHD calculations using the ELITE code indicate that I-mode pedestals are strongly stable to edge peeling-ballooning instabilities. Likewise, numerical modeling of the KBM turbulence onset, as well as scalings of the pedestal width with poloidal beta, indicates that I-mode pedestals are not limited by KBM turbulence-both features identified with the trigger for large ELMs, consistent with the observed suppression of large ELMs in I-mode. [ABSTRACT FROM AUTHOR]

X-ray spectra of n = 2 to 1 transitions in hydrogen-like Ca19 +, helium-like Ca18 + and nearby satellites have been obtained from Alcator C-Mod tokamak plasmas using a spatially imaging high resolution x-ray spectrometer system. For Ca19 +, the intensity ratio of Lyα2 (1s 1S1/2–2p 2P1/2) to Lyα1 (1s 1S1/2–2p 2P3/2) was found to be ∼0.531 ± 0.005 over a range of plasma parameters, which is somewhat greater than the ratio of the statistical weights of the upper n = 2 levels, 1/2. This difference is mainly due to interaction with the 2S1/2 fine structure sub-level. Experimental results are compared to calculations from COLRAD, a collisional-radiative modelling code, and good agreement is shown. For Ca18 +, the intensity ratio of the dielectronic satellite k (1s22p 2P1/2–1s2p22D3/2) to the resonance line w (1s21S0–1s2p 1P1) is sensitive only to the electron temperature. The observed brightness ratio scaling with Te is in good agreement with the calculated ratio of the respective dielectronic recombination to the collisional excitation rates. [ABSTRACT FROM AUTHOR]

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.8MW, the interior (r/a < 0.5) rotation increments (on a time scale of order the current relaxation time) in the counter-current direction if ne(1020 m-3) > q95/11.5, and in the co-current direction if ne(1020 m-3) < q95/11.5. All discharges with co-current rotation changes have q0 > 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Experiments on Alcator C-Mod from 2009-2012 have examined the rampup, flattop and rampdown phases of the proposed ITER 15 MA baseline scenario. Rampup studies show ICRF heating can significantly reduce the V-s requirement, and that an H-mode late in the ramp can reduce this further. ICRF modifications to li in L-mode are minimal, although the Te profile is peaked relative to ohmic in the plasma centre, and reduces the sawtooth onset times. Flattop plasmas targeting ITER baseline parameters have been sustained for 20τE or 8 - 13τCR, but only reach H98 ~ 0.6 at n/nGr = 0.85, rising to 0.9 at n/nGr = 0.65. Rampdown studies show H-modes can be routinely sustained with ICRF power injection, avoiding an OH coil over-current associated with the H-L transition. In addition, faster current rampdowns are preferred to avoid an over-current when an H-L transition ultimately does occur. In the H-mode rampdown the density is found to drop with Ip, preserving the n/nGr ratio, so long as ICRF power is injected. [ABSTRACT FROM AUTHOR]

The first measurements of long wavelength (kyps < 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 (p ~ 0.8) where the changes in electron temperature fluctuations are measured. In contrast, deeper in the core (p < 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. [ABSTRACT FROM AUTHOR]

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 fluctuation level ratio (Te/Te)/(ñe/ne), but underprediction of electron heat flux. We find that changes in momentum transport (rotation profiles changing from peaked to hollow) do not correlate with changes in particle transport, and also do not correlate with changes in linear mode dominance, e.g., Ion Temperature Gradient versus Trapped Electron Mode. The new C-Mod results suggest that the drives for momentum transport differ from drives for heat and particle transport. The experimental results are inconsistent with present quasilinear models, and the strong sensitivity of core rotation to density remains unexplained. [ABSTRACT FROM AUTHOR]

Measurements of poloidal variation, ñz/z>, in high-Z impurity density have been made using photodiode arrays sensitive to vacuum ultraviolet and soft x-ray emission in Alcator C-Mod plasmas. In/out asymmetries in the range of -0.2z,cos/z><0.3 are observed for r/a<0.8, and accumulation on both the high-field side, nz,cos<0, and low-field side, nz,cos>0, of a flux surface is found to be well described by a combination of centrifugal, poloidal electric field, and ion-impurity friction effects. Up/down asymmetries, -0.05z,sin/z><0.10, are observed over 0.5z,sin>0 corresponding to accumulation opposite the ion ∇B drift direction. Measurements of the up/down asymmetry of molybdenum are found to disagree with predictions from recent neoclassical theory in the trace limit, nzZ2/ni<1. Non-trace levels of impurities are expected to modify the main-ion poloidal flow and thus change friction-driven impurity density asymmetries and impurity poloidal rotation, vθ,z. Artificially modifying main-ion flow in parallel transport simulations is shown to impact both ñz/z> and vθ,z, but simultaneous agreement between measured and predicted up/down and in/out asymmetry as well as impurity poloidal rotation is not possible for these C-Mod data. This link between poloidal flow and poloidal impurity density variation outlines a more stringent test for parallel neoclassical transport theory than has previously been performed. Measurement and computational techniques specific to the study of poloidal impurity asymmetry physics are discussed as well. [ABSTRACT FROM AUTHOR]

We demonstrate that existing theories are insufficient to explain up/down asymmetries of argon x-ray emission in Alcator C-Mod ohmic plasmas. Instead of the poloidal variation, ñz/z>, being of order the inverse aspect ratio, ∈, and scaling linearly with Btñe/Ip², it is observed over 0.8 < r/a < 1.0 to be of order unity and exhibits a threshold behaviour between 3.5 < Btñe/Ip < 4.0 (T1020 m-3 MA-1). The transition from a poloidally symmetric to asymmetric impurity distribution is shown to occur at densities just below those that trigger a reversal of the core toroidal rotation direction, thought to be linked to the transition between the linear and saturated ohmic confinement regimes. A possible drive is discussed by which anomalous radial transport might sustain the impurity density asymmetry as the ratio of the perpendicular to parallel equilibration times, τ⊥,z/τ||,z, approaches unity. This explanation requires a strong up/down asymmetry in radial flux which, while not observable on C-Mod, has been measured in TEXT and Tore Supra ohmic plasmas. [ABSTRACT FROM AUTHOR]

Several seemingly unrelated effects in Alcator C-Mod ohmic L-mode plasmas are shown to be closely connected: non-local heat transport, core toroidal rotation reversals, energy confinement saturation and up/down impurity density asymmetries. These phenomena all abruptly transform at a critical value of the collisionality. At low densities in the linear ohmic confinement regime, with collisionality v* ≤ 0.35 (evaluated inside of the q = 3/2 surface), heat transport exhibits non-local behaviour, core toroidal rotation is directed co-current, edge impurity density profiles are up/down symmetric and a turbulent feature in core density fluctuations with kθ up to 15 cm-1 (kθρs ~ 1) is present. At high density/collisionality with saturated ohmic confinement, electron thermal transport is diffusive, core rotation is in the counter-current direction, edge impurity density profiles are up/down asymmetric and the high kθ turbulent feature is absent. The rotation reversal stagnation point (just inside of the q = 3/2 surface) coincides with the non-local electron temperature profile inversion radius. All of these observations suggest a possible unification in a model with trapped electron mode prevalence at low collisionality and ion temperature gradient mode domination at high collisionality. [ABSTRACT FROM AUTHOR]

The physical mechanisms that cause tokamak plasmas to rotate toroidally without external momentum input are of considerable interest to the plasma physics community. This paper documents a substantial change in both the magnitude of the core-rotation frequency, -1 < ω(r/a = 0) < +10 kHz, and the sign of rotation shear at mid-radius, u' = -R² dω/dr/νth,i, which varies in the range -0.6 < u' < +0.8 in response to very small changes in the electron density. In 0.8 MA, 5.4 T Alcator C-Mod L-mode plasmas using 1.2MW of on-axis ion-cyclotron resonance heating, plasmas with line-averaged densities in the range 1.0 < ne < 1.2 × 1020 m-3 exhibit a transition from a peaked intrinsic rotation profile to one that is hollow. Gradient scale lengths of the temperature and density profiles, the drive for plasma turbulence thought to play a role in intrinsic rotation, are indistinguishable within experimental uncertainties between the plasmas, and linear stability analysis using GYRO shows the plasmas to be in the ion temperature gradient-dominated turbulence regime. The impact of changes in the rotation profile in response to minor changes under target plasma conditions is discussed in relation to established analysis techniques and cross-machine rotation scaling studies, with comparisons made with existing ASDEX-Upgrade work on intrinsic rotation shear. [ABSTRACT FROM AUTHOR]

Transitions from the L-mode regime to the I-mode regime, with an energy transport barrier, and to the H-mode regime with both an energy and particle transport barrier are studied on the Alcator C-Mod tokamak. Steady I-mode plasmas have been produced over a wide range of plasma field (3-6 T), current (0.8-1.35 MA), density and shaping in the unfavourable ion B×B configuration. The power threshold for the L-I transition is higher than scalings for the L-H transition with favourable drift, and increases with plasma current as well as density. Threshold conditions for the I-H transition are more variable. In some conditions I-mode is maintained up to the maximum available ICRF power of 5MW, nearly a factor of two above the L-I threshold, giving a robust operating window. Edge Te at the L-I transition is in the range 250-450 eV, over a range of current and density, about a factor of two higher than with favourable drift, while at the I-H transition it can be much higher (up to 1.1 keV) but is again widely variable. Heat pulses due to sawteeth may play a role in transitions. Controlled I-L back transitions indicate little power hysteresis. [ABSTRACT FROM AUTHOR]

This research describes advancements in the spectral analysis and error propagation techniques associated with x-ray imaging crystal spectroscopy (XICS) that have enabled this diagnostic to be used to accurately constrain particle, momentum, and heat transport studies in a tokamak for the first time. Doppler tomography techniques have been extended to include propagation of statistical uncertainty due to photon noise, the effect of non-uniform instrumental broadening as well as flux surface variations in impurity density. These methods have been deployed as a suite of modeling and analysis tools, written in interactive data language (IDL) and designed for general use on tokamaks. Its application to the Alcator C-Mod XICS is discussed, along with novel spectral and spatial calibration techniques. Example ion temperature and radial electric field profiles from recent I-mode plasmas are shown, and the impact of poloidally asymmetric impurity density and natural line broadening is discussed in the context of the planned ITER x-ray crystal spectrometer. [ABSTRACT FROM AUTHOR]

H-mode power thresholds (Pth) are examined in Alcator C-Mod lower single null plasmas with standard vertical-plate and slot-type divertor configuration (BT = 5.4T, Ip = 0.9 MA, q95 sim;4.0 and = (1.3-1.6) × 1020 m-3). The required power to access H-mode is found to be significantly reduced when the plasma is operated with a slot divertor. The lowest Pth achieved is ∼0.7MW, about 1/3 of the normal level and 40% of the multi-machine scaling law prediction. Pth does not explicitly correlate with plasma shaping parameters (X-point locations, elongation, triangularities, etc), but somehow shows an overall inverse correlation with outer divertor leg length (the poloidal distance between X-point and outer strike point) and/or the outer divertor connection length. This implies that another, still unknown parameter hidden in the vicinity of the separatrix outer leg might be responsible for the observed difference in H-mode access power. Despite the large differences in Pth,the edge electron temperature (Te) and density (ne) profiles prior to L-H transition are not strongly affected by the divertor configuration. Values of Te,95 (i.e. Te at 95% normalized poloidal flux surface) at the L-H transition are similar for the two divertor configurations [ABSTRACT FROM AUTHOR]

This paper reports on newexperiments at AlcatorC-Modthat were performed in order to investigate the long-standing, unresolved discrepancy between Thomson scattering (TS) and electron cyclotron emission (ECE) measurements of electron temperature in high temperature tokamak plasmas. Ion cyclotron range of frequency (ICRF) heating is used to produce high temperature conditions where the type of TS-ECE discrepancy observed in the past at JET and TFTR should become observable. At Alcator C-Mod, plasmas with Te(0) up to 8 keV are obtained using ion cyclotron resonance heating (ICRH), ICRF mode conversion heating and a combination of the two heating methods in order to explore the hypothesis that the presence of ICRH-generated fast ions may be related to the discrepancy. In all high temperature cases, the TS and ECE measurements of electron temperature agree to within experimental uncertainties. We find no evidence for the type of discrepancy reported at JET and TFTR. These results show that the TS-ECE discrepancy does not depend on high temperatures alone and also that the presence of ICRH-generated fast ions is insufficient to cause the TS-ECE discrepancy. [ABSTRACT FROM AUTHOR]