On frequency estimations in phase fitted variational integrators for the general N-body problem

In the present work, the advantages of high order variational integrator methods are combined with phase lag properties for the numerical integration of the general N-body problem. Expressing the action integral at any intermediate points along the curve segment using a discrete Lagrangian that depends only on the end points of the interval, high order integrators can be obtained by defining the discrete Lagrangian in any time segment as a weighted sum on intermediate points, whose expressions for positions and velocities use Galerkin interpolation techniques. When oscillatory behavior is taken into account, the methods derived use trigonometric interpolation functions that depend on a frequency, which needs to be estimated. For that, using phase lag analysis, a new way to derive methods has been developed, that uses frequency estimation for each body at every time step. Results on special cases of the N-body problem show more stable orbits and less energy error when compared with the linear interpolation scheme. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)