Speaker
Daniel Cornelis van Vugt
Description
See the full Abstract at http://ocs.ciemat.es/EPS2018ABS/pdf/P1.1049.pdf
Coupled nonlinear MHD-particle simulations for ITER with the
JOREK+particle-tracking code
D.C. van Vugt1,2 , G.T.A. Huijsmans1,3 , M.Hoelzl4 , N.J. Lopes Cardozo1 , A. Loarte2
1 Eindhoven University of Technology, Eindhoven, The Netherlands
2
ITER Organization, 13067 St. Paul Lez Durance Cedex, France
3 CEA Cadarache, IRFM, 13108 St. Paul Lez Durance Cedex, France
4 IPP Garching, Boltzmannstraße 2, 85748 Garching bei München, Germany
There are many effects in tokamak physics which cannot be described with solely a magne-
tohydrodynamic (MHD) code or with particle tracking methods. These include heavy impurity
transport, behaviour of fast ions and neutral particles, radiation modelling, and impurity pro-
duction by sputtering, especially in a non-stationary background plasma and with feedback to
this plasma. Their consistent modelling is very important to reproduce plasma behaviour, par-
ticularly for future tokamaks such as ITER. For instance the behaviour of W (production by
sputtering, transport, ion/recombination balance, radiation emission, etc.) can be modelled in a
stationary or time-varying MHD plasma background by coupling to a particle tracer model.
In this paper we present the extension of the non-linear MHD code JOREK [1] with a particle
tracking code. This code follows particles with a kinetic 6D full-orbit or 5D guiding-center ap-
proximation in the JOREK mesh. Additionally the code contains modules to calculate the ioni-
sation and recombination probabilities of atoms/ions in the plasma, as well as the corresponding
radiated power emitted, with the rates for the atomic processes from ADAS data. Particle colli-
sions (e.g. between W and background DT ions) are modelled with the binary collision model
(BCM) [2], which has been found to reproduce key impurity transport mechanisms such as the
thermal force and the corresponding temperature screening effects. Sputtering sources are im-
plemented using the Eckstein formulation. The particle density and radiation is finally projected
onto the JOREK finite element representation.
The JOREK+particle-tracking code can be used either one-way (evaluation of the conse-
quence of the JOREK plasma behavior on particles) or through a coupled run, where the pro-
jected quantities from the particle distributions are used as source terms in the MHD equations
in JOREK. Such simulations are required, for example, when the modelled W radiation levels
are high and can decrease the plasma temperature, when fast particles affect MHD stability, etc.
Examples of both one-way and coupled simulations will be shown in the paper.
References
[1] G.T.A. Huysmans and O. Czarny. In: Nuclear Fusion 47.7 (2007). DOI: 10.1088/0029-5515/47/7/016.
[2] Y. Homma and A. Hatayama. In: Journal of Computational Physics 231.8 (2012). DOI: 10.1016/j.jcp.
2011.12.037.
