HEPP-Seminar 2020

Raum: Zoom

Hybrid driftkinetic-kinetic implementations and simulations for uniform magnetized space plasma

HEPP Colloquium
Kinetic numerical simulations, applied to study local heating in the solar wind, are computationally expensive due to the different evolution scales involved in the dynamics. Therefore, simplified models, such as hybrid fluid-kinetic and gyrokinetic, are widely employed. Gyrokinetics is missing waves with frequencies above the cyclotron frequencies of the species involved and can be applied only to cases of strong magnetization and where the magnetic moment is conserved. Instead, the hybrid-fluid model is missing electron kinetic effects, which can be important even at ion scales (Told et al., New J. Phys. 2016).Consequently, we are working on a new computationally lighter hybrid model, composed of kinetic ions and gyrokinetic electrons. As a first step, we are considering driftkinetic electrons and uniform magnetic field. The distribution functions are evolved through semilagrangian schemes separately for ions and electrons and coupled through the field, computed using a domain decomposition method and an iterative scheduled relaxation method. We have perfomed turbulence simulations and analysed them in detail to compare them with linear results and with previous works (Tatsuno et al., Phys. Rev. Lett. 2009). [mehr]

MHD simulations of ELM cycles in ASDEX Upgrade

HEPP Colloquium
Edge Localized Modes (ELMs) in tokamaks cause severe concern for future devices like ITER. Large ELMs lead to an expulsion of hot plasma from the edge of the confined region to the tokamak plasma facing components in 0.1–1 milliseconds repetitively every 10–100 milliseconds.Simulations of single ELM crashes with the non-linear 3D magnetohydrodynamic (MHD) code JOREK [GTA Huysmans and O Czarny, NF 47 7 2007] have been validated qualitatively and quantitatively showing good agreement against experimentally observed ELM crashes. Such simulations start with unstable plasma equilibria. To become predictive the entire ELM cycle needs to be simulated. Here, we present simulations of ELM cycles in ASDEX Upgrade and thorough comparisons against experimental measurements. The difficulties related to simulating ELM cycles, how they were overcome with JOREK, and further steps necessary for a better and more comprehensive understanding of ELM dynamics will be discussed at length. [mehr]
Wendelstein 7-X (W7-X) is planned to work at high plasma densities aiming at detached steady-state operation for improved plasma confinement. For plasma densities beyond 1.2×10^20 m-3, electron cyclotron emission (ECE) from the optically thick X2-mode (120-160 GHz) is in the cut-off. Hence, the electron temperature profiles cannot be accessed from ECE for overdense plasmas, but X3-mode is still available to be investigated. A Martin-Puplett interferometer was commissioned in operational phase OP1.2b of W7-X to scan the higher harmonics of ECE for a broad spectral range of 50-500 GHz. The experimental results indicate that X3-mode (180-220 GHz) is optically thick enough to be explored for its diagnostic capabilities as high-density access to electron temperature. The forward modeling of experimental results is done in the bayesian Minerva modeling framework. The results of forward calculations show that X3-mode measurements can be used to provide the electron temperature profiles for overdense plasmas. [mehr]

PICLS: a gyrokinetic full-f particle-in-cell code for open field line simulations

HEPP Colloquium
While in recent years, gyrokinetic simulations have become the workhorse for theoretical turbulence and transport studies in the plasma core, their application to the edge and scrape-off layer (SOL) region presents significant challenges. The ``full-f" code PICLS has been developed, to in particular study the SOL region with its steep density and temperature gradients as well as large fluctuation amplitudes. PICLS is based on an electrostatic full-f model with a linearized field equation and uses kinetic electrons. The electrostatic potential is calculated via the polarization equation, with the help of B-spline finite-elements for the charge deposition and the field solver. In this talk, we will introduce the PICLS model and show our results of applying it to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) for the non-collisional and collisional case. Our current progress on extending PICLS towards three spatial domains, will be presented and key features for the 3D extension, such as field solver and particle pusher, will be shown. [mehr]

Deep Learning for Tokamak Plasma Confinement Mode Classification

