Recent ASDEX upgrade research in support of ITER and DEMO

Recent experiments on the ASDEX Upgrade tokamak aim at improving the physics base for ITER and DEMO to aid the machine design and prepare efficient operation. Type I edge localized mode (ELM) mitigation using resonant magnetic perturbations (RMPs) has been shown at low pedestal collisionality (nu(*)(ped) < 0.4). In contrast to the previous high nu* regime, suppression only occurs in a narrow RMP spectral window, indicating a resonant process, and a concomitant confinement drop is observed due to a reduction of pedestal top density and electron temperature. Strong evidence is found for the ion heat flux to be the decisive element for the L-H power threshold. A physics based scaling of the density at which the minimum P-LH occurs indicates that ITER could take advantage of it to initiate H-mode at lower density than that of the final Q = 10 operational point. Core density fluctuation measurements resolved in radius and wave number show that an increase of R/L-Te introduced by off-axis electron cyclotron resonance heating (ECRH) mainly increases the large scale fluctuations. The radial variation of the fluctuation level is in agreement with simulations using the GENE code. Fast particles are shown to undergo classical slowing down in the absence of large scale magnetohydrodynamic (MHD) events and for low heating power, but show signs of anomalous radial redistribution at large heating power, consistent with a broadened off-axis neutral beam current drive current profile under these conditions. Neoclassical tearing mode (NTM) suppression experiments using electron cyclotron current drive (ECCD) with feedback controlled deposition have allowed to test several control strategies for ITER, including automated control of (3,2) and (2,1) NTMs during a single discharge. Disruption mitigation studies using massive gas injection (MGI) can show an increased fuelling efficiency with high field side injection, but a saturation of the fuelling efficiency is observed at high injected mass as needed for runaway electron suppression. Large locked modes can significantly decrease the fuelling efficiency and increase the asymmetry of radiated power during MGI mitigation. Concerning power exhaust, the partially detached ITER divertor scenario has been demonstrated at P-sep/R = 10 MW m(-1) in ASDEX Upgrade, with a peak time averaged target load around 5 MW m(-2), well consistent with the component limits for ITER. Developing this towards DEMO, full detachment was achieved at P-sep/R = 7 MW m(-1) and stationary discharges with core radiation fraction of the order of DEMO requirements (70% instead of the 30% needed for ITER) were demonstrated. Finally, it remains difficult to establish the standard ITER Q = 10 scenario at low q(95) = 3 in the all-tungsten (all-W) ASDEX Upgrade due to the observed poor confinement at low beta(N). This is mainly due to a degraded pedestal performance and hence investigations at shifting the operational point to higher beta(N) by lowering the current have been started. At higher q(95), pedestal performance can be recovered by seeding N-2 as well as CD4, which is interpreted as improved pedestal stability due to the decrease of bootstrap current with increasing Z(eff). Concerning advanced scenarios, the upgrade of ECRH power has allowed experiments with central ctr-ECCD to modify the q-profile in improved H-mode scenarios, showing an increase in confinement at still good MHD stability with flat elevated q-profiles at values between 1.5 and 2

Role of fast ion pressure in the isotope effect in JET L-mode plasmas

This paper presents results of JET ITER-like wall L-mode experiments in hydrogen and deuterium (D) plasmas, dedicated to the study of the isotope dependence of ion heat transport by determination of the ion critical gradient and stiffness by varying the ion cyclotron resonance heating power deposition. When no strong role of fast ions in the plasma core is expected, the main difference between the two isotope plasmas is determined by the plasma edge and the core behavior is consistent with a gyro-Bohm scaling. When the heating power (and the fast ion pressure) is increased, in addition to the difference in the edge region, also the plasma core shows substantial changes. The stabilization of ion heat transport by fast ions, clearly visible in D plasmas, appears to be weaker in H plasmas, resulting in a higher ion heat flux in H with apparent anti-gyro-Bohm mass scaling. The difference is found to be caused by the different fast ion pressure between H and D plasmas, related to the heating power settings and to the different fast ion slowing down time, and is completely accounted for in non-linear gyrokinetic simulations. The application of the TGLF quasi-linear model to this set of data is also discussed

Long-lived coupled peeling ballooning modes preceding ELMs on JET

