EUROfusion-theory and advanced simulation coordination (E-TASC): programme and the role of high performance computing

The paper is a written summary of an overview oral presentation given at the 1st Spanish Fusion HPC Workshop that took place on the 27th November 2020 as an online event. Given that over the next few years ITER will move to its operation phase and the European-DEMO design will be significantly advanced, the EUROfusion consortium has initiated a coordination effort in theory and advanced simulation to address some of the challenges of the fusion research in Horizon EUROPE (2021-2027), i.e. the next EU Framework Programme for Research and Technological Development. This initiative has been called E-TASC that stands for EUROfusion-Theory and Advanced Simulation Coordination. The general and guiding principles of E-TASC are summarized in the paper. In addition, an overview of the scientific results obtained in a pilot phase (2019-2020) of E-TASC are provided while highlighting the importance of the required progress in computational methods and HPC techniques. In the initial phase, five pilot theory and simulation tasks were initiated: 1. Towards a validated predictive capability of the L-H transition and pedestal physics; 2. Electron runaway in tokamak disruptions in the presence of massive material injection; 3. Fast code for the calculation of neoclassical toroidal viscosity in stellarators and tokamaks; 4. Development of a neutral gas kinetics modular code; 5. European edge and boundary code for reactor-relevant devices. In this paper we report on recent progress made by each of these projects.

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The Generalized Perturbed Equilibrium Code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. This work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.