Construction and analysis of guiding center distributions for tokamak plasmas with ambient radial electric field

TL;DR

This paper extends the VisualStart code to include a static radial electric field, enabling detailed construction of guiding center distributions in tokamak plasmas, with applications to stability analysis and resonance phenomena.

Contribution

It introduces a method to incorporate ambient radial electric fields into guiding center orbit calculations, improving accuracy and efficiency in tokamak plasma modeling.

Findings

Midplane-based coordinates are effective for sampling GC orbit space with $E_r$.

Radial electric field modulates parallel acceleration and affects trapped-passing boundaries.

Electric frequency shifts depend on the reference point, impacting resonance analysis.

Abstract

The contribution of a time-independent toroidally-symmetric radial electric field $E_r$ is implemented in VisualStart [Comp. Phys. Comm. 275 (2022) 108305; arXiv:2111.08224], a code whose purposes include the construction of guiding center (GC) drift orbit databases for the study of plasma instabilities in tokamaks. $E_r$ is important for the thermal part of the velocity distribution and for fast particle resonances in the kHz frequency range. KSTAR, JT-60U and ITER tokamak cases are used as working examples to test our methods and discuss practical issues connected with $E_r$. Two points are worth noting: First, the GC orbit space is sampled in the magnetic midplane as before, and we find that in the presence of $E_r$, midplane-based coordinates are not only equivalent but superior to conventional constants of motion, allowing to attain high numerical accuracy and efficiency with a relatively simple mesh. Second, the periodic parallel acceleration and deceleration of GCs via the mirror force is modulated by $E_r$. Although this parallel electric acceleration averages to zero during a poloidal transit (or bounce) period, it has important consequences, one being the known shift of the trapped-passing boundary. Another consequence is that electric frequency shifts depend on the chosen reference point, so that some care is required when evaluating the $E_r$-dependence of transit frequencies for resonance analyses.