Engineering Papers⌕ Search

Engineering topics

Greenwald, M.

Publications and source records attributed to Greenwald, M..

At least 37 records · Page 2

Alcator C-Mod: research in support of ITER and steps beyond

This paper presents an overview of recent highlights from research on Alcator C-Mod. Significant progress has been made across all research areas over the last two years, with particular emphasis on divertor physics and power handling, plasma–material interaction studies, edge localized mode-suppressed pedestal dynamics, core transport and turbulence, and RF heating and current drive utilizing ion cyclotron and lower hybrid tools. Specific results of particular relevance to ITER include: inner wall SOL transport studies that have led, together with results from other experiments, to the change of the detailed shape of the inner wall in ITER; runaway electron studies showing that the critical electric field required for runaway generation is much higher than predicted from collisional theory; core tungsten impurity transport studies reveal that tungsten accumulation is naturally avoided in typical C-Mod conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multispecies density peaking in gyrokinetic turbulence simulations of low collisionality Alcator C-Mod plasmas

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/a = 0.44 to that at r/a = 0.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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved confinement in high-density H-modes via modification of the plasma boundary with lower hybrid wavesa)

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 imply that edge LHRF can be used as an actuator to increase energy confinement via modification of boundary quantities. H 98 -factor increases of up to ~35% (e.g., H 98 from 0.75 to 1.0) are seen when moderate amounts of LH power (P LH /P tot ~0.15) are applied to H-modes of densities $\overline{n}$ e ~3 × 10 20 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., H 98 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–300 kHz range. Additionally, there is 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-resonant destabilization of (1/1) internal kink mode by suprathermal electron pressure

New experimental observations are reported on the structure and dynamics of short-lived periodic (1, 1) “fishbone”-like oscillations that appear during radio frequency heating and current-drive experiments in tokamak plasmas. For the first time, measurements can directly relate changes in the high energy electrons to the mode onset, saturation, and damping. In the relatively high collisionality of Alcator C-Mod with lower hybrid current drive, the instability appears to be destabilized by the non-resonant suprathermal electron pressure—rather than by wave-particle resonance, rotates toroidally with the plasma and grows independently of the (1, 1) sawtooth crash driven by the thermal plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ADX: a high field, high power density, advanced divertor and RF tokamak

The MIT Plasma Science and Fusion Center and collaborators are proposing a high-performance Advanced Divertor and RF tokamak eXperiment (ADX) – a tokamak specifically designed to address critical gaps in the world fusion research program on the pathway to next step devices: fusion nuclear science facility (FNSF), fusion pilot plant (FPP), and/or demonstration power plant (DEMO). This high field (≥ 6.5 tesla, 1.5 MA), high power density facility (P/S ~ 1.5 MW/m2) will test innovative divertor ideas, including an ‘X-point target divertor’ concept, at the required performance parameters – reactor-level boundary plasma pressures, magnetic field strengths and parallel heat flux densities entering into the divertor region – while simultaneously producing high performance core plasma conditions that are prototypical of a reactor: equilibrated electrons and ions, regimes with low or no torque, and no fueling from external heating and current drive systems. Equally important, the experimental platform will test innovative concepts for lower hybrid current drive (LHCD) and ion-cyclotron range of frequency (ICRF) actuators with the unprecedented ability to deploy launch structures both on the lowmagnetic- field side and the high-magnetic-field side – the latter being a location where energetic plasma-material interactions can be controlled and favorable RF wave physics leads to efficient current drive, current profile control, heating and flow drive. This triple combination – advanced divertors, advanced RF actuators, reactor-prototypical core plasma conditions – will enable ADX to explore enhanced core confinement physics, such as made possible by reversed central shear, using only the types of external drive systems that are considered viable for a fusion power plant. Such an integrated demonstration of high-performance core-divertor operation with steady state sustainment would pave the way toward an attractive pilot plant, as envisioned in the ARC concept (Affordable, Robust, Compact) [B. N. Sorbom, et al., submitted to Fusion Engineering Design, 2014] that makes use of high-temperature superconductor technology – a high-field (9.25 tesla) tokamak the size of the Joint European Torus that produces 270 MW of net electricity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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. Furthermore, 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlation ECE diagnostic in Alcator C-Mod

