Optimising TGLF for a Q=10 Burning Spherical Tokamak

TGLF transport model predictions have been assessed in the vicinity of a theoretical high β burning plasma spherical tokamak at Q=10. Linear micro-stability calculations from TGLF have been compared on a surface at mid-radius with the gyrokinetic code GS2. Differences between TGLF and GS2 spectra can be characterised by the RMS difference in growth rates, σγ. We find considerable improvement in the quality of TGLF growth rate spectrum can be achieved by increasing the number of parallel basis functions and by tuning the TGLF parameter used in the model for trapped particles, θtrap

Gyrokinetic analysis and simulation of pedestals, to identify the culprits for energy losses using fingerprints

Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experiments, find that each mode type produces characteristic ratios of transport in the various channels: density, heat and impurities. This, together with experimental observations of transport in some channel, or, of the relative size of the driving sources of channels, can identify or determine the dominant modes causing energy transport. In multiple ELMy H-mode cases that are examined, these fingerprints indicate that MHD-like modes are apparently not the dominant agent of energy transport; rather, this role is played by Micro-Tearing Modes (MTM) and Electron Temperature Gradient (ETG) modes, and in addition, possibly Ion Temperature Gradient (ITG)/Trapped Electron Modes (ITG/TEM) on JET. MHD-like modes may dominate the electron particle losses. Fluctuation frequency can also be an important means of identification, and is often closely related to the transport fingerprint. The analytical arguments unify and explain previously disparate experimental observations on multiple devices, including DIII-D, JET and ASDEX-U, and detailed simulations of two DIII-D ETBs also demonstrate and corroborate this

Self-consistent simulation of plasma scenarios for ITER using a combination of 1.5D transport codes and free-boundary equilibrium codes

Self-consistent transport simulation of ITER scenarios is a very important tool for the exploration of the operational space and for scenario optimisation. It also provides an assessment of the compatibility of developed scenarios (which include fast transient events) with machine constraints, in particular with the poloidal field (PF) coil system, heating and current drive (H&CD), fuelling and particle and energy exhaust systems. This paper discusses results of predictive modelling of all reference ITER scenarios and variants using two suite of linked transport and equilibrium codes. The first suite consisting of the 1.5D core/2D SOL code JINTRAC [1] and the free boundary equilibrium evolution code CREATE-NL [2,3], was mainly used to simulate the inductive D-T reference Scenario-2 with fusion gain Q=10 and its variants in H, D and He (including ITER scenarios with reduced current and toroidal field). The second suite of codes was used mainly for the modelling of hybrid and steady state ITER scenarios. It combines the 1.5D core transport code CRONOS [4] and the free boundary equilibrium evolution code DINA-CH [5].Comment: 23 pages, 18 figure

Modelling of 3D fields due to ferritic inserts and test blanket modules in toroidal geometry at ITER

Computations in toroidal geometry are systematically performed for the plasma response to 3D magnetic perturbations produced by ferritic inserts (FIs) and test blanket modules (TBMs) for four ITER plasma scenarios: the 15 MA baseline, the 12.5 MA hybrid, the 9 MA steady state, and the 7.5 MA half-field helium plasma. Due to the broad toroidal spectrum of the FI and TBM fields, the plasma response for all the n = 1-6 field components are computed and compared. The plasma response is found to be weak for the high-n (n > 4) components. The response is not globally sensitive to the toroidal plasma flow speed, as long as the latter is not reduced by an order of magnitude. This is essentially due to the strong screening effect occurring at a finite flow, as predicted for ITER plasmas. The ITER error field correction coils (EFCC) are used to compensate the n = 1 field errors produced by FIs and TBMs for the baseline scenario for the purpose of avoiding mode locking. It is found that the middle row of the EFCC, with a suitable toroidal phase for the coil current, can provide the best correction of these field errors, according to various optimisation criteria. On the other hand, even without correction, it is predicted that these n = 1 field errors will not cause substantial flow damping for the 15 MA baseline scenario