SOLPS-ITER simulations to study the impact of aspect ratio on edge fueling neutrals in tokamaks
Not Available
Engineering topics
Publications and source records attributed to Fitzpatrick, Richard.
Not Available
The tearing mode stability of an inverse aspect-ratio expanded tokamak plasma equilibrium of general shape is investigated using asymptotic matching techniques. Particular emphasis is placed on the conservation of toroidal electromagnetic angular momentum. The TJ code, which is a specific implementation of the results of the investigation, is described.
The critical n=1, n=2, and n=3 error-field amplitudes needed to trigger error-field penetration in ITER and SPARC, steady-state, ohmically heated plasmas are calculated using a standard asymptotic matching approach. The calculation incorporates plasma impurities, trapped particles, the bootstrap current, and neoclassical poloidal flow-damping. The energy confinement time is specified by the neo-Alcator scaling law in the low-density linear ohmic confinement (LOC) regime, and by the ITER-89P L-mode scaling law in the high-density saturated ohmic confinement regime (SOC) regime. The response of the plasma in the inner region is calculated using a linearized version of the four-field model. At the normal operating electron number density, diamagnetic levels of rotation are found to be sufficient to protect ITER and SPARC ohmically heated plasmas from m=2/n=1 error-field penetration. On other hand, SPARC, and especially ITER, ohmically heated plasmas may be vulnerable to n>1 error-field penetration. ITER and SPARC ohmically heated plasmas are also slightly more susceptible to error-field penetration when the electron fluid at the rational surface rotates in the ion diamagnetic direction, rather than the electron diamagnetic direction. At electron number densities that are sufficiently low that the plasma lies in the LOC confinement regime, the error-field penetration threshold increases with increasing density. However, as soon as the electron number density becomes large enough that the plasma enters the SOC regime, the increase in the error-field penetration threshold with increasing density levels off.
A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.
According to a recent paper [Hu et al ., Phys. Plasmas 26, 120702 (2019)], mode penetration at the top of the pedestal is a necessary and sufficient condition for the suppression of edge localized modes (ELMs) by a resonant magnetic perturbation (RMPs) in an H-mode tokamak discharge. This paper employs asymptotic matching theory to model a particular DIII-D discharge in which ELMs were suppressed by an externally generated, static, n ¼ 2, RMP whose amplitude was modulated at a frequency of 1 Hz. It is demonstrated that the response of the plasma to the applied RMP, in the immediate vicinities of the rational (i.e., resonant) surfaces, is governed by nonlinear, rather than by linear, physics. This is the case because the magnetic island widths associated with driven reconnection exceed the linear layer widths, even in cases where driven reconnection is strongly suppressed by plasma rotation. The natural frequency at a given rational surface (i.e., the helical frequency at which the locally resonant component of the RMP would need to propagate in order to maximize driven reconnection) is found to be offset from the local E B frame in the ion diamagnetic direction. The size of the offset is mostly determined by neoclassical poloidal rotation. Finally, the predictions of a fully nonlinear plasma response model are found to be broadly consistent with the DIII-D experimental data.