HEPP Colloquium
During a discharge at the TCV tokamak, the plasma be classified as varying between Low (L), High (H) and, in some cases, a temporary (intermediate) mode, called Dithering (D). In addition, while the plasma is in H mode, Edge Localized Modes (ELMs) can occur. The ability to accurately, and automatically, detect changes between these states, and ELMs, is considered important for future tokamak operation. However, it is difficult to design a traditional rulebased system that can accurately account for all the possible reasons behind these phenomena.The alternative is to use an approach whereby data generated in fusion experiments is used by algorithms which can, by themselves, learn the underlying rules that explain these events. Deep learning algorithms are exactly suited for this task. By feeding them with enough data, these models can automatically extract any existing correlations that allow for accurately detecting plasma confinement states and ELMs. In this work, we present a series of different deep learning algorithms for this task, namely, convolutional neural networks, recurrent neural networks, and sequence to sequence encoder-decoder models. The algorithms presented differ from each other with regards to the assumptions made regarding the data that they process, their architectures, and their capacity to accurately carry out the classification task. We will show, in particular, that a sequence to sequence model can achieve the best results, while also allowing for explicit incorporation of domain knowledge into the classification task. [mehr]

Fast-ion confinement studies in W7-X using active Balmer-alpha spectroscopy

HEPP Colloquium
Investigation of fast-ion transport in fusion plasmas plays a central role as good fast-particle confinement is essential for burning plasmas. For its recent experimental campaign, the W7-X stellarator was equipped with two neutral hydrogen beam injectors (NBI) which provide fast-particles of our interest and neutral hydrogen particles which make active Balmer-alpha spectroscopy possible.Due to the complex shape of the measured spectra forward modelling is required for its interpretation. It is done with a code called FIDASIM which takes into account different assumed kinetic plasma profiles, fast-ion densities and information of the beam- and observation geometry.It was found that most measured spectral components are well reproducible with FIDASIM but the observed active emission, coming from the beam neutral - confined fast-ion charge-exchange reaction (FIDA) cannot fully explain the measured intensities. This suggests that fast-ions in the plasma edge region could interact with the cold neutral population from the plasma vessel, causing additional passive FIDA emission. This needs to be understood in order to address the question of edge charge-exchange fast-ion losses and to infer information on the edge neutral density. [mehr]

Carbon distribution and transport in ECRH and NBI heated plasmas with Charge Exchange Spectroscopy on W7-X

HEPP Colloquium
Impurity transport plays a crucial role in the optimization of fusion plasmas, as impurities affect the plasma radiation and can cause power losses. If neoclassical effects dominate the transport, strong impurity accumulation is predicted in the plasma core. According to simulations, neoclassically dominated impurity transport is a possibility in the optimized stellarator Wendelstein 7-X (W7-X) plasma. To quantify impurity confinement in W7-X, carbon concentration profiles are investigated and used with the impurity transport modeling code STRAHL to determine the transport coefficients (diffusivity and radial convective velocity). The results are compared with neoclassical predictions in order to assess the anomalous contribution. The profiles are derived from the Charge Exchange Recombination Spectroscopy (CXRS) diagnostic that observes the Neutral Beam Injection (NBI) which is well-suited for determining spatially resolved profiles of fully-stripped low-Z impurities. This work concentrates on carbon, the main intrinsic impurity in W7-X. Different configurations, densities and heating scenarios with different NBI and ECRH power ratios are explored. Of particular interest are discharges with pure NBI heating phases or with very low ECRH power, where indications of unusually high impurity confinement times have been observed. [mehr]

Gyrokinetic investigation of the damping channels of Alfvén modes in ASDEX Upgrade

HEPP Colloquium
The linear destabilization and nonlinear saturation of energetic-particle driven Alfvénic instabilities in tokamaks strongly depend on the damping channels. In this work, the collisionless damping mechanisms of Alfvénic modes are investigated within a gyrokinetic framework, by means of global simulations with the particle-in-cell code ORB5, and compared with the eigenvalue code LIGKA and reduced models. In particular, the continuum damping and the Landau damping (of ions and electrons) are considered. The electron Landau damping is found to be dominant on the ion Landau damping for experimentally relevant cases. As an application, the linear and nonlinear dynamics of toroidicity induced Alfvén eigenmodes and energetic-particle driven modes in ASDEX Upgrade is investigated theoretically and compared with experimental measurements. [mehr]

