Superosculating intermediate orbits: Fourth-and fifth-order tangencies to trajectories of perturbed motion

The theory of superosculating intermediate orbits previously suggested by the author is developed. A new class of orbits with a fourth-order tangency to the actual trajectory of a celestial body at the initial time is constructed. Orbits with a fifth-order tangency have been constructed for the first time. The motion in the constructed orbits is represented as a combination of two motions: the motion of a fictitious attracting center with a variable mass and the motion relative to this center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden—Mestschersky problem. The variation in the mass of the fictitious center obeys Mestschersky’s first and combined laws. The new orbits represent more accurately the actual motion in the initial segment of the trajectory than an osculating Keplerian orbit and other existing analogues. Encke’s generalized methods of special perturbations in which the constructed intermediate orbits are used as reference orbits are presented. Numerical simulations using the approximations of the motions of Asteroid Toutatis and Comet P/Honda—Mrkos—Pajdu?áková as examples confirm that the constructed orbits are highly efficient. Their application is particularly beneficial in investigating strongly perturbed motion.     

Osculating and Superosculating Intermediate Orbits and their Applications

A method of construction of intermediate orbits for approximating the real motion of celestial bodies in the initial part of trajectory is proposed. The method is based on introducing a fictitious attracting centre with a time-variable gravitational parameter. The variation of thisparameter is assumed to obey the Eddington–Jeans mass-variationlaw. New classes of orbits having first-, second-, and third-order tangency to the perturbed trajectory at the initial instant of time are constructed. For planar motion, the tangency increases by one or two orders. The constructed intermediate orbits approximate the perturbed motion better than the osculating Keplerian orbit and analogous orbits of otherauthors. The applications of the orbits constructed in Encke's methodfor special perturbations and in the procedure for predicting themotion in which the perturbed trajectory is represented by a sequenceof short arcs of the intermediate orbits are suggested.The use of the constructed orbits is especially advantageous in the investigation of motion under the action of large perturbations.     

The Determination of an Intermediate Perturbed Orbit from Two Position Vectors

Based on the theory of intermediate orbits developed earlier by the author of this paper, a new approach is proposed to the solution of the problem of finding the orbit of a celestial body with the use of two position vectors of this body and the corresponding time interval. This approach makes it possible to take into account the main part of perturbations. The orbit is constructed, the motion along which is a combination of two motions: the uniform motion along a straight line of a fictitious attracting center, whose mass varies according to the first Meshchersky law, and the motion around this center. The latter is described by the equations of the Gylden–Meshchersky problem. The parameters of the constructed orbit are chosen so that their limiting values at any reference epoch determine a superosculating intermediate orbit with third-order tangency. The accuracy of approximation of the perturbed motion by the orbits calculated by the classical Gauss method and the new method is illustrated by an example of the motion of the unusual minor planet 1566 Icarus. Comparison of the results obtained shows that the new method has obvious advantages over the Gauss method. These advantages are especially prominent in cases where the angular distances between the reference positions are small.     

Determination of an intermediate perturbed orbit from three position vectors

We suggest a new approach to solving the problem of finding the orbit of a celestial body from its three spatial position vectors and the corresponding times. It allows most of the perturbations in the motion of a celestial body to be taken into account. The approach is based on the theory of intermediate orbits that we developed previously. We construct the orbit the motion along which is a combination of two motions: the motion of a fictitious attracting center whose mass varies according to Mestschersky’s first law and the motion relative to the fictitious center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden-Mestschersky problem. The constructed orbit has such parameters that their limiting values at any reference epoch define a superosculating intermediate orbit with a fourth-order tangency. We have performed a numerical analysis to estimate the accuracy of approximating the perturbed motion of two minor planets, 145 Adeona and 4179 Toutatis, by the orbits computed using two-position procedures (the classical Gauss method and the method that we suggested previously), a three-position procedure based on the Herrick-Gibbs equation, and the new method. Comparison of the results obtained suggests that the latter method has an advantage.     

On the accuracy of approximation of motion of a small celestial body by intermediate perturbed orbits calculated on the basis of three position vectors and three observations