In some JET discharges, type-I edge localised modes (ELMs) are preceded by a class of low-frequency oscillations (Perez et al 2004 Nucl. Fusion 44 609). While in many cases the ELM is triggered during the growth phase of this oscillation, it is also observed that this type of oscillation can saturate and last for several tens of ms until an ELM occurs. In order to identify the nature of these modes, a wide pre-ELM oscillation database, including detailed pedestal profile information, has been assembled and analysed in terms of MHD stability parameters. The existence domain of these pre-ELM oscillations and the statistical distribution of toroidal mode numbers (n) up to n = 16 have been mapped in ballooning alpha (alpha(ball)) and either edge current density (J(edge)) or pedestal collisionality (nu(ee,ped)*) coordinates and compared to linear MHD stability predictions. The pre-ELM oscillations are reliably observed when the J/alpha ratio is high enough for the pedestal to access the coupled peeling-ballooning (PB) domain (aka stability nose). Conversely, when the pedestal is found to be in or near the high-n ballooning domain (which is at low J/alpha ratio), ELMs are usually triggered promptly, i.e. with no detectable pre-ELM oscillations, or with pre-ELM oscillations only observable on ECE whose n appears to be too high to be resolved by the magnetics. Individual discharges can sometimes exhibit a fairly wide range of pre-ELM mode numbers, but for a wider database, the statistical n-number domains are found to be well ordered along the J - alpha stability boundary and behave as expected from PB theory: the higher the J/alpha ratio, the lower the mode's measured n tends to be. Within the measurement uncertainties, the measured n is usually found to be compatible with the most unstable n predicted by the linear stability code MISHKA1. These results confirm the earlier hypothesis that these modes are coupled peeling-ballooning modes, and extend and generalise to higher-mode numbers the work by Huysmans et al (1998 Nucl. Fusion 38 179), who identified the lowest n modes as pure external kink modes. Since the destabilisation of PB modes is widely accepted to give rise to ELMs, the mode saturation and delayed ELM triggering that is sometimes observed is rather unexpected. Possibilities to reconcile the extended lifetime of these modes with current ELM models are briefly discussed, but will require further investigation

Diagnostic of fast-ion energy spectra and densities in magnetized plasmas

The measurement of the energy spectra and densities of alpha-particles and other fast ions are part of the ITER measurement requirements, highlighting the importance of energy-resolved energetic-particle measurements for the mission of ITER. However, it has been found in recent years that the velocity-space interrogation regions of the foreseen energetic-particle diagnostics do not allow these measurements directly. We will demonstrate this for gamma-ray spectroscopy (GRS), collective Thomson scattering (CTS), neutron emission spectroscopy and fast-ion D-alpha spectroscopy by invoking energy and momentum conservation in each case, highlighting analogies and differences between the different diagnostic velocity-space sensitivities. Nevertheless, energy spectra and densities can be inferred by velocity-space tomography which we demonstrate using measurements at JET and ASDEX Upgrade. The measured energy spectra agree well with corresponding simulations. At ITER, alpha-particle energy spectra and densities can be inferred for energies larger than 1.7 MeV by velocity-space tomography based on GRS and CTS. Further, assuming isotropy of the alpha-particles in velocity space, their energy spectra and densities can be inferred by 1D inversion of spectral single-detector measurements down to about 300 keV by CTS. The alpha-particle density can also be found by fitting a model to the CTS measurements assuming the alpha-particle distribution to be an isotropic slowing-down distribution

EDGE2D-EIRENE simulations of the influence of isotope effects and anomalous transport coefficients on near scrape-off layer radial electric field

