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

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM

Theoretical interpretation of frequency sweeping observations in the Mega-Amp Spherical Tokamak

Frequency sweeping (chirping) of high frequency magnetohydrodynamic modes is widely observed in tokamak plasmas. In this paper observations of chirping in neutral-beam-heated plasmas in the Mega-Amp Spherical Tokamak (MAST) [A. Sykes, R. J. Akers, L. C. Appel et al., Nucl. Fusion 41, 1423 (2001)] are considered, and it is shown that these may be interpreted using the Berk–Breizman augmentation of the Vlasov–Maxwell equations. This model includes an energetic particle source: it leads not only to a single chirp but also to a series of bursting events. This repetitious behavior is characteristic of the chirping seen in experiments such as MAST. The similarity between features in velocity space and features in frequency space reinforces the theory that hole-clump pair formation is responsible for the observed frequency sweeping

The Effect of Energetic Particles on Resistive Wall Mode Stability in MAST

Resistive wall mode (RWM) stability limits have been probed by MHD spectroscopy and numerical modelling. MAST plasmas have operated up to βN = 5.7, well above the predicted ideal kink no-wall limit or measured resonant field amplification limits due to a combination of rotation and kinetic damping. By varying the density, both the rotation and the fast ion distribution function have been changed dramatically. Detailed drift-kinetic modelling shows that whilst the contribution of energetic beam ions to RWM damping does increase at sufficiently high plasma rotation as to allow resonance with the fast ion precession frequency, the thermal ion damping always dominates over the fast ion contribution

Modelling resonant field amplification due to low-n peeling modes in JET

The MHD code MARS-F is used to model low-n, low-frequency, large-amplitude resonant field amplification peaks observed in JET low-pressure plasmas. The resonant response of a marginally stable, n = 1 ideal peeling mode is offered as a candidate to explain the experimental observation. It is found that, unlike the response of a stable resistive wall mode, the peeling mode response is not sensitive to the plasma rotation, nor to the kinetic effects

Modelling of plasma response to RMP fields in MAST and ITER

The resonant magnetic perturbation (RMP) fields, including the plasma response, are computed within a linear, full toroidal, single-fluid resistive magnetohydrodynamic (MHD) model, and under realistic plasma conditions for MAST and ITER. The response field is found to be considerably reduced, compared with the vacuum field produced by the magnetic perturbation coils. This field reduction relies strongly on the screening effect from the toroidal plasma rotation. Computations also quantify three-dimensional (3D) distortions of the plasma surface, caused by RMP fields. A correlation is found between the computed mode structures, the plasma surface displacement and the observed density pump-out effect in MAST experiments. Generally, the density pump-out tends to occur when the surface displacement peaks near the X-points