Magnetic flux pumping in 3D nonlinear magnetohydrodynamic simulations

A self-regulating magnetic flux pumping mechanism in tokamaks that maintains the core safety factor at $q\approx 1$, thus preventing sawteeth, is analyzed in nonlinear 3D magnetohydrodynamic simulations using the M3D-C$^1$ code. In these simulations, the most important mechanism responsible for the flux pumping is that a saturated $(m=1,n=1)$ quasi-interchange instability generates an effective negative loop voltage in the plasma center via a dynamo effect. It is shown that sawtoothing is prevented in the simulations if $\beta$ is sufficiently high to provide the necessary drive for the $(m=1,n=1)$ instability that generates the dynamo loop voltage. The necessary amount of dynamo loop voltage is determined by the tendency of the current density profile to centrally peak which, in our simulations, is controlled by the peakedness of the applied heat source profile.Comment: submitted to Physics of Plasmas (23 pages, 15 Figures

Coupling JOREK and STARWALL for Non-linear Resistive-wall Simulations

The implementation of a resistive-wall extension to the non-linear MHD-code JOREK via a coupling to the vacuum-field code STARWALL is presented along with first applications and benchmark results. Also, non-linear saturation in the presence of a resistive wall is demonstrated. After completion of the ongoing verification process, this code extension will allow to perform non-linear simulations of MHD instabilities in the presence of three-dimensional resistive walls with holes for limited and X-point plasmas.Comment: Contribution for "Theory Of Fusion Plasmas, Joint Varenna - Lausanne International Workshop, Villa Monastero, Varenna, Italy (27.-31.8.2012)", accepted for publication in Journal of Physics Conference Serie

Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents

The dynamics of large scale plasma instabilities can strongly be influenced by the mutual interaction with currents flowing in conducting vessel structures. Especially eddy currents caused by time-varying magnetic perturbations and halo currents flowing directly from the plasma into the walls are important. The relevance of a resistive wall model is directly evident for Resistive Wall Modes (RWMs) or Vertical Displacement Events (VDEs). However, also the linear and non-linear properties of most other large-scale instabilities may be influenced significantly by the interaction with currents in conducting structures near the plasma. The understanding of halo currents arising during disruptions and VDEs, which are a serious concern for ITER as they may lead to strong asymmetric forces on vessel structures, could also benefit strongly from these non-linear modeling capabilities. Modeling the plasma dynamics and its interaction with wall currents requires solving the magneto-hydrodynamic (MHD) equations in realistic toroidal X-point geometry consistently coupled with a model for the vacuum region and the resistive conducting structures. With this in mind, the non-linear finite element MHD code JOREK has been coupled with the resistive wall code STARWALL, which allows to include the effects of eddy currents in 3D conducting structures in non-linear MHD simulations. This article summarizes the capabilities of the coupled JOREK-STARWALL system and presents benchmark results as well as first applications to non-linear simulations of RWMs, VDEs, disruptions triggered by massive gas injection, and Quiescent H-Mode. As an outlook, the perspectives for extending the model to halo currents are described.Comment: Proceeding paper for Theory of Fusion Plasmas (Joint Varenna-Lausanne International Workshop), Varenna, Italy (September 1-5, 2014); accepted for publication in: to Journal of Physics: Conference Serie