Electron Transport Stiffness and Heat Pulse Propagation on DIII-D by C.C. Petty, J.C. DeBoo, C.H. Holland, 1 S.P. Smith, A.E. White, 2 K.H. Burrell, J.C. Hillesheim 3 and T.C. Luce 1 University of California San Diego 2 Massachusetts Institute of Technology 3 University of California Los Angeles Presented at the 54 th American Physical Society Division of Plasma Physics Meeting Providence, Rhode Island October 29 November 2, 212 1
Fundamental Behavior of Drift Wave Turbulent Transport is Tested Using Heat Pulse Propagation Carefully constructed experiment directly probes diffusive transport to test key predicted behaviors: Instability threshold in T e Electron transport stiffness Off-axis ECH varies electron heat flux to scan T e over large range Moved one gyrotron from outside to inside on shot-to-shot basis Modulated one gyrotron (outside) for measurement of heat pulse propagation 2
Modulation in Electron Temperature Profile is Fitted to Linearized Energy Conservation Equation Fourier-transformed first-order equation: "D HP # 2 T e + V HP # T & 1 e + ( $ + i 3 ' HP 2 % ) + T * e = S e Multiple harmonics of T e oscillations are simultaneously fit to determine D HP, V HP, τ HP 3
Heat Pulse Propagation is a Good Test of Transport Stiffness Models This talk focuses on the relation between the heat pulse and power balance diffusivities: D HP = " 1 #Q PB = # n e #$T e #$T e where ( D PB $T e ) Q PB = "n e D PB #T e + n e V PB T e Comparing the incremental (D HP ) to equilibrium diffusivity (D PB ) relates stiffness to the shape of the heat flux curve: Q/n D PB D HP ( ) ( ) S = DHP D = " ln Q/n PB " ln #T T 4
The Heat Pulse Diffusivity at ρ =.6 Rapidly Increases for - T e > 3.2 kev/m Critical Gradient Threshold? D HP (m 2 /s) 25 2 15 1 5 1 2 3 4 5 Key analysis step is to determine the power balance diffusivity by numerical integration of the measured heat pulse diffusivity: "T e D PB = 1 # D HP d "T e "T e ( ) This yields the purely diffusive portion of the equilibrium heat flux 5
Diffusive Heat Flux Falls Short of Total Heat Flux from Power Balance Indicating Something is Missing 6 Diffusive heat flux is Q e /n e (1-19 Wm) 5 4 3 2 1 Q diff = "n e D PB #T e #T e $ ( ) = n e D HP d "#T e The difference between the heat fluxes can be reconciled by a nonzero power balance convective velocity 1 2 3 4 5 V PB = QPB n e T e " 1 T e #T e $ D HP d ("#T e ) 6
The Power Balance Diffusivity Increases Rapidly Above - T crit While Convection is Mainly Outwards "T e D PB = 1 # D HP d "T e "T e ( ) V PB = QPB n e T e " 1 T e #T e $ D HP d ("#T e ) D PB (m 2 /s) V PB (m/s) T e (kev) 7
Measured Stiffness Factor Jumps Up ~4 Times When - T e /T e Exceeds Critical Value S = DHP D = " ln(qpb /n e #V PB T e ) PB " ln($t e ) D HP /D PB 8
Measured Stiffness Factor Jumps Up ~4 Times When - T e /T e Exceeds Critical Value S = DHP D = " ln(qpb /n e #V PB T e ) PB " ln($t e ) D HP /D PB 9
Measured Stiffness Factor Jumps Up ~4 Times When - T e /T e Exceeds Critical Value S = DHP D = " ln(qpb /n e #V PB T e ) PB " ln($t e ) D HP /D PB D HP /D PB T e /T e (m -1 ) 1
Critical Gradient and Stiffness Factor from Nonlinear GYRO Simulations Agree With Experiment Comparison of ECH + co- NBI case GYRO stiffness factor determined using total heat flux S = " lnq e " ln#t e Stiffness 5 4 3 2 1 2 3 4 5 11
The Minimal Critical Gradient Transport Model by Garbet Can Be Tested Against DIII-D Data Simple transport model that preserves some basic properties of turbulent transport Main hypothesis is gyrobohm-like turbulent transport that is switched on above a threshold crit = R T crit /T D & = " s q 1.5 (# R$T D gb ' T #% crit ) & + H (# R$T * ' T #% crit ) + + " * q 1.5! Note that there is no convective term χ s and χ are dimensionless coefficients to be fitted to data χ s is the effective stiffness factor X. Garbet et al. 24 Plasma Phys. Control. Fusion 46 1351 12
Diffusion Coefficients Exhibit Expected Behavior for Transport Switched On Above a Critical Gradient D HP (m 2 /s) model D PB (m 2 /s) D HP /D PB ECH only: L -1 crit = 3.7 m -1 χ s =.8 χ = 1.4 13
TGLF Transport Model has Similar Effective Stiffness Factor as Experiment But Predicts a Lower Critical Gradient Comparison of ECH-only case TGLF diffusion coefficient determined using total heat flux Q e D TGLF = " n e #T e Experiment: L -1 crit = 3.7 m -1 χ s =.8 TGLF: L crit -1 = 3.1 3.4 m -1 χ s =.73 D PB (m 2 /s) 8 6 4 2 1 2 3 4 5 14
Conclusions We have developed a new method to look directly at diffusive behavior by combining heat pulse propagation and power balance analysis In L-mode plasmas with off-axis ECH, a critical value of T e is observed, above which there is a sudden increase (~4 ) in the electron transport stiffness Predicted electron transport stiffness and critical T e from GYRO and TGLF are in good agreement ( 1%) with experiment Will extend this study to H-mode plasmas 15