Abstract:A variational symplectic integrator for the guiding center motion of charged particles in general magnetic fields is developed for long time simulation studies of magnetized plasmas [1]. Instead of discretizing the differential equations of the guiding center motion, the action of the guiding center motion is discretized and minimized to obtain the iteration rules for advancing the dynamics. Standard integrators only guarantee the error to be small in each time-step. The errors at different time-steps often accumulate coherently, and result in a large error over a large number of time-steps. The variational symplectic integrator conserves exactly a discrete Lagrangian symplectic structure, and has better numerical properties over long integration time, compared with standard integrators, such as the standard and variable time-step 4th order Runge-Kutta methods. The symplectic integrator conserves the symplectic structure exactly, and guarantees that the energy error is bounded by a small number for all the time-steps. To construct the symplectic algorithm, it is necessary to adopt the variational approach because standard symplectic integrators are only valid for canonical Hamiltonian systems, and the guiding-center dynamics in general magnetic field does not possess a (global) canonical symplectic structure. Numerical examples with more than 25 million timesteps are given to demonstrate the superiority of the variational symplectic integrator. This significant improvement in long term simulation capability of gyrokinetics is a direct, otherwise-impossible result of the geometric formulation of the gyrokinetic theory using the modern language of differential geometry [2, 3].

Variational Integrators for the Gravitational N-Body Problem

High-order geometric integrators for the variational Gaussian approximation

Surrogate Lagrangians for Variational Integrators: High Order Convergence with Low Order Schemes

Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories

Neural Symplectic Integrator with Hamiltonian Inductive Bias for the Gravitational $N$-body Problem

Energy-preserving Variational Integrators for Forced Lagrangian Systems

Slimplectic Integrators: Variational Integrators for General Nonconservative Systems

Variational Symplectic Integrator for the Guiding Center Motion of Charged Particles for Long Time Simulations in General Magnetic Fields

A Tutorial on Variational Integrators

Variational Time Integration Approach for Smoothed Particle Hydrodynamics Simulation of Fluids

New Way to Construct High Order Hamiltonian Variational Integrators

Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler-Lagrange equations

High Order Numerical Integrators for Relativistic Charged Particle Tracking

Collocation in Conjunction with the Local Variational Iteration Method

Multi-Point Hermite Methods for the N-Body Problem

Variational Symplectic Integrator for Long-Time Simulations of the Guiding-Center Motion of Charged Particles in General Magnetic Fields.

A Symplectic Integrator for Hill's Equations

Multiple timestep reversible N -body integrators for close encounters in planetary systems

High-order conservative and accurately dissipative numerical integrators via auxiliary variables

Variational Discretizations for Hamiltonian Systems

A New Variational Integrator for Constrained Mechanical System Dynamics