Rapid Fourier space solution of linear partial differential equations in toroidal magnetic confinement geometries

Fluctuating quantities in magnetic confinement geometries often inherit a strong anisotropy along the field lines. One technique for describing these structures is the use of a certain set of Fourier components on the tori of nested flux surfaces. We describe an implementation of this approach for solving partial differential equations, like Poisson's equation, where a different set of Fourier components may be chosen on each surface according to the changing safety factor profile. Allowing the resolved components to change to follow the anisotropy significantly reduces the total number of degrees of freedom in the description. This can permit large gains in computational performance. We describe, in particular, how this approach can be applied to rapidly solve the gyrokinetic Poisson equation in a particle code, ORB5 (Jolliet et al., (2007) [5]), with a regular (non-field-aligned) mesh. (C) 2009 Published by Elsevier B.V.


Published in:
Computer Physics Communications, 181, 4, 715
Year:
2009
Publisher:
Elsevier
Keywords:
Laboratories:
SPC
CRPP




 Record created 2010-07-02, last modified 2018-01-28

External links:
Download fulltextURL
Download fulltextURL
Download fulltextn/a
Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)