V.A. Izzo, 1 N. Commaux, 2 N.W. Eidietis, 3 R.S. Granetz, 4 E. - PowerPoint PPT Presentation

The Role of MHD in 3D Aspects of Massive Gas Injection V.A. Izzo, 1 N. Commaux, 2 N.W. Eidietis, 3 R.S. Granetz, 4 E. Hollmann, 1 G. Huijsmans, 7 D.A. Humphreys, 3 C.J. Lasnier, 3 M. 25 th IAEA Fusion Lehnen, 7 A. Loarte, 7 R.A. Moyer, 1 Energy Conference P.B. Parks, 3 C. Paz-Soldan, 5 R. Raman, 6 D. Shiraki, 2 E.J. Strait 3 16 October 2014 1 UCSD, 2 ORNL, 3 GA, 4 MIT, 5 ORISE, TH/4-1 6 UW, 7 ITER IO

Massive Gas Injection is a leading candidate for disruption mitigation on ITER In the event that a disruption is unavoidable, the goal of massive gas injection (MGI) shutdown is to radiate plasma stored energy in order to: 1) Avoid conduction of large heat loads to the divertor during the thermal quench (TQ), and … 2) Appropriately tailor the current quench (CQ) time to avoid large vessel forces

Goal of massive gas injection is to isotropically radiate plasma stored energy # of valves & location(s) MGI valve Radiation toroidal peaking factor (TPF)

NIMROD modeling finds a more complicated relationship # of valves & location(s) MGI ? ? valve MHD Impurity transport Heat flux Radiation toroidal peaking factor (TPF)

Outline # of valves & location(s) PART I. Key 3D MGI ? ? Physics of Massive valve Gas Injection MHD Impurity transport PART II. DIII-D TPF Heat flux Predictions & Comparison with Measurements Radiation toroidal peaking PART III. ITER TPF factor (TPF) Predictions

PART I. Key 3D Physics of Massive Gas Injection NIMROD 4-stage MGI shutdown Early TQ CQ Pre-TQ Late TQ Ne MGI at t=0 MHD None m/n >1 m=1/n=1 Particle Ne plume expansion || to B Radial mixing Transport Fast ||  B r V ·  T Heat Slow  conduction conduction convection Transport

PART I. Key 3D Physics of Massive Gas Injection NIMROD 4-stage MGI shutdown NIMROD predictions Early TQ Pre-TQ CQ Late TQ concerning the role of the n=1 mode have been Ne MGI tested experimentally at t=0 NIMROD finds asymmetric impurity spreading for off-midplane injection MHD None m/n >1 m=1/n=1 Particle Ne plume expansion || to B Radial mixing Transport NIMROD multi-valve MGI Fast ||  B r V ·  T Heat Slow  conduction conduction convection simulations reveal implications of Transport both effects for optimum valve positioning

PART I. Key 3D Physics of Massive Gas Injection NIMROD 4-stage MGI shutdown NIMROD predictions Early TQ Pre-TQ CQ Late TQ concerning the role of the n=1 mode have been Ne MGI tested experimentally at t=0 NIMROD finds asymmetric impurity spreading for off-midplane injection MHD None m/n >1 m=1/n=1 Particle Ne plume expansion || to B Radial mixing Transport NIMROD multi-valve MGI Fast ||  B r V ·  T Heat Slow  conduction conduction convection simulations reveal implications of Transport both effects for optimum valve positioning

NIMROD simulations produced two predictions regarding the role of the 1/1 in an MGI TQ* 1) 1/1 phase determines location of toroidal radiation peaking due to asymmetric convected heat flux 180 º 2) Absent other 0 º m/n=1/1 MGI asymmetries, 1/1 phase is anti-aligned with gas jet *IZZO, V.A., Phys. Plasmas 20 (2013) 056107.

DIII-D experiments: Initial n=1 phase corresponds to NIMROD prediction NIMROD MGI predicted location n=1 phase

DIII-D experiments: n=1 phase at TQ influenced by rotation, error fields Rotation and error field effects (not in simulations) also determine final mode phase at TQ

Experiments verify: the phase of the n=1 mode (relative to the gas jet) affects asymmetry TQ W rad asymmetry vs. applied n=1 phase 0.1 Radiated energy asymmetry  DIII-D experiments: 0.0 Changing phase of applied n=1 fields -0.1 changes measured radiation asymmetry -0.2 during TQ -0.3

PART I. Key 3D Physics of Massive Gas Injection NIMROD 4-stage MGI shutdown NIMROD predictions Early TQ Pre-TQ CQ Late TQ concerning the role of the n=1 mode have been Ne MGI tested experimentally at t=0 NIMROD finds asymmetric impurity spreading for off-midplane injection MHD None m/n >1 m=1/n=1 Particle Ne plume expansion || to B Radial mixing Transport NIMROD multi-valve MGI Fast ||  B r V ·  T Heat Slow  conduction conduction convection simulations reveal implications of Transport both effects for optimum valve positioning

Injected Ne plume spreads along B-field in one direction toroidally  toward HFS poloidally Contours/isosurface of MGI15U ionized Ne density 2.25 ms MGI135L 0.25 ms MGI135L MGI15U MGI15U

