GA A27948 FIRST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION

Similar documents
GA A27806 TURBULENCE BEHAVIOR AND TRANSPORT RESPONSE APPROACHING BURNING PLASMA RELEVANT PARAMETERS

GA A27418 THE ROLE OF ZONAL FLOWS AND PREDATOR- PREY OSCILLATIONS IN THE FORMATION OF CORE AND EDGE TRANSPORT BARRIERS

GA A THERMAL ION ORBIT LOSS AND RADIAL ELECTRIC FIELD IN DIII-D by J.S. degrassie, J.A. BOEDO, B.A. GRIERSON, and R.J.

Relating the L-H Power Threshold Scaling to Edge Turbulence Dynamics

The role of zonal flows and predator-prey oscillations in triggering

GA A27805 EXPANDING THE PHYSICS BASIS OF THE BASELINE Q=10 SCENRAIO TOWARD ITER CONDITIONS

GA A26866 REDUCED ELECTRON THERMAL TRANSPORT IN LOW COLLISIONALITY H-MODE PLASMAS IN DIII-D AND THE IMPORTANCE OF SMALL-SCALE TURBULENCE

GA A27933 TURBULENCE BEHAVIOR AND TRANSPORT RESPONSE APPROACHING BURNING PLASMA RELEVANT PARAMETERS

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT

Multi-scale turbulence, electron transport, and Zonal Flows in DIII-D

Observation of Reduced Core Electron Temperature Fluctuations and Intermediate Wavenumber Density Fluctuations in H- and QH-mode Plasmas

GA A26474 SYNERGY IN TWO-FREQUENCY FAST WAVE CYCLOTRON HARMONIC ABSORPTION IN DIII-D

GA A26057 DEMONSTRATION OF ITER OPERATIONAL SCENARIOS ON DIII-D

Particle transport results from collisionality scans and perturbative experiments on DIII-D

GA A25853 FAST ION REDISTRIBUTION AND IMPLICATIONS FOR THE HYBRID REGIME

STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK

Reduced Electron Thermal Transport in Low Collisionality H-mode Plasmas in DIII-D and the Importance of Small-scale Turbulence

GA A27290 CALCULATION OF IMPURITY POLOIDAL ROTATION FROM MEASURED POLOIDAL ASYMMETRIES IN THE TOROIDAL ROTATION OF A TOKAMAK PLASMA

GA A23114 DEPENDENCE OF HEAT AND PARTICLE TRANSPORT ON THE RATIO OF THE ION AND ELECTRON TEMPERATURES

GA A27235 EULERIAN SIMULATIONS OF NEOCLASSICAL FLOWS AND TRANSPORT IN THE TOKAMAK PLASMA EDGE AND OUTER CORE

GA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIII-D

RWM FEEDBACK STABILIZATION IN DIII D: EXPERIMENT-THEORY COMPARISONS AND IMPLICATIONS FOR ITER

Transport Improvement Near Low Order Rational q Surfaces in DIII D

PREDICTIVE MODELING OF PLASMA HALO EVOLUTION IN POST-THERMAL QUENCH DISRUPTING PLASMAS

GA A22443 STUDY OF H MODE THRESHOLD CONDITIONS IN DIII D

GA A23713 RECENT ECCD EXPERIMENTAL STUDIES ON DIII D

GA A22571 REDUCTION OF TOROIDAL ROTATION BY FAST WAVE POWER IN DIII D

GA A26785 GIANT SAWTEETH IN DIII-D AND THE QUASI-INTERCHANGE MODE

Scaling of divertor heat flux profile widths in DIII-D

GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIII-D

Role of Flow Shear in Enhanced Core Confinement

FAR SCRAPE-OFF LAYER AND NEAR WALL PLASMA STUDIES IN DIII D

GA A24016 PHYSICS OF OFF-AXIS ELECTRON CYCLOTRON CURRENT DRIVE

TARGET PLATE CONDITIONS DURING STOCHASTIC BOUNDARY OPERATION ON DIII D

IMPACT OF EDGE CURRENT DENSITY AND PRESSURE GRADIENT ON THE STABILITY OF DIII-D HIGH PERFORMANCE DISCHARGES

GA A26565 LIMITS TO H-MODE PEDESTAL PRESSURE GRADIENT IN DIII-D

GA A23977 DETAILED MEASUREMENTS OF ECCD EFFICIENCY ON DIII D FOR COMPARISON WITH THEORY

UNDERSTANDING TRANSPORT THROUGH DIMENSIONLESS PARAMETER SCALING EXPERIMENTS

Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas

Spatial, temporal and spectral structure of the turbulence-flow interaction at the L-H transition

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER

CORE TURBULENCE AND TRANSPORT REDUCTION IN DIII-D DISCHARGES WITH WEAK OR NEGATIVE MAGNETIC SHEAR

