Direct Comparison of Gyrokinetic and Fluid Scrape-Off Layer Simulations

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.