A fast solver for the gyrokinetic field equation with adiabatic electrons

JOURNAL OF COMPUTATIONAL PHYSICS

卷:231

Article

Borchardt, M.¹  Kleiber, R.¹  Hackbusch, W.² 

[1] EURATOM, Max Planck Inst Plasmaphys, D-17491 Greifswald, Germany

[2] Max Planck Inst Math Nat Wissensch, D-04103 Leipzig, Germany

关键词: Fast solver;    Gyrokinetic field equation;   

DOI  :  10.1016/j.jcp.2012.06.001

来源: Elsevier

【 摘 要 】

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. (C) 2012 Elsevier Inc. All rights reserved.

【 授权许可】

Free