Dynamics of the Disruption Halo Current Toroidal Asymmetry in NSTX

This paper describes the dynamics of disruption halo current non-axisymmetries in the lower divertor of the National Spherical Torus Experiment [M. Ono, et al. Nuclear Fusion 40, 557 (2000)]. While. The halo currents typically have a strongly asymmetric structure where they enter the divertor floor, and this asymmetry has been observed to complete up to 7 toroidal revolutions over the duration of the halo current pulse. However, the rotation speed and toroidal extend of the asymmetry can vary significantly during the pulse. The rotation speed, halo current pulse duration, and total number of revolutions tend to be smaller in cases with large halo currents. The halo current pattern is observed to become toroidally symmetric at the end of the halo current pulse. It is proposed that this symmeterization is due to the loss of most or all of the closed field line geometry in the final phase of the vertical displacement event

A multi-machine scaling of halo current rotation

Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machinespecific minimum current quench time, τCQ; (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, fh . These data reveal that the normalized rotation duration, trot/τCQ, and the average rotation velocity, vh , are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above fh = 20 Hz is to be expected. More importantly, depending on the projected value of τCQ in ITER, substantial rotation could also occur in the resonant frequency range of 6–20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out.EURATOM 633053RCUK Energy Programme EP/ I501045Princeton University DE-AC02-09CH1146

Exploration of the equilibrium operating space for NSTX-Upgrade

This paper explores a range of high-performance equilibrium scenarios available in the NSTX-Upgrade device [J.E. Menard, submitted for publication to Nuclear Fusion]. NSTX-Upgrade is a substantial upgrade to the existing NSTX device [M. Ono, et al., Nuclear Fusion 40, 557 (2000)], with significantly higher toroidal field and solenoid capabilities, and three additional neutral beam sources with significantly larger current drive efficiency. Equilibria are computed with freeboundary TRANSP, allowing a self consistent calculation of the non-inductive current drive sources, the plasma equilibrium, and poloidal field coil current, using the realistic device geometry. The thermal profiles are taken from a variety of existing NSTX discharges, and different assumptions for the thermal confinement scalings are utilized. The no-wall and idealwall n=1 stability limits are computed with the DCON code. The central and minimum safety factors are quite sensitive to many parameters: they generally increases with large outer plasmawall gaps and higher density, but can have either trend with the confinement enhancement factor. In scenarios with strong central beam current drive, the inclusion of non-classical fast ion diffusion raises qmin, decreases the pressure peaking, and generally improves the global stability, at the expense of a reduction in the non-inductive current drive fraction; cases with less beam current drive are largely insensitive to additional fast ion diffusion. The non-inductive current level is quite sensitive to the underlying confinement and profile assumptions. For instance, for BT=1.0 T and Pinj=12.6 MW, the non-inductive current level varies from 875 kA with ITER-98y,2 thermal confinement scaling and narrow thermal profiles to 1325 kA for an ST specific scaling expression and broad profiles. This sensitivity should facilitate the determination of the correct scaling of transport with current and field to use for future fully non-inductive ST devices. Scenarios are presented which can be sustained for 8-10 seconds, or (20-30)τCR, at βN=3.8-4.5, facilitating, for instance, the study of disruption avoidance for very long pulse. Scenarios have been documented which can operate with βT~25% and equilibrated qmin>1. The value of qmin can be controlled at either fixed non-inductive fraction of 100% or fixed plasma current, by varying which beam sources are used, opening the possibility for feedback qmin control. In terms of quantities like collisionality, neutron emission, non-inductive fraction, or stored energy, these scenarios represent a significant performance extension compared to NSTX and other present spherical torii

Notes on a paper of Mess

These notes are a companion to the article "Lorentz spacetimes of constant curvature" by Geoffrey Mess, which was first written in 1990 but never published. Mess' paper will appear together with these notes in a forthcoming issue of Geometriae Dedicata.Comment: 26 page

Experimental investigation and validation of neutral beam current drive for ITER through ITPA Joint Experiments

Joint experiments investigating the off-axis neutral beam current drive (NBCD) capability to be utilized for advanced operation scenario development in ITER were conducted in four tokamaks (ASDEX Upgrade (AUG), DIII-D, JT-60U and MAST) through the international tokamak physics activity (ITPA). The following results were obtained in the joint experiments, where the toroidal field, B t, covered 0.4-3.7 T, the plasma current, Ip, 0.5-1.2 MA, and the beam energy, Eb, 65-350 keV. A current profile broadened by off-axis NBCD was observed in MAST. In DIII-D and JT-60U, the NB driven current profile has been evaluated using motional Stark effect diagnostics and good agreement between the measured and calculated NB driven current profile was observed. In AUG (at low δ ∼ 0.2) and DIII-D, introduction of a fast-ion diffusion coefficient of Db ∼ 0.3-0.5 m2 s-1 in the calculation gave better agreement at high heating power (5 MW and 7.2 MW, respectively), suggesting anomalous transport of fast ions by turbulence. It was found through these ITPA joint experiments that NBCD related physics quantities reasonably agree with calculations (with Db = 0-0.5 m2 s-1) in all devices when there is no magnetohydrodynamic (MHD) activity except ELMs. Proximity of measured off-axis beam driven current to the corresponding calculation with Db = 0 has been discussed for ITER in terms of a theoretically predicted scaling of fast-ion diffusion that depends on Eb/Te for electrostatic turbulence or βt for electromagnetic turbulence. © 2011 IAEA, Vienna

Dynamics of the Disruption Halo Current Toroidal Asymmetry in NSTX

This paper describes the dynamics of disruption halo current non-axisymmetries in the lower divertor of the National Spherical Torus Experiment [M. Ono, et al. Nuclear Fusion 40, 557 (2000)]. While. The halo currents typically have a strongly asymmetric structure where they enter the divertor floor, and this asymmetry has been observed to complete up to 7 toroidal revolutions over the duration of the halo current pulse. However, the rotation speed and toroidal extend of the asymmetry can vary significantly during the pulse. The rotation speed, halo current pulse duration, and total number of revolutions tend to be smaller in cases with large halo currents. The halo current pattern is observed to become toroidally symmetric at the end of the halo current pulse. It is proposed that this symmeterization is due to the loss of most or all of the closed field line geometry in the final phase of the vertical displacement event

Characterization of the plasma current quench during disruptions in the National Spherical Torus Experiment

A detailed analysis of the plasma current quench in the National Spherical Torus Experiment [M.Ono, et al Nuclear Fusion 40, 557 (2000)] is presented. The fastest current quenches are fit better by a linear waveform than an exponential one. Area-normalized current quench times down to .4 msec/m2 have been observed, compared to the minimum of 1.7 msec/m2 recommendation based on conventional aspect ratio tokamaks; as noted in previous ITPA studies, the difference can be explained by the reduced self-inductance at low aspect ratio and high-elongation. The maximum instantaneous dIp/dt is often many times larger than the mean quench rate, and the plasma current before the disruption is often substantially less than the flat-top value. The poloidal field time-derivative during the disruption, which is directly responsible for driving eddy currents, has been recorded at various locations around the vessel. The Ip quench rate, plasma motion, and magnetic geometry all play important roles in determining the rate of poloidal field change

A Reduced Resistive Wall Mode Kinetic Stability Model for Disruption Forecasting

Kinetic modification of ideal stability theory from stabilizing resonances of mode-particle interaction has had success in explaining resistive wall mode (RWM) stability limits in tokamaks. With the goal of real-time stability forecasting, a reduced kinetic stability model has been implemented in the new Disruption Event Characterization and Forecasting (DECAF) code, which has been written to analyze disruptions in tokamaks. The reduced model incorporates parameterized models for ideal limits on beta, a ratio of plasma pressure to magnetic pressure, which are shown to be in good agreement with DCON code calculations. Increased beta between these ideal limits causes a shift in the unstable region of delta W_K space, where delta W_K is the change in potential energy due to kinetic effects that is solved for by the reduced model, such that it is possible for plasmas to be unstable at intermediate beta but stable at higher beta. Gaussian functions for delta W_K are defined as functions of E cross B frequency and collisionality, with parameters reflecting the experience of the National Spherical Torus Experiment (NSTX). The reduced model was tested on a database of discharges from NSTX and experimentally stable and unstable discharges were separated noticeably on a stability map in E cross B frequency, collisionality space. The reduced model only failed to predict an unstable RWM in 15.6% of cases with an experimentally unstable RWM and performed well on predicting stability for experimentally stable discharges as well.Readme, digital data file