Kosmas O, Leyendecker S (2018)
Publication Language: English
Publication Type: Journal article, Review article
Publication year: 2018
Pages Range: 1-21
In this work, we present a new derivation of higher order variational integration methods that exploit the phase lag properties for numerical integrations of systems with oscillatory solutions. More specifically, for the derivation of these integrators, the action integral along any curve segment is defined using a discrete Lagrangian that depends on the endpoints of the segment and on a number of intermediate points of interpolation. High order integrators are then obtained by writing down the discrete Lagrangian at any time interval as a weighted sum of the Lagrangians corresponding to a set of the chosen intermediate points. The respective positions and velocities are interpolated using trigonometric functions. The methods derived this way depend on a frequency, which in general needs to be accurately estimated. The new methods, which improve the phase lag characteristics by re-estimating the frequency at every time step, are presented and tested on the general N-body problem as numerical examples.
Kosmas, O., & Leyendecker, S. (2018). Variational integrators for orbital problems using frequency estimation. Advances in Computational Mathematics, 1-21. https://dx.doi.org/10.1007/s10444-018-9603-y
Kosmas, Odysseas, and Sigrid Leyendecker. "Variational integrators for orbital problems using frequency estimation." Advances in Computational Mathematics (2018): 1-21.