Frederic Imbeaux

A unified set of conventions and normalizations for the inputs and outputs of delta-f flux-tube gyrokinetic simulations has been developed within the ITER Integrated Modelling and Analysis Suite (IMAS). This new standard aims at facilitating benchmarks and data exchange between different code users, increasing the modularity of frameworks manipulating inputs and outputs from gyrokinetic simulations and serving as the basis for databases storing simulations from multiple gyrokinetic codes.

The WEST tokamak is equipped with a superconducting toroidal magnetic field system, a multi-megawatt radiofrequency auxiliary power system, and an actively cooled ITER-grade tungsten divertor. As such, it is well adapted to explore experimental aspects related to the long pulse operation of next-step devices. Supported by predict-first integrated modeling, bespoke scenario development has allowed zero-loop voltage pulses to be achieved. The resulting discharges, with plasma currents in the range I p ∼ 0.22 − 0.28 MA exclusively sustained by the lower hybrid current drive (LHCD) system as an auxiliary power source, have achieved durations in excess of 22 min and injected/extracted energies up to 2.61 GJ. Plasma performance is characterized by ranges of poloidal beta β p ∼ 1.6 − 2.0 , normalized toroidal beta β N ∼ 0.6 − 0.9 and confinement factor H 96 L = 1 − 1.3 . Mild MHD activity, identified as resulting from the interaction of 3/1 and 4/1 tearing modes, is occasionally present, depending on the LHCD antenna combination used. This article describes the predict-first approach thathas been employed in the context of this long-pulse scenario development endeavor. The main achievements and the physics analyses performed are reviewed, including post-experiment integrated modeling aspects. Prospects for further long-pulse developments are drawn.

Bridging the gap to next-step devices – while saving valuable time and resources – requires more than semi-empirical models, which struggle to predict plasma behavior in unexplored parameter regimes. Instead, validated simulation tools are essential, leveraging high-fidelity exascale computing and multi-fidelity models, including AI-based surrogates. To address these challenges, the “EUROfusion Theory and Advanced Simulation Coordination (E-TASC)” initiative was launched in 2021. It includes 15 TSVV (Theory, Simulation, Verification, and Validation) projects supported by five Advanced Computing Hubs. This “team of teams” has made substantial progress toward developing digital twins of fusion systems, with key scientific achievements to be presented in the paper. This includes the following topical areas: core performance in burning plasmas, magnetohydrodynamic transients, L-H transitions and ELM-free regimes, plasma exhaust, and plasma-wall interactions – with applications to both tokamaks and stellarators.

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

This paper is a written summary of an overview oral presentation given at the 1st Spanish Fusion High Performance Computer (HPC) Workshop that took place on the 27 November 2020 as an online event. Given that over the next few years ITER ²⁴ 24 ITER (‘The Way’ in Latin) is the world’s largest tokamak under construction in the south of France: a magnetic fusion device that has been designed to prove the feasibility of fusion as a large-scale and carbon-free source of energy (https://www.iter.org/). 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, which stands for EUROfusion-Theory and Advanced Simulation Coordination. The general and guiding principles of E-TASC are summarized in this paper. In addition, an overview of the scientific results obtained in the 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: towards a validated predictive capability of the low to high transition and pedestal physics; runaway electrons in tokamak disruptions in the presence of massive material injection; fast code for the calculation of neoclassical toroidal viscosity in stellarators and tokamaks; development of a neutral gas kinetics modular code; European edge and boundary code for reactor-relevant devices. In this paper, we report on recent progress made by each of these projects.

Self-consistent full-size turbulent-transport simulations of the divertor and SOL of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted Ohmic L-mode discharges on the TCV tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from 5 diagnostic systems -- giving a total of 45 1- and 2-D comparison observables in two toroidal magnetic field directions. The dataset is used to validate three flux-driven 3D fluid-turbulence models: GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP -- both in terms of profile shape and absolute magnitude -- while a poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

Integrating several important technological elements required for long pulse operation in magnetic fusion devices, the Tore Supra tokamak routinely addresses the physics and technology issues related to this endeavor and, as a result, contributes essential information on critical issues for ITER. During the last experimental campaign, components of the radiofrequency system including an ITER relevant launcher (passive active multijunction (PAM)) and continuous wave/3.7?GHz klystrons, have been extensively qualified, and then used to develop steady state scenarios in which the lower hybrid (LH), ion cyclotron (IC) and electron cyclotron (EC) systems have been combined in fully stationary shots (duration ?150?s, injected power up to ?8?MW, injected/extracted energy up to ?1?GJ). Injection of LH power in the 5.0?6.0?MW range has extended the domain of accessible plasma parameters to higher densities and non-inductive currents. These discharges exhibit steady electron internal transport barriers (ITBs). We report here on various issues relevant to the steady state operation of future devices, ranging from operational aspects and limitations related to the achievement of long pulses in a fully actively cooled fusion device (e.g. overheating due to fast particle losses), to more fundamental plasma physics topics. The latter include a beneficial influence of IC resonance heating on the magnetohydrodynamic (MHD) stability in these discharges, which has been studied in detail. Another interesting observation is the appearance of oscillations of the central temperature with typical periods of the order of one to several seconds, caused by a nonlinear interplay between LH deposition, MHD activity and bootstrap current in the presence of an ITB.