Correlation ECE (CECE) is a diagnostic technique that allows measurement of small amplitude electron temperature, Te, 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 Te 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]. Te 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 Te fluctuations may be applied to both ohmic and RF heated L-mode plasmas. In I-mode plasmas, we observe the reduction of core Te 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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, yet 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 demonstrate 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. Moreover, we 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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. Lastly, a summary of important achievements and discoveries are included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synergistic cross-scale coupling of turbulence in a tokamak plasma

For the first time, nonlinear gyrokinetic simulations spanning both the ion and electron spatio-temporal scales have been performed with realistic electron mass ratio ((m D /m e ) 1/2 = 60.0), realistic geometry, and all experimental inputs, demonstrating the coexistence and synergy of ion ( ( k θ ρ s ∼ O ( 1.0 ) ) ) and electron-scale ( ( k θ ρ e ∼ O ( 1.0 ) ) ) turbulence in the core of a tokamak plasma. All multi-scale simulations utilized the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] to study the coupling of ion and electron-scale turbulence in the core (r/a = 0.6) of an Alcator C-Mod L-mode discharge shown previously to exhibit an under-prediction of the electron heat flux when using simulations only including ion-scale turbulence. Electron-scale turbulence is found to play a dominant role in setting the electron heat flux level and radially elongated (k r ≪ k θ ) “streamers” are found to coexist with ion-scale eddies in experimental plasma conditions. Inclusion of electron-scale turbulence in these simulations is found to increase both ion and electron heat flux levels by enhancing the transport at the ion-scale while also driving electron heat flux at sub-ρ i scales. The combined increases in the low and high-k driven electron heat flux may explain previously observed discrepancies between simulated and experimental electron heat fluxes and indicates a complex interaction of short and long wavelength turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impurity transport, turbulence transitions and intrinsic rotation in Alcator C-Mod plasmas

Linear and nonlinear gyrokinetic simulations are used to probe turbulent impurity transport in intrinsically rotating tokamak plasmas. For this simulation-based study, experimental input parameters are taken from a pair of ICRF heated Alcator C-Mod discharges exhibiting a change in the sign of the normalized toroidal rotation gradient at mid-radius (i.e. a change from hollow to peaked intrinsic rotation profiles). The simulations show that there is no change in the peaking of the calcium impurity between the plasmas with peaked and hollow rotation profiles, suggesting that the impurity transport and the shape of the rotation do not always change together. Furthermore, near mid-radius, r/a = 0.5 (normalized midplane minor radius), the linear and nonlinear gyrokinetic simulations exhibit no evidence of a transition from ion temperature gradient (ITG) to trapped electron mode dominance when the intrinsic rotation profile changes from peaked to hollow. Finally, extensive nonlinear sensitivity analysis is performed, and there is no change in the ITG critical gradient or in the stiffness of ion heat transport with the change in the intrinsic toroidal rotation profile shape, which suggests that the shape of the rotation profile is not dominated by the ITG onset in these cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantify Plasma Response to Non-Axisymmetric (3D) Magnetic Fields in Tokamaks, Final Report for FES (Fusion Energy Sciences) FY2014 Joint Research Target

The goal of the 2014 Joint Research Target (JRT) has been to conduct experiments and analysis to investigate and quantify the response of tokamak plasmas to non-axisymmetric (3D) magnetic fields. Although tokamaks are conceptually axisymmetric devices, small asymmetries often result from inaccuracies in the manufacture and assembly of the magnet coils, or from nearby magnetized objects. In addition, non-axisymmetric fields may be deliberately applied for various purposes. Even at small amplitudes of order 10 -4 of the main axisymmetric field, such “3D” fields can have profound impacts on the plasma performance. The effects are often detrimental (reduction of stabilizing plasma rotation, degradation of energy confinement, localized heat flux to the divertor, or excitation of instabilities) but may in some case be beneficial (maintenance of rotation, or suppression of instabilities). In general, the magnetic response of the plasma alters the 3D field, so that the magnetic field configuration within the plasma is not simply the sum of the external 3D field and the original axisymmetric field. Typically the plasma response consists of a mixture of local screening of the external field by currents induced at resonant surfaces in the plasma, and amplification of the external field by stable kink modes. Thus, validated magnetohydrodynamic (MHD) models of the plasma response to 3D fields are crucial to the interpretation of existing experiments and the prediction of plasma performance in future devices. The non-axisymmetric coil sets available at each facility allow well-controlled studies of the response to external 3D fields. The work performed in support of the 2014 Joint Research Target has included joint modeling and analysis of existing experimental data, and collaboration on new experiments designed to address the goals of the JRT. A major focus of the work was validation of numerical models through quantitative comparison to experimental data, in order to gain confidence in their ability to predict the performance of ITER and other burning plasmas. An important element of this effort was also to verify experimentally the limits of validity of the models. The rest of this overview will summarize the key results of the work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

