For the first time: Realistic simulation of plasma edge instabilities in tokamaks

Credit: CC0 Public Domain

Edge Localized Modes, ELMs for short, are one of the disturbances of the plasma confinement that are caused by the interaction between the charged plasma particles and the confining magnetic field cage. During ELM events, the edge plasma loses its confinement for a short time and periodically throws plasma particles and energy outwards onto the vessel walls. Typically, one tenth of the total energy content can thus be ejected abruptly. While the present generation of medium-sized fusion devices can cope with this, large devices such as ITER or a future power plant would not be able to withstand this strain.

Experimental methods to attenuate, suppress or avoid ELMs have already been successfully developed in current devices (see PI 3/2020). After extensive previous work, it has now been possible for the first time by means of computational simulations to identify the trigger responsible for the explosive onset of these edge instabilities and to reconstruct the course of several ELM cycles—in good agreement with experimentally observed values. A publication accepted in the Nuclear Fusion explains this important prerequisite for predicting and avoiding ELM instabilities in future fusion devices.

The ELM instability builds up after a quiet phase of about 5 to 20 milliseconds—depending on the external conditions—until in half a millisecond between 5 and 15 percent of the energy stored in the plasma is flung onto the walls. Then the equilibrium is restored until the next ELM eruption follows.

The plasma theorists around first author Andres Cathey of IPP, who come from several laboratories of the European fusion program EUROfusion, were able to describe and explain the complex behind this phenomenon in detail: as a non-linear interplay between destabilizing effects—the steep rise in plasma pressure at the plasma edge and the increase in current density—and the stabilizing plasma flow. If the heating power fed into the plasma is changed in the simulation, the calculated result shows the same effect on the repetition rate of the ELMs, i.e. the frequency, as an increase of the heating power in a experiment at ASDEX Upgrade tokamak: experiment and simulation are in Agreement.

Computer simulation of ELM crashes at the lower plasma edge of ASDEX Upgrade fusion device. The video shows the evolution of the plasma pressure over several ELM cycles during 35 milliseconds. Credit: Max Planck Institute for Plasma Physics / Andres Cathey, Matthias Hoelzl

Although the processes take place in a very short time, their simulation requires a great deal of computing effort. This is because the must resolve into small calculation steps both the short ELM crash and the long development phase between two ELMs—a calculation problem that could only be solved with one of the fastest supercomputers currently available.

For the simulations the JOREK code was used, a non-linear code for the calculation of tokamak plasmas in realistic geometry, which is being developed in European and international cooperation with strong contributions from IPP.

More information: Andres Cathey et al, Non-linear extended MHD simulations of type-I edge localised mode cycles in ASDEX Upgrade and their underlying triggering mechanism, Nuclear Fusion (2020). DOI: 10.1088/1741-4326/abbc87

Provided by Max Planck Society

Citation: For the first time: Realistic simulation of plasma edge instabilities in tokamaks (2020, October 22) retrieved 28 May 2024 from
This document is subject to copyright. Apart from any fair dealing for the purpose of private study or research, no part may be reproduced without the written permission. The content is provided for information purposes only.

Explore further

Promising computer simulations for stellarator plasmas


Feedback to editors