J.C. Wesley

Control and dissipation of runaway electron beams created during rapid shutdown experiments in DIII-D

DIII-D experiments on rapid shutdown runaway electron (RE) beams have improved the understanding of the processes involved in RE beam control and dissipation. Improvements in RE beam feedback control have enabled stable confinement of RE beams out to the volt-second limit of the ohmic coil, as well as enabling a ramp down to zero current. Spectroscopic studies of the RE beam have shown that neutrals tend to be excluded from the RE beam centre. Measurements of the RE energy distribution function indicate a broad distribution with mean energy of order several MeV and peak energies of order 30?40?MeV. The distribution function appears more skewed towards low energies than expected from avalanche theory. The RE pitch angle appears fairly directed (????0.2) at high energies and more isotropic at lower energies (??<?100?keV). Collisional dissipation of RE beam current has been studied by massive gas injection of different impurities into RE beams; the equilibrium assimilation of these injected impurities appears to be reasonably well described by radial pressure balance between neutrals and ions. RE current dissipation following massive impurity injection is shown to be more rapid than expected from avalanche theory?this anomalous dissipation may be linked to enhanced radial diffusion caused by the significant quantity of high-Z impurities (typically argon) in the plasma. The final loss of RE beams to the wall has been studied: it was found that conversion of magnetic to kinetic energy is small for RE loss times smaller than the background plasma ohmic decay time of order 1?2?ms.

Runaway electron confinement modelling for rapid shutdown scenarios in DIII-D, Alcator C-Mod and ITER

MHD simulations of rapid shutdown scenarios by massive particle injection in DIII-D, Alcator C-Mod and ITER are performed in order to study runaway electron (RE) transport during mitigated disruptions. The simulations include a RE confinement model using drift-orbit calculations for test particles. A comparison of limited and diverted plasma shapes is studied in DIII-D simulations, and improved confinement in the limited shape is found due to both spatial localization and reduced toroidal spectrum in the nonlinear MHD activity. C-Mod simulations compare shutdown scenarios in which impurity (Ar) fuelling is concentrated in the edge versus the core, and the confinement of REs in the core is maintained until the onset of the m = 1/n = 1 mode, which is delayed in the case of edge deposition, relative to core deposition. But, the overall RE loss fraction is 100% regardless of Ar fuelling profile. A comparison of simulations across the three devices points to a trend of increased RE confinement with increasing device size, wherein all REs are lost in C-Mod, all are confined in ITER, and a partial loss is observed in DIII-D. This trend is related to a reduction in the fluctuating field amplitude near the plasma edge during the thermal-quench-induced MHD activity. The result bodes poorly for RE mitigation strategies in ITER that rely on MHD deconfinement of REs.

Demonstration of rapid shutdown using large shattered deuterium pellet injection in DIII-D

A severe consequence of a disruption on large tokamaks such as ITER could be the generation of multi-megaelectronvolt electron beams that could damage the vacuum vessel and the structures of the machine if they hit the wall unmitigated. The mitigation of runaway electron beams is thus a key requirement for reliable operation of ITER. In order to achieve reliable disruption mitigation, a new fast shutdown technique has been developed: the injection of a large shattered cryogenic pellet in the plasma, which is expected to increase the electron density up to levels where the beam generation processes are mitigated by collisional losses. This technique has been implemented and tested for the first time ever on DIII-D. The first tests show evidence of an almost instantaneous deposition of more than 260 Pa m3 of deuterium deep in the core. Record local densities during the thermal quench were observed for each injection with a very high reliability. Pellet mass and plasma energy content scans show an improvement of the assimilation of the particles for higher plasma energy and larger pellet mass.

Experiments in DIII-D toward achieving rapid shutdown with runaway electron suppression

Experiments have been performed in the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] toward understanding runaway electron formation and amplification during rapid discharge shutdown, as well as toward achieving complete collisional suppression of these runaway electrons via massive delivery of impurities. Runaway acceleration and amplification appear to be well explained using the zero-dimensional (0D) current quench toroidal electric field. 0D or even one-dimensional modeling using a Dreicer seed term, however, appears to be too small to explain the initial runaway seed formation. Up to 15% of the line-average electron density required for complete runaway suppression has been achieved in the middle of the current quench using optimized massive gas injection with multiple small gas valves firing simultaneously. The novel rapid shutdown techniques of massive shattered pellet injection and shell pellet injection have been demonstrated for the first time. Experiments using external magnetic perturbations to deconfine runaways have shown promising preliminary results.

Effect of applied toroidal electric field on the growth/decay of plateau-phase runaway electron currents in DIII-D

Large relativistic runaway electron currents (0.1?0.5?MA) persisting for ~100?ms are created in the DIII-D tokamak during rapid discharge shut down caused by argon pellet injection. Slow upward and downward ramps in runaway currents were found in response to externally applied loop voltages. Comparison between the observed current growth/decay rate and the rate expected from the knock-on avalanche mechanism suggests that classical collisional dissipation of runaways alone cannot account for the measured growth/damping rates. It appears that a fairly constant anomalous dissipation rate of order 10?s?1 exists, possibly stemming from radial transport or direct orbit losses to the vessel walls, although the possibility of an apparent loss due to current profile shrinking cannot be ruled out at present.

