Global electromagnetic gyrokinetic simulations of internal transport barriers in reversed-shear tokamaks

TL;DR

This study uses electromagnetic gyrokinetic simulations to explore the formation of internal transport barriers in reversed-shear tokamaks, highlighting the importance of electron dynamics and specific q-profile conditions.

Contribution

It demonstrates the critical role of electron dynamics and q-profile rationality in ITB formation through comprehensive flux-driven gyrokinetic simulations.

Findings

Electron dynamics induce zonal flows and profile corrugation.

ITB forms at specific q-profile minima close to rational values.

Transport barrier development depends on q-profile and turbulence interactions.

Abstract

This work aims at improving our understanding of the conditions enabling the development of an Internal transport barriers (ITB), using a more comprehensive physical model, including low-$\beta$ electromagnetic flux-driven simulations. Our key findings are that electron dynamics is crucial for ITB formation even in an ITG scenario and that having $q_{\text{min}}$ close to a lowest order rational value (2 in our simulations) to allow for eddies self-interaction is a necessary ingredient. Electron dynamics has two critical effects. First, it leads to a structure formation characterized by strong zonal flows shearing rate, reduction of turbulence and profile corrugation. Second, it leads to zonal current sheets that result in a broadening of the minimum-q region, qualitatively consistent with the flux-tube simulations of Vol\v{c}okas et al. [1]. Flux-driven simulations performed with $q_{\text{min}}=2$ reveal the development of the transport barrier in the ion channel, forming at inner and outer radial positions with respect to the $q_{\text{min}}$ position. The ITB formation in flux-driven setup is not recovered if $q_{\text{min}} = 2.03$. Additionally, a simulation at higher $\rho^*$ indicates that the extent of the flattened region of the q-profile due to turbulent self-interaction does not change proportionally to $\rho^*$ or to $\rho_i$, but somewhere in between. On the other hand, the input power required to achieve similar on-axis temperatures appears to exhibit almost GyroBohm scaling (for the two considered $\rho^*$ values). Furthermore, considering an initial q-profile with $q_{\text{min}} = 2.01$, flux-driven simulations show that partial self-interaction can evolve to complete self-interaction. This occurs due to turbulent-driven zonal currents that lower and flatten the q-profile down to $q_{\text{min}} = 2.0$, in line with what is reported in Vol\v{c}okas et al.[1].