Work

An electromagnetic model in SOLEDGE3X for edge plasma turbulence simulations in tokamaks

Year: 2025

Type: article

Abstract: Even small fluctuations in the magnetic field are known to impact edge plasma turbulence and transport properties in magnetic confinement fusion devices. Magnetic induction modifies the parallel electric field, and as such, it impacts the parallel current in Ohm's law. In contrast, magnetic fluctuations can induce geometrical and topological changes in B, leading to parallel transport across the equilibrium magnetic surfaces.This paper presents the new drift-reduced electromagnetic model implemented in SOLEDGE3X [Bufferand et al., Nuc. Fus. 2021]. Based on a domain decomposition, a specific numerical scheme is proposed using conservative second-order finite volumes associated with a semi-implicit time advancement. The coupling between the parallel current [[EQUATION]] and the parallel electromagnetic potential [[EQUATION]] in Ohm's and Ampere's laws is treated using a new toroidally and poloidally staggered grid. While adding [[EQUATION]] doubles the size of the vorticity operator to be inverted, numerical tests show that the inclusion of the finite mass of the electrons in the new model acts as an upper limit on the parallel diffusion coefficient, thus improving the conditioning of the matrix when the parallel resistivity [[EQUATION]] approaches zero.The changes in the algorithm are verified using the method of manufactured solutions. The implementation is first validated with respect to theoretical linear stability results, recovering the transition from Alfvén to thermal electron waves as the perpendicular wavenumber increased. Moreover, comparisons with available results in the literature show an overall agreement on a single blob propagation in a limited slab geometry. The first 3D simulations of plasma edge turbulence in TCV demonstrate the capability of the new solver to handle realistic tokamak configurations. The explicit-implicit time integration scheme enables one to compare electrostatic and electromagnetic effects using the same solver, with or without electron inertia and magnetic flutter. This ability opens the way to a better understanding of the impact of the various electrostatic and electromagnetic mechanisms on the transport and turbulence properties in realistic tokamak configurations.

Source: Journal of Computational Physics

DOI: 10.1016/j.jcp.2025.114052