Intermediate perturbed orbits, which were proposed earlier by the first author and are calculated based on three position vectors and three measurements of angular coordinates of a small celestial body, are examined. Provided that the reference time interval encompassing the measurements is short, these orbits are close in the accuracy of approximation of actual motion to an orbit with fourth-order tangency. The shorter the reference time interval is, the better is the approximation. The laws of variation of the errors of methods for constructing such intermediate orbits with the length of the reference time interval are formulated. According to these laws, the rate of convergence of the methods to an exact solution in the process of shortening of the reference time interval is, in general, three orders of magnitude higher than that of conventional methods relying on an unperturbed Keplerian orbit. The considered orbits are among the most accurate of their class that is defined by the order of tangency. The obtained theoretical results are verified by numerical experiments on determining the orbit of 99942 Apophis.     

Computation of Partial Derivatives of Perturbed Motion Using the Arcs of Osculating and Superosculating Orbits

The methods for analytical determination of partial derivatives of the current parameters of motion with respect to their initial values are described. The methods take into account principal perturbations and are based on the use of the osculating and superosculating intermediate orbits constructed earlier by the author. These orbits ensure the first-, second-, and third-order contact to the real trajectory at the initial time. The solution for parameters of the intermediate motion and partial derivatives of these parameters is given in a universal closed form. The partial derivatives on long time intervals are computed using a step-by-step procedure combined with the Encke method of special perturbations, in which the intermediate orbits are used as the reference. The numerical results show that the new approach can be efficiently used for solving the problem of differential correction of orbits of asteroids and comets on the basis of observational data.     

Short-Period Comets with a High Tisserand Constant: III. Kinematics of Low-Velocity Encounters

Within the framework of a pair two-body problem (Sun–Jupiter, Sun–comet), the kinematics of the encounter of a minor body with a planet is investigated. The notion of points of low-velocity tangency of the orbits of the comet and Jupiter, as well as the point of Jovicentric velocity and the low-velocity tangent section of a cometary orbit, is introduced. The conditions and definitions of low-velocity and high-velocity encounters are proposed. The systems of inequalities relating the aand eparameters, which make it possible to single out those comets that are likely to be objects with low-velocity encounters, are presented. The regions of orbits that have low-velocity tangent sections, i.e., regions of low-velocity tangency of orbits, are singled out on the (a, e) plane. These regions agree well with the corresponding parameters of the orbits of real comets whose evolution contains low-velocity encounters with Jupiter.     

Matrix perturbation methods using regularized coordinates

Matrix methods for computing perturbations of non-linear perturbed systems, as formulated by Alexeev, involve an expression for the full solution of the first variational equations of the system evaluated about a reference orbit. These cannot be immediately applied to a regularized system of equations where perturbations about Keplerian motion are considered since the solution of the variational equations of regularized Keplerian motion does not in general correspond to the solution of the variational equations of the unregularized equations. But, as Kustaanheimo and Stiefel have pointed out, the regularized equations of Keplerian motion should be excellent for the initiation of a perturbation theory since they are linear in form. This paper describes a method for applying Alexeev's theorem to a regularized system where full advantage is taken of the basic linear form of the unperturbed equations.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969.     

The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

Approximate Analytical Theory of Satellite Orbit Prediction Presented in Spherical Coordinates Frame

A comparative review of analytic theories for the motion of Earth satellites in quasi-circular orbits written in the spherical coordinate frame is presented. The theory of motion is developed for satellites in quasi-circular and quasi-equatorial orbits subjected to geopotential, luni-solar and solar radiation pressure force perturbations. The intermediate orbit is Keplerian and the equations of motion are solved by the Lyapunov–Poincaré small parameter method. Both resonant and non-resonant cases are considered. The results can be useful for the development of a complete theory of weakly eccentric orbits.     

The equations of relative motion in the orbital reference frame

