Electrostatic gyrokinetic simulations in Wendelstein 7-X geometry: benchmark between the codes stella and GENE

TL;DR

This paper presents a detailed benchmark comparison of two gyrokinetic simulation codes, stella and GENE, in Wendelstein 7-X geometry, focusing on electrostatic turbulence, stability, and heat flux calculations to validate and improve stellarator modeling.

Contribution

It provides the first comprehensive benchmark between stella and GENE in W7-X geometry, aiding code validation and development for stellarator turbulence studies.

Findings

Consistent stability analysis results between codes.

Agreement in linear zonal flow response calculations.

Comparable ITG-driven heat flux predictions.

Abstract

The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attemps to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode (TEM)} stability analyses, computation of linear zonal flow responses and calculation of ITG-driven heat fluxes.