Below midplane jet spreads in the opposite toroidal direction, also toward HFS Contours/isosurface of ionized Ne density MGI135L MGI15U 2.25 ms MGI15U 0.25 ms MGI135L MGI135L

PART I. Key 3D Physics of Massive Gas Injection NIMROD 4-stage MGI shutdown NIMROD predictions Early TQ Pre-TQ CQ Late TQ concerning the role of the n=1 mode have been Ne MGI tested experimentally at t=0 NIMROD finds asymmetric impurity spreading for off-midplane injection MHD None m/n >1 m=1/n=1 Particle Ne plume expansion || to B Radial mixing Transport NIMROD multi-valve MGI Fast ||  B r V ·  T Heat Slow  conduction conduction convection simulations reveal implications of Transport both effects for optimum valve positioning

NIMROD: I p direction affects direction of impurity spreading MGI15U NORMAL HELICITY MGI135L Ionized Ne density contours/isosurface MGI15U REVERSED HELICITY MGI135L

Relative spacing of gas valves affects interaction with 1/1 mode 135 º 15 º Temperature contours Radiated power and n=1 amplitude Time (ms)

MGI15U and MGI135L will tend to drive the same 1/1 mode phase 135 º 15 º Gas jets are separated by Temperature 120 º poloidally contours and toroidally Radiated power and n=1 amplitude Time (ms)

Simulation with both gas jets drives same mode phase as single jet 135 º 15 º Normal Helicity Temperature contours Radiated power and n=1 amplitude Time (ms)

Heat flux due to 1/1 convection is simultaneously away from both jets 135 º 15 º 1/1 convection also mixes impurities Temperature inward radially contours at both locations Radiated power and n=1 amplitude Time (ms)

In reversed helicity, spacing of two jets no longer coheres with 1/1 symmetry 135 º 15 º Reversed Helicity Temperature contours Radiated power and n=1 amplitude Time (ms)

Interaction of 1/1 mode with each of the two impurity plumes is very different 135 º 15 º No coherent 1/1 mode can interact Temperature with both contours jets in the same way Radiated power and n=1 amplitude Time (ms)

PART II. NIMROD asymmetry predictions and comparison with DIII-D measurements  DIII-D has two fast radiated power measurements  Both jets are closer to Prad90 Diagnostic locations Radiated Energy Prad90 TPF = Max(Wrad)/Mean(Wrad) Clearly, asymmetry calculated from 2 measurement locations is an Toroidal angle approximation…

NIMROD predicts improved symmetry when both DIII-D jets are used Pre-TQ TQ ONLY MGI135L ONLY MGI135L ONLY MGI15U ONLY MGI15U BOTH BOTH All cases in normal helicity

NIMROD predicts improved symmetry when both DIII-D jets are used Pre-TQ TQ ONLY MGI135L ONLY MGI135L ONLY MGI15U ONLY MGI15U BOTH BOTH All cases in normal helicity

DIII-D finds little or no variation in the asymmetry for one vs two gas jets DIII-D measured Radiated energy asymmetry asymmetry Pre-TQ TQ CQ Asymmetry calculated MGI135L MGI15U from 90 and ONLY ONLY BOTH 210 degree detectors t MGI135L – t MGI15U (ms)

NIMROD synthetic asymmetry diagnostic largely reproduces missing trend in DIII-D data NIMROD 2-point DIII-D “TPF” measured Radiated energy asymmetry asymmetry Pre-TQ Comparison of asymmetry TQ CQ using only information from 90 and MGI135L MGI15U 210 degrees ONLY ONLY BOTH t MGI135L – t MGI15U (ms)

NIMROD: 2- point “TPF” does not capture real trend in TPF NIMROD real TPF Pre-TQ TQ ONLY MGI15U ONLY MGI135L ONLY MGI135L ONLY MGI15U BOTH BOTH NIMROD “synthetic 2- point TPF”

NIMROD: reversing helicity increases TQ TPF with 2 jets Pre-TQ ONLY MGI15U TQ ONLY MGI135L ONLY MGI135L ONLY MGI15U BOTH BOTH Reversed Helicity Case

Part III. ITER simulations use three upper ports allocated for TQ mitigation part of DMS Normalized Ne injection rate Fraction of plenum injected Total particle injection rate vs. time based on FLUENT calculations  Assumes 1 m delivery tube: unrealistically short!

3-valves and 1-valve have same TPF, different TQ durations • Slight decrease in TPF during • Single valve has higher pre-TQ with 3 valves maximum Prad • Virtually no change in TPF • Three valve has longer TQ during TQ duration TPF Number of valves Time (ms)

NIMROD modeling provides new physics insights into MGI with single or multiple gas valves NIMROD predicts that DIII-D 2-valve configuration reduces TPF, but increased diagnostic resolution is needed to capture trend, validate model On ITER, 3 upper valve configuration is not found to reduce TPF compared to single upper valve during TQ  Single jet TPF during the thermal quench is not very severe in DIII-D or ITER