Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices

Publikation: Beitrag in einer FachzeitschriftArtikelForschungBegutachtung

Abstract

We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The technique combines time-stepping in canonical coordinates with quadrature in non-canonical coordinates and is applicable in systems where a global transformation to canonical coordinates is known explicitly but its inverse is not. A fully implicit class of symplectic Runge-Kutta schemes has recently been introduced and applied to guiding-center motion by Zhang et al. (2014) [9]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance performance, in particular avoiding evaluation of non-canonical variables at full time steps. For application in toroidal plasma confinement configurations with nested magnetic flux surfaces a global canonicalization of coordinates for the guiding-center Lagrangian by a spatial transform is presented that allows for pre-computation of the required map in a parallel algorithm in the case of time-independent magnetic field geometry. Guiding-center orbits are studied in stationary magnetic equilibrium fields of an axisymmetric tokamak and a realistic three-dimensional stellarator configuration. Superior long-term properties of symplectic methods are demonstrated in comparison to a conventional adaptive Runge-Kutta scheme. Finally statistics of fast fusion alpha particle losses over their slowing-down time are computed in the stellarator field on a representative sample, reaching a speed-up of the symplectic Euler scheme by more than a factor three compared to usual Runge-Kutta schemes while keeping the same statistical accuracy and linear scaling with the number of computing threads.
Originalspracheenglisch
Aufsatznummer109065
Seitenumfang23
FachzeitschriftJournal of computational physics
Jahrgang403
DOIs
PublikationsstatusVeröffentlicht - Feb 2020

Fingerprint

Symplectic Integration
Plasma confinement
quadratures
Quadrature
Orbits
Orbit
Runge-Kutta Schemes
orbits
Hamiltonians
Alpha particles
Magnetic flux
Parallel algorithms
stellarators
Fusion reactions
Statistics
Magnetic fields
Plasmas
Plasma
Geometry
Symplectic Methods

Fields of Expertise

  • Information, Communication & Computing

Kooperationen

  • NAWI Graz

Dies zitieren

@article{1c88cbc1117f4cb8a0e138ae6679747a,
title = "Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices",
abstract = "We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The technique combines time-stepping in canonical coordinates with quadrature in non-canonical coordinates and is applicable in systems where a global transformation to canonical coordinates is known explicitly but its inverse is not. A fully implicit class of symplectic Runge-Kutta schemes has recently been introduced and applied to guiding-center motion by Zhang et al. (2014) [9]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance performance, in particular avoiding evaluation of non-canonical variables at full time steps. For application in toroidal plasma confinement configurations with nested magnetic flux surfaces a global canonicalization of coordinates for the guiding-center Lagrangian by a spatial transform is presented that allows for pre-computation of the required map in a parallel algorithm in the case of time-independent magnetic field geometry. Guiding-center orbits are studied in stationary magnetic equilibrium fields of an axisymmetric tokamak and a realistic three-dimensional stellarator configuration. Superior long-term properties of symplectic methods are demonstrated in comparison to a conventional adaptive Runge-Kutta scheme. Finally statistics of fast fusion alpha particle losses over their slowing-down time are computed in the stellarator field on a representative sample, reaching a speed-up of the symplectic Euler scheme by more than a factor three compared to usual Runge-Kutta schemes while keeping the same statistical accuracy and linear scaling with the number of computing threads.",
author = "Albert, {Christopher G.} and Kasilov, {Sergei V.} and Winfried Kernbichler",
year = "2020",
month = "2",
doi = "10.1016/j.jcp.2019.109065",
language = "English",
volume = "403",
journal = "Journal of computational physics",
issn = "0021-9991",
publisher = "Academic Press",

}

TY - JOUR

T1 - Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices

AU - Albert, Christopher G.

AU - Kasilov, Sergei V.

AU - Kernbichler, Winfried

PY - 2020/2

Y1 - 2020/2

N2 - We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The technique combines time-stepping in canonical coordinates with quadrature in non-canonical coordinates and is applicable in systems where a global transformation to canonical coordinates is known explicitly but its inverse is not. A fully implicit class of symplectic Runge-Kutta schemes has recently been introduced and applied to guiding-center motion by Zhang et al. (2014) [9]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance performance, in particular avoiding evaluation of non-canonical variables at full time steps. For application in toroidal plasma confinement configurations with nested magnetic flux surfaces a global canonicalization of coordinates for the guiding-center Lagrangian by a spatial transform is presented that allows for pre-computation of the required map in a parallel algorithm in the case of time-independent magnetic field geometry. Guiding-center orbits are studied in stationary magnetic equilibrium fields of an axisymmetric tokamak and a realistic three-dimensional stellarator configuration. Superior long-term properties of symplectic methods are demonstrated in comparison to a conventional adaptive Runge-Kutta scheme. Finally statistics of fast fusion alpha particle losses over their slowing-down time are computed in the stellarator field on a representative sample, reaching a speed-up of the symplectic Euler scheme by more than a factor three compared to usual Runge-Kutta schemes while keeping the same statistical accuracy and linear scaling with the number of computing threads.

AB - We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The technique combines time-stepping in canonical coordinates with quadrature in non-canonical coordinates and is applicable in systems where a global transformation to canonical coordinates is known explicitly but its inverse is not. A fully implicit class of symplectic Runge-Kutta schemes has recently been introduced and applied to guiding-center motion by Zhang et al. (2014) [9]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance performance, in particular avoiding evaluation of non-canonical variables at full time steps. For application in toroidal plasma confinement configurations with nested magnetic flux surfaces a global canonicalization of coordinates for the guiding-center Lagrangian by a spatial transform is presented that allows for pre-computation of the required map in a parallel algorithm in the case of time-independent magnetic field geometry. Guiding-center orbits are studied in stationary magnetic equilibrium fields of an axisymmetric tokamak and a realistic three-dimensional stellarator configuration. Superior long-term properties of symplectic methods are demonstrated in comparison to a conventional adaptive Runge-Kutta scheme. Finally statistics of fast fusion alpha particle losses over their slowing-down time are computed in the stellarator field on a representative sample, reaching a speed-up of the symplectic Euler scheme by more than a factor three compared to usual Runge-Kutta schemes while keeping the same statistical accuracy and linear scaling with the number of computing threads.

UR - https://doi.org/10.1016/j.jcp.2019.109065

U2 - 10.1016/j.jcp.2019.109065

DO - 10.1016/j.jcp.2019.109065

M3 - Article

VL - 403

JO - Journal of computational physics

JF - Journal of computational physics

SN - 0021-9991

M1 - 109065

ER -