Analysis of equilibrium and turbulent fluxes across the separatrix in a gyrokinetic simulation

TL;DR

This paper investigates the particle and heat fluxes across the separatrix in a realistic tokamak simulation, revealing turbulence's dominant role in SOL width and the complex interplay of drifts, flows, and potential structures.

Contribution

It provides a detailed gyrokinetic analysis of turbulence and fluxes at the separatrix, incorporating realistic geometry and effects often neglected in simplified models.

Findings

Turbulent electron flux primarily determines SOL width.

Edge turbulence is regulated by $ abla B$-drifts and ion X-point losses.

Shear flows suppress turbulence and influence potential structures.

Abstract

The SOL width is a parameter of paramount importance in modern tokamaks as it controls the power density deposited at the divertor plates, critical for plasma-facing material survivability. An understanding of the parameters controlling it has consequently long been sought (Connor et al. 1999 NF 39 2). Prior to Chang et al.(2017 NF 57 11), studies of the tokamak edge have been mostly confined to reduced fluid models and simplified geometries, leaving out important pieces of physics. Here, we analyze the results of a DIII-D simulation performed with the full-f gyrokinetic code XGC1 which includes both turbulence and neoclassical effects in realistic divertor geometry. More specifically, we calculate the particle and heat ExB fluxes along the separatrix, discriminating between equilibrium and turbulent contributions. We find that the density SOL width is impacted almost exclusively by the turbulent electron flux. In this simulation, the level of edge turbulence is regulated by a mechanism we are only beginning to understand: $\nabla B$-drifts and ion X-point losses at the top and bottom of the machine, along with ion banana orbits at the low field side (LFS), result in a complex poloidal potential structure at the separatrix which is the cause of the ExB drift pattern that we observe. Turbulence is being suppressed by the shear flows that this potential generates. At the same time, turbulence, along with increased edge collisionality and electron inertia, can influence the shape of the potential structure by making the electrons non-adiabatic. Moreover, being the only means through which the electrons can lose confinement, it needs to be in a balance with the original direct ion orbit losses to maintain charge neutrality.