Abstract:Symplectic integrators offer many advantages for the numerical solution of Hamiltonian differential equations, including bounded energy error and the preservation of invariant sets. Two of the central Hamiltonian systems encountered in plasma physics — the flow of magnetic field lines and the guiding center motion of magnetized charged particles — resist symplectic integration by conventional means because the dynamics are most naturally formulated in non-canonical coordinates, i.e., coordinates lacking the familiar (q, p) partitioning. Recent efforts made progress toward non-canonical symplectic integration of these systems by appealing to the variational integration framework; however, those integrators were multistep methods and later found to be numerically unstable due to parasitic mode instabilities. This work eliminates the multistep character and, therefore, the parasitic mode instabilities via an adaptation of the variational integration formalism that we deem “degenerate variational integration”. Both the magnetic field line and guiding center Lagrangians are degenerate in the sense that their resultant Euler-Lagrange equations are systems of first-order ODEs. We show that retaining the same degree of degeneracy when constructing a discrete Lagrangian yields one-step variational integrators preserving a non-canonical symplectic structure on the original Hamiltonian phase space. The advantages of the new algorithms are demonstrated via numerical examples, demonstrating superior stability compared to existing variational integrators for these systems and superior qualitative behavior compared to non-conservative algorithms.

Preservation of Quadratic Invariants by Semiexplicit Symplectic Integrators for Non-separable Hamiltonian Systems

Energy and Quadratic Invariants Preserving Methods for Hamiltonian Systems with Holonomic Constraints

Structure-preserving quantum algorithms for linear and nonlinear Hamiltonian systems

High-order linearly implicit exponential integrators conserving quadratic invariants with application to scalar auxiliary variable approach

Preservation of Stability Properties Near Fixed Points of Linear Hamiltonian Systems by Symplectic Integrators

A new type of energy-preserving integrators for quasi-bi-Hamiltonian systems

Conservative Multiquadric Quasi-Interpolation Method for Hamiltonian Wave Equations

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

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

Efficient Energy-Preserving Integrators for Oscillatory Hamiltonian Systems

High-order, large time-stepping integrators for scalar hyperbolic conservation laws

Novel semi-explicit symplectic schemes for nonseparable stochastic Hamiltonian systems

Quinpi: Integrating stiff hyperbolic systems with implicit high order finite volume schemes

Variational symplectic diagonally implicit Runge-Kutta methods for isospectral systems

Volume-preserving Exponential Integrators and Their Applications.

Volume-preserving Exponential Integrators

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

High Order Mimetic Symplectic Methods For Hamiltonian Systems

Volume Preservation by Runge-Kutta Methods

Linearly Implicit Global Energy Preserving Reduced-order Models for Cubic Hamiltonian Systems

On the preservation of second integrals by Runge-Kutta methods