Quasi-geometric integration of guiding-center orbits in piecewise linear toroidal fields
- 1 December 2020
- journal article
- research article
- Published by AIP Publishing in Physics of Plasmas
- Vol. 27 (12), 122508
- https://doi.org/10.1063/5.0022117
Abstract
A numerical integration method for guiding-center orbits of charged particles in toroidal fusion devices with three-dimensional field geometry is described. Here, high order interpolation of electromagnetic fields in space is replaced by a special linear interpolation, leading to locally linear Hamiltonian equations of motion with piecewise constant coefficients. This approach reduces computational effort and noise sensitivity, while the conservation of total energy, magnetic moment and phase space volume is retained. The underlying formulation treats motion in piecewise linear fields exactly and, thus, preserves the non-canonical symplectic form. The algorithm itself is only quasi-geometric due to a series expansion in the orbit parameter. For practical purposes, an expansion to the fourth order retains geometric properties down to computer accuracy in typical examples. When applied to collisionless guiding-center orbits in an axisymmetric tokamak and a realistic three-dimensional stellarator configuration, the method demonstrates stable long-term orbit dynamics conserving invariants. In Monte Carlo evaluation of transport coefficients, the computational efficiency of quasi-geometric integration is an order of magnitude higher than with a standard fourth order Runge-Kutta integrator.Funding Information
- OeAD-GmbH (UA 04/2017)
- Helmholtz Association (ZT-I-0010)
- H2020 Euratom (633053)
This publication has 42 references indexed in Scilit:
- VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particlesComputer Physics Communications, 2014
- Neoclassical Toroidal Viscosity Calculations in Tokamaks Using aMonte Carlo Simulation and Their VerificationsPhysical Review Letters, 2011
- Variance reduction in computations of neoclassical transport in stellarators using a δf methodPhysics of Plasmas, 2008
- Study of Neoclassical Transport in LHD Plasmas by Applying the DCOM/NNW Neoclassical Transport DatabasePlasma and Fusion Research, 2008
- A global collisionless PIC code in magnetic coordinatesComputer Physics Communications, 2007
- VENUS+δf: A Bootstrap Current Calculation Module for 3-D ConfigurationsFusion Science and Technology, 2006
- Particle Simulation of the Neoclassical PlasmasJournal of Computational Physics, 2001
- Monte Carlo estimation of neoclassical transport for the TJ-II stellaratorPhysics of Plasmas, 2001
- Monte Carlo computations of neoclassical transportPhysics of Fluids, 1988
- Monte Carlo evaluation of transport coefficientsPhysics of Fluids, 1981