Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.

Hamilton–Pontryagin spectral-collocation methods for the orbit propagation

According to the discrete Hamilton–Pontryagin variational principle,we construct a class of variational integrators in the real vector spaces and extend to the Lie groups for the left-trivialized Lagrangian mechanical systems by employing the spectral-collocation method to discretize the corresponding Lagrangian and kinematic constraints.The constructed framework can be transformed easily to the well-known symplectic partitioned Runge–Kutta methods and the higher order symplectic partitioned Lie Group methods by choosing same interpolation nodes and quadrature points.Two numerical experiments about the orbit propagation of Kepler two-body system and the rigid-body flow propagation of a free rigid body are conducted,respectively.The simulating results reveal that the constructed update schemes can possess simultaneously the excellent exponent convergence rates of spectral methods and the attractive long-term structure-preserving properties of geometric numerical algorithms.

Symplectic orbit propagation based on Deprit’s radial intermediary

The constantly challenging requirements for orbit prediction have opened the need for better onboard propagation tools.Runge-Kutta(RK)integrators have been widely used for this purpose;however RK integrators are not symplectic,which means that RK integrators may lead to incorrect global behavior and degraded accuracy.Emanating from Deprit’s radial intermediary,obtained by the elimination of the parallax transformation,we present the development of symplectic integrators of different orders for spacecraft orbit propagation.Through a set of numerical simulations,it is shown that these integrators are more accurate and substantially faster than Runge-Kutta-based methods.Moreover,it is also shown that the proposed integrators are more accurate than analytic propagation algorithms based on Deprit’s radial intermediary solution,and even other previously-developed symplectic integrators.