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].

Multi-Symplectic Variational Integrators for the Gross-Pitaevskii Equations in Bec

Multi-symplectic Variational Integrators for Nonlinear Schrödinger Equations with Variable Coefficients

An efficient and conservative compact finite difference scheme for the coupled Gross-Pitaevskii equations describing spin-1 Bose-Einstein condensate.

Nonstandard Finite Difference Variational Integrators for Multisymplectic Pdes

Multisymplectic Variational Integrators for Nonsmooth Lagrangian Continuum Mechanics

Nonstandard Finite Difference Variational Integrators for Nonlinear Schrödinger Equation with Variable Coefficients

High-order geometric integrators for the variational Gaussian approximation

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

Numerical solution of the Gross–Pitaevskii equation for Bose–Einstein condensation

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

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

i-SPin 2: An integrator for general spin-s Gross-Pitaevskii systems

A conservative splitting high-order finite difference method for coupled Gross–Pitaevskii equations in 2D

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

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

Stationary solutions and dynamical evolution of the nonlocal Gross–Pitaevskii equation in spin-2 Bose–Einstein condensates

Multisymplectic Variational Integrators and Space/time Symplecticity

Variational Symplectic Algorithm for Guiding Center Dynamics and Its Application in Tokamak Geometry

Finite volumes for the Gross-Pitaevskii equation

On the Exact Solutions for a Type of Gross-Pitaevskill Equation