Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
An Isogeometric Analysis approach for the study of the gyrokinetic quasi-neutrality equation
Journal of Computational Physics
| Hi-index | 31.45 |
Describing turbulence and microinstabilities in fusion devices is often modelled with the gyrokinetic equation. During the time evolution of the distribution function a field equation for the electrostatic potential needs to be solved. In the case of adiabatic electrons it contains a flux-surface-average term resulting in an integro-differential equation. Its numerical solution is time and memory intensive for three-dimensional configurations. Here a new algorithm is presented which only requires the numerical inversion of a simpler differential operator and a subsequent addition of a correction term. This new procedure is as fast as solving the equation without the surface average.