Projected Variational Integrators for Degenerate Lagrangian Systems

Abstract

We propose and compare several projection methods applied to variational integrators for degenerate Lagrangian systems, whose Lagrangian is of the form $L=\vartheta (q)\cdot \dot{q}-H(q)$ and thus linear in velocities. While previous methods for such systems only work reliably in the case of $\vartheta$ being a linear function of $q$, our methods are long-time stable also for systems where $\vartheta$ is a nonlinear function of $q$. We analyse the properties of the resulting algorithms, in particular with respect to the conservation of energy, momentum maps and symplecticity. In numerical experiments, we verify the favourable properties of the projected integrators and demonstrate their excellent long-time fidelity. In particular, we consider a two-dimensional Lotka-Volterra system, planar point vortices with position-dependent circulation and guiding centre dynamics.

67 Pages

Cited by in corpus (8)

• Symplectic integration with non-canonical quadrature for guiding-center orbits in magnetic confinement devices
• Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories
• Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics
• Field theory and structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence
• Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds
• Learning of discrete models of variational PDEs from data
• Time-symmetry, symplecticity and stability of Euler-Maclaurin and Lanczos-Dyche integration
• Charge-conserving, variational particle-in-cell method for the drift-kinetic Vlasov-Maxwell system