Magnetohydrodynamic simulations of massive gas injection into Alcator C-Mod and DIII-D plasmas

Disruption mitigation experiments using massive gas injection (MGI) on Alcator C-Mod [Hutchinson et al., Phys. Plasmas 1, 1511 (1994)] and DIII-D [Luxon and Davis, Fusion Technol. 8, 441 (1985)] have shown that magnetohydrodynamics (MHD) plays an important role. The three-dimensional MHD code NIMROD [Sovinec et al., J. Comput. Phys. 195, 355 (2004)] has been extended to include atomic physics taken from the KPRAD code to perform simulations of MGI. Considerable benchmarking of the code has been done against Alcator C-Mod for neon and helium gas jet experiments. The code successfully captures the qualitative sequence of events observed in MGI experiments up to the end of the thermal quench. Neon jet simulations also show quantitative agreement with the experimental thermal quench onset time. For helium gas jets, we show that a small percent boron density can significantly alter the results even in the presence of a helium jet with three orders of magnitude higher density. The thermal quench onset time is c...

Control of post-disruption runaway electron beams in DIII-Da)

Recent experiments in the DIII-D tokamak have demonstrated real-time control and dissipation of post-disruption runaway electron (RE) beams. In the event that disruption avoidance, control, and mitigation schemes fail to avoid or suppress RE generation, active control of the RE beam may be an important line of defense to prevent the rapid, localized deposition of RE beam energy onto vulnerable vessel sections. During and immediately after the current quench, excessive radial compression of the runaway beams is avoided by a combination of techniques, improving the likelihood of the beams surviving this dynamic period without a fast termination. Once stabilized, the runaway beams are held in a steady state (out to the ohmic flux limit) with the application of active plasma current and position controls. Beam interaction with the vessel wall is minimized by avoiding distinct thresholds for enhanced wall interaction at small and large radii, corresponding to inner wall and outer limiter interaction, respectiv...

Neoclassical islands, -limits, error fields and ELMs in reactor scale tokamaks

An assessment is presented of the impact of recent magnetohydrodynamic research results on performance projections for reactor scale tokamaks as exemplified by the ITER Final Design Report (ITER/FDR) facility. For nominal ELMy H mode operation, the presence and amplitude of neoclassical tearing modes governs the achievable β value. Recent work finds that the scaling of β at which such modes onset agrees well with a polarization drift model, with the consequence that, with reasonable assumptions regarding seed island width, the mode onset β will be lower in reactor scale tokamaks than in contemporary devices. Confinement degradation by such modes, on the other hand, depends on relative saturated island size which is governed principally by β and secondarily by ν* effects on bootstrap current density. Relative saturated island size should be comparable in present and reactor devices. DT ITER demonstration discharges in JET exhibited no confinement degradation at the planned ITER operating value of βN = 2.2. Theory indicates that electron cyclotron current drive can either stabilize these modes or appreciably reduce saturated island size. Turning to operation in candidate steady state, reverse shear, high bootstrap fraction configurations, wall stabilization of external kink modes is effective while the plasma is rotating but (so far) rotation has not been maintained. Recent error field observations in JET imply an error field size scaling that leads to a projection that the ITER/FDR facility will be somewhat more tolerant to error fields than thought previously. ICRF experiments on JET and Alcator C-Mod indicate that plasmas heated by central energetic particles have benign ELMs compared with the usual type 1 ELM of NBI heated discharges.

Growth and decay of runaway electrons above the critical electric field under quiescent conditions

Extremely low density operation free of error field penetration supports the excitation of trace-level quiescent runaway electron (RE) populations during the flat-top of DIII-D Ohmic discharges. Operation in the quiescent regime allows accurate measurement of all key parameters important to RE excitation, including the internal broadband magnetic fluctuation level. RE onset is characterized and found to be consistent with primary (Dreicer) generation rates. Impurity-free collisional suppression of the RE population is investigated by stepping the late-time main-ion density, until RE decay is observed. The transition from growth to decay is found to occur 3–5 times above the theoretical critical electric field for avalanche growth and is thus indicative of anomalous RE loss. This suggests that suppression of tokamak RE avalanches can be achieved at lower density than previously expected, though extrapolation requires predictive understanding of the RE loss mechanism and magnitude.

Measurements of hard x-ray emission from runaway electrons in DIII-D

The spatial distribution of runaway electron (RE) strikes to the wall during argon pellet-initiated rapid shutdown of diverted and limited plasma shapes in DIII-D is studied using a new array of hard x-ray (HXR) scintillators. Two plasma configurations were investigated: an elongated diverted H-mode and a low-elongation limited L-mode. HXR emission from MeV level REs generated during the argon pellet injection is observed during the thermal quench (TQ) in diverted discharges from REs lost into the divertor. In limiter discharges, this prompt TQ loss is reduced, suggesting improved TQ confinement of REs in this configuration. During the plateau phase when the plasma current is carried by REs, toroidally symmetric HXR emission from remaining confined REs is seen. Transient HXR bursts during this RE current plateau suggest the presence of a small level of wall losses due to the presence of an unidentified instability. Eventually, an abrupt final loss of the remaining RE current occurs. This final loss HXR emission shows a strong toroidal peaking and a consistent spatiotemporal evolution that suggests the development of a kink instability.