proposed [6]. This scaling is found for single null divertor configurations with the VB

Plasma Response Control Using Advanced Feedback Techniques

GA A27398 COMPARISON OF DEUTERIUM TOROIDAL AND POLOIDAL ROTATION TO NEOCLASSICAL THEORY

Excitation of Alfvén eigenmodes with sub-alfvénic neutral beam ions in JET and DIII-D plasmas

GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELM-SUPPRESSED REGIMES, AND ITER PREDICTIONS

GA A23428 IMPROVEMENT OF CORE BARRIERS WITH ECH AND COUNTER-NBI IN DIII D

Characterization of the Perpendicular Rotation Velocity of the Turbulence by Reflectometry in the Stellarator TJ-II

GA A26732 DEUTERIUM VELOCITY AND TEMPERATURE MEASUREMENTS ON THE DIII-D TOKAMAK

GA A22993 EFFECTS OF PLASMA SHAPE AND PROFILES ON EDGE STABILITY IN DIII D

GA A26741 SCINTILLATOR-BASED DIAGNOSTIC FOR FAST ION LOSS MEASUREMENTS ON DIII-D

GA A26686 FAST ION EFFECTS DURING TEST BLANKET MODULE SIMULATION EXPERIMENTS IN DIII-D

Validation of Theoretical Models of Intrinsic Torque in DIII-D and Projection to ITER by Dimensionless Scaling

Predicting the Rotation Profile in ITER

Variation of Turbulence and Transport with the Te/Ti Ratio in H-Mode Plasmas

GA A22689 SLOW LINER FUSION

Nonlinear Gyrokinetic Simulations of Ion Turbulence in Impurity Seeded and High Density Toroidal Plasmas

Simulation of Double-Null Divertor Plasmas with the UEDGE Code

GA A24166 SUPER-INTENSE QUASI-NEUTRAL PROTON BEAMS INTERACTING WITH PLASMA: A NUMERICAL INVESTIGATION

Generation of Plasma Rotation in a Tokamak by Ion-Cyclotron Absorption of Fast Alfven Waves

GA A25410 DISRUPTION CHARACTERIZATION AND DATABASE ACTIVITIES FOR ITER

EFFECT OF PLASMA FLOWS ON TURBULENT TRANSPORT AND MHD STABILITY*

GA A26123 PARTICLE, HEAT, AND SHEATH POWER TRANSMISSION FACTOR PROFILES DURING ELM SUPPRESSION EXPERIMENTS ON DIII-D

PLASMA MASS DENSITY, SPECIES MIX AND FLUCTUATION DIAGNOSTICS USING FAST ALFVEN WAVE

Scaling of divertor heat flux profile widths in DIII-D

STATIONARY, HIGH BOOTSTRAP FRACTION PLASMAS IN DIII-D WITHOUT INDUCTIVE CURRENT CONTROL

Thermal ion orbit loss and radial electric field in DIII-D

UCLA. Broadband Magnetic and Density Fluctuation Evolution Prior to First ELM in DIII-D Edge Pedestal. Presented by G. Wang a, In collaboration with

GA A23698 ELECTRON CYCLOTRON WAVE EXPERIMENTS ON DIII D

Plasma Fusion Center Massachusetts Institute of Technology Cambridge, MA Burrell, K.H. General Atomics PO Box San Diego, CA

The I-phase and its relation to other phenomena at AUG

GA A26119 MEASUREMENTS AND SIMULATIONS OF SCRAPE-OFF LAYER FLOWS IN THE DIII-D TOKAMAK

GA A26247 EFFECT OF RESONANT AND NONRESONANT MAGNETIC BRAKING ON ERROR FIELD TOLERANCE IN HIGH BETA PLASMAS

Controlling H-Mode Particle Transport with Modulated Electron Heating in DIII-D and Alcator C-Mod via TEM Turbulence

W.M. Solomon 1. Presented at the 54th Annual Meeting of the APS Division of Plasma Physics Providence, RI October 29-November 2, 2012

GA A27849 APPLICATION OF ELECTRON CYCLOTRON HEATING TO THE STUDY OF TRANSPORT IN ITER BASELINE SCENARIO-LIKE DISCHARGES IN DIII-D

Modeling of MHD Equilibria and Current Profile Evolution during the ERS Mode in TFTR

GA A25592 STABILITY AND DYNAMICS OF THE EDGE PEDESTAL IN THE LOW COLLISIONALITY REGIME: PHYSICS MECHANISMS FOR STEADY-STATE ELM-FREE OPERATION

GA A27437 SCALING OF THE DIVERTOR HEAT FLUX WIDTH IN THE DIII-D TOKAMAK

GA A27444 PROBING RESISTIVE WALL MODE STABILITY USING OFF-AXIS NBI

QTYUIOP LOCAL ANALYSIS OF CONFINEMENT AND TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL.

