The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.Peer ReviewedPostprint (published version

Comparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations

Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. This paper presents non-linear simulations of spontaneous type-I ELMs and pellet-triggered ELMs in ASDEX Upgrade performed with the extended MHD code JOREK. A thorough comparison of the non-linear dynamics of these events is provided. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in A Cathey et al (2020 Nuclear Fusion 60 124007). Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to accurately capture the balance between stabilising and destabilising terms for the edge instabilities. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n¿=¿1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and a narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluence is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.Peer ReviewedPostprint (published version

Raffinement de maillage adaptatif pour la simulation numérique des instabilités MHD dans les tokamaks : le code JOREK

The purpose of this paper is to illustrate both validity and advantages of the implementation of the adaptive mesh raffinement strategy in the recent version of the 3D non-linear MHD code JOREK which uses a technique based on the bicubic Bezier surfaces developed in the paper of Czarny-Huijsmans. We describe the physcal model and establish a refinement criteria. Then, we also present the numerical results of adaptive mesh raffinement simulation for the a tearing instability test case and to the test case of injection mechanism of a small pellet of frozen hydrogen into a tokamak

Numerical study of tearing mode seeding in tokamak X-point plasma

A detailed understanding of island seeding is crucial to avoid (N)TMs and their negative consequences like confinement degradation and disruptions. In the present work, we investigate the growth of 2/1 islands in response to magnetic perturbations. Although we use externally applied perturbations produced by resonant magnetic perturbation (RMP) coils for this study, results are directly transferable to island seeding by other MHD instabilities creating a resonant magnetic field component at the rational surface. Experimental results for 2/1 island penetration from ASDEX Upgrade are presented extending previous studies. Simulations are based on an ASDEX Upgrade L-mode discharge with low collisionality and active RMP coils. Our numerical studies are performed with the 3D, two fluid, non-linear MHD code JOREK. All three phases of mode seeding observed in the experiment are also seen in the simulations: first a weak response phase characterized by large perpendicular electron flow velocities followed by a fast growth of the magnetic island size accompanied by a reduction of the perpendicular electron velocity, and finally the saturation to a fully formed island state with perpendicular electron velocity close to zero. Thresholds for mode penetration are observed in the plasma rotation as well as in the RMP coil current. A hysteresis of the island size and electron perpendicular velocity is observed between the ramping up and down of the RMP amplitude consistent with an analytically predicted bifurcation. The transition from dominant kink/bending to tearing parity during the penetration is investigated

MHD simulations of formation, sustainment and loss of Quiescent H-mode in the all-tungsten ASDEX Upgrade

Periodic edge localized modes (ELMs) are the non-linear consequences of pressure-gradient-driven ballooning modes and current-driven peeling modes becoming unstable in the pedestal region of high confinement fusion plasmas. In future tokamaks like ITER, large ELMs are foreseen to severely affect the lifetime of wall components as they transiently deposit large amounts of heat onto a narrow region at the divertor targets. Several strategies exist for avoidance, suppression, or mitigation of these instabilities, such as the naturally ELM-free quiescent H-mode (QH-mode). In the present article, an ASDEX Upgrade equilibrium that features a QH-mode is investigated through non-linear extended MHD simulations covering the dynamics over tens of milliseconds. The equilibrium is close to the ideal peeling limit and non-linearly develops saturated modes at the edge of the plasma. A dominant toroidal mode number of $n=1$ is found, for which the characteristic features of the edge harmonic oscillation are recovered. The saturated modes contribute to heat and particle transport preventing pedestal build-up to the ELM triggering threshold. The non-linear dynamics of the mode, in particular its interaction with the evolution of the edge safety factor is studied, which suggest a possible new saturation mechanism for the QH-mode. The simulations show good qualitative and quantitative agreement to experiments in AUG. In particular, the processes leading to the termination of QH-mode above a density threshold is studied, which results in the transition into an ELM regime. In the vicinity of this threshold, limit cycle oscillations are observed.Comment: Revised version with modifications from review process include

ELMs in High Recycling ITER Plasmas

International audienc

Challenges in MHD modelling in fusion devices

International audienc

MHD stability in X-point geometry:simulation of ELMs

A non-linear MHD code, named JOREK, is under development with the aim of studying the non-linear evolution of the MHD instabilities thought to be responsible for edge localized modes (ELMs): external kink (peeling) and medium-n ballooning modes. The full toroidal X-point geometry is taken into account including the separatrix, open and closed field lines. Analysis of the influence of the separatrix shows a strong stabilization of the ideal and resistive MHD external kink/peeling modes. One instability remains unstable in the presence of the X-point, characterized by a combination of a tearing and a peeling mode. The so-called peeling-tearing mode shows a much weaker dependence on the edge q. Non-linearly the n = 1 peeling-tearing mode saturates at a constant amplitude yielding a mostly kink-like perturbation of the boundary with an island-like structure close to the X-point. The non-linear evolution of a medium-n ballooning mode shows the formation of density filaments. The density filaments are sheared off from the main plasma by an n = 0 flow non-linearly induced by the Maxwell stress. The amplitude of the ballooning mode is limited by this n = 0 flow and multiple (in time) density filaments can develop to bring the plasma below the stability boundary