The analysis of relative motion of two spacecraft in Earth-bound orbits is usually carried out on the basis of simplifying assumptions. In particular, the reference spacecraft is assumed to follow a circular orbit, in which case the equations of relative motion are governed by the well-known Hill–Clohessy–Wiltshire equations. Circular motion is not, however, a solution when the Earth’s flattening is accounted for, except for equatorial orbits, where in any case the acceleration term is not Newtonian. Several attempts have been made to account for the \(J_2\) effects, either by ingeniously taking advantage of their differential effects, or by cleverly introducing ad-hoc terms in the equations of motion on the basis of geometrical analysis of the \(J_2\) perturbing effects. Analysis of relative motion about an unperturbed elliptical orbit is the next step in complexity. Relative motion about a \(J_2\)-perturbed elliptic reference trajectory is clearly a challenging problem, which has received little attention. All these problems are based on either the Hill–Clohessy–Wiltshire equations for circular reference motion, or the de Vries/Tschauner–Hempel equations for elliptical reference motion, which are both approximate versions of the exact equations of relative motion. The main difference between the exact and approximate forms of these equations consists in the expression for the angular velocity and the angular acceleration of the rotating reference frame with respect to an inertial reference frame. The rotating reference frame is invariably taken as the local orbital frame, i.e., the RTN frame generated by the radial, the transverse, and the normal directions along the primary spacecraft orbit. Some authors have tried to account for the non-constant nature of the angular velocity vector, but have limited their correction to a mean motion value consistent with the \(J_2\) perturbation terms. However, the angular velocity vector is also affected in direction, which causes precession of the node and the argument of perigee, i.e., of the entire orbital plane. Here we provide a derivation of the exact equations of relative motion by expressing the angular velocity of the RTN frame in terms of the state vector of the reference spacecraft. As such, these equations are completely general, in the sense that the orbit of the reference spacecraft need only be known through its ephemeris, and therefore subject to any force field whatever. It is also shown that these equations reduce to either the Hill–Clohessy–Wiltshire, or the Tschauner–Hempel equations, depending on the level of approximation. The explicit form of the equations of relative motion with respect to a \(J_2\)-perturbed reference orbit is also introduced.     

Earth–Moon Separation Problem In The Motion Of Near Earth Asteroids

The analysis of application of two dynamical models (``Earth–Moon' and``barycentre' model) in the motion of Near Earth Asteroids was performed. Mainaim was the quantitative estimation of the influence of lunar perturbations on the motionof NEA. Additionally, basic tests of application of numerical methods weremade (RMVS3 and B–S methods). The orbits of 1083 Apollo–Aten–Amor and 7selected AAA objects were adopted as test particles in numerical integrationof the motion. The comparison between results obtained by both dynamicalmodels is discussed in detail. In specific cases, the application of the``Earth–Moon' dynamical model is very important and cannot be neglected incomputations of orbits.     

Analytical method for perturbed frozen orbit around an asteroid in highly inhomogeneous gravitational fields: a first approach

This article provides a method for finding initial conditions for perturbed frozen orbits around inhomogeneous fast rotating asteroids. These orbits can be used as reference trajectories in missions that require close inspection of any rigid body. The generalized perturbative procedure followed exploits the analytical methods of relegation of the argument of node and Delaunay normalisation to arbitrary order. These analytical methods are extremely powerful but highly computational. The gravitational potential of the heterogeneous body is firstly stated, in polar-nodal coordinates, which takes into account the coefficients of the spherical harmonics up to an arbitrary order. Through the relegation of the argument of node and the Delaunay normalization, a series of canonical transformations of coordinates is found, which reduces the Hamiltonian describing the system to a integrable, two degrees of freedom Hamiltonian plus a truncated reminder of higher order. Setting eccentricity, argument of pericenter and inclination of the orbit of the truncated system to be constant, initial conditions are found, which evolve into frozen orbits for the truncated system. Using the same initial conditions yields perturbed frozen orbits for the full system, whose perturbation decreases with the consideration of arbitrary homologic equations in the relegation and normalization procedures. Such procedure can be automated for the first homologic equation up to the consideration of any arbitrary number of spherical harmonics coefficients. The project has been developed in collaboration with the European Space Agency (ESA).     

Numerical stabilization of the differential equations of Keplerian motion

