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

DOI10.1016/J.JCP.2019.109065MaRDI QIDQ2222984zbMATHOpenAlexWikidataFDO

Authors Christopher G. Albert, Sergei V. Kasilov, Winfried Kernbichler

Publication date 28 January 2021

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1903.06885

zbMATH Keywords

plasmanumerical integrationHamiltonian systemssymplectic integrationmagnetic confinementguiding-center dynamics

Mathematics Subject Classification ID

Numerical methods for Hamiltonian systems including symplectic integrators (65P10) Applications to the sciences (65Z05)

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 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., Phys. Plasmas 21, 32504 (2014)]. Here a generalization of this approach with emphasis on semi-implicit partitioned schemes is described together with methods to enhance their 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.

Recommendations

• Symplectic integration of magnetic systems
• Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
• Hamiltonian theory of guiding-center motion
• scientific article; zbMATH DE number 4215292
• Symplectic integration of Hamiltonian systems

Cites work

• scientific article; zbMATH DE number 1745051 (Why is no real title available?)
• scientific article; zbMATH DE number 1762480 (Why is no real title available?)
• A family of embedded Runge-Kutta formulae
• An explicit symplectic integration scheme for plasma simulations
• Discrete mechanics and variational integrators
• Geometric integration for particle accelerators
• Gyrokinetics from variational averaging: existence and error bounds
• Hamiltonian guiding center drift orbit calculation for plasmas of arbitrary cross section
• Hamiltonian theory of guiding-center motion
• On the Numerical Integration of Ordinary Differential Equations by Symmetric Composition Methods
• Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria

Cited in

(4)

• Poisson integrators based on splitting method for Poisson systems
• Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics
• Symplectic Gaussian process regression of maps in Hamiltonian systems
• VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

This page was built for publication: Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222984)