Drift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer

GA A23403 GAS PUFF FUELED H MODE DISCHARGES WITH GOOD ENERGY CONFINEMENT ABOVE THE GREENWALD DENSITY LIMIT ON DIII D

Characteristics of the H-mode H and Extrapolation to ITER

Triggering Mechanisms for Transport Barriers

GA A22863 PLASMA PRESSURE AND FLOWS DURING DIVERTOR DETACHMENT

Investigation of Intrinsic Rotation Dependencies in Alcator C-Mod

Measurements of Core Electron Temperature Fluctuations in DIII-D with Comparisons to Density Fluctuations and Nonlinear GYRO Simulations

Two-Screen Method for Determining Electron Beam Energy and Deflection from Laser Wakefield Acceleration

DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH

Observation of Neo-Classical Ion Pinch in the Electric Tokamak*

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Overview of equilibrium reconstruction on DIII-D using new measurements from an expanded motional Stark effect diagnostic

Electron Transport Stiffness and Heat Pulse Propagation on DIII-D

Edge Zonal Flows and Blob Propagation in Alcator C-Mod P5.073 EPS 2011

PROGRESS IN STEADY-STATE SCENARIO DEVELOPMENT IN THE DIII-D TOKAMAK

Effect of ion B drift direction on density fluctuation poloidal flow and flow shear

Effect of Resonant and Non-resonant Magnetic Braking on Error Field Tolerance in High Beta Plasmas

GA A25229 FINAL REPORT ON THE OFES ELM MILESTONE FOR FY05

Transcription:

GA A27948 FIRST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION By L. SCHMITZ, B.A. GRIERSON, L. ZENG, J.A. BOEDO, T.L. RHODES, Z. YAN, G.R. McKEE, P. GOHIL, D. ELDON, L. BARDOCZI, C. CHRYSTAL, P.H. DIAMOND, W.A. PEEBLES, G.R. TYNAN, R.J. GROEBNER, K.H. BURRELL, C.C. PETTY, E.J. DOYLE and G. WANG SEPTEMBER 2014

DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