A stabilization of the classical equations of two-body motion is offered. It is characterized by the use of the regularizing independent variable (eccentric anomaly) and by the addition of a control-term to the differential equations. This method is related to the KS-theory (Stiefel, 1970) which performed for the first time a stabilization of the Kepler motion. But in contrast to the KS-theory our method does not transform the coordinates of the particle. As far as the theory of stability and the numerical experiments are concerned we restrict ourselves to thepure Kepler motion. But, of course, the stabilizing devices will also improve the accuracy of the computation of perturbed orbits. We list, therefore, also the equations of the perturbed motion.     

Orbital evolution dynamics of two satellites in encounter phase using multiple scales expansion

An approximate solution of the encounter problem of two small satellites describing initially elliptical orbits around a massive oblate primary is obtained. The equations of motion of the center of mass of the two masses are developed in the most general form without any restrictions on the orbital elements. The method of multiple scales which seeks a solution whose behavior depends on several time scales is used. To overcome the singularity the equations of motion are transformed to the Struble variables. An analytical second order theory of the evolution dynamics is obtained. A MATHEMATICA program is constructed. The evolution dynamics of the orbital parameters between the perturbed and the unperturbed cases are plotted. The effect of changing eccentricity and changing inclination on the orbital parameters are highlighted.     

The stability of periodic orbits in the three-body problem

A method is developed to study the stability of periodic motions of the three-body problem in a rotating frame of reference, based on the notion of surface of section. The method is linear and involves the computation of a 4×4 variational matrix by integrating numerically the differential equations for time intervals of the order of a period. Several properties of this matrix are proved and also it is shown that for a symmetric periodic motion it can be computed by integrating for half the period only.This linear stability analysis is used to study the stability of a family of periodic motions of three bodies with equal masses, in a rotating frame of reference. This family represents motion such that two bodies revolve around each other and the third body revolves around this binary system in the same direction to a distance which varies along the members of the family. It was found that a large part of the family, corresponding to the case where the distance of the third body from the binary system is larger than the dimensions of the binary system, represents stable motion. The nonlinear effects to the linear stability analysis are studied by computing the intersections of several perturbed orbits with the surface of sectiony ₃=0. In some cases more than 1000 intersections are computed. These numerical results indicate that linear stability implies stability to all orders, and this is true for quite large perturbations.     

Oblateneess and drag effects on the motion of satellites in the set of Eulerian redundant parameters

In this paper, the perturbed motion of the artificial satellites under the effects of the earth's oblateness and atmospheric drag will be established. These equations will be expressed in terms of the Eulerian redundant parameters. Applications of the method for the perturbed motion are illustrated by numerical examples for some test cases of the orbits.     

Accurate computation of highly eccentric satellite orbits

For computing highly eccentric (e0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen.     

Studies in the application of recurrence relations to special perturbation methods

Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the Encke and Cowell methods of integration. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relationships. It is shown that for fairly highly eccentric orbits in which the perturbing force is less than one thousandth of the two-body force the Encke method achieves a considerable saving in machine time. This is also true for almost circular orbits when low or moderate accuracy is required. When very high accuracy is required, however, the Cowell method is faster unless the perturbing force is less than 10^–6 of the two-body force. There is little difference in the accuracy of the two methods, the Cowell method being slightly more accurate when a low or moderate accuracy criterion is imposed.     

Relative motion of satellites exploiting the super-integrability of Kepler’s problem

This paper builds upon the work of Palmer and Imre exploring the relative motion of satellites on neighbouring Keplerian orbits. We make use of a general geometrical setting from Hamiltonian systems theory to obtain analytical solutions of the variational Kepler equations in an Earth centred inertial coordinate frame in terms of the relevant conserved quantities: relative energy, relative angular momentum and the relative eccentricity vector. The paper extends the work on relative satellite motion by providing solutions about any elliptic, parabolic or hyperbolic reference trajectory, including the zero angular momentum case. The geometrical framework assists the design of complex formation flying trajectories. This is demonstrated by the construction of a tetrahedral formation, described through the relevant conserved quantities, for which the satellites are on highly eccentric orbits around the Sun to visit the Kuiper belt.