Integrated modelling and multiscale gyrokinetic validation study of ETG turbulence in a JET hybrid H-mode scenario

Previous studies with first-principle-based integrated modelling suggested that ETG turbulence may lead to an anti-GyroBohm isotope scaling in JET high-performance hybrid H-mode scenarios. A dedicated comparison study against higher-fidelity turbulence modelling invalidates this claim. Ion-scale turbulence with magnetic field perturbations included, can match the power balance fluxes within temperature gradient error margins. Multiscale gyrokinetic simulations from two distinct codes produce no significant ETG heat flux, demonstrating that simple rules-of-thumb are insufficient criteria for its onset

Impact of W Events and Dust on JET-ILW Operation

ABSTRACT The occurrence of transient impurity events (TIE) leading to intense radiation spikes in JET plasma discharges has been studied since the installation of the ITER-like Wall (ILW). To generate the observed average increase in radiated power of 1.5MW, a 100µm-radius sphere of solid W dust would be required. The drop in plasma energy caused by W-TIEs is fully recovered in 90% of all cases, only 1% inducing a longer term loss in plasma energy which sometimes leads to the shutdown of plasma operation. TIEs are correlated with disruptions and with measurements of the dust mobilized by disruptions using the high resolution Thomson scattering (HRTS) diagnostic. The dust characteristics giving rise to TIEs have been studied using the dust transport code DTOKS and the 1D impurity transport code STRAH

Role of NBI fuelling in contributing to density peaking between the ICRH and NBI identity plasmas on JET

Density peaking has been studied between an ICRH and NBI identity plasma in JET. The comparison shows that 8 MW of NBI heating/fueling increases the density peaking by a factor of two, being R/L (n) = 0.45 for the ICRH pulse and R/L (n) = 0.93 for the NBI one averaged radially over rho (tor) = 0.4, 0.8. The dimensionless profiles of q, rho *, upsilon *, beta (n) and T (i)/T (e) approximate to 1 were matched within 5% difference except in the central part of the plasma (rho (tor) < 0.3). The difference in the curvature pinch (same q-profile) and thermo-pinch (T (i) = T (e)) between the ICRH and NBI discharges is virtually zero. Both the gyro-kinetic simulations and integrated modelling strongly support the experimental result where the NBI fuelling is the main contributor to the density peaking for this identity pair. It is to be noted here that the integrated modeling does not reproduce the measured electron density profiles, but approximately reproduces the difference in the density profiles between the ICRH and NBI discharge. Based on these modelling results and the analyses, the differences between the two pulses in impurities, fast ions (FIs), toroidal rotation and radiation do not cause any such changes in the background transport that would invalidate the experimental result where the NBI fuelling is the main contributor to the density peaking. This result of R/L (n) increasing by a factor of 2 per 8 MW of NBI power is valid for the ion temperature gradient dominated low power H-mode plasmas. However, some of the physics processes influencing particle transport, like rotation, turbulence and FI content scale with power, and therefore, the simple scaling on the role of the NBI fuelling in JET is not necessarily the same under higher power conditions or in larger devices

High Z neoclassical transport: Application and limitation of analytical formulae for modelling JET experimental parameters

Heavy impurities, such as tungsten (W), can exhibit strongly poloidally asymmetric density profiles in rotating or radio frequency heated plasmas. In the metallic environment of JET, the poloidal asymmetry of tungsten enhances its neoclassical transport up to an order of magnitude, so that neoclassical convection dominates over turbulent transport in the core. Accounting for asymmetries in neoclassical transport is hence necessary in the integrated modeling framework. The neoclassical drift kinetic code, NEO [E. Belli and J. Candy, Plasma Phys. Controlled Fusion P50, 095010 (2008)], includes the impact of poloidal asymmetries on W transport. However, the computational cost required to run NEO slows down significantly integrated modeling. A previous analytical formulation to describe heavy impurity neoclassical transport in the presence of poloidal asymmetries in specific collisional regimes [C. Angioni and P. Helander, Plasma Phys. Controlled Fusion 56, 124001 (2014)] is compared in this work to numerical results from NEO. Within the domain of validity of the formula, the factor for reducing the temperature screening due to poloidal asymmetries had to be empirically adjusted. After adjustment, the modified formula can reproduce NEO results outside of its definition domain, with some limitations: When main ions are in the banana regime, the formula reproduces NEO results whatever the collisionality regime of impurities, provided that the poloidal asymmetry is not too large. However, for very strong poloidal asymmetries, agreement requires impurities in the Pfirsch-Schlüter regime. Within the JETTO integrated transport code, the analytical formula combined with the poloidally symmetric neoclassical code NCLASS [W. A. Houlberg et al., Phys. Plasmas 4, 3230 (1997)] predicts the same tungsten profile as NEO in certain cases, while saving a factor of one thousand in computer time, which can be useful in scoping studies. The parametric dependencies of the temperature screening reduction due to poloidal asymmetries would need to be better characterised for this faster model to be extended to a more general applicability.</p