Geometry and Kinetics of Astrophysical Plasmas: A gyrokinetic approach

HEPP Colloquium
In the context of astrophysical plasmas, various methods are used in order to study problems such as dissipation of turbulent energy and magnetic reconnection. The use of fluid models allow us to understand macroscopic phenomena, but lacks the dynamics of kinetic physics. On another hand, kinetic models usually consume an enormous amount of computing time. The use of reduced models such as gyrokinetics are foreseen to bridge the gap between the fluid and kinetic approaches. In the present work, we aim to investigate the use of gyrokinetics in two different scenarios. Firstly we are going to consistently derive a hybrid hamiltonian field theoretical system, based on the lagragian formulation of a symplectic two-form. In this system, ions are treated fully kinetically, and electrons gyrokinetically. With this model, we wish to develop a cost effective kinetic computational framework. The second aspect of the present work addresses a well known problem in space physics, namely magnetic reconnection with guide field. We start with a gyrokinetic analysis using the code GENE. Firstly we analyze the dynamics of the parallel electric field and reconnection rate on the X point, and proceed with benchmarking GENE with a fully kinetic PIC code. [mehr]

Gyrokinetic investigation of the ASDEX Upgrade I-mode pedestal

HEPP Colloquium
Characterizing pedestal turbulence in the tokamak I-mode is a crucial step in understanding how particle and heat transport decouple during I-mode operation. This work models an ASDEX Upgrade I-mode discharge for the first time via linear and nonlinear gyrokinetic simulations with the GENE code. L-mode and I-mode regimes at two different pedestal locations are investigated. A microtearing mode which is not apparent in initial value linear L-mode simulations is found to dominate in I-mode simulations at both radial positions, and ion-scale instabilities are characterized for all four scenarios linearly. Computed nonlinear heat flux values approach experimental measurements with nominal input parameters in three of the four cases, and heat transport is found to be dominated by ion-scale electrostatic turbulence. Electrostatic potential oscillation frequencies, as well as potential-temperature and potential-density crossphases are compared linearly and nonlinearly, and agreement is found at wavenumber ranges corresponding with peaks in the simulated heat flux spectra at one radial position for L-mode and I-mode. [mehr]

Runaway electron modeling in massive material injection scenarios in ASDEX Upgrade

HEPP Colloquium
In current-carrying fusion devices, the conversion of a large fraction of the plasma current to runaway electrons (RE) following the sudden loss of thermal energy poses a threat to the integrity of the plasma vessel. However, RE formation may be suppressed by massive material injection (MMI); a concept presently being investigated experimentally across various devices. To complement extrapolation to future devices, RE model development and validation is mandatory. In this talk, we present the 1.5D transport toolset ASTRA-STRAHL for the simulation of RE generation during MMI scenarios. The toolset is then applied in first-time integrated simulations of MMI, background plasma response, and RE generation in ASDEX Upgrade (AUG). Employing state-of-the-art models for RE generation considering the impact of partially ionized impurities, the evolution of key plasma parameters (plasma current decay, line integrated electron density, etc.) is calculated well in agreement with experimental observations of AUG. Considering instead commonly used formulae without the impact of these impurities, simulations cannot capture experimental trends, thus demonstrating the importance of these kinetic effects on RE generation. [mehr]

Quantitative investigation of the neutron production in ASDEX Upgrade

HEPP Colloquium
The neutron production in ASDEX Upgrade(AUG) neutral beam injection (NBI) heated discharges is dominated by beam-target fusion reactions. Hence, the neutron rate and energy distributions are footprints of the fast ion distribution. This motivates the achievement of reliable neutron measurements and good agreement between experiment and theoretical calculations. However, comparisons at AUG between the experimental neutron rate and the one predicted by TRANSP show systematic discrepancies. Potential reason for this is the delicate absolute calibration of the neutron detectors. Therefore, a different calibration technique was performed, allowing for longer calibration time, better statistics and thus less uncertainty. A toy train carrying a radioactive source (238Pu/B) over two poloidal positions on the equatorial plane shows a clearly periodical neutron rate on the epithermal 3He neutron detector. The calibration results are compared to a Monte Carlo (MC) transport simulation using the Serpent code. [mehr]
Zur Redakteursansicht