Assessing the Impact of Alpha Particles on Thermal Confinement in JET D-T Plasmas through Global GENE-Tango Simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge - despite this being the shot with the highest fusion output - is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization

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

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 ITER24 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

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L-H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’ q-profile) were established on AUG and characterised concerning their plasma performance

Observation of magnetic islands in tokamak plasmas during the suppression of edge-localized modes

In tokamaks, a leading platform for fusion energy, periodic filamentary plasma eruptions known as edge-localized modes occur in plasmas with high-energy confinement and steep pressure profiles at the plasma edge. These edge-localized modes could damage the tokamak wall but can be suppressed using small three-dimensional magnetic perturbations. Here we demonstrate that these magnetic perturbations can change the magnetic topology just inside the steep gradient region of the plasma edge. We identify signatures of a magnetic island, and their observation is linked to the suppression of edge-localized modes. We compare high-resolution measurements of perturbed magnetic surfaces with predictions from ideal magnetohydrodynamic theory where the magnetic topology is preserved. Although ideal magnetohydrodynamics adequately describes the measurements in plasmas exhibiting edge-localized modes, it proves insufficient for plasmas where these modes are suppressed. Nonlinear resistive magnetohydrodynamic modelling supports this observation. Our study experimentally confirms the predicted role of magnetic islands in inhibiting the occurrence of edge-localized modes. This will be beneficial for physics-based predictions in future fusion devices to control these modes

Characterization of the I-phase regime at TCV

The I-phase is an H-mode confinement regime of tokamaks characterized by limit cycle oscillations, the so-called LCOs or bursts. These bursts are the manifestation of a periodic flattening of the plasma edge pressure profile. The profile flattening is caused by increased radial transport, driven by a high-frequency plasma edge mode that periodically appears. This short-living mode is intrinsically connected to each burst. It vanishes once the profiles are fully flattened, and it reestablishes during profile recovery once critical gradients are reached and a new cycle begins. In this paper, we describe for the first time the unambiguous presence of the I-phase at the tokamak `a configuration variable (TCV). As the I-phase confinement regime is found in the parameter regime between the L-mode and the fully developed H-mode, it is often confused with dithers between H-mode and L-mode. Therefore, we are highlighting the differences between these two phenomena. Furthermore, we show the two-dimensional dynamics of the I-phase mode and bursts and the associated filamentary transport, enabled by the outstanding capabilities of the 2D TCV Gas Puff Imaging diagnostics

Peripheral temperature gradient screening of high-Z impurities in optimised 'hybrid' scenario H-mode plasmas in JET-ILW

Screening of high-Z (W) impurities from the confined plasma by the temperature gradient at the plasma periphery of fusion-grade H-mode plasmas has been demonstrated in the JET-ILW (ITER-like wall) tokamak. Through careful optimisation of the hybrid-scenario, deuterium plasmas with sufficient heating power (greater than or similar to 32 MW), high enough ion temperature gradients at the H-mode pedestal top can be achieved for the collisional, neo-classical convection of the W impurities to be directed outwards, expelling them from the confined plasma. Measurements of the W impurity fluxes between and during edge-localised modes (ELMs) based on fast bolometry measurements show that in such plasmas there is a net efflux (loss) between ELMs but that ELMs often allow some W back into the confined plasma. Provided steady, high-power heating is maintained, this mechanism allows such plasmas to sustain high performance, with an average D-D neutron rate of similar to 3.2 x 10(16) s(-1) over a period of similar to 3 s, after an initial overshoot (equivalent to a D-T fusion power of similar to 9.4 MW), without an uncontrolled rise in W impurity radiation, giving added confidence that impurity screening by the pedestal may also occur in ITER, as has previously been predicted (Dux et al 2017 Nucl. Mater. Energy 12 28-35)

Testing a prediction model for the H-mode density pedestal against JET-ILW pedestals

The neutral ionisation model proposed by Groebner et al (2002 Phys. Plasmas 9 2134) to determine the plasma density profile in the H-mode pedestal, is extended to include charge exchange processes in the pedestal stimulated by the ideas of Mahdavi et al (2003 Phys. Plasmas 10 3984). The model is then tested against JET H-mode pedestal data, both in a 'standalone' version using experimental temperature profiles and also by incorporating it in the Europed version of EPED. The model is able to predict the density pedestal over a wide range of conditions with good accuracy. It is also able to predict the experimentally observed isotope effect on the density pedestal that eludes simpler neutral ionization models

Predictive JET current ramp-up modelling using QuaLiKiz-neural-network

This work applies the coupled JINTRAC and QuaLiKiz-neural-network (QLKNN) model on the ohmic current ramp-up phase of a JET D discharge. The chosen scenario exhibits a hollow T-e profile attributed to core impurity accumulation, which is observed to worsen with the increasing fuel ion mass from D to T. A dynamic D simulation was validated, evolving j, n(e), T-e, T-i, n(Be), n(Ni), and n(W) for 7.25 s along with self-consistent equilibrium calculations, and was consequently extended to simulate a pure T plasma in a predict-first exercise. The light impurity (Be) accounted for Z(eff) while the heavy impurities (Ni, W) accounted for Prad. This study reveals the role of transport on the Te hollowing, which originates from the isotope effect on the electron-ion energy exchange affecting T-i. This exercise successfully affirmed isotopic trends from previous H experiments and provided engineering targets used to recreate the D q-profile in T experiments, demonstrating the potential of neural network surrogates for fast routine analysis and discharge design. However, discrepancies were found between the impurity transport behaviour of QuaLiKiz and QLKNN, which lead to notable T-e hollowing differences. Further investigation into the turbulent component of heavy impurity transport is recommended

First-Principles Density Limit Scaling in Tokamaks Based on Edge Turbulent Transport and Implications for ITER

A first-principles scaling law, based on turbulent transport considerations, and a multimachine database of density limit discharges from the ASDEX Upgrade, JET, and TCV tokamaks, show that the increase of the boundary turbulent transport with the plasma collisionality sets the maximum density achievable in tokamaks. This scaling law shows a strong dependence on the heating power, therefore predicting for ITER a significantly larger safety margin than the Greenwald empirical scaling [Greenwald et al., Nucl. Fusion, 28, 2199 (1988)] in case of unintentional high-to-low confinement transition