Flow Measurements in the DIII-D Divertor Flow Measurements in the Divertor Region of DIII-D and Plasma Characterization using a Reciprocating Probe J. Boedo, R. Lehmer, R. Moyer, J. Watkins, D. Hill, T. Evans, A. Leonard, C. Lasnier, R. Maingi, M. Schaffer QTYUIOP
Motivation Flows in the divertor region are important because they are involved in: Power flow to the target plates. Impurity transport. Particle transport. Particle and heat removal from the divertor area. Researching the private region and its boundary with the divertor leg and the X-point is important because: It affects divertor design, baffled and non-baffled. The interaction of the core/divertor leg/private region affects all of the above zones but a lot of the underlying physics is not yet clearly understood.
Results Measurements show that a the plasma flows at Mach 1 over most of the outer leg in a radiative divertor, these flows carry most of the heat to the plates by convection. Reverse flow has been observed at the separatrix between the private region and the divertor leg, consistent with measured radial electric field. Convected power deposition at the plate for attached a detached divertor is compared to that measured with an infrared camera. Steep gradients in plasma potential have been observed at the transition between the private region and the inner and outer legs. These electric fields should drive poloidal flows.
Attached vs Detached Plasmas
Geometry of the measurement Outer Leg Probe The divertor has been explored by moving the plasma so that the various regions of interest can be probed.
Flow Measurements in the DIII-D Divertor Mach Number 1.60 1.40 1.20 1.00 0.80 0.60 Mach Number vs Height from Floor Discharge 93541 Attached Discharge 93543 Dettached Outer Leg A radiating divertor is able to generate high Mach # flows over a large section of the divertor. 0.40 0.20 0.00 0 2 4 6 8 10 12 14 Z(cm) These flows can convect heat in situations when conduction is not efficient.
Density Profiles Density vs Height from Divertor Floor Discharge 93541 --Attached Discharge 93543 --Dettached Outer Leg N e x10 12 (cm -3 ) 140 120 100 80 60 40 20 The dettached plasma features high density (1E14 cm-3 vs 2E13 cm-3 ) which is fairly constant along the leg. 0 0 2 4 6 8 10 12 14 Z (cm)
Temperature profiles 30 25 20 Temperature vs Height from the Divertor Floor Discharge 93541 --Attached Discharge 92543 --Dettached Outer Leg The temperature drops due to the enhanced radiation losses. T e (ev) 15 10 5 0 0 2 4 6 8 10 12 14 Z (cm)
e q = κt 5 / 2 5 1 2 + nv [ ( T + T ) + mv + I ] II Parallel energy transport Parallel energy transport is classically given by: e Conduction 1 2 3 dt ds II 2 e i 2 II Convection In detached plasmas conduction is poor due to flat temperature profiles. Convection (Terms 1, 2, 3) is dominant. In detached plasmas, the parallel velocity is high 2-4E6 cm/s and the density is very high (1E14 cm-3), thereby heat removal is performed by term 1, UNLESS Te is very low so the ionization energy I 0 (~13.6eV) becomes significant and term 3 dominates. 0
Convected heat flux measurement Heat Flux to Divertor Floor from IR camera Heat Flux (W/cm 2 ) 160 140 120 100 80 60 40 20 Private Region Discharge 93543 --Dettached Discharge 93541 --Attached Separatrix from EFIT Divertor Leg Convected heat flux is 80% of the total for dettached plasmas. Conected heat flux is 30% of the total for attached plasmas 0 100 120 140 160 180 200 Radius (cm) Probe
UEDGE calculations have shown that flow reversal is possible at the outer strike point and along the separatrix towards the X point. We have observed this behaviour experimentally.
Geometry of the measurement Private Region to Outer Leg In this configuration, the private region is sampled first and the separatrix is crossed at 5 cm
Radial electric field is measured across the separatrix Vfloat (V) 40 20 0-20 -40-60 Floating Potential (V) Private Region Discharge 94008 t= 3440ms -80 0 0 2 4 6 8 10 12 14 Z (cm) Plasma Potential (V) Outer Leg from EFIT Outer Leg 120 100 80 60 40 20 Private Region to Outer Leg Vplasma (V) A radial electric Field of 20 V/cm is observed at the boundary between the private region and the divertor leg. The electric field induces poloidal flows of 2.0E5 cm/s towards the X point.
Reversed flow is measured along the separatrix 2.00 1.50 Discharge 94008 t=3440ms Mach Number Flow Velocity (cm/s) Outer Leg from EFIT 4 10 6 Private Region to Outer Leg Mach 1.00 0.50 0.00-0.50-1.00 2 10 6 0 10 0-2 10 6-4 10 6 Vflow(cm/s) The Mach probe shows that flow reverses at the boundary between the private region and the outer leg. -1.50 0 2 4 6 8 10 12 14 Z(cm)
Pressure increases across the separatrix 350 Discharge 94008 t=3440 ms Density Temperature Outer Leg from EFIT 50 Private Region to Outer Leg 300 n e x10 12 (cm -3 ) 250 200 150 100 50 40 30 20 10 T e (ev) Pressure increases towards the boundary between the private region and the outer leg and then stays constant 0 0 0 2 4 6 8 10 12 14 Z(cm)
Plasma flow drags the impurities (carbon) towards the X-point. Thereby, understanding the flows has implications for impurity transport and core contamination. Reverse flow is needed to explain the measured carbon radiation distribution at the divertor (T. Evans)
Attached divertor plasma
Geometry for inner leg and private region measurements
Mach number at inner leg towards the floor Mach # 0.0-0.2-0.4-0.6-0.8 Mach # and Flow Velocity vs Height from Floor Discharge 93637 Private Region Mach Vflow(cm/s) Separatrix from EFIT Inner Divertor Leg -1.5 10 6-2 10 6-2.5 10 6-3 10 6-3.5 10 6-4 10 6-4.5 10 6-1.0-5 10 6 0.00 2.00 4.00 6.00 8.00 10.00 12.00 14.00 Z (cm) Vflow (cm/s) The plasma flows towards the floor as expected and no flow reversal (but a flow reduction) is seen at the inner separatrix for this discharge
Temperature and Density profiles T e and N e vs Height from the Floor N e x 10 12 (cm -3 ) 500 400 300 200 100 Private Region Density Discharge 93637 Inner Divertor Leg Temperature 50 40 30 20 10 Te (ev) The temperature and density gradients in the boundary between the private region and the inner leg are very large. 0 0 0 2 4 6 8 10 12 14 Z (cm)
Plasma potential shows transition from private region to inner leg V float (V) Floating and Plasma Potential vs Height from the Floor Discharge 93637 V V float plasma 40 160 30 140 20 120 10 100 0 80-10 60-20 40 Vplasma (V) Large gradient in the plasma potential marks the transition from private region to the inner leg. Poloidal flow towards the floor is driven. -30 20 0 2 4 6 8 10 12 14 Z (cm)
A basic treatment and equations for the interpretation of Mach probe data are found in K.S. Chung and Hutchinson Phys Rev A Vol 38, No 9 (1988). The sound speed is given by: c s = (T e + T is )/m i For normalized drift velocities, defined as; V s = ZT e m u = v V s The ratio of upstream and downstream current densities is: R J up J down = e Ku d where K=1.9, 1.7, 1.3 for Ti=0.2Te, 1.0Te, 2.0Te. In our case we will consider Ti=Te and therefore K=1.7.
3D carbon fiber center support rod Mach pin 1D carbon fiber BN insert pyrolytic washer stack Tip Areas: Mach Up: 0.019 cm2 Mach Do: 0.021 cm2 Double: 0.032 cm2 pyrolytic graphite washer stack with BN interlayers inside 3 B t 4 7 1 divertor floor 6 stack compression spring Boron Nitride spring cover Insulating mechanical interface fast drive tube outside height =0.030 0.030 0.027 Double Probe #1 Double Probe #2 Vf pin 3 pin 7 pin 6 pin height end radius /60 rod radius 1 0.0344 0.03725 0.03573 3 0.030 0.0342 0.03432 4 0.3255 0.414 0.03593 6 0.0270 0.0413/0.0405 0.03582 7 0.030 0.0393/0.392 0.03820 8/18/95
There are many discharges without useful upstream Mach pin data, so we can use the envelope of the double probe current signal, which is equivalent to a full pin Isat. The total current collected by a fully exposed pin can be written as a sum of the downstream and upstream current. I envelope double = I up + I down The downstream current can be writen in terms of the Mach tip signal as: I down = I Mach A d A Mach = J Mach A Mach A d A Mach thus the current densities are: = J Mach A dp J double 2 A d = J up A d + J Mach A d Finally we obtain: J up = 2 J double J down which implies R 2 J double J down J down = e Ku d and the Mach number can be writen as: M = (1 / K )ln( 2 J double J down J down )
The area of the top of the probe is (round with 1mm radius) 3.1E-6 m2. One can estimate the radial fluxes using Γ r = D r r n For a D of 1 m2/sec (for L mode) and 0.1-0.3 m2/sec (H-mode) and density gradients of 1E19 and 1E20 m-4 respectively we get a current of: 0.0004-0.024 ma Another alternative is to use the turbulent fluxes directly. The turbulent flux is of the order of 15E20 10^20 /m2/sec for L-mode (worst case), producing a current of: I = q Ä Γ r A r = q Ä Γ r πr 2 or 0.072 ma, to be compared to a minimum of 40mA to 4 A measured currents. We therefore conclude that the effect of the radial currents on the pin caps are negligible
Does the probe perturb the plasma? hp=0.6 mm h=1.77 mm Γ Γ Γ = D n3dl II Γ II Γ II = 0.5nc s d 2 Γ L II = c s d2 6D For DIII-D Lpar=1 cm!!! As long as Lfloor>Lpar, life is good!
Conclusions For detached or semi-detached divertor Plasma flows at Mach 1 over an extended region in a radiative divertor. Heat removal is then predominantly done by convection over a large volume. The detached divertor features a continuosly accelerating plasma towards the plate, in agreement with the classical model. For attached divertor Flow reversal at the separatrix between the private region and the divertor leg has been observed in highly radiating plasmas