EDGE2D-EIRENE (the 'code') simulations show that radial electric field, Er, in the near scrape-off layer (SOL) of tokamaks can have large variations leading to a strong local E x B shear greatly exceeding that in the core region. This was pointed out in simulations of JET plasmas with varying divertor geometry, where the magnetic configuration with larger predicted near SOL E-r was found to have lower H-mode power threshold, suggesting that turbulence suppression in the SOL by local E. x. B shear can be a player in the L-H transition physics (Delabie et al 2015 42nd EPS Conf. on Plasma Physics (Lisbon, Portugal, 22-26 June 2015) paper O3.113 (http://ocs.ciemat.es/EPS2015PAP/pdf/O3.113.pdf), Chankin et al 2017 Nucl. Mater. Energy 12 273). Further code modeling of JET plasmas by changing hydrogen isotopes (H-D-T) showed that the magnitude of the near SOL E-r is lower in H cases in which the H-mode threshold power is higher (Chankin et al 2017 Plasma Phys. Control. Fusion 59 045012). From the experiment it is also known that hydrogen plasmas have poorer particle and energy confinement than deuterium plasmas, consistent with the code simulation results showing larger particle diffusion coefficients at the plasma edge, including SOL, in hydrogen plasmas (Maggi et al 2018 Plasma Phys. Control. Fusion 60 014045). All these experimental observations and code results support the hypothesis that the near SOL E x B shear can have an impact on the plasma confinement. The present work analyzes neutral ionization patterns of JET plasmas with different hydrogen isotopes in L-mode cases with fixed input power and gas puffing rate, and its impact on target electron temperature, T-e, and SOL E-r. The possibility of a self-feeding mechanism for the increase in the SOL E-r via the interplay between poloidal E x B drift and target T-e is discussed. It is also shown that reducing anomalous turbulent transport coefficients, particle diffusion and electron and ion heat conductivities, leads to higher peak target T-e and larger E-r, suggesting the possibility of a positive feedback loop, under an implicitly made assumption that the E x B shear in the SOL is capable of suppressing turbulence

A power-balance model of the density limit in fusion plasmas: application to the L-mode tokamak

A power-balance model, with radiation losses from impurities and neutrals, gives a unified description of the density limit (DL) of the stellarator, the L-mode tokamak, and the reversed field pinch (RFP). The model predicts a Sudo-like scaling for the stellarator, a Greenwald- like scaling, alpha I-p(8/9), for the RFP and the ohmic tokamak, a mixed scaling, alpha (PIp4/9)-I-4/9, for the additionally heated L-mode tokamak. In a previous paper (Zanca et al 2017 Nucl. Fusion 57 056010) the model was compared with ohmic tokamak, RFP and stellarator experiments. Here, we address the issue of the DL dependence on heating power in the L-mode tokamak. Experimental data from high-density disrupted L-mode discharges performed at JET, as well as in other machines, arc taken as a term of comparison. The model fits the observed maximum densities better than the pure Greenwald limit

A new mechanism for increasing density peaking in tokamaks: improvement of the inward particle pinch with edge E x B shearing

Developing successful tokamak operation scenarios, as well as confident extrapolation of present-day knowledge requires a rigorous understanding of plasma turbulence, which largely determines the quality of the confinement. In particular, accurate particle transport predictions are essential due to the strong dependence of fusion power or bootstrap current on the particle density details. Here, gyrokinetic turbulence simulations are performed with physics inputs taken from a JET power scan, for which a relatively weak degradation of energy confinement and a significant density peaking is obtained with increasing input power. This way physics parameters that lead to such increase in the density peaking shall be elucidated. While well-known candidates, such as the collisionality, previously found in other studies are also recovered in this study, it is furthermore found that edge E x B shearing may adopt a crucial role by enhancing the inward pinch. These results may indicate that a plasma with rotational shear could develop a stronger density peaking as compared to a non-rotating one, because its inward convection is increased compared to the outward diffusive particle flux as long as this rotation has a significant on E x B flow shear stabilization. The possibly significant implications for future devices, which will exhibit much less torque compared to present day experiments, are discussed

Ion cyclotron resonance heating scenarios for DEMO

The present paper offers an overview of the potential of ion cyclotron resonance heating (ICRH) or radio frequency heating for the DEMO machine. It is found that various suitable heating schemes are available. Similar to ITER and in view of the limited bandwidth of about 10 MHz that can be achieved to ensure optimal functioning of the launcher, it is proposed to make core second harmonic tritium heating the key ion heating scheme, assisted by fundamental cyclotron heating He-3 in the early phase of the discharge; for the present design of DEMO-with a static magnetic field strength of B-o = 5.855 T-that places the T and 3He layers in the core for f = 60 MHz and suggests centering the bandwidth around that main operating frequency. In line with earlier studies for hot, dense plasmas in large-size magnetic confinement machines, it is shown that good single pass absorption is achieved but that the size as well as the operating density and temperature of the machine cause the electrons to absorb a non-negligible fraction of the power away from the core when core ion heating is aimed at. Current drive and alternative heating options are briefly discussed and a dedicated computation is done for the traveling wave antenna, proposed for DEMO in view of its compatibility with substantial antenna-plasma distances. The various tasks that ICRH can fulfill are briefly listed. Finally, the impact of transport and the sensitivity of the obtained results to changes in the machine parameters is commented on