-
The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
Authors:
M Hoelzl,
GTA Huijsmans,
SJP Pamela,
M Becoulet,
E Nardon,
FJ Artola,
B Nkonga,
CV Atanasiu,
V Bandaru,
A Bhole,
D Bonfiglio,
A Cathey,
O Czarny,
A Dvornova,
T Feher,
A Fil,
E Franck,
S Futatani,
M Gruca,
H Guillard,
JW Haverkort,
I Holod,
D Hu,
SK Kim,
SQ Korving
, et al. (28 additional authors not shown)
Abstract:
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and phys…
▽ More
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. 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 and current quench triggered by massive gas injection (MGI) and shattered pellet injection (SPI), runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events (VDEs). Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring thermal quenches triggered by locked modes, and radiative collapses are being studied.
△ Less
Submitted 21 April, 2021; v1 submitted 18 November, 2020;
originally announced November 2020.
-
Simulation of the electromagnetic wall response during Vertical Displacement Events (VDE) in ITER tokamak
Authors:
Cǎlin V. Atanasiu,
Leonid E. Zakharov,
Karl Lackner,
Matthias Hoelzl
Abstract:
The key basis for tokamak plasma disruption modeling is to understand how currents flow to the plasma facing surfaces during plasma disruption events. In ITER tokamak, the occurrence of a limited number of major disruptions will definitively damage the chamber with no possibility to restore the device. In the current exchange plasma-wall-plasma, according to the Helmholtz decomposition theorem, ou…
▽ More
The key basis for tokamak plasma disruption modeling is to understand how currents flow to the plasma facing surfaces during plasma disruption events. In ITER tokamak, the occurrence of a limited number of major disruptions will definitively damage the chamber with no possibility to restore the device. In the current exchange plasma-wall-plasma, according to the Helmholtz decomposition theorem, our surface current density in the conducting shell - the unknown of our problem - being a vector field twice continuously differentiable in 3D, has been splited into two components: an irrotational (curl-free) vector field and a solenoidal (divergence-free) vector field. Developing a weak formulation form and minimizing the correspondent energy functionals in a Finite Element approach, we have obtained the space and time distribution of the surface currents. We verified successfully our numerical simulation with an analytical solution with pure homogeneous Neumann B.C. and satisfying the necessary existence condition. By considering the iron core presence in JET tokamak, we have split the magnetization currents - the unknowns in some integral equations - into two components, the first producing a magnetic field in the iron region only and the second producing a magnetic field in the vacuum, obtaining thus a better evaluation of the influence of the iron core on the plasma equilibrium. To reduce the influence of the singularities appearing during the surface currents determination in multiply connected domains (L-shaped domains) we have used a conformal transformation method.
△ Less
Submitted 24 October, 2018;
originally announced October 2018.
-
Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents
Authors:
M. Hoelzl,
G. T. A. Huijsmans,
P. Merkel,
C. Atanasiu,
K. Lackner,
E. Nardon,
K. Aleynikova,
F. Liu,
E. Strumberger,
R. McAdams,
I. Chapman,
A. Fil
Abstract:
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)…
▽ More
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.
△ Less
Submitted 23 September, 2014; v1 submitted 27 August, 2014;
originally announced August 2014.