Connections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod I. Cziegler J.L. Terry, B. LaBombard, J.W. Hughes MIT - Plasma Science and Fusion Center th 19 Plasma Surface Interactions San Diego, CA - May 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 Diagnostics Lá Array ionization profiles: ionization source Fluxes are then estimated from a simple 1D balance [1] What is increasing the flux so dramatically? [1] B. LaBombard et al., Nucl. Fusion 40 (2000) 2041 10 20 m -3 10 20 m -2 s -1 1.0 0.1 L-mode profiles Near SOL Far SOL -20-10 0 10 20 Distance into SOL(mm) 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 n e /n G 0.43 0.37 0.28 0.23 0.17
Observations: Enhanced D-Alpha H-mode QC mode related to cross-field transport (m 2 /sec) 0.04 0.03 0.02 0.01 D eff 200 150 50 250 EDA H-mode spectrogram QC mode 0 QC mode D = Ã/ n eff Ä e 16 /m 2 ) D not meant to imply diffusive transport: eff à is measured in the same way as before[2] [2] J. L. Terry et al., Nucl. Fusion 45 (2005) 1321 0.96 0.98 1.00 1.02 1.04 1.06 1.08 Time (s) The Quasi-Coherent Mode is an EM edge oscillation Characteristic to EDA H-modes Frequency range ~ 50-150 khz Propagates in the EDD in lab frame Provides edge relaxation instead of large ELM s
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
Full radial profiles of turbulence spectra in wavenumber and frequency space - routinely recorded Ohmic L-modes poloidal arrays + time history: S( k í)=s( k, í)/s( í) dispersions are very nearly linear: k (í)> v = í/ k è 300 250 200 150 50 è è ñ = -1.1cm ñ = 0.1cm ñ = 1.2cm EDGE ph IDD : EDD IDD : EDD IDD : EDD v ph è SEPARATRIX REGION S(k,í) è 300 250 200 150 50 300 250 200 150 50-6 -4-2 0 2 4 6 k pol (cm -1 ) SOL shot 1091014033 S(í)=S (í)+s (í) Relative fluctuation power 10 0 10-1 10-2 10-3 10-4 10-5 10-6 k<0 k>0 EDD propagating part 1 10 0 edge: dominated by electron diam. drift (EDD) prop. feature SOL: dominated by blobs in the ion diam. drift (IDD) direction shot 1091014006 t=1.148-1.178s -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 )
Velocities in edge and far SOL are what we expect for drift wave turbulence and blobs Measured pol. velocity (km/s) 6 4 2 0-2 I. Cziegler et al, Phys. Plasmas 065091 (2010) B =6.8T t I =1.0MA p -4-2 -1 0 1 2 Distance into SOL (cm) IDD : EDD Poloidal velocity in the far SOL where blobs are seen matches ExB velocity estimated from probe measurements EDD propagation in the edge is close to the electron diamagnetic flow velocity estimated from TS profile 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 1 10 0 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 1 10 0 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 1 10 0 ã = -1.3±0.30 1 ã = -4.6±0.15 2 ã 1 c í ã 2 crit through k (í) gives k è 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 - ã 2 indicates forward cascade
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 1 10 low intensity broad plateau 1 10 med intensity pl. diminishing 1 10 ã = -4.6±0.15 2 forward transfer (features breaking up) is reproducible in the range of known theories - Itoh interchange: ã = -4 - Wakatani-Hasegawa (collisional drift): ã = -3 (Kolmogorov-Kraichnan 2D) - Hasegawa (1979): ã = -6 break-in-slope spectral position develops a peak at low densities, indicating hindered inverse dynamics k ñ s 0.20 0.15 0.10 10-4 10-5 10-6 10-7 10-8 10-9 high intensity no plateau 1 10 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
Evidence that the EDD propagating edge turbulence is a major factor in edge particle transport and may be related to the 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 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 better 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 0. 20-3 n e (10 m ) 200 150 L-mode ELMfree EDA 200 150 200 150 shot 1091014020 0.010 50 50 50 0.001 10-2 n e 10-3 10-4 10-5 n e n e n e IDD, SOL EDD, edge L-mode ELM free EDA H-mode 1.10 1.15 1.20 1.25 1.30 Time (s) Relative fluctuation power 10 0 10-1 10-2 10-3 10-4 10-5 10-6 -6-4 -2 0 2 4 6 k pol (cm -1 ) EDD propagating part k ñ =0.11 s -6-4 -2 0 2 4 6 k pol (cm -1 ) shot 1091014006 t=1.22-1.25s 1 10 0-6 -4-2 0 2 4 6 k pol (cm -1 ) 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
Conclusions Found evidence that the edge turbulence is responsible for cross field particle transport, ie creating blobs, which hit the walls - located in the last 2 cm inside the separatrix, the same location the H-mode pedestal forms - turbulence power in goes up rapidly as n e/n G and the perpendicular particle flux increases - in ELMfree H-mode EDD edge, rather than IDD SOL turbulence normalized power is much reduced The possible energy input scale is k ñ s~0.1 showing likely connection to QCM (resistive ballooning drive) and possibly a drift-interchange mixed regime (Scott) Turbulence power may be better organized by grad-p source â indicating further connection to a p The spectral distribution responds sensitively to n e/n G, becoming more developed as n e/ng increases (more inverse cascade dynamics) - Why? What theory predicts the right spectra and what consequences does the right model have for transport?