Plasma Science and Fusion Center Turbulence and transport studies in ALCATOR C Mod using Phase Contrast Imaging (PCI) Diagnos@cs and Comparison with TRANSP and Nonlinear Global GYRO Miklos Porkolab (in collabora/on with Liang Lin, E. Edlund, J. Dorris, C. Rost and the Alcator team) MIT Plasma Science and Fusion Center (PSFC) with special thanks to Ron Waltz and Jeff Candy, General Atomics References: L. Lin, M. Porkolab,E. M. Edlund, et al, Phys. Plasmas 16, 012502 (2009); also, in Plasma Phys. Contr. Fusion 51, 065006 (2009) Presented at EFTSOMP2010, June 29, 2010, Dublin, Ireland
Well known experimental results from Alcator C (Greenwald, Parker, 1984), as well as other tokamaks up to the present @me, indicate that Ohmic confinement follows a linear scaling with density, followed by a saturated confinement regime in L mode; Why? Gyrokine@cs? The results imply that in the linear confinement (neo Alcator) regime the electron thermal diffusivity should decrease as the density increases.
Turbulence measurements in Alcator C Mod plasmas as a func@on of density show good agreement with the GYRO in ohmic plasmas at frequencies above 80 KHz with phase veloci@es dominantly in the ion diamagne@c direc@on
Turbulence measured with Phase Contrast Imaging (PCI) at high frequencies in H mode plasmas (high densi@es) shows good agreement with GYRO nonlinear simula@ons with phase veloci@es in the ion diamagne@c direc@on
( Experimental results using TRANSP analysis)
Measured transport coefficients at low densi@es significantly deviate form the GYRO code predic@on and cannot be resolved even with varia@on of density and temperature profiles
As density increases the dri` speed decreases to c s and confinement improves Do current driven dri- waves play a role? Theorists say no!
Summary of Turbulence and Transport Study
SUMMARY PHASE CONTRAST IMAGING IS A POWERFUL DIAGNOSTIC TOOL TO MEASURE WAVES, INSTABILITIES AND TURBULENCE IN TOKAMAK PLASMAS THE COST IS VERY REASONABLE ( 150 200 K$) THE SHORTCOMING OF THE DIAGNOSTIC IS LACK OF GOOD LOCALIZATION ALONG THE LASER BEAM
Phase Contrast Imaging (PCI) in Alcator C Mod
Edge Turbulence on Alcator C-Mod I. Cziegler J.L. Terry, B. LaBombard, J.W. Hughes MIT - Plasma Science and Fusion Center rd 3 EFTSOMP Dublin, Ireland, 28-29 June 2010
Background and motivations: Edge-SOL transport schematic Open field line region, Scrape-Off-Layer limiter Closed flux surfaces, plasma edge Material surfaces
Background and motivations: Edge-SOL transport schematic limiter In the SOL the turbulence is manifested as filaments (blobs) Blobs propagate rapidly radially, carrying plasma to the main walls Dominate transport 0 LCFS Z (cm) -2-4 1 cm APD spot size -6 1080108019 @ 0.939933 s 86 88 90 92 major radius (cm)
Background and motivations: Edge-SOL transport schematic limiter In the SOL the turbulence is manifested as filaments (blobs) Blobs propagate rapidly radially, carrying plasma to the main walls Dominate transport through SOL Source of SOL plasma Wavelike edge fluct. Why is SOL transport the way it is? 0 LCFS Z (cm) -2-4 1 cm APD spot size -6 1080108019 @ 0.939933 s 86 88 90 92 major radius (cm)
How does cross-field particle transport work? Observations in Ohmic L-modes 1) As central density is raised, density profiles flatten in the SOL 2) Fluxes across the LCFS increase exponentially with density 10 20 m -2 s -1 100 10 Cross-Field Particle Flux 1mm outside LCFS 1 0.5 1.0 1.5 2.0 2.5 Line-Averaged Density (10 20 m -3 ) factor of ~20 10 20 m -3 1.0 0.1 L-mode profiles Near SOL Far SOL -20-10 0 10 20 (mm) Distance into SOL n e /n G 0.43 0.37 0.28 0.23 0.17 [1] Search for the underlying cause of this increase of transport as n grows. e
Observations: Enhanced D-Alpha H-mode QC mode leads to cross-field transport 200 150 100 50 Frequency (khz)250 EDA H-mode spectrogram QC mode (m 2 /sec) 0.04 0.03 0.02 0.01 D eff 0.96 0.98 1.00 1.02 1.04 1.06 1.08 Time (s) In EDA H-mode QCM provides edge relaxation instead of large ELM s The Quasi-Coherent Mode is an EM edge oscillation Frequency range ~ 50-150 khz Propagates in the EDD in lab 0 QC mode D = Ã/ n eff Ä e [2] 16 /m 2 ) D not meant to imply diffusive transport: eff
Experimental setup Gas-Puff-Imaging diagnostics Side View Top View? gas puff injects? neutral D2, sensitive to ne, Te limiter Toroidal view of Outboard Edge? small toroidal extent (~5cm) allows localization? 90 channels cover ~ 5cm x 5cm nozzle? views coupled to APD arrays, sampled @ 2MHz? limiter gas puff
The edge of Ohmic L-mode plasmas is dominated by EDD while the SOL is dominated by IDD propagating features Frequency (khz) poloidal arrays + time history: S(k è í)=s(k è,í)/s(í) Shows dispersion and 300 250 dominant frequencies 200 150 Dispersions are found to be 100 very nearly linear: 50 k è(í)> v ph=2ð í/ k è ñ = -1.1cm ñ = 0.1cm ñ = 1.2cm 300 250 200 150 100 EDGE IDD : EDD IDD : EDD IDD : EDD v ph Frequency (khz) SEPARATRIX REGION Frequency (khz) S(k,í) è 300 250 200 150 100 Frequency (khz) -6-4 -2 0 2 4 6 k pol (cm -1 ) SOL S(í)=S (í)+s (í) k<0 k>0 Shows power spectra Characteristic scales can be seen shot 1091014033 EDGE SOL 50 50-6 -4-2 0 2 4 6 k pol (cm -1 ) -6-4 -2 0 2 4 6 k pol (cm -1 ) -6-4 -2 0 2 4 6 k pol (cm -1 )
EDD velocity in edge is close to ù /k pol * far SOL IDD velocity matches ExB (expected for blobs) Measured pol. velocity (km/s) 6 4 2 0-2 -4-2 -1 0 1 2 Distance into SOL (cm) EDD propagation in the edge is close to the electron diamagnetic flow velocity estimated from TS profile Poloidal velocity in the far SOL where blobs are seen matches ExB velocity estimated from probe measurements B =6.8T t I =1.0MA p Overlap region between the two features [3] IDD : EDD Normalized total power Normalized total power 0.20 0.15 0.10 0.05 0.00-2 -1 0 1 2 Distance into SOL (cm) 0.20 0.15 0.10 0.05 e. prop. feature I. prop. feature n e /n G 0.16 0.40 0.00-2 -1 0 1 2 Distance into SOL (cm)
Wavenumber filtered spectra on underlying dynamics S k>0 (í) can yield information Relative fluctuation power 10 0 10-1 10-2 10-3 10-4 10-5 10-6 EDD propagating part -3-2 -1 0 10 10 10 10 k è ñ s shot 1091014006 t=1.148-1.178s
Wavenumber filtered spectra on underlying dynamics S k>0 (í) can yield information Relative fluctuation power 10 0 10-1 10-2 10-3 10-4 10-5 10-6 EDD propagating part ã 1-3 -2-1 0 10 10 10 10 k è ñ s shot 1091014006 t=1.148-1.178s Spectra are often clear power laws Spectral indices contain information about the spectral transfer dynamics similar to Kolmogorov cascades
Wavenumber filtered spectra on underlying dynamics S k>0 (í) can yield information Relative fluctuation power 10 0 10-1 10-2 10-3 10-4 10-5 10-6 EDD propagating part -3-2 -1 0 10 10 10 10 ã = -1.3±0.40 1 ã = -4.6±0.30 2 ã 1 k è ñ s ã 2 shot 1091014006 t=1.148-1.178s Spectra are often clear power laws Spectral indices contain information about the spectral transfer dynamics similar to Kolmogorov cascades c Break-in-slope í may indicate: - dissipation scale - scale of energy input In the latter case: - ã 1 indicates inverse cascade - indicates forward cascade ã 2
The total spectral power and the distribution in L-mode show a strong dependence on n e /n G n e/n G = 0.15 n e/n G = 0.16 n e /n G = 0.31-3 n e /n G = 0.45 10 10-2 Total relative fluctuation power 10-4 10-5 10-6 10-7 10-8 10-9 clear peak at high freq part -3-2 -1 0 10 10 10 10 k ñ ã = -4.6±0.30 2 low intensity broad plateau med intensity pl. diminishing forward transfer (features breaking up) is reproducible - Wakatani-Hasegawa (collisional drift)[4]: ã = -3 (Kolmogorov-Kraichnan 2D) - Hasegawa-Mima (theory)[5], Fyfe (sim)[6], Xia (exp)[7]: ã = -6 - Chen, Fowler (dimensional)[8]: ã = -5 - Chen (exp): ã = -4.6 break-in-slope spectral position develops a peak hindered inverse dynamics? increased dissipation? k ñ s 0.20 0.15 0.10 10-4 10-5 10-6 ã 2 ã 2 ã 2 ã 2-3 -2-1 0-3 -2-1 0 10-7 10-8 10-9 high intensity no plateau -3-2 -1 0 10 10 10 10 10 10 10 10 10 10 10 10 è s k è ñ s k è ñ s k è ñ s c k ñ s~0.1 0.05 0.15 0.20 0.25 0.30 n e /n G B t = 4.6T 5.4T 6.8T
Bispectral analysis: estimation of three-wave interactions and spectral transfer Bicoherence spectra? strength of phase coupling? no information on energy flow, coupling coefficients Spectral transfer function can be estimated from two separate times:
Bispectral analysis results: energy input at the same scale for H-mode and L-mode stronger transfer towards large scales in L-mode From S(k í) we can move to spatial separation rather than time using è k (í)> v =2ðí/k, and the appropriate Galilean transformation for the velocities. è ph è viewing chans 0.8 0.6 0.4 0.2 0.0-0.2 large errors } Äz } S(k í) è -0.4 1 10 100 1000 Frequency (khz) EDD propagating part large errors time average instead of ensemble average 10 0 10-1 10-2 10-3 10-4 10-5 10-6 k ñ =0.11 s shot 1091014006 t=1.22-1.25s 1 10 100 1000 Frequency (khz)
Total fluctuation power increases with Greenwald fraction Total relative fluctuation power n e/n G = 0.15 10-4 10-5 10-6 10-7 10-8 10-9 n e/n G = 0.16 n e /n G = 0.31-3 n e /n G = 0.45-3 -2-1 0-3 -2-1 0-3 -2-1 0 10 10 10 10 10 10 10 10 10 10 10 10 n e/n G total power 10 10-2 10-4 10-5 10-6 10-7 10-8 10-9 -3-2 -1 0 10 10 10 10
Evidence that the EDD propagating edge turbulence is a major factor in edge particle transport and may be related to the robust tokamak (Greenwald) density limit EDGE REGION SOL densities Total relative fluctuation power (x10-5 ) 20 15 10 5 EDD propagating part 1 0 20 m -3 1.0 0.1 L-mode profiles Near SOL Far SOL -20-10 0 10 20 Distance into SOL(mm) 10 20 m -2 s -1 100 10 Cross-Field Particle Flux 1mm outside LCFS n e /n G 0.43 0.37 0.28 0.23 0.17 factor of ~20 0 0.15 0.20 0.25 0.30 0.35 0.40 0.45 n /n e G 1 0.5 1.0 1.5 2.0 2.5 Line-Averaged Density (10 20 m -3 )
The physical variable â p seems to parametrize the turbulence more sensitively than the Greenwald fraction Total relative fluctuation power (x10-5 ) 20 15 10 5 0 EDGE REGION EDD propagating part I. Cziegler et al, PoP 065091 (2010) 0.15 0.20 0.25 0.30 0.35 0.40 0.45 n /n e G 25 20 15 10 5 âp- normalized pressure better to the free energy source p and perhaps an interchangeballooning type drive? 0 0.05 0.10 0.15 0.20 0.25 â p
Significant reduction in EDD turb. power at L-H transition and suggested connection to the QC mode 3.0 2.5 2.0 1.5 1.0 0.5 0.0 10-2 n e 10-3 20-3 n e (10 m ) n e EDD, edge Frequency (khz) 200 150 100 50 L-mode ELMfree EDA Frequency (khz) 200 150 100 50 Frequency (khz) 200 150 100 50 shot 1091014020 10-4 10-5 0.100 0.010 0.001 n e n e L-mode ELM free IDD, SOL EDA H-mode 1.10 1.15 1.20 1.25 1.30 Time (s) 10 0 10-1 10-2 10-3 10-4 10-5 10-6 0 2 4 6 k (cm -1 ) pol EDD propagating part k ñ =0.11 s shot 1091014006 t=1.22-1.25s 1 10 100 1000 Frequency (khz) 0 2 4 6 0 2 4 6 k (cm -1 ) pol k (cm -1 ) pol Edge turbulence drops at L-H transition No change in norm. SOL fluct. level Further evidence to transport-relevance of EDD edge turbulence critical wavenumber is at the same physical scale as QCM
Fluctuations in the far SOL show no change in character or relative amplitude 3.0 2.5 2.0 1.5 1.0 0.5 0.0 10-2 n e 10-3 20-3 n e (10 m ) n e EDD, edge Frequency (khz) 200 150 100 50 L-mode ELMfree EDA Frequency (khz) 200 150 100 50 Frequency (khz) 200 150 100 50 shot 1091014020 10-4 10-5 0.100 0.010 0.001 n e n e L-mode ELM free IDD, SOL EDA H-mode 1.10 1.15 1.20 1.25 1.30 Time (s) 10 0 10-1 10-2 10-3 10-4 10-5 10-6 0 2 4 6 k (cm -1 ) pol EDD propagating part k ñ =0.11 s shot 1091014006 t=1.22-1.25s 1 10 100 1000 Frequency (khz) 0 2 4 6 0 2 4 6 k (cm -1 ) pol k (cm -1 ) pol 10 0 10-1 10-2 10-3 10-4 10-5 IDD propagating part SOL 1 10 100 Frequency (khz)
Summary Drift-wave-like edge turbulence overlaps with IDD features at the SOL boundary Short wavelength spectra ã = -4.6 quite reproducible but not explained 2 Found evidence that the edge turbulence may be responsible for cross field particle transport, ie creating blobs The energy input scale is k ñ s~0.1 and shows likely connection to QCM (resistive ballooning drive) through spectral transfer Sensitive dependence on n /n (or e G â p - indication of grad-p source?) What theory predicts the right spectra and what consequences does the right model have for transport?
References [1] B. LaBombard et al., Nucl. Fusion 40 (2000) 2041 [2] J. L. Terry et al., Nucl. Fusion 45 (2005) 1321 [3] I. Cziegler et al, Phys. Plasmas 065091 (2010) [4] A. Hasegawa, M. Wakatani, Phys. Rev. Lett. 50, 9 (1983) 682 [5] A. Hasegawa et al., Phys. Fluids 22 (1979) 11 [6] D. Fyfe, D. Montgomery, Phys. Fluids 22 (1978) 2 [7] H. Xia, M.G. Shats, Phys. Rev. Lett. 91(15) 155001 (2003) [8] F. Chen, Phys. Rev. Lett. 15 381 (1965)