Algorithms and optimizations for global non-linear hybrid fluid-kinetic finite element stellarator simulations

TL;DR

This paper introduces a novel, high-fidelity computational framework for simulating stellarator plasmas using hybrid fluid-kinetic models, overcoming complex 3D coupling challenges with efficient algorithms.

Contribution

It presents a globally coupled projection scheme within the JOREK finite element framework, enabling accurate and scalable simulation of stellarator plasma dynamics.

Findings

Achieved spectral convergence in realistic stellarator geometries.

Demonstrated computational efficiency with FFT and R-Tree spatial indexing.

Validated the framework's accuracy and scalability.

Abstract

Predictive modeling of stellarator plasmas is crucial for advancing nuclear fusion energy, yet it faces unique computational difficulties. One of the main challenges is accurately simulating the dynamics of specific particle species that are not well captured by fluid models, which necessitates the use of hybrid fluid-kinetic models. The non-axisymmetric geometry of stellarators fundamentally couples the toroidal Fourier modes, in contrast to what happens in tokamaks, requiring different numerical and computational treatment. This work presents a novel, globally coupled projection scheme inside the JOREK finite element framework. The approach ensures a self-consistent and physically accurate transfer of kinetic markers to the fluid grid, effectively handling the complex 3D mesh by constructing and solving a unified linear system that encompasses all toroidal harmonics simultaneously. To manage the computational complexity of this coupling, the construction of the system's matrix is significantly accelerated using the Fast Fourier Transform (FFT). The efficient localization of millions of particles is made possible by implementing a 3D R-Tree spatial index, which supports this projection and ensures computational tractability at scale. On realistic Wendelstein 7-X stellarator geometries, the fidelity of the framework is rigorously shown. In sharp contrast to the uncoupled approaches' poor performance, quantitative convergence tests verify that the coupled scheme attains the theoretically anticipated spectral convergence. This study offers a crucial capability for the predictive analysis and optimization of next-generation stellarator designs by developing a validated, high-fidelity computational tool.