VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

D. Pfefferlé; W. A. Cooper; J. P. Graves; C. Misev

doi:10.1016/j.cpc.2014.08.007

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

D. Pfefferlé, W. A. Cooper, J. P. Graves, C. Misev

Research output: Contribution to journal › Article › peer-review

33 Citations (Scopus)

Abstract

Curvilinear guiding-centre drift and full-orbit equations of motion are presented as implemented in the VENUS-LEVIS code. A dedicated interpolation scheme based on Fourier reconstruction in the toroidal and poloidal directions and cubic spline in the radial direction of flux coordinate systems is detailed. This interpolation method exactly preserves the order of the RK4 integrating scheme which is crucial for the investigation of fast particle trajectories in 3D magnetic structures such as helical saturated tokamak plasma states, stellarator geometry and resonant magnetic perturbations (RMP). The initialisation of particles with respect to the guiding-centre is discussed. Two approaches to implement RMPs in orbit simulations are presented, one where the vacuum field is added to the 2D equilibrium, creating islands and stochastic regions, the other considering 3D nested flux-surfaces equilibrium including the RMPs.

Original language	English
Pages (from-to)	3127-3140
Number of pages	14
Journal	Computer Physics Communications
Volume	185
Issue number	12
DOIs	https://doi.org/10.1016/j.cpc.2014.08.007
Publication status	Published - 1 Dec 2014
Externally published	Yes

Access to Document

10.1016/j.cpc.2014.08.007

Cite this

@article{954a814f762e45009b4782678579ad5d,

title = "VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles",

abstract = "Curvilinear guiding-centre drift and full-orbit equations of motion are presented as implemented in the VENUS-LEVIS code. A dedicated interpolation scheme based on Fourier reconstruction in the toroidal and poloidal directions and cubic spline in the radial direction of flux coordinate systems is detailed. This interpolation method exactly preserves the order of the RK4 integrating scheme which is crucial for the investigation of fast particle trajectories in 3D magnetic structures such as helical saturated tokamak plasma states, stellarator geometry and resonant magnetic perturbations (RMP). The initialisation of particles with respect to the guiding-centre is discussed. Two approaches to implement RMPs in orbit simulations are presented, one where the vacuum field is added to the 2D equilibrium, creating islands and stochastic regions, the other considering 3D nested flux-surfaces equilibrium including the RMPs.",

keywords = "Curvilinear coordinates, Fast ion transport, Magnetic equilibrium representation, Orbit simulations, Tokamak plasmas",

author = "D. Pfefferl{\'e} and Cooper, {W. A.} and Graves, {J. P.} and C. Misev",

year = "2014",

month = dec,

day = "1",

doi = "10.1016/j.cpc.2014.08.007",

language = "English",

volume = "185",

pages = "3127--3140",

journal = "Computer Physics Communications",

issn = "0010-4655",

publisher = "Elsevier",

number = "12",

}

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles. / Pfefferlé, D.; Cooper, W. A.; Graves, J. P. et al.
In: Computer Physics Communications, Vol. 185, No. 12, 01.12.2014, p. 3127-3140.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

AU - Pfefferlé, D.

AU - Cooper, W. A.

AU - Graves, J. P.

AU - Misev, C.

PY - 2014/12/1

Y1 - 2014/12/1

N2 - Curvilinear guiding-centre drift and full-orbit equations of motion are presented as implemented in the VENUS-LEVIS code. A dedicated interpolation scheme based on Fourier reconstruction in the toroidal and poloidal directions and cubic spline in the radial direction of flux coordinate systems is detailed. This interpolation method exactly preserves the order of the RK4 integrating scheme which is crucial for the investigation of fast particle trajectories in 3D magnetic structures such as helical saturated tokamak plasma states, stellarator geometry and resonant magnetic perturbations (RMP). The initialisation of particles with respect to the guiding-centre is discussed. Two approaches to implement RMPs in orbit simulations are presented, one where the vacuum field is added to the 2D equilibrium, creating islands and stochastic regions, the other considering 3D nested flux-surfaces equilibrium including the RMPs.

AB - Curvilinear guiding-centre drift and full-orbit equations of motion are presented as implemented in the VENUS-LEVIS code. A dedicated interpolation scheme based on Fourier reconstruction in the toroidal and poloidal directions and cubic spline in the radial direction of flux coordinate systems is detailed. This interpolation method exactly preserves the order of the RK4 integrating scheme which is crucial for the investigation of fast particle trajectories in 3D magnetic structures such as helical saturated tokamak plasma states, stellarator geometry and resonant magnetic perturbations (RMP). The initialisation of particles with respect to the guiding-centre is discussed. Two approaches to implement RMPs in orbit simulations are presented, one where the vacuum field is added to the 2D equilibrium, creating islands and stochastic regions, the other considering 3D nested flux-surfaces equilibrium including the RMPs.

KW - Curvilinear coordinates

KW - Fast ion transport

KW - Magnetic equilibrium representation

KW - Orbit simulations

KW - Tokamak plasmas

UR - http://www.scopus.com/inward/record.url?scp=85027939426&partnerID=8YFLogxK

U2 - 10.1016/j.cpc.2014.08.007

DO - 10.1016/j.cpc.2014.08.007

M3 - Article

AN - SCOPUS:85027939426

SN - 0010-4655

VL - 185

SP - 3127

EP - 3140

JO - Computer Physics Communications

JF - Computer Physics Communications

IS - 12

ER -

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

Abstract

Access to Document

Other files and links

Fingerprint

Identification of an Optimized Heating and Fast Ion Generation Scheme for the Wendelstein 7-X Stellarator

The VENUS-LEVIS orbit solver

Cite this