Validation of Theoretical Models of Intrinsic Torque in DIII-D and Projection to ITER by Dimensionless Scaling by B.A. Grierson1, C. Chrystal2, W.X. Wang1, J.A. Boedo3, J.S. degrassie2, W.M. Solomon2, G.M. Staebler2, D.J. Battaglia1, T. Tala4, A. Salmi4 S-P. Pehkonen4 J. Ferreira5, C. Giroud6 and JET Contributors 1 PPPL, 2GA, 3UCSD, 4VTT, 5IST, 6CCFE Presented at the 26th IAEA Fusion Energy Conference Kyoto, Japan October 17-22, 216 B.A. Grierson / IAEA FEC / Oct 216
Experiments on DIII-D have Advanced the Predictive Capability for Intrinsic Rotation in ITER Toroidal rotation in ITER is expected to be strongly influenced by intrinsic sources Dimensionless Scaling ITER Norm. Intrinsic Torque Empirical dimensionless scalings provide extrapolation from multimachine databases First-principles based predictions require momentum diffusion, pinch and residual fluxes plus boundary condition DIII-D.1 JET.1.1 * 1 5-5 Ω ϕ D Intrinsic Rotation (krad/s) GTS 159396.2225-1..2.4.6.8 2
Prediction of Toroidal Rotation Required to Assess Stability, Transport, ELM Suppression for ITER 1-1 Peaking 2 Impurity Ω (krad/s) 153523.37 Hollowing Notch DIII-D ITER Baseline Electron Heating 2..2.4.6.8 1.2 N Can we predict overall angular momentum? Can we predict the core gradients and structure? What is the role of the boundary condition? Overall rotation for MHD stability, W transport Core rotation gradients for E r shear and confinement improvement Rotation and E r in the pedestal for RMP ELM suppression 1 and access to QH-mode 1 C. Paz-Soldan EX/1-2 2 C. Holland TH/6-1 3
Outline Dimensionless scaling of intrinsic torque Testing physics based models of core intrinsic rotation profile Dependence of Edge Flow on Geometry 4
Joint Experiment 1 with JET and ASDEX-U Tests * Scaling of Intrinsic Torque * for DIII-D, JET, ITER.1,.4,.2 ( 2.5, 2.) Dimensionless parameter scan used to match conditions 2 1 T. Tala (EPS 216) 2 Petty, Phys. Plasmas 15 851 (29) 3 W.M. Solomon, Phys. Plasmas 17 5618 (21). 5
Joint Experiment 1 with JET and ASDEX-U Tests * Scaling of Intrinsic Torque * for DIII-D, JET, ITER.1,.4,.2 ( 2.5, 2.) Dimensionless parameter scan used to match conditions 2 Torque perturbation method used to extract intrinsic torque 3 1 T. Tala (EPS 216) 2 Petty, Phys. Plasmas 15 851 (29) 3 W.M. Solomon, Phys. Plasmas 17 5618 (21). 6
Dimensionless Match Obtained on DIII-D for Factor of 1.4 Variation in * ELMy H-mode conditions Low * 2.2 T (Un-Scaled) High * 1.4 T (scaled) Ip scaled with BT, n~bt 4/3, T~BT 2/3 same βn, q95 3. 2.5 2. 1.5.5. 3 2 1 8 T e (kev) T i (kev) Good confinement H98(y,2) ~..2.4.6.8 6 4 2 n e (1 19 m -3 ) Low * (16141.322) High * Scaled (16177.322) 7
Intrinsic Torque Dominantly in Outer Region and Increases When Reducing * Higher torque at low * projects favorably to ITER Shape consistent with previous studies of intrinsic torque 1 r r intdv/t*1.1 Low * (16141) High * (16177) Region of sufficient measurement certainty Measurement cross-validated with zero rotation technique 2 Norm. Int. Torque..2.4.6.8 intdv (Nm).8.6.4 Int. Torque (Nm) Torque Pertub. (16168) Zero Rotation (161112) 1 W.M. Solomon, Nucl. Fusion 51 731 (211) 2 W.M. Solomon, Nucl. Fusion 49 855 (29)..2. B.A. B.A. Grierson / IAEA / IAEA FEC FEC / October / 216 212..2.4.6 EX/P3-19.8 8
Projecting the DIII-D Intrinsic Torque to ITER Yields Approximately 45 Nm of Intrinsic Torque More than NBI Torque ITER has 33 Nm of available neutral beam torque Norm. Int. Torque Intrinsic torque is near available torque from NBI int dv Projection to ITER scaled using T i to dimensionalize Predicted 1 toroidal rotation from τ NBI +τ int Ω 12 krad/s 1% of Alfven speed, marginal for RWM stability.1 ITER = DIII =.75 1.49±.79 =.85 * T i,diii D D ITER * DIII D.1 1.5 T i,it ER 1 C. Chrystal Phys. Plasmas (submitted) 9
Multi-machine Intrinsic Torque Scaling Confirms DIII-D Results and Normalization Similar experiments executed on JET and ASDEX-U 1 ITER Norm. Intrinsic Torque Scaling with ion temperature best organizes data set 1,2 Projection to ITER with multiple parameters in progress 2.1.1.1 * We Project Relatively Low Toroidal Rotation Finite performance improvement for ITER 2 Details of the rotation profile matter DIII-D JET 1 T. Tala (EPS 216) 2 C. Chrystal (APS 216) 1
Outline Dimensionless scaling of intrinsic torque Testing physics based models of core intrinsic rotation profile Dependence of Edge Flow on Geometry 11
Electron Heating Causes a Rotation Reversal and Global NL GTS Simulation Captures both Shape and Magnitude Rotation reversal occurs because of turbulence residual stress 1 5 Ω ϕ D (krad/s) GTS 159396.2225 Validation of residual stress required to predict ITER rotation Simulations indicate low-k ITG with zonal flow E B and global effects responsible 1-5 -1 1 Ω D..2.4.6.8 5-5 -1 ϕ (krad/s) ECH 159396.128 159396.2225.5 MW MW 1 W.X. Wang TH/P3-12 12
Rotation Profile Reversal Correlates with Onset of Temperature Profile Resilience and Confinement Degradation Density ne 2.5-3. 1 19 m -3 and no NBI promotes Te~Ti T e (kev) 2. 1.5 ECH Power.5.5 Direct electron heating raises both temperatures. 2. 1.5 Between.5 MW power degradation sets in T i (kev).5...2.4.6.8 1.2 N 13
Rotation Profile Reversal Correlates with Onset of Temperature Profile Resilience and Confinement Degradation Density ne 2.5-3. 1 19 m -3 and no NBI promotes Te~Ti T e (kev) 2. 1.5 ECH Power.5.5 Direct electron heating raises both temperatures. 2. 1.5 Between.5 MW power degradation sets in T i (kev).5...2.4.6.8 1.2 N 14
Rotation Profile Reversal Correlates with Onset of Temperature Profile Resilience and Confinement Degradation Density ne 2.5-3. 1 19 m -3 and no NBI promotes Te~Ti T e (kev) 2. 1.5 ECH Power 1.7.5.5 Direct electron heating raises both temperatures. 2. 1.5 Between.5 MW power degradation sets in T i (kev).5...2.4.6.8 1.2 N 15
Rotation Profile Reversal Correlates with Onset of Temperature Profile Resilience and Confinement Degradation Density ne 2.5-3. 1 19 m -3 and no NBI promotes Te~Ti T e (kev) 2. 1.5 ECH Power 2.2 1.7.5.5 Direct electron heating raises both temperatures. 2. 1.5 Between.5 MW power degradation sets in T i (kev).5...2.4.6.8 1.2 N 16
Rotation Profile from Main-ion CER Show Rotation Reversal When Power Degradation Sets In On DIII-D direct measurements of the mainion (D) toroidal flow available 1 Other machines use impurities Deuterium carries energy, momentum fluxes Location where rotation gradient changes indicates non-diffusive flux Not necessarily the sign (+/-) of rotation velocity 1 5-5 Ω ϕ D (krad/s).5-1..2.4.6.8 1 Grierson et. al. Phys. Plasmas (214) 17
Rotation Profile from Main-ion CER Show Rotation Reversal When Power Degradation Sets In On DIII-D direct measurements of the mainion (D) toroidal flow available 1 Other machines use impurities Deuterium carries energy, momentum fluxes Location where rotation gradient changes indicates non-diffusive flux Not necessarily the sign (+/-) of rotation velocity 1 5-5 Ω ϕ D (krad/s).5-1..2.4.6.8 1 Grierson et. al. Phys. Plasmas (214) 18
Rotation Profile from Main-ion CER Show Rotation Reversal When Power Degradation Sets In On DIII-D direct measurements of the mainion (D) toroidal flow available 1 Other machines use impurities Deuterium carries energy, momentum fluxes Location where rotation gradient changes indicates non-diffusive flux Not necessarily the sign (+/-) of rotation velocity 1 5-5 -1 Ω ϕ D (krad/s).5 1.7..2.4.6.8 1 Grierson et. al. Phys. Plasmas (214) 19
Rotation Profile from Main-ion CER Show Rotation Reversal When Power Degradation Sets In On DIII-D direct measurements of the mainion (D) toroidal flow available 1 Other machines use impurities Deuterium carries energy, momentum fluxes Location where rotation gradient changes indicates non-diffusive flux Not necessarily the sign (+/-) of rotation velocity 1 5-5 -1 Ω ϕ D (krad/s) Boundary Effect.5 1.7 2.2..2.4.6.8 1 Grierson et. al. Phys. Plasmas (214) 2
At Onset of Power Degradation Plasma Becomes Linearly Unstable to ITG at Radius where Rotation Gradient Increases Direct electron heating raises both Te, Ti with clear increase of a/lti at mid-radius 1 5-5 Ω ϕ D TGLF 1 indicates excitation of ITG 2 turbulence (krad/s).5-1..2.4.6.8 1.4 1.2.8.6.4.2. 1.4 1.2.8.6.4 T e (kev) T i (kev) 159396.128 159396.2225.2...2.4.6.8 1.2 1 G.M. Staebler et. al. Phys. Plasmas 12 (25) B.A. B.A. Grierson / IAEA / IAEA FEC FEC / October / 216 212 2 C.L. Retting, et. al. Phys. Plasmas EX/P3 8 -(21) 19 21
At Onset of Power Degradation Plasma Becomes Linearly Unstable to ITG at Radius where Rotation Gradient Increases Direct electron heating raises both Te, Ti with clear increase of a/lti at mid-radius 1 5-5 Ω ϕ D TGLF 1 indicates excitation of ITG 2 turbulence (krad/s).5-1..2.4.6.8 1.4 1.2.8.6.4.2. 1.4 1.2.8.6.4 T e (kev) T i (kev) 159396.128 159396.2225 Largest Increase in a/lti -.3 Growth Rate γ(a/cs ).2....2.4.6.8 1.2..2.4.6.8 k θ s 1 G.M. Staebler et. al. Phys. Plasmas 12 (25) B.A. B.A. Grierson / IAEA / IAEA FEC FEC / October / 216 212 2 C.L. Retting, et. al. Phys. Plasmas EX/P3 8 -(21) 19.5. -.5.2.1 Frequency ω(a/c s ) electron direction ion direction 159396.128 159396.2225 =.5 Mean-field E B (CER) γ E B 22
Global Nonlinear Gyrokinetic Simulation Shows Rotation Profile is Balance of Residual Stress and Diffusion Momentum flux ' = decomposed by series of simulations 1 Three simulations produce Π resid., V p, χ φ m i n i hr 2 r i ' d d V p ' 6 Residual Stress 4 Π ϕ resid + Resid. Diffusive Flux Π ϕ - Π ϕ resid Residual stress balanced by momentum diffusion 2 Momentum pinch small Spatial integration of Π φ + B.C. predicts intrinsic rotation profile 2-2 -4.4.5.6.7.8 1 W.X. Wang Phys. Plasmas 13 (26) 2 W.X. Wang APS (216) 23
Prediction of Intrinsic Rotation Profile from Electrostatic Turbulence Matches Both Shape and Magnitude of Experiment Prandtl number used to relate energy and momentum flux Rotation profile and χ φ is not available in ab. initio. prediction Experiment and theory Pr = χ φ /χ i.7 Qualitative shape and quantitative magnitude in agreement with experiment 1 5-5 Ω ϕ D (krad/s) GTS 159396.2225-1..2.4.6.8 Global Nonlinear GTS Simulation can Predict Intrinsic Rotation Profile from Electrostatic Fluctuation-Induced Residual Stress 24
Outline Dimensionless scaling of intrinsic torque Testing physics based models of core intrinsic rotation profile Dependence of Edge Flow on Geometry 25
Edge Velocity Layer Exhibits Novel Dependece on Boundary Shape - Maximized for ITER Configuration Edge velocity w/pinch possible source of intrinsic co-current angular momentum I p =+/-.96 MA B T =-2. T 5.. +I p V (km/s) Inverting plasma shape test orbit-loss mechanism Find V always co-ip and maximized for LSN, favorable B 1 Same as ITER configuration 15791 -I p -5..95 5 1.1 1.15 1.2 ψ n 1 Boedo et. al. Phys. Plasmas accepted (216) 26
Edge Velocity Layer Exhibits Novel Dependece on Boundary Shape - Maximized for ITER Configuration Edge velocity w/pinch possible source of intrinsic co-current angular momentum I p =+/-.96 MA B T =-2. T 5.. +I p V (km/s) Inverting plasma shape test orbit-loss mechanism Find V always co-ip and maximized for LSN, favorable B 1 Same as ITER configuration 157911 -I p -5..95 5 1.1 1.15 1.2 ψ n 1 Boedo et. al. Phys. Plasmas accepted (216) 27
Edge Velocity Layer Exhibits Novel Dependece on Boundary Shape - Maximized for ITER Configuration Edge velocity w/pinch possible source of intrinsic co-current angular momentum I p =+.96 MA B T =-2. T 5.. +Ip LSN (15791) +Ip USN (157911) -Ip USN (15798) V (km/s) Inverting plasma shape test orbit-loss mechanism Find V always co-ip and maximized for LSN, favorable B 1 Same as ITER configuration Core rotation correlates with edge rotation layer 1,2 V CER (km/s) ψ n =.8 3 2 1 15791 157911 -Ip LSN (15796) -5..95 5 1.1 1.15 1.2 ECH OH ψ n -5 2 4 6 8 V Mach (km/s) ψ n = 1 Boedo et. al. Phys. Plasmas accepted (216) 2 S.R.Haksey APS-DPP (216) 28
Rotation Experiments on DIII-D Are Producing Scalings and First-principles-based Validation of Momentum Transport for ITER Scaling of increased intrinsic torque towards ITER * observed on DIII-D and projects twice as much torque from NBI only Successful capture of hollow intrinsic rotation demonstrates ability to predict rotation profile self-organization Edge rotation where intrinsic torque is maximized correlates with core rotation 1 5-5 Ω ϕ D ITER DIII-D.1 JET.1.1 * (krad/s) GTS Norm. Intrinsic Torque 159396.2225-1..2.4.6.8 B.A. Grierson / IAEA FEC / Oct 216
Bonus Slides B.A. Grierson / IAEA FEC / Oct 216
Predicted Intrinsic Rotation Improves ITER Performance via Direct E B and Indirect Multi-Channel Effect (Density Peaking) C. Chrystal APS-DPP (216) Invited Talk 31
How Does an Intrinsic Rotation Gradient Appear? Need a mechanism to break the toroidal symmetry Ω mean E B shear, ZF E B shear, up/down asymmetry, profile and turbulence intensity variation Intrinsic torque generated by residual stress Π resid. Low-k turbulence in the presence of symmetry breaking Πresid. η=- Π 1 B.C. Both diffusion and pinch can balance residual stress here Ω adjusts with local rotation gradient creating Π diff. until total Π = Πdiff. Residual stress and momentum diffusion balance 32
Balance of Residual Stress and Momentum Diffusion Responsible for Steady-State Experimental Rotation Profile Qualitative balance residual stress and momentum diffusion producing a rotation hollowing" is realized in GTS simulation But quantitative, firstprinciples prediction of rotation profile requires additional information Time (v th /L Ti ) 7 6 5 4 3 2.4.5.6.7.8 Πresid. Π resid. /Π norm 15 1 5 5 1 Cannot use experimental rotation profile to derive diffusive flux, needs Pr η=- Π Πdiff. Residual stress and momentum diffusion balance 33