Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userComparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations
165
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. 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. Nuclear Fusion 60, 124007 (2020)]. Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to capture the balance between stabilising and destabilising terms for the edge instabilities accurately. 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 narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluency is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.
rate research
Read More
Pellet ELM triggering is a well established scheme for decreasing the time between two successive ELM crashes below its natural value. Reliable ELM pacing has been demonstrated experimentally in several devices increasing the ELM frequency considerably. However, it was also shown that the frequency cannot be increased arbitrarily due to a so-called lag-time. During this time after a preceding natural or triggered ELM crash, neither a natural ELM crash occurs nor the triggering of an ELM crash by pellet injection is possible. For this article, pellet ELM triggering simulations are advanced beyond previous studies in two ways. Firstly, realistic ExB and diamagnetic background flows are included. And secondly, the pellet is injected at different stages of the pedestal build-up. This allows to recover the lag-time for the first time in simulations and investigate it in detail. A series of non-linear extended MHD simulations is performed to investigate the plasma dynamics resulting from an injection at different time points during the pedestal build-up. The experimentally observed lag-time is qualitatively reproduced well. In particular, a sharp transition is observed between the regime where no ELMs can be triggered and the regime where pellet injection causes an ELM crash. Via variations of pellet parameters and injection time, the two regimes are studied and compared in detail revealing pronounced differences in the non-linear dynamics. The toroidal mode spectrum is significantly broader when an ELM crash is triggered enhancing the stochasticity and therefore also the losses of thermal energy along magnetic field lines. In the heat fluxes to the divertor targets, pronounced toroidal asymmetries are observed. In case of high injection velocities leading to deep penetration, also the excitation of core modes like the $2/1$ neoclassical tearing mode is observed.
A triggering mechanism responsible for the explosive onset of edge localised modes (ELMs) in fusion plasmas is identified by performing, for the first time, non-linear magnetohydrodynamic simulations of repetitive type-I ELMs. Briefly prior to the ELM crash, destabilising and stabilising terms are affected at different timescales by an increasingly ergodic magnetic field caused by non-linear interactions between the axisymmetric background plasma and growing non-axisymmetric perturbations. The separation of timescales prompts the explosive, i.e. faster than exponential, growth of an ELM crash which lasts ${sim}$ 500 ${mu}$s. The duration and size of the simulated ELM crashes compare qualitatively well with type-I ELMs in ASDEX Upgrade. As expected for type-I ELMs, a direct proportionality between the heating power in the simulations and the ELM repetition frequency is obtained. The simulations presented here are a major step forward towards predictive modelling of ELMs and of the assessment of mitigation techniques in ITER and other future tokamaks.
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.
The radiation response and the MHD destabilization during the thermal quench after a mixed species Shattered Pellet Injection (SPI) with impurity species neon and argon are investigated via 3D non-linear MHD simulation using the JOREK code. Both the $n=0$ global current profile contraction and the local helical cooling at each rational surface caused by the pellet fragments are found to be responsible for MHD destabilization after the injection. Significant current driven mode growth is observed as the fragments cross low order rational surfaces, resulting in rapidly inward propagating stochastic magnetic field, ultimately causing the core temperature collapse. The Thermal Quench (TQ) is triggered as the fragments arrive on the $q=1$ or $q=2$ surface depending on the exact $q$ profile and thus mode structure. When injecting from a single toroidal location, strong radiation asymmetry is found before and during the TQ as a result of the unrelaxed impurity density profile along the field line and asymmetric outward heat flux. Such asymmetry gradually relaxes over the course of the TQ, and is entirely eliminated by the end of it. Simulation results indicate that the aforementioned asymmetric radiation behavior could be significantly mitigated by injection from toroidally opposite locations, provided that the time delay between the two injectors is shorter than $1ms$. It is also found that the MHD response are sensitive to the relative timing and injection configuration in these multiple injection cases.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
