Distance in the space of energetically bounded Keplerian orbits

In this paper we derive an explicit, analytic formula for the geodesic distance between two points in the space of bounded Keplerian orbits in a particular topology. The specific topology we use is that of a cone passing through the direct product of two spheres. The two spheres constitute the configuration manifold for the space of bounded orbits of constant energy. We scale these spheres by a factor equal to the semi-major axis of the orbit, forming a linear cone. This five-dimensional manifold inherits a Riemannian metric, which is induced from the Euclidean metric on \mathbbR⁶{\mathbb{R}^6}, the space in which it is embedded. We derive an explicit formula for the geodesic distance between any two points in this space, each point representing a physical, gravitationally bound Keplerian orbit. Finally we derive an expression for the Riemannian metric that we used in terms of classical orbital elements, which may be thought of as local coordinates on our configuration manifold.     

Metrics in Keplerian orbits quotient spaces

Quotient spaces of Keplerian orbits are important instruments for the modelling of orbit samples of celestial bodies on a large time span. We suppose that variations of the orbital eccentricities, inclinations and semi-major axes remain sufficiently small, while arbitrary perturbations are allowed for the arguments of pericentres or longitudes of the nodes, or both. The distance between orbits or their images in quotient spaces serves as a numerical criterion for such problems of Celestial Mechanics as search for common origin of meteoroid streams, comets, and asteroids, asteroid families identification, and others. In this paper, we consider quotient sets of the non-rectilinear Keplerian orbits space \(\mathbb H\). Their elements are identified irrespective of the values of pericentre arguments or node longitudes. We prove that distance functions on the quotient sets, introduced in Kholshevnikov et al. (Mon Not R Astron Soc 462:2275–2283, 2016), satisfy metric space axioms and discuss theoretical and practical importance of this result. Isometric embeddings of the quotient spaces into \(\mathbb R^n\), and a space of compact subsets of \(\mathbb H\) with Hausdorff metric are constructed. The Euclidean representations of the orbits spaces find its applications in a problem of orbit averaging and computational algorithms specific to Euclidean space. We also explore completions of \(\mathbb H\) and its quotient spaces with respect to corresponding metrics and establish a relation between elements of the extended spaces and rectilinear trajectories. Distance between an orbit and subsets of elliptic and hyperbolic orbits is calculated. This quantity provides an upper bound for the metric value in a problem of close orbits identification. Finally the invariance of the equivalence relations in \(\mathbb H\) under coordinates change is discussed.     

Orbit Determination of Binary Asteroids

In addition to the detection of an asteroid moon or a binary asteroid, the knowledge of the satellite’s true orbit is of high importance to derive fundamental physical parameters of the binary system such as its mass and to shed light on its possible formation history and dynamical evolution (prograde/retrograde orbit, large/small eccentricity or inclination, etc.). A new methodology for preliminary orbit determination of binary asteroids – and visual binaries in general – is proposed. It is based on Thiele–Innes method combined with a ‘trial and error’ Monte-Carlo technique. This method provides the full set of solutions (bundle of orbits, with the 7 orbital elements) even for a reduced number of observations. The mass is a direct by-product of this orbit determination, from which one can next infer the bulk-density and porosity. In addition to the bundle of orbits, the method provides the marginal probability densities of the foreseen parameters. Such error analysis – since it avoids linear approximation – can be of importance for the prediction of the satellite’s position in the plane-of-sky during future stellar occultations or subsequent observations, but also for the analysis of the orbit’s secular evolution. After briefly describing the method, we present the algorithm and its application to some practical cases, with particular emphasis on asteroids binaries and applications on orbital evolution.     

Lagrange’s planetary equations for the motion of electrostatically charged spacecraft

A spacecraft that generates an electrostatic charge on its surface in a planetary magnetic field will be subject to a perturbative Lorentz force. Active modulation of the surface charge can take advantage of this electromagnetic perturbation to modify or to do work on the spacecraft’s orbit. Lagrange’s planetary equations are derived using the Lorentz force as the perturbation on a Keplerian orbit, incorporating orbital inclination and true anomaly for the first time for an electrostatically charged vehicle. The planetary equations reveal that orbital inclination is a second-order effect on the perturbation, explaining results found in earlier studies through numerical integration. All of the orbital elements are coupled, but the coupling notably does not depend on the magnitude of the electrostatic charge or on the strength of the magnetic field. Analytical expressions that characterize this coupling are tested with a propellantless escape example at Jupiter. A closed-form solution exists that constrains the set of equatorial orbits for which planetary escape is possible, and a sufficient condition is identified for escape from inclined orbits. The analytical solutions agree with results from the numerically integrated equations of motion to within a fraction of a percent.     

