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

DOI10.1016/J.JCP.2019.109065zbMATH Open1453.65436arXiv1903.06885OpenAlexW2981783246WikidataQ126981457 ScholiaQ126981457MaRDI QIDQ2222984FDOQ2222984

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)

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.

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

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

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)

Cites Work

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

Cited In (4)

• Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics
• Poisson integrators based on splitting method for Poisson systems
• 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)