GA A27948 FIRST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION By L. SCHMITZ,* B.A. GRIERSON, L. ZENG,* J.A. BOEDO, T.L. RHODES,* Z. YAN, G.R. McKEE, P. GOHIL, D. ELDON, L. BARDOCZI,* C. CHRYSTAL, P.H. DIAMOND, W.A. PEEBLES,* G.R. TYNAN, R.J. GROEBNER, K.H. BURRELL, C.C. PETTY, E.J. DOYLE* and G. WANG* This is a preprint of a paper to be presented at the Twenty-Fifth IAEA Fusion Energy Conf., October 13-18, 2014 in Saint Petersburg, Russia and published in the Proceedings. *University of California Los Angeles, Los Angeles, California. Princeton Plasma Physics Laboratory, Princeton, New Jersey. University of California San Diego, La Jolla, California. University of Wisconsin-Madison, Madison, Wisconsin. Work supported by the U.S. Department of Energy under DE-FG02-08ER54984, DE-AC02-09CH11466, DE-FG02-07ER54917, DE-FG02-89ER53296, and DE-FG02-08ER54999 GENERAL ATOMICS PROJECT 30200 SEPTEMBER 2014

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION First Direct Evidence of Turbulence-Driven Ion Flow Triggering the L- to H-Mode Transition L. Schmitz, et al. EX/11-4 L. Schmitz, 1 B.A. Grierson, 2 L. Zeng, 1 J.A. Boedo, 3 T.L. Rhodes, 1 Z. Yan, 4 G.R. McKee, 4 P. Gohil, 5 D. Eldon, 3 L. Bardoczi, 1 C. Chrystal, 3 P.H. Diamond, 3 W.A. Peebles, 1 G.R. Tynan, 3 R.J. Groebner, 5 K.H. Burrell, 5 C.C. Petty, 5 E.J. Doyle, 1 and G. Wang 1 1 University of California Los Angeles, PO Box 957099, Los Angeles, CA 90095-7099, USA 2 Princeton Plasma Physics Laboratory, PO Box 451, Princeton, NJ 08543-0451, USA 3 University of California San Diego, 9500 Gilman Dr., La Jolla, CA 92093-0417, USA 4 University of Wisconsin-Madison, 1500 Engineering Dr., Madison, WI 53706, USA 5 General Atomics, PO Box 85608, San Diego, CA 92186-5608, USA email: lschmitz@ucla.edu Abstract Simultaneous measurements of main ion flow, E B flow, and turbulence level ñ/n inside the separatrix (LCFS) show for the first time that the initial turbulence collapse preceding the L-H transition is due to turbulence-driven E B flow and ion flow in the ion diamagnetic direction, opposing the pressure-gradientdriven equilibrium E B flow in the L-mode phase. Low to high confinement (L-H) transitions characterized by limit cycle oscillations (LCO) allow probing the trigger dynamics and synergy of turbulence-driven meso-scale flows, and pressure-gradient driven flows with high spatio-temporal resolution. A density scan indicates that the seed flow shear at the L-mode to LCO transition is lowest near the power threshold minimum. Causality of shear-flow generation has been established: early during LCO, the E B shearing rate leads the ion pressure gradient response; as the LCO evolves, the edge pressure gradient and ion diamagnetic flow increase substantially, and the shearing rate lags the ion pressure gradient. Pressure-gradient-driven shear flow then becomes sufficiently large to secure the final LCO-H-mode transition. A two-predator, one-prey model, retaining arbitrary polarity of turbulence-driven flow, captures essential aspects of the transition dynamics, including the magnitude and direction of the driven poloidal main ion flow. 1. Introduction Developing a physics-based model of the L-H transition is critical for confidently extrapolating the auxiliary heating requirements for ITER from the existing empirical L-H transition power threshold scaling. For the first time, it is shown here that the initial (transient) turbulence collapse preceding the L-H transition is caused by turbulence-generated poloidal ion flow and E B flow opposing the equilibrium (L-mode) edge plasma E B flow. Subsequently, the edge pressure gradient P i and the pressure-gradient-driven E B flow shear increase periodically, eventually sustaining turbulence suppression and H-mode confinement. Evidence from several recent experiments [1 7] has pointed towards a synergistic role of turbulence-driven flows [Zonal Flows (ZFs)] and pressure-gradient-driven flows in the trigger and evolution of the L-H transition. The Reynolds stress, and the nonlinear energy transfer between the turbulence spectrum and the turbulence-driven flow, have been evaluated recently for specific plasma parameters [8 10]. Near power threshold, the transition dynamics is substantially expanded/slowed via limit cycle oscillations (LCO) [4 7] between the turbulence amplitude and the fluid E B velocity, allowing profile and flow measurements to be made with unprecedented spatial and temporal resolution. L-H transitions preceded by LCO have been obtained in DIII-D near marginal input power P th (0.75 P/P th 2) across a range of plasma currents (0.6 I p 1.5 MA), edge safety factor (3.8 q 95 8), and plasma density (1x10 19 m -3 n 5.5x10 19 m -3 ). In the experiments reported here we have measured evolution of the turbulence (density fluctuation level), the total E B flow velocity, the ion pressure gradient, and the poloidal main ion mean flow, across the high and low density branch of the L-H transition power threshold scaling. GENERAL ATOMICS REPORT GA-A27948 1

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING L. Schmitz, et al. THE L- TO H-MODE TRANSITION 2. Experimental Evidence for Turbulence-Driven E B and Ion Flow Figure 1(a,b) show the time evolution of the density fluctuation level ñ, measured via Doppler backscattering (DBS) [11,12] at a normalized wave-number k θ ρ s ~0.5 (measured with a half-width wave-number resolution Δk θ /k θ ~0.4), and the total E B velocity determined from the Doppler shift of the backscattered signal, in a lower single null (LSN) diverted deuterium plasma (L-mode density n = 3x10 19 m -3, I p = 0.75 MA, B~1.8 T). The data shown in Fig. 1(a-d) is obtained 3 cm above the tokamak midplane, outboard of the minimum in the radial electric field (in the outer E B shear layer). The radius and poloidal wave-number probed by DBS are extracted via GENRAY ray tracing, using density radial profiles reconstructed with high time resolution (25 µs) from profile reflectometry. Figure 1(c) shows the v i B component of the E B velocity, calculated from the main ion momentum balance equation v E B = ( Z i en i B) n i T i r + ( v i B) B by subtracting the measured pressure gradient term from the total E B velocity (where Z i, n i, T i, and v i are the main ion charge number, density, temperature and fluid velocity). P i is approximated using the electron density from profile reflectometry (as Z eff ~1.6 and n ic <<n i ) and the carbon ion temperature from charge exchange recombination (CER) (assuming T ic ~T i ). Fluctuation suppression is first observed in the outer shear layer at t=t 0 during a positive peak of the v i B term [Fig. 1(c)]. A local E B flow reversal is also observed at this time [Fig. 1(d)], leading to concomitant large positive and negative shearing rates within a narrow (2 cm) radial region, characteristic of turbulencedriven (meso-scale) flows. Figure 1(e) shows that the (negative) shearing rate in the inner shear layer (inboard of the minimum) does not initially increase significantly as fluctuations are first suppressed, but starts oscillating at the LCO frequency and eventually becomes periodically more negative as the LCO evolves. Hence the outer shear layer is initially most significant for fluctuation suppression and initiating the LCO. Evidence for a dipolar meso-scale E B flow structure in the LCO edge electric field layer is presented in Fig. 2, which shows a contour plot of the radial electric field evolution across the L-mode LCO transition (obtained from the measured Doppler shift of the DBS signal), along with the density fluctuation level from DBS. In L-mode, the radial electric field exhibits a shallow well structure with E B flow in the electron diamagnetic direction ñ/n (au) vexb (km/s) (v i xb)/b (km/s) E B (10 5 rad/s) ExB (10 5 rad/s) 0.8 0.4 0-2 -4-6 2 0-4 8 4 0 2 0-2 -4 L-Mode 153624 R=Rs-1.4cm V Di V De R=R s -1.1 cm R=R s -1.8 cm R=Rs-2.6cm 2433 2434 2435 2436 Time (ms) t 0 LCO (a) (b) (c) (d) (e) Fig. 1. Time evolution of (a) density fluctuation level; (b) E B velocity; (c) main ion v B contribution to the E B velocity; (d) shearing rate ω ExB at two radii inside the LCFS; (e) shearing rate in the inner shear layer. Time of initial turbulence suppression is indicated by a grey bar. R s is the LCFS radius. 2 GENERAL ATOMICS REPORT GA-A27948

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION L. Schmitz, et al. [Fig. 2(a)], compensating the ion diamagnetic flow due to the edge ion pressure gradient. As the LCO is triggered, the electric field plot shows periodic positive (ion-diamagnetic direction) E B flow near the LCFS. The negative (electron-diamagnetic) excursions of ve B in the center of the electric field well (around R~2.26 m) become progressively larger and are clearly phase-shifted with respect to the positive peaks at/outside the LCFS. Furthermore, in the inboard shear layer, around R~2.23 m, the positive (ion-diamagnetic direction) peaks are delayed with respect to the positive peaks at the LCFS/outer shear layer, consistent with radial inward propagation of the flow structure [6]. The observed radial flow structure is also qualitatively consistent with the previously reported Reynolds stress measure-ments near the center of the electric field well [8], and the radial stress profile inferred from beam emission spectroscopy Fig. 2. (a) Contour plot of E B velocity across the L-mode LCO transition; (b) contour plot (BES) velocimetry measurements [13]. of density fluctuation level. The LCFS position Comparing Figs 2(a) and 2(b), it is clearly seen is indicated. Start of turbulence suppression in that fluctuations are initially suppressed inside each LCO cycle is indicated via white dashed lines. the LCFS concomitantly with a short (~200 µs ) positive excursion of ve B [as shown in Fig. 1(b) and Fig. 2(a) after 2434.5 ms as a green-blue feature spanning the outboard and inboard shear layers]. The positive flow extends across the center of the electric field well as fluctuations are suppressed, counteracting the L-mode E B flow, and leading to a transient shear reversal. We show now that the observed positive flow transients are consistent with turbulencedriven poloidal ion flow and E B flow, and inconsistent with turbulence suppression via diamagnetic (profile) shear in the early LCO Fig. 3. Cross correlation coefficient of density phase. The cross-correlations between turbufluctuation envelope ñ and poloidal main ion lence envelope, main ion flow, and pressureflow velocity vi, at the L-LCO transition and gradient driven flow, and their detailed spatio1.5 ms into the LCO phase. The cross correlation of the total E B velocity with vi, temporal evolution, have been measured in a 1.5 ms into the LCO phase, is also shown. helium plasma with dominant electron 19-3 cyclotron heating (ECH) ( n ~ 3.5x10 m, Ip ~ 0.8 MA, B=1.8 T). Figure 3 shows that the main ion poloidal mean flow viθ (obtained via main ion CER) lags ñ after the L-mode-LCO transition and is nearly in phase with vexb. Both results are consistent with turbulence-driven poloidal flow/zf dominating ve B early during the LCO phase. Weak anti-correlation of ñ, and vexb with the toroidal ion flow is found in the early LCO, and the measured modulation of the toroidal flow velocity is a fraction (~15%θ θ GENERAL ATOMICS REPORT GA-A27948 3

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING L. Schmitz, et al. THE L- TO H-MODE TRANSITION 30%) of the co-current toroidal velocity resulting from co-neutral beam injection (NBI) (which reaches ~35 km/s near the outboard LCFS as measured here via main ion (CER). Together with the results presented below in Fig. 5, this indicates that the poloidal velocity modulation is dominant early during the LCO and that the toroidal ion velocity modulation in the LCO is a fraction of the mean toroidal velocity resulting from the co-nbi beam torque. 3. Causality of Shear-Flow Generation The causality of shear flow generation has been investigated via correlation analysis of the E B flow-shearing rate and the ion pressure gradient/diamagnetic flow. Figure 4(a) shows that the total E B shearing rate ω E B leads the ion pressure gradient P i early in the LCO (t o + 0.7 ms), establishing clearly that the E B flow shear modulation is not caused by P i changes, but that the v θ component is dominant. Later in the LCO (t o + 5.2 ms), and closer to the final H-mode transition (t H = t 0 +13.5 ms), the periodic reduction in edge turbulence and edge transport enables a gradual increase (and periodic modulation) of the edge pressure gradient and the ion diamagnetic flow shear becomes dominant such that ω E B then lags P i. In Fig. 4(b) the time evolution of the correlation delay between the shearing rate and the ion pressure gradient is shown. The correlation delay becomes negative ~3 ms after the LCO starts, indicating that the profile shear becomes dominant fairly early in this case. -0.3-0.2-0.1 0 0.1 0.2 0.3 t(ms) Fig. 4. (a) Cross-correlation coefficient between the E B shearing rate w ExB and the ion pressure gradient; shown are data for three times after the L-mode - LCO transition; (b) measured correlation delay vs elapsed time after LCO start. Figure 5(a-c) show the time evolution of the density gradient (used here as proxy for the ion pressure gradient), the E B flow and the density fluctuation level. It has been confirmed that the ion pressure gradient oscillations are in phase with the density gradient oscillations [13]. An expanded view [Fig. 5(d,e)] shows that P i increases periodically, generating strong negative E B flow (blue bar). The data suggest that later in the LCO phase, the turbulence growth rate (and n ) increase due to the increasing P i within each cycle, in turn driving radial transport and reducing P i. The turbulence then collapses due to depletion of turbulence energy [10]. Reduced radial transport then allows P i to rise again, transiently maintaining fluctuation suppression via profile shear. During the final phase of the LCO the pressure gradient (diamagnetic flow) dominates the mean flow E B shearing rate, which becomes sufficiently large to sustain fluctuation suppression and secure the LCO-H-mode transition. A heuristic dynamical model has been used to describe the L-LCO-H-mode transition, using the predator-prey Eqns (1) and (2) for the evolution of ñ and the turbulence-driven flow v ZF [14,15] along with an equation describing diamagnetic flow ( v p ) evolution (3), and an equation for the meann poloidal ion flow (4). This model is based on earlier theoretical work describing two-predator flow-turbulence interactions at the L-H transition [15,16]. In contrast to previous work, however, opposite polarity of P i driven and fluctuation-driven E B/poloidal ion flow is retained here: C( ExB, p i ) Correlation Delay t (ms) 0.8 0.6 0.4 0.2 0-0.2-0.4 0.10 0 (a) -0.10 (b) LCO 140426 t0+0.7 ms t0+5.2 ms t0+13.5 ms 0 4 8 12 16 t-t 0 (ms) H 4 GENERAL ATOMICS REPORT GA-A27948

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION L. Schmitz, et al. Fig. 5. (a) Contour plot of density gradient time evolution, (b) density fluctuation level ñ/n, (c) E B velocity in the outer shear layer. The position of the LCFS is indicated; (d) expanded time slice of density gradient; (e) expanded view of ñ/n and v E B ; a time window of increased density gradient is indicated. ñ / t = γ ñ c 1 v ZF ñ c 2 v P c 3 ñ 2 (1) v ZF / t = c 4 ñ 2 v ZF γ ZF v ZF (2) v P / t = Q + c 5 ñ 2 v P + c 6 v P (3) v iθ / t = c 7 ñ 2 c 8 v iθ c 9 v. P (4) Here, Q is the (slowly increasing) radial heat flux driving the system, ñ is the turbulence amplitude, and the constants c 1 -c 3 describe, in order, turbulence depletion via turbulencedriven mesoscale shear flow, turbulence depletion via diamagnetically driven E B shear flow, and the nonlinear turbulence saturation rate. c 4 ñ 2 and γ ZF describe the turbulent flow drive via the Reynolds stress radial gradient, and the Zonal Flow damping rate; c 5 ñ 2 and c 6 describe the turbulent and neoclassical radial diffusion rates; and c 7 ñ 2, c 8, and c 9 describe, in order, poloidal mean ion flow drive via the Reynolds stress, and neoclassical poloidal flow damping and gradient drive [4]. Toroidal flow is neglected here for simplicity. The initial modeling has been carried out using 0D equations, although a more advanced 1D model has also been developed and is needed in order to address the spatial evolution of the edge flow layer and the radial flow reversal described earlier. 1D modeling results will be discussed in a separate paper. Typical 0D modeling results applicable to the radial region at or just inside the LCFS are presented in Fig. 6 The LCO is triggered via an increase in the turbulence level ñ (green) that drives concomitant meso-scale flow v ZF [Fig. 6(b)] and mean poloidal ion flow v iθ [Fig. 6(a)] in the positive (ion-diamagnetic) direction. In practice, the LCO can be triggered via an increase in turbulence amplitude leading to an increase in the energy transfer GENERAL ATOMICS REPORT GA-A27948 5

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING L. Schmitz, et al. THE L- TO H-MODE TRANSITION Fig. 6. 0D predator-prey modeling results: (a) Time evolution of the turbulence level ñ, the E B velocity component associated with the mean poloidal ion flow v i θbφ/b, and the total E B velocity, across the L- mode-lco-h-mode transition, (b) turbulence level ñ, turbulence-driven flow v ZF, and diamagnetic component of the E B velocity; (c) expanded view and phase delay of turbulence-driven flow, and total E B flow with respect to the turbulence amplitude; (d) phase delay of mean poloidal ion flow and diamagnetic flow with respect to ñ. rate from the turbulence spectrum into the low-frequency flow [8 10]. An increase in turbulence level can occur for example due to a step in auxiliary heating power, a radial turbulence avalanche event, an internal sawtooth crash and concomitant radial heat pulse. The trigger criterion for the LCO is that the turbulent energy transfer rate via the Reynolds stress exceeds the ambient turbulence decorrelation rate and the low frequency poloidal flowdamping rate, which will then cause the turbulence amplitude to collapse and the ZF and v θ to increase markedly. The 0D modeling results shown here are in agreement with and illustrate several important experimental observations. (i) The initial turbulence quench occurs during a positive v iθ and v ExB transient; the poloidal ion flow shows a small phase lag with respect to ñ (~25º in the modeling results, similar to the experimental observation (Fig. 3). (ii) The driven mean ion poloidal flow excursions (blue) oppose the diamagnetic component of the E B velocity, v dia. (iii) Positive viθ (blue) excursions lag ñ with a ~25º similar to that observed experimentally (Fig. 3); the phase delay of the positive Zonal E B flow component with respect to ñ, is 90º in the LCO as expected for a predator-prey relationship and observed experimentally in Fig. 1(a,b)]. (iv) The LCO decreases in frequency (and eventually ceases) due to increasing diamagnetic flow (orange), and pressure-gradient driven shear - / r(v dia ) and concomitantly increasing delay of the turbulence recovery; increasing pressure-gradientdriven shear is clearly observed to secure the final transition to H-mode in the modeling results. (v) Pressure gradient oscillations during the LCO are out of phase with ñ [Fig. 6(d)], indicating that the pressure gradient increases after turbulence suppression in each LCO cycle, again in agreement with experimental observations. 6 GENERAL ATOMICS REPORT GA-A27948

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING THE L- TO H-MODE TRANSITION L. Schmitz, et al. 4. Role of Seed Flow Shear The ultimate goal of a detailed physics-based L-H transition model is to interpret the observed power threshold scaling and allow threshold predictions for future burning plasmas. The energy transfer from the ambient L-mode turbulence into directional mesoscale flows is the important quantity that couples the LCO trigger dynamics to the power threshold. This energy transfer has been described by a rate equation for the (low frequency) flow energy W (where v θ v r is the Reynolds stress and v θ is the poloidal ion flow [10]): W t = v θ v r v θ r γ ZF W r v θ v r v θ. (5) If the flow damping and the third term on the right hand side (describing radial spreading of the low frequency flow energy) are neglected, the turbulence-flow energy transfer is proportional to the Reynolds stress and the radial shear of the poloidal ion velocity. A density dependence of the turbulent energy transfer could result from a density dependence of the Reynolds stress, or the (seed) flow-shearing rate. In addition to macroscopic turbulence levels, we are examining the maximum E B flow-shearing rate and the flow shearing rates due to the diamagnetic and the vxb term in the radial ion momentum balance in the outer shear layer [the latter evaluated from ω E B via subtracting the diamagnetic (profile) shear]. Preliminary results for the different contributions to the L-mode shearing rate are shown in Fig. 7, using a ~10 ms time average just before the L-mode-LCO transition. Data for D,H, and He plasmas are included. Both the (total) E B shearing rate and its vxb component are seen to increase towards low and high plasma densities [Fig. 7(a)], with minimum values close to the density where the minimum transition threshold power is observed (<n> min ~ 2-3x10 19 m -3 ). In contrast, the diamagnetic L-mode shearing rate peaks near <n> min. If the Fig. 7 (a) L-mode E B flow shearing rates and ω ExB -ω dia just before the L-mode-LCO transition; ~1 cm inside the LCFS, vs lineaveraged plasma density; shown are data for D, H, and He plasmas; (b) diamagnetic (profile) shearing rates. flux-surface averaged toroidal velocity shear [17] is small, the vxb shearing rate reflects primarily the poloidal velocity shear. The density dependence of the diamagnetic shear is less pronounced [Fig. 7(b)]; no or only a weak density dependence of the latter had also been observed previously in the ASDEX-U and JET tokamaks [18,19]. Hence the data indicates that the poloidal velocity seed flow shear may scale similarly with density as the L-H transition power threshold, and that enhanced poloidal velocity seed flow shear may be required to initiate the transition at low density. GENERAL ATOMICS REPORT GA-A27948 7

FAST DIRECT EVIDENCE OF TURBULENCE-DRIVEN ION FLOW TRIGGERING L. Schmitz, et al. THE L- TO H-MODE TRANSITION 5. Summary We provide direct evidence of the role of turbulence-driven ZF/main ion poloidal flow in triggering the LCO, and initiating the transition to H-mode confinement. The synergy of turbulence-driven and pressure-gradient-driven flows in accessing H-mode confinement has been confirmed, with pressure-gradient driven flows dominating late in the LCO phase and securing the final LCO-H-mode transition. Causality of turbulence-driven flow in shear flow generation at the L-mode-LCO transition has been established. The observation of a turbulence-driven poloidal main ion flow in the ion diamagnetic direction also explains the direction of the ñ, v ExB limit cycle observed in the outboard shear layer in DIII-D [13] and recently in JFT-2M [20] and in HL-2A under certain conditions [21], which is consistent only with positive (ion diamagnetic direction) flow drive. It should be noted that positive flow transients preceding fluctuation suppression have also been observed in regular fast L-H transitions in DIII-D using the same diagnostic tools. These results confirm the role of turbulence-driven flow in initiating the L-H transition as observed in regular (fast) L-H transitions in DIII-D [9] and C-Mod [10]), and preceding the subsequent pressure gradient/diamagnetic shear flow evolution leading to H-mode confinement [10]. A predatorprey model incorporating turbulence-driven flow, pressure gradient evolution and diamagnetic flow shear as well as mean ion poloidal flow evolution, describes qualitatively the essential experimental observations, including the phasing of turbulence amplitude, main ion poloidal flow, the pressure gradient oscillation and the total E B flow. To investigate the role of seed flow shear, the scaling of the E B, diamagnetic, and poloidal flow shearing rates with density have been measured. Preliminary results indicate that the L-mode poloidal seed flow shear just before the L-mode-LCO transition varies with density in a similar fashion than the L-H transition threshold power, and may play a t role in the observed density scaling of the L-H power threshold. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, using the DIII-D National Fusion Facility, a DOE Office of Science user facility, under Awards DE-FG02-08ER54984, DE-AC02-09CH11466, DE-FG02-07ER54917, DE-FG02-89ER53296, DE-FG02-08ER54999, and DE-FC02-04ER54698. References [1] G.R. McKee, et al., Nucl. Fusion 49, 115016 (2009). [2] S.J. Zweben, R.J.Maqueda, R. Hager, et al., Phys. Plasmas 17, 102502 (2010). [3] G.D. Conway, C. Angioni, F. Ryter, et al., Phys. Rev. Lett. 106, 065001 (2011). [4] G.S. Xu, et al., Phys. Rev. Lett. 107, 125001 (2011). [5] T. Estrada, C. Hidalgo, T. Happel, et al., Phys. Rev. Lett. 107, 245004 (2011). [6] L. Schmitz, L. Zeng, T.L. Rhodes, et al., Phys. Rev. Lett. 108, 155002-5 (2012). [7] J. Cheng, J.Q. Dong, K. Itoh, et al., Phys. Rev. Lett. 110, 265002 (2013). [8] G. Tynan, et al., Nucl. Fusion 53, 073053 (2012). [9] Z. Yan, G.R. McKee, R. Fonck, et al., Phys. Rev. Lett. 112, 125002 (2014). [10] I. Cziegler, G.R. Tynan, P.H. Diamond, et al., Plasma Phys. Control. Fusion 56, 075013 (2014). [11] W.A.Peebles, T.L. Rhodes, J.C. Hillesheim, et al., Rev. Sci. Instrum. 81, 10D902 (2010). [12] J.C. Hillesheim, et al., Rev. Sci. Instrum. 81, 10D907 (2010). [13] Z. Yan, G.R. McKee, private communication [14] L. Schmitz et al., Nucl. Fusion 54, 073012 (2014). [15] E.J. Kim and P.H. Diamond, Phys. Rev. Lett. 90, 185006 (2003). [16] K. Miki and P.H. Diamond, Phys. Plasmas 19, 092306 (2012). [17] J. Boedo, E. Belli, E. Hollmann, et al., Phys. Plasmas 18, 032510 (2011). [18] P. Sauter, C. Puetterich, F. Ryter, et al., Nucl. Fusion 52, 012001 (2012). [19] C. Maggi, E. Dalabie, T.M. Biewer, et al. Nucl. Fusion 54, 023007 (2014) [20] T. Kobayashi, et al., Phys. Rev. Lett. 111, 035002 (2013) 8 GENERAL ATOMICS REPORT GA-A27948