Explicit Symplectic Integrator for Rotating Satellites

An explicit symplectic integrator is constructed for the problem of a rotating planetary satellite on a Keplerian orbit. The spin vector is fixed perpendicularly to the orbital plane. The integrator is constructed according to the Wisdom-Holman approach: the Hamiltonian is separated in two parts so that one of them is multiplied by a small parameter. The parameter depends on the satellite’s shape or the eccentricity of its orbit. The leading part of the Hamiltonian for small eccentricity orbits is similar to the simple pendulum and hence integrable; the perturbation does not depend on angular momentum which implies a trivial ‘kick’ solution. In spite of the necessity to evaluate elliptic function at each step, the explicit symplectic integrator proves to be quite efficient. This revised version was published online in July 2006 with corrections to the Cover Date.     

Secular orbit variation due to solar radiation effects: a detailed model for BYORP

A detailed derivation of the effect of solar radiation pressure on the orbit of a body about a primary orbiting the Sun is given. The result is a set of secular equations that can be used for long-term predictions of changes in the orbit. Solar radiation pressure is modeled as a Fourier series in the body’s rotation state, where the coefficients are based on the shape and radiation properties of the body as parameters. In this work, the assumption is made that the body is in a synchronous orbit about the primary and rotates at a constant rate. This model is used to write explicit variational equations of the energy, eccentricity vector, and angular momentum vector for an orbiting body. Given that the effect of the solar radiation pressure and the orbit are periodic functions, they are readily averaged over an orbit. Furthermore, the equations can be averaged again over the orbit of the primary about the Sun to give secular equations for long-term prediction. This methodology is applied to both circular and elliptical orbits, and the full equations for secular changes to the orbit in both cases are presented. These results can be applied to natural systems, such as the binary asteroid system 1999 KW4, to predict their evolution due to the Binary YORP effect, or to artificial Earth orbiting, nadir-pointing satellites to enable more precise models for their orbital evolution.     

Elliptical multi-sun-synchronous orbits for Mars exploration

The multi-sun-synchronous orbits allow cycles of observation of the same area in which solar illumination repetitively changes according to the value of the orbit elements and returns to the initial condition after a temporal interval multiple of the repetition of observation. This paper generalizes the concept of multi-sun-synchronous orbits, whose classical sun-synchronous orbits represent particular solutions, taking into consideration the elliptical case. The feasibility of using this typology of orbits, referred to as elliptical periodic multi-sun-synchronous orbits, has been investigated for the exploration of Mars and particular solutions have been selected. Such solutions considerably reduce the manoeuvre of velocity variation at the end of the interplanetary transfer with respect to the case of a target circular orbit around Mars. They are based on the use of quasi-critical inclinations in order to minimize the apsidal line motion and thus reduce orbit maintenance costs. Moreover, in the case of high eccentricities, the argument of pericentre may be set in order to obtain, around the apocentre, a condition of quasi-synchronism with the planet (the footprint of the probe on the surface presents a small shift with respect to a fixed point on the Martian surface). The low altitude of pericentre allows observation of the planet at a higher spatial resolution, while the orbit arc around the apocentre may be used to observe Mars with a wide spatial coverage in quasi-stationary conditions. This latter characteristic is useful for analysing atmospheric and meteorological phenomena and it allows for most of the orbital period a link between a rover on the surface of Mars and a probe orbiting around the planet.     

引力常数随时间变化对双星轨道演变的长期效应(椭圆轨道情形)

给出了以偏近点角为自变量的变引力常数的摄动方程组的解.解包括轨道半长轴的长期和周期变化项,其他轨道根数在一阶解中无长期项,只有周期项.近星点经度和平经度在二阶解中显示长期项变化.给出了由于引力常数变化对双星轨道演变情况的数值估计,对结果做了讨论并给出结论.     

On constructing the general Earth’s rotation theory

