D. R. Hatch

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the SOL is highly collisional, an assumption that is valid in many present day experiments but is questionable in the high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the scrape off layer (SOL) for parameters and geometry representative of the Spherical Tokamak for Energy Production (STEP) fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional (2D) axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS's form of fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in reactor-like regimes the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

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.