Numerical and diagnostic extensions to the gyrokinetic δf code stella
Abstract
Understanding turbulence in magnetised plasmas requires robust numerical tools capable of capturing complex free-energy transfer processes across scales. We demonstrate the utility of free energy as a diagnostic for both physical turbulence characteristics and numerical stability in gyrokinetic simulations with the stella code. Applied to the W7-X high-mirror configuration, free-energy diagnostics reveal characteristic power-law behaviour in the turbulent cascade, act as sensitive indicators of numerical instabilities undetected by conventional diagnostics, and expose localised structures in real space not apparent in the electrostatic potential alone. To support these analyses, we present a systematic comparison of explicit time-integration schemes, showing that the optimal choice depends critically on the regime: SSP22 achieves the lowest computational cost in linear ITG simulations, while SSP33 is most efficient in the nonlinear saturated state, with nonlinear simulations requiring two to three orders of magnitude more computational cost than linear cases. These results establish free-energy diagnostics as a valuable tool for interpreting plasma turbulence and provide practical guidance for scheme selection in stella simulations.
Related
- Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE
- ORB5X: a performance-portable global electromagnetic gyrokinetic PIC code in C++/Kokkos built using Agentic AI
- A Hybrid Simulation Code for Hall Thrusters and its Sensitivity to Numerical Parameters
- Design of the IMAS interface data structure for flux-tube gyrokinetic simulations
- Gyrokinetic Theory of Linear Gravitational Flute Interchanges with Flow Shear Stabilization
BibTeX
@misc{stroteich2026numerical,
title = {{Numerical and diagnostic extensions to the gyrokinetic $δf$ code stella}},
author = {Sven Stroteich and Georgia Acton and Michael Barnes and Pavlos Xanthopoulos and Ralf Schneider},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.02958},
eprint = {2610.02958},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.02958},
}
RIS
TY - UNPB
TI - Numerical and diagnostic extensions to the gyrokinetic δf code stella
AU - Sven Stroteich
AU - Georgia Acton
AU - Michael Barnes
AU - Pavlos Xanthopoulos
AU - Ralf Schneider
PY - 2026
DA - 2026/10/02
PB - arXiv
DO - 10.48550/arxiv.2610.02958
UR - https://arxiv.org/abs/2610.02958
AB - Understanding turbulence in magnetised plasmas requires robust numerical tools capable of capturing complex free-energy transfer processes across scales. We demonstrate the utility of free energy as a diagnostic for both physical turbulence characteristics and numerical stability in gyrokinetic simulations with the stella code. Applied to the W7-X high-mirror configuration, free-energy diagnostics reveal characteristic power-law behaviour in the turbulent cascade, act as sensitive indicators of numerical instabilities undetected by conventional diagnostics, and expose localised structures in real space not apparent in the electrostatic potential alone. To support these analyses, we present a systematic comparison of explicit time-integration schemes, showing that the optimal choice depends critically on the regime: SSP22 achieves the lowest computational cost in linear ITG simulations, while SSP33 is most efficient in the nonlinear saturated state, with nonlinear simulations requiring two to three orders of magnitude more computational cost than linear cases. These results establish free-energy diagnostics as a valuable tool for interpreting plasma turbulence and provide practical guidance for scheme selection in stella simulations.
ER -