In the present paper the equations of the orbital motion of the major planets and the Moon and the equations of the three–axial rigid Earth’s rotation in Euler parameters are reduced to the secular system describing the evolution of the planetary and lunar orbits (independent of the Earth’s rotation) and the evolution of the Earth’s rotation (depending on the planetary and lunar evolution). Hence, the theory of the Earth’s rotation can be presented by means of the series in powers of the evolutionary variables with quasi-periodic coefficients with respect to the planetary–lunar mean longitudes. This form of the Earth’s rotation problem is compatible with the general planetary theory involving the separation of the short–period and long–period variables and avoiding the appearance of the non–physical secular terms.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

On a two-body problem with periodically changing equivalent gravitational parameter

Assigning to the equivalent gravitational parameter of a two-body dynamic system, a periodic change of a small amplitude B and arbitrary frequency and phase, the behaviour of an elliptic-type orbit is studied. The first order (in B) perturbations of the orbital elements are determined by using Delaunay's canonical variables. According to the value of the ratio between oscillation frequency and dynamic frequency, three cases (non-resonant (NR), quasi-resonant (QR), and resonant (R) ones) are pointed out. The solution of motion equations shows that only in the QR and R cases there are elements (argument of pericentre and mean anomaly) affected by secular perturbations. The solutions are valid over prediction times of order of pericentre and mean anomaly) affected by secular perturbations. The solutions are valid over prediction times of order B⁻¹ in the NR case and B^−1/2 in the QR and R cases.     

Gravitational Radiation Damping and Evolution of the Orbit of Compact Binary Stars (Solution by the Second Perturbation Method)

The influence of the gravitational radiation damping on the evolution of the orbital elements of compact binary stars is examined by using the method of perturbation. The perturbation equations with the true anomaly as an independent variable are given. This effect results in both the secular and periodic variation of the semi-major axis, the eccentricity, the mean longitude at the epoch and the mean longitude. However, the longitude of periastron exhibits no secular variation, but only periodic variation. The effect of secular variation of the orbit would lead to collapse of the system of binary stars. The deduced formulae are applied to the calculation of secular variation of the orbital elements for three compact binary stars: PSR19 13 + 16, PSR J0737-3039 and M33X-7. The results obtained are discussed.     

Orbital evolution of giant comet-like objects

We investigated by numerical integrations the long-term orbital evolution of four giant comets or comet-like objects. They are Chiron, P/Schwassmann-Wachmann 1 (SW1), Hidalgo, and 1992AD (5145), and their orbits were traced for 100–200 thousand years (kyr) toward both the past and the future. For each object, 13 orbits were calculated, one for the nominal orbital elements and other 12 with slightly modified elements based on the rms residual of the orbit determination and on the number of observations. As past studies indicate, their orbital evolution is found to be very chaotic, and thus can be described only in terms of probability. Plots of the semi-major axis (a) and perihelion distance (q) of the objects treated here seem to cross each other frequently, suggesting a possibility of their common evolutionary paths. About a half of all the calculated orbits showedq- ora-decreasing evolution. This indicates that, at least on the time scale in question, the giant comet-like objects are possibly on a dynamical track that can lead to capture from the outer solar system. We could hardly find the orbits with perihelia far outside the orbit of Saturn (q>15 AU). This is perhaps because the evolution of the orbits beyond Saturn is so slow that substantial orbital changes do not take place within 100–200 kyr.     

Position and velocity perturbations for the determination of geopotential from space geodetic measurements

Although space geodetic observing systems have been advanced recently to such a revolutionary level that low Earth Orbiting (LEO) satellites can now be tracked almost continuously and at the unprecedented high accuracy, none of the three basic methods for mapping the Earth’s gravity field, namely, Kaula linear perturbation, the numerical integration method and the orbit energy-based method, could meet the demand of these challenging data. Some theoretical effort has been made in order to establish comparable mathematical modellings for these measurements, notably by Mayer-Gürr et al. (J Geod 78:462–480, 2005). Although the numerical integration method has been routinely used to produce models of the Earth’s gravity field, for example, from recent satellite gravity missions CHAMP and GRACE, the modelling error of the method increases with the increase of the length of an arc. In order to best exploit the almost continuity and unprecedented high accuracy provided by modern space observing technology for the determination of the Earth’s gravity field, we propose using measured orbits as approximate values and derive the corresponding coordinate and velocity perturbations. The perturbations derived are quasi-linear, linear and of second-order approximation. Unlike conventional perturbation techniques which are only valid in the vicinity of reference mean values, our coordinate and velocity perturbations are mathematically valid uniformly through a whole orbital arc of any length. In particular, the derived coordinate and velocity perturbations are free of singularity due to the critical inclination and resonance inherent in the solution of artificial satellite motion by using various types of orbital elements. We then transform the coordinate and velocity perturbations into those of the six Keplerian orbital elements. For completeness, we also briefly outline how to use the derived coordinate and velocity perturbations to establish observation equations of space geodetic measurements for the determination of geopotential.     

