Transport Characteristics and Modelling of ST40 Hot Ion Plasmas

TL;DR

This study investigates turbulent transport in ST40 hot ion plasmas, comparing models and performing predictive simulations to understand confinement and transport phenomena in spherical tokamaks.

Contribution

It provides a comprehensive comparison of gyro-kinetic and gyro-fluid models and introduces a novel reduced SOL model for predictive plasma transport simulations.

Findings

Core ion transport is near neoclassical levels due to turbulence suppression.

Edge transport remains dominated by anomalous processes.

Predictive simulations show good agreement with experimental data.

Abstract

In this paper, the turbulent transport properties of ST40 hot ion plasmas are examined and fully predictive time evolving modelling of a hot ion plasma pulse was performed. Understanding turbulent transport on spherical tokamaks (STs) is challenging due to their unique geometry characteristics. ST40 hot ion plasmas are typically unstable to ion scale Trapped Electron Modes (TEMs) and Ubiquitous Modes (UMs), driven from the kinetic response of trapped particles and passing ions, and electron scale Electron Temperature Gradient Modes (ETGs) at the edge of the plasma. A comparison between the linear unstable modes of the gyro-kinetic code GS2 and the gyro-fluid code TGLF showed that both models agree to a satisfactory level. However, some discrepancy was observed at the core of the plasma where a large fraction of beams ions exists, and electromagnetic effects are potentially important. Turbulent fluxes were also observed to be somewhat overpredicted with TGLF. The core heat ion transport is observed to be close to neoclassical levels due to turbulence suppression from high rotation and fast ion stabilisation, while the edge region is dominated by anomalous transport in both ions and electrons. As a result, enhanced energy confinement is observed in those plasmas driven by the reduced turbulent core region and the confined beam ions. Fully predictive simulations using the ASTRA transport solver coupled with SPIDER, NUBEAM, NCLASS and TGLF together with a novel reduced scrape of layer (SOL) model for the simulation of the last closed flux surface (LCFS) boundary conditions was attempted. Agreement in global quantities but also kinetic profiles between the predictive and interpretative modelling as well as experimental measurements was observed.