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.     

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.     

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.     

New Numerical/Analytical Methods for Solving Equations of Orbital Motion

New methods are proposed for solving equations of motion of celestial bodies. The methods are based on the use of superosculating orbits with second- and third-order tangency to the trajectory of the real motion of a body. The construction of these orbits is related to the concept of a fictitious attracting center, whose mass varies in accordance with the first Meshchersky law. In the original reference methods, the perturbed trajectory is represented by a sequence of small arcs of superosculating orbits. The order of accuracy of the reference methods coincides with the order of tangency of the superosculating orbit used in calculations. Using Runge's rule and Richardson's extrapolation scheme leads to the methods of higher order. The efficiency of the new methods in comparison with the numerical integration of equations of motion based on the well-known fourth- and seventh-order Runge–Kutta–Fehlberg methods is illustrated by examples of the calculation of perturbed orbits of some asteroids.     

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.     

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.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.     

A Method for the Determination of an Intermediate Orbit from Three Positions of the Small Body on the Celestial Sphere

A new method is suggested for finding the preliminary orbit from three complete measurements of the angular coordinates of a celestial body developed by analogy with the classic Lagrange–Gauss method. The proposed method uses the intermediate orbit that we had constructed in an earlier paper based on two position vectors and the corresponding time interval. This intermediate orbit allows for most of the perturbations in the motion of the body. Using the orbital motion of asteroid 1566 Icarus as an example, we compare the results obtained by applying the classic and the new method. The comparison shows the new method to be highly efficient for studying perturbed motion. It is especially efficient if applied to high-precision observational data covering short orbital arcs.     

Explicit Symplectic Integrator for Highly Eccentric Orbits

An explicit symplectic integrator is constructed for perturbed elliptic orbits of an arbitrary eccentricity. The perturbation should be Hamiltonian, but it may depend on time explicitly. The main feature of the integrator is the use of KS variables in the ten-dimensional extended phase space. As an example of its application the motion of an Earth satellite under the action of the planet's oblateness and of lunar perturbations is studied. The results confirm the superiority of the method over a classical Wisdom–Holman algorithm in both accuracy and computation time. This revised version was published online in July 2006 with corrections to the Cover Date.     

Optimal transfers between unstable periodic orbits using invariant manifolds

This paper presents a method to construct optimal transfers between unstable periodic orbits of differing energies using invariant manifolds. The transfers constructed in this method asymptotically depart the initial orbit on a trajectory contained within the unstable manifold of the initial orbit and later, asymptotically arrive at the final orbit on a trajectory contained within the stable manifold of the final orbit. Primer vector theory is applied to a transfer to determine the optimal maneuvers required to create the bridging trajectory that connects the unstable and stable manifold trajectories. Transfers are constructed between unstable periodic orbits in the Sun–Earth, Earth–Moon, and Jupiter-Europa three-body systems. Multiple solutions are found between the same initial and final orbits, where certain solutions retrace interior portions of the trajectory. All transfers created satisfy the conditions for optimality. The costs of transfers constructed using manifolds are compared to the costs of transfers constructed without the use of manifolds. In all cases, the total cost of the transfer is significantly lower when invariant manifolds are used in the transfer construction. In many cases, the transfers that employ invariant manifolds are three times more efficient, in terms of fuel expenditure, than the transfer that do not. The decrease in transfer cost is accompanied by an increase in transfer time of flight.     

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.     

A new method for the determination of an intermediate orbit using three positions of a small body on the celestial sphere

We propose a new method for the determination of the preliminary orbit of a small celestial body using three pairs of its angular coordinates in three moments of time. The method is based on the use of the intermediate orbit we constructed earlier using three position vectors and the corresponding time moments. This intermediate orbit accounts for the main part of the perturbations of the motion of the body under study. We compare the results obtained by the classical Lagrange-Gauss method, Herrick-Gibbs method, generalized Herrick-Gibbs method, and the new method by the examples of the determination of the orbit of the small planet 1566 Icarus. The comparison showed that the new method is a highly efficient tool for the study of perturbed motion. It is especially efficient when applied to high-precision observational data covering short arcs of the orbit.     

Element sets for high-order Poincaré mapping of perturbed Keplerian motion

The propagation and Poincaré mapping of perturbed Keplerian motion is a key topic in Celestial Mechanics and Astrodynamics, e.g., to study the stability of orbits or design bounded relative trajectories. The high-order transfer map (HOTM) method enables efficient mapping of perturbed Keplerian orbits using the high-order Taylor expansion of a Poincaré or stroboscopic map. The HOTM is only accurate close to the expansion point and therefore the number of revolutions for which the map is accurate tends to be limited. The proper selection of coordinates is of key importance for improving the performance of the HOTM method. In this paper, we investigate the use of different element sets for expressing the high-order map in order to find the coordinates that perform best in terms of accuracy. A new set of elements is introduced that enables extremely accurate mapping of the state, even for high eccentricities and higher-order zonal perturbations. Finally, the high-order map is shown to be very useful for the determination and study of fixed points and center manifolds of Poincaré maps.     

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.     

Planetary orbital equations in externally-perturbed systems: position and velocity-dependent forces

The increasing number and variety of extrasolar planets illustrates the importance of characterizing planetary perturbations. Planetary orbits are typically described by physically intuitive orbital elements. Here, we explicitly express the equations of motion of the unaveraged perturbed two-body problem in terms of planetary orbital elements by using a generalized form of Gauss’ equations. We consider a varied set of position and velocity-dependent perturbations, and also derive relevant specific cases of the equations: when they are averaged over fast variables (the “adiabatic” approximation), and in the prograde and retrograde planar cases. In each instance, we delineate the properties of the equations. As brief demonstrations of potential applications, we consider the effect of Galactic tides. We measure the effect on the widest-known exoplanet orbit, Sedna-like objects, and distant scattered disk objects, particularly with regard to where the adiabatic approximation breaks down. The Mathematica code which can help derive the equations of motion for a user-defined perturbation is freely available upon request.     

A keplerian method to estimate perturbations in the restricted three-body problem

We develop a new and fast method to estimate perturbations by a planet on cometary orbits. This method allows us to identify accurately the cases of large perturbations in a set of fictitious orbits. Hence, it can be used in constructing perturbation samples for Monte Carlo simulations in order to maximize the amount of information. Furthermore, the estimated perturbations are found to yield a good approximation to the real perturbation sample. This is shown by a comparison of the perturbations obtained by the new estimator with the results of numerical integration of regularized equations of motion for the same orbits in the same dynamical model: the three-dimensional elliptic restricted three-body problem (Sun-Jupiter-comet).     

A two-level perturbation method for connecting unstable periodic orbits with low fuel cost and short time of flight: application to a lunar observation mission

The proposed method connects two unstable periodic orbits by employing trajectories of their associated invariant manifolds that are perturbed in two levels. A first level of velocity perturbations is applied on the trajectories of the discretized manifolds at the points where they approach the nominal unstable periodic orbit in order to accelerate them. A second level of structured velocity perturbations is applied to trajectories that have already been subjected to first level perturbations in order to approximately meet the necessary conditions for a low \(\varDelta \text {V}\) transfer. Due to this two-level perturbation approach, the number of the trajectories obtained is significantly larger compared with approaches that employ traditional invariant manifolds. For this reason, the problem of connecting two unstable periodic orbits through perturbed trajectories of their manifolds is transformed into an equivalent discrete optimization problem that is solved with a very low computational complexity algorithm that is proposed in this paper. Finally, the method is applied to a lunar observation mission of practical interest and is found to perform considerably better in terms of \(\varDelta \text {V}\) cost and time of flight when compared with previous techniques applied to the same project.     

Some particularities in motion of retrograde satellites of Jupiter

Some peculiarities in the motion of retrograde satellites of Jupiter have been investigated. The intermediate orbits were obtained by approximated solution of differential equations before transformation by the Zeipel's method. These orbits are non-keplerian ellipses. For their construction the secular motion of nodes, perijoves, and essential periodic perturbations were taken into account.The eccentricities and inclinations of all the retrograte satellites change in a large range. The motion may happen in a region, which is located very near to the limit cases of our theory. For some satellites the sign of the constant, which characterizes the type of orbit, librating or circular, may change. In some cases the value of this constant may be close to zero. Then the motion of the longitude of perijove will reduce the speed and in some moment the circular orbit may change its direction.