The effect of the radiation pressure on the orbital evolution of geosynchronous objects

The effect of the radiation pressure and Poynting-Robertson effect on the evolution of the orbits of geosynchronous satellites is studied, depending on their area to mass ratio. The qualitative changes of the orbital evolution caused by these disturbances are considered. The reflection coefficient of the satellite’s surface was assumed to be 1.44. In the vicinity of the stable point with the longitude of 75° the exit from the libration resonance mode was registered when the area to mass ratio value changed from 5.9 to 6.0 m²/kg; in the vicinity of the unstable point at 345° with the area to mass ratio of 1.4 it occurred at 1.5 m²/kg. Re-entry to Earth occurs at values of the area to mass ratio above 32.2 m²/kg, and hyperbolic exit from the low-Earth orbit occurs at values of the area to mass ratio over 5267 m²/kg. At high values of the area to mass ratio, slopes of initially equatorial orbits can reach 49°. It is shown that due to the Poynting-Robertson effect the secular decrease in the semimajor axis of orbit in libration resonance region is 3–4 orders of magnitude less than outside of it.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

Secular influence of change in the heliocentric gravitation constant <Emphasis Type="Italic">GM</Emphasis><Subscript>⊙</Subscript> on evolution of orbits of Meteor Streams

The Secular influence of the change in the heliocentric gravitational constant on the evolution of orbits of Meteor Streams is examined by using the method of celestial mechanics with variable mass and variable gravitational constant. The change in the heliocentric gravitational constant includes the combined changes in the sun’s mass and gravitational constant obtained from the modern observation of planets and spacecraft. The perturbation equations are solved by expanding series with mean anomaly. The solutions of the secular and periodic variation of orbital elements are derived. The theoretical results for the secular variables of the semi-major axes, solar distances at perihelion and orbital periods are given for three Meteor Streams: Dracorids, Quadrantids, and Ursids. The numerical results are shown in Table 2. The discussion and conclusion are drawn.     

A numerical-analytical method for studying the orbital evolution of distant planetary satellites

We describe an approximate numerical-analytical method for calculating the perturbations of the elements of distant satellite orbits. The model for the motion of a distant satellite includes the solar attraction and the eccentricity and ecliptic inclination of the orbit of the central planet. In addition, we take into account the variations in planetary orbital elements with time due to secular perturbations. Our work is based on Zeipel’s method for constructing the canonical transformations that relate osculating satellite orbital elements to the mean ones. The corresponding transformation of the Hamiltonian is used to construct an evolution system of equations for mean elements. The numerical solution of this system free from rapidly oscillating functions and the inverse transformation from the mean to osculating elements allows the evolution of distant satellite orbits to be studied on long time scales on the order of several hundred or thousand satellite orbital periods.     

On the Milankovitch orbital elements for perturbed Keplerian motion

We consider sets of natural vectorial orbital elements of the Milankovitch type for perturbed Keplerian motion. These elements are closely related to the two vectorial first integrals of the unperturbed two-body problem; namely, the angular momentum vector and the Laplace–Runge–Lenz vector. After a detailed historical discussion of the origin and development of such elements, nonsingular equations for the time variations of these sets of elements under perturbations are established, both in Lagrangian and Gaussian form. After averaging, a compact, elegant, and symmetrical form of secular Milankovitch-like equations is obtained, which reminds of the structure of canonical systems of equations in Hamiltonian mechanics. As an application of this vectorial formulation, we analyze the motion of an object orbiting about a planet (idealized as a point mass moving in a heliocentric elliptical orbit) and subject to solar radiation pressure acceleration (obeying an inverse-square law). We show that the corresponding secular problem is integrable and we give an explicit closed-form solution.