Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSelf-consistent modeling of runaway electron generation in massive gas injection scenarios in ASDEX Upgrade
69
0
0.0
(
0
)
Authors
O. Linder
Ask ChatGPT about the research
No Arabic abstract
We present the first successful simulation of a induced disruption in ASDEX Upgrade from massive material injection (MMI) up to established runaway electron (RE) beam, thus covering pre-thermal quench, thermal quench and current quench (CQ) of the discharge. For future high-current fusion devices such as ITER, the successful suppression of REs through MMI is of critical importance to ensure the structural integrity of the vessel. To computationally study the interplay between MMI, background plasma response, and RE generation, a toolkit based on the 1.5D transport code coupling ASTRA-STRAHL is developed. Electron runaway is described by state-of-the-art reduced kinetic models in the presence of partially ionized impurities. Applied to argon MMI in ASDEX Upgrade discharge #33108, key plasma parameters measured experimentally, such as temporal evolution of the line averaged electron density, plasma current decay rate and post-CQ RE current, are well reproduced by the simulation presented. Impurity ions are transported into the central plasma by the combined effect of neoclassical processes and additional effects prescribed inside the $q = 2$ rational surface to explain experimental time scales. Thus, a thermal collapse is induced through strong impurity radiation, giving rise to a substantial RE population as observed experimentally.
rate research
Read More
The formation of a substantial post-disruption runaway electron current in ASDEX Upgrade material injection experiments is determined by avalanche multiplication of a small seed population of runaway electrons. For the investigation of these scenarios, the runaway electron description of the coupled 1.5D transport solvers ASTRA-STRAHL is amended by a fluid-model describing electron runaway caused by the hot-tail mechanism. Applied in simulations of combined background plasma evolution, material injection, and runaway electron generation in ASDEX Upgrade discharge #33108, both the Dreicer and hot-tail mechanism for electron runaway produce only $sim$ 3$~$kA of runaway current. In colder plasmas with core electron temperatures $T_mathrm{e,c}$ below 9$~$keV, the post-disruption runaway current is predicted to be insensitive to the initial temperature, in agreement with experimental observations. Yet in hotter plasmas with $T_mathrm{e,c} > 10~mathrm{keV}$, hot-tail runaway can be increased by up to an order of magnitude, contributing considerably to the total post-disruption runaway current. In ASDEX Upgrade high temperature runaway experiments, however, no runaway current is observed at the end of the disruption, despite favourable conditions for both primary and secondary runaway.
First simulations of deuterium shattered pellet injection (SPI) into an ASDEX Upgrade H-Mode plasma with the JOREK MHD code are presented. Resistivity is increased by one order of magnitude in most simulations to reduce computational costs and allow for extensive parameter scans. The effect of various physical parameters onto MHD activity and thermal quench (TQ) dynamics is studied and the influence of MHD onto ablation is shown. TQs are obtained quickly after injection in most simulations with a typical duration of 100 microseconds, which slows down at lower resistivity. Although the n=1 magnetic perturbation dominates in the simulations, toroidal harmonics up to n=10 contribute to stochastization and stochastic transport in the plasma core. The post-TQ density profile remains hollow for a few hundred microseconds. However, when flux surfaces re-form around the magnetic axis, the density has become monotonic again suggesting a beneficial behaviour for runaway electron avoidance/mitigation. With $10^{21}$ atoms injected, the TQ is typically incomplete and triggered when the shards reach the q=2 rational surface. At a larger number of injected atoms, the TQ can set in even before the shards reach this surface. For low field side injection considered here, repeated formation of outward convection cells is observed in the ablation region reducing material assimilation. Responsible is a sudden rise of pressure in the high density cloud when the stochastic region expands further releasing heat from the hot core. After the TQ, strong sheared poloidal rotation is created by Maxwell stress, which contributes to re-formation of flux surfaces.
Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quench, together with an avalanche runaway tail which has an exponentially decreasing energy spectrum. We find that, although the avalanche component carries the vast majority of the current, it is the high-energy seed remnant that dominates synchrotron emission. With insights from the fluid-kinetic simulations, an analytic model for the evolution of the runaway seed component is developed and used to reconstruct the radial density profile of the RE beam. The analysis shows that the observed change of the synchrotron pattern from circular to crescent shape is caused by a rapid redistribution of the radial profile of the runaway density.
The plasma response to Resonant Magnetic Perturbations (RMPs) in ASDEX Upgrade is modeled with the non-linear resistive MHD code JOREK, using input profiles that match those of the experiments as closely as possible. The RMP configuration for which Edge Localized Modes are best mitigated in experiments is related to the largest edge kink response observed near the X-point in modeling. On the edge resonant surfaces $q = m/n$, the coupling between the m + 2 kink component and the m resonant component is found to induce the amplification of the resonant magnetic perturbation. The ergodicity and the 3D-displacement near the X-point induced by the resonant amplification can only partly explain the density pumpout observed in experiments.
Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concern. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