External excitation of a short-wavelength fluctuation in the Alcator C-Mod edge plasma and its relationship to the quasi-coherent mode

A novel “Shoelace” antenna has been used to inductively excite a short-wavelength edge fluctuation in a tokamak boundary layer for the first time. The principal design parameters, k{sub ⊥}=1.5±0.1 cm{sup −1} and 45<f<300 kHz, match the Quasi-Coherent Mode (QCM, k{sub ⊥}∼1.5 cm{sup −1}, f∼50−150 kHz) in Alcator C-Mod, responsible for exhausting impurities in the steady-state, ELM-free Enhanced D{sub α} H-mode. In H-mode, whether or not there is a QCM, the antenna drives coherent, field-aligned perturbations in density, n{sup ~}{sub e}, and field, B{sup ~}{sub θ}, which are guided by field lines, propagate in the electron diamagnetic drift direction, and exhibit a weakly damped (γ/ω{sub 0}∼5%−10%) resonance near the natural QCM frequency. This result is significant, offering the possibility that externally driven modes may be used to enhance particle transport. In L-mode, the antenna drives only a non-resonant B{sup ~}{sub θ} response. The facts that the driven mode has the same wave number and propagation direction as the QCM, and is resonant at the QCM frequency, suggest the antenna may couple to this mode, which we have shown elsewhere to be predominantly drift-mode-like [B. LaBombard et al., Phys. Plasmas 21, 056108 (2014)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New insights on boundary plasma turbulence and the quasi-coherent mode in Alcator C-Mod using a Mirror Langmuir Probe

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/cm{sup 2}. 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-scale gyrokinetic simulation of Alcator C-Mod tokamak discharges

Alcator C-Mod tokamak discharges have been studied with nonlinear gyrokinetic simulation simultaneously spanning both ion and electron spatiotemporal scales. These multi-scale simulations utilized the gyrokinetic model implemented by GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] and the approximation of reduced electron mass (μ = (m{sub D}/m{sub e}){sup .5} = 20.0) to qualitatively study a pair of Alcator C-Mod discharges: a low-power discharge, previously demonstrated (using realistic mass, ion-scale simulation) to display an under-prediction of the electron heat flux and a high-power discharge displaying agreement with both ion and electron heat flux channels [N. T. Howard et al., Nucl. Fusion 53, 123011 (2013)]. These multi-scale simulations demonstrate the importance of electron-scale turbulence in the core of conventional tokamak discharges and suggest it is a viable candidate for explaining the observed under-prediction of electron heat flux. In this paper, we investigate the coupling of turbulence at the ion (k{sub θ}ρ{sub s}∼O(1.0)) and electron (k{sub θ}ρ{sub e}∼O(1.0)) scales for experimental plasma conditions both exhibiting strong (high-power) and marginally stable (low-power) low-k (k{sub θ}ρ{sub s} < 1.0) turbulence. It is found that reduced mass simulation of the plasma exhibiting marginally stable low-k turbulence fails to provide even qualitative insight into the turbulence present in the realistic plasma conditions. In contrast, multi-scale simulation of the plasma condition exhibiting strong turbulence provides valuable insight into the coupling of the ion and electron scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel transport studies of high-Z impurities in the core of Alcator C-Mod plasmas

Measurements of poloidal variation, ñ{sub 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.2 <0.3 are observed for r/a<0.8, and accumulation on both the high-field side, n{sub z,cos}<0, and low-field side, n{sub z,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.05 <0.10, are observed over 0.5 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, n{sub z}Z{sup 2}/n{sub i}≪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{sub θ,z}. Artificially modifying main-ion flow in parallel transport simulations is shown to impact both ñ{sub z}/ and v{sub θ,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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