Comparison of Gyrokinetic and Fluid Axisymmetric Simulations of a Prospective Spherical Tokamak Pilot Plant Scrape-Off Layer
POSTER
Abstract
Typically, fluid simulations are used for designing divertors for power handling and helium exhaust, and to study the dynamics of the Scrape-Off Layer (SOL) in tokamaks. These simulations are only valid if the SOL is highly collisional, but this assumption can break down in a reactor SOL, even in a high recycling regime. We find that kinetic effects can be significant in the upstream SOL in some reactor-relevant cases.
We conducted SOL simulations with parameters representative of a spherical tokamak pilot plant using a gyrokinetic code, Gkeyll, and fluid code, SOLPS. Simulations exhibit significant ion trapping in the upstream SOL when upstream parameters typical of a high recycling regime are used.
Of particular interest is confinement of impurities to the divertor region. 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 violated. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions (Kotschenreuther 2023).
We conducted SOL simulations with parameters representative of a spherical tokamak pilot plant using a gyrokinetic code, Gkeyll, and fluid code, SOLPS. Simulations exhibit significant ion trapping in the upstream SOL when upstream parameters typical of a high recycling regime are used.
Of particular interest is confinement of impurities to the divertor region. 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 violated. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions (Kotschenreuther 2023).
*Funding: CEDAC SciDAC (Center for Computational Evaluation and Design of Actuators for Core-edge Integration) Grant Number 26399299, Exofusion (partially funded by STEP), and other DOE sources
Presenters
-
Akash Shukla
- The University of Texas at Austin