The Characteristics and Related Problems of the Orbits Around the Earth-Moon Libration Points

Due to the specific dynamics, the probes located at the halo orbits or Lissajous orbits around the Earth-Moon collinear libration point L₁ or L₂ are always studied in the synodic system to understand their trajectories. In fact, they are also orbiting the Earth in a distant Keplerian ellipse. Because of their intrinsic orbital instability, in the orbit prediction the initial errors propagate more prominently than those of the normal orbiting satellites, this requires special attention in the orbit design, maneuver, and control. Despite of all this, they are similar to the normal orbiting satellites in orbit determination and hardly require other special attentions. In this paper, the quantitative results of error propagation under the unstable dynamics, together with the theoretical analysis are presented. The results of precise orbit determination and short-arc orbit predictions are also shown, and compared with the results from the Beijing Aerospace Control Center.     

Optimal periodic relative orbit and rectilinear relative orbits with eccentric reference orbits

The problem of two-body linearized periodic relative orbits with eccentric reference orbits is studied in this paper. The periodic relative orbit in the target-orbital coordinate system can be used in fly-around and formation-flying orbit design. Based on the closed-form solutions to the Tschauner–Hempel equations, the initial condition for periodic relative orbits is obtained. Then the minimum-fuel periodic-orbit condition with a single impulse is analytically derived for given initial position and velocity vectors. When considering the initial coasting time, the impulse position of the global minimum-fuel periodic orbit is proved to be near to the perigee of the target and can be obtained by numerical optimization algorithms. Moreover, the condition for a special periodic orbit, i.e., the rectilinear relative orbit in the target-orbital frame, is obtained. Numerical simulations are used to demonstrate the efficacy of the method, and show the geometry of the periodic relative orbit and the rectilinear relative orbit.     

Low-eccentricity motion of asteroids near the 2/1 Jovian resonance

(903) Nealley moves on an orbit of low eccentricity with a mean motion that is slightly larger than the 2/1 value of resonance. This orbit and some related fictious orbits are studied by numerical integrations of the four-body problem Sun-Jupiter-Saturn-asteroid over an interval of 110000 yr. The author's experience on related cases of resonance allows a study of the variation of suitably defined orbital parameters. The long-term evolution of the orbits is compared with earlier predictions. Some of the librating orbits are temporarily captured in a secondary resonance that refers to three-dimensional motion and is demonstrated by a special example.     

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies in introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The linear stability of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.     

月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

An exact solution to determination of an open orbit

We present an exact solution of the equations for orbit determination of a two body system in a hyperbolic or parabolic motion. In solving this problem, we extend the method employed by Asada, Akasaka and Kasai (AAK) for a binary system in an elliptic orbit. The solutions applicable to each of elliptic, hyperbolic and parabolic orbits are obtained by the new approach, and they are all expressed in an explicit form, remarkably, only in terms of elementary functions. We show also that the solutions for an open orbit are recovered by making a suitable transformation of the AAK solution for an elliptic case.     

A new family of periodic orbits for the restricted problem

A new family of periodic orbits of the three-dimensional restricted three-body problem which continue off from a consecutive collision orbit are numerically studied. Their behavior for varying energy is unexpected. In particular, associated with our system is a countable set of resonant energy values and each time the energy passes through one of them the periodic orbit forms a loop by self-intersection. Any number of loops can form by this process and the resulting orbits take on an interesting appearance.     

Libration of arguments of circumbinary-planet orbits at resonance

The paper refers to fictitious resonant orbits of planet type that surround both components of a binary system. In case of 16 studied examples a suitable choice of the starting values leads to a process of libration of special angular arguments and to an evolution with an at least temporary stay of the planet in the resonant orbit. The ratio of the periods of revolution of the binary and a planet is equal to 1:5. Eight orbits depend on the ratio 1:5 of the masses of the binary components, but two other ratios appear as well. The basis of this study is the planar, elliptic or circular restricted problem of three bodies, but remarks at the end of the text refer to a four-body problem.     

Spacecraft Orbit Determination with The B-spline Approximation Method

It is known that the dynamical orbit determination is the most common way to get the precise orbits of spacecraft. However, it is hard to build up the precise dynamical model of spacecraft sometimes. In order to solve this problem, the technique of the orbit determination with the B-spline approximation method based on the theory of function approximation is presented in this article. In order to verify the effectiveness of this method, simulative orbit determinations in the cases of LEO (Low Earth Orbit), MEO (Medium Earth Orbit), and HEO (Highly Eccentric Orbit) satellites are performed, and it is shown that this method has a reliable accuracy and stable solution. The approach can be performed in both the conventional celestial coordinate system and the conventional terrestrial coordinate system. The spacecraft's position and velocity can be calculated directly with the B-spline approximation method, it needs not to integrate the dynamical equations, nor to calculate the state transfer matrix, thus the burden of calculations in the orbit determination is reduced substantially relative to the dynamical orbit determination method. The technique not only has a certain theoretical significance, but also can serve as a conventional algorithm in the spacecraft orbit determination.     

Properties of linearized mappings associated with periodic orbits in the three-body problem

In this paper three results on the linearized mapping associated with the plane three body problem near a periodic orbit are established. It is first shown that linear stability of such an orbit is independent of initial position on the orbit and of coordinate system. Second, the relation of Hénon connecting the rates of change of rotation angle and period on an isoenergetic family of periodic orbits is proved, together with a similar relation for families of orbits closing exactly in a rotating coordinate system. Finally, a condition for a critical orbit is given which is applicable to any family of periodic orbits.     

The Symmetrical One-dimensional Newtonian Four-body Problem: A Numerical Investigation

Numerical simulations of the one-dimensional Newtonian four-body problem have been conducted for the special case in which the bodies are distributed symmetrically about the centre of mass. Simulations show a great similarity between this problem and the one-dimensional Newtonian three-body problem. As in that problem the orbits can be divided into three different categories which form well-defined regions on a Poincaré section: there is a region of quasiperiodic orbits about a Schubart-like periodic orbit, there is a region of fast-scattering encounters and in between these two regions there is a chaotic scattering region. The Schubart-like periodic orbit's stability to perturbation is studied. It is apparently stable in one-dimension but is unstable in three-dimensions.This revised version was published online in October 2005 with corrections to the Cover Date.     

The three-dimensional general three-body problem: Determination of periodic orbits

The three-dimensional general three-body problem is formulated suitably for the numerical determination of periodic orbits either directly or by continuation from the three-dimensional periodic orbits of the restricted problem. The symmetry properties of the equations of motion are established and the algorithms for the numerical determination of families of periodic orbits are outlined. A normalization scheme based on the concept of the invariable plane is introduced to simplify the process. All three types of symmetric orbit, as well as the general type of asymmetric orrbit, are considered. Many threedimmensional p periodic orbits are given.     

Exchange orbits: an interesting case of co-orbital motion

In this investigation we treat a special configuration of two celestial bodies in 1:1 mean motion resonance namely the so-called exchange orbits. There exist—at least—theoretically—two different types: the exchange-a orbits and the exchange-e orbits. The first one is the following: two celestial bodies are in orbit around a central body with almost the same semi-major axes on circular orbits. Because of the relatively small differences in semi-major axes they meet from time to time and exchange their semi-major axes. The inner one then moves outside the other planet and vice versa. The second configuration one is the following: two planets are moving on nearly the same orbit with respect to the semi-major axes, one on a circular orbit and the other one on an eccentric one. During their dynamical evolution they change the characteristics of the orbit, the circular one becomes an elliptic one whereas the elliptic one changes its shape to a circle. This ‘game’ repeats periodically. In this new study we extend the numerical computations for both of these exchange orbits to the three dimensional case and in another extension treat also the problem when these orbits are perturbed from a fourth body. Our results in form of graphs show quite well that for a large variety of initial conditions both configurations are stable and stay in these exchange orbits.     

On the effect of the eccentricity of a planetary orbit on the stability of satellite orbits

The effect of the eccentricity of a planet’s orbit on the stability of the orbits of its satellites is studied. The model used is the elliptic Hill case of the planar restricted three-body problem. The linear stability of all the known families of periodic orbits of the problem is computed. No stable orbits are found, the majority of them possessing one or two pairs of real eigenvalues of the monodromy matrix, while a part of a family with complex instability is found. Two families of periodic orbits, bifurcating from the Lagrangian points L₁, L₂ of the corresponding circular case are found analytically. These orbits are very unstable and the determination of their stability coefficients is not accurate, so we compute the largest Liapunov exponent in their vicinity. In all cases these exponents are positive, indicating the existence of chaotic motions     

Lambert problem solution in the hill model of motion

The goal of this paper is obtaining a solution of the Lambert problem in the restricted three-body problem described by the Hill equations. This solution is based on the use of pre determinate reference orbits of different types giving the first guess and defining the sought-for transfer type. A mathematical procedure giving the Lambert problem solution is described. This procedure provides step-by-step transformation of the reference orbit to the sought-for transfer orbit. Numerical examples of the procedure application to the transfers in the Sun–Earth system are considered. These examples include transfer between two specified positions in a given time, a periodic orbit design, a halo orbit design, halo-to-halo transfers, LEO-to-halo transfer, analysis of a family of the halo-to-halo transfer orbits. The proposed method of the Lambert problem solution can be used for the two-point boundary value problem solution in any model of motion if a set of typical reference orbits can be found.     

Orbit determination of space debris: admissible regions

The main problem in the orbit determination of the space debris population orbiting our planet is identifying which separate sets of data belong to the same physical object. The observations of a given object during a passage above an observing station are collectively called a Too Short Arc (TSA): data from a TSA cannot allow for a complete determination of an orbit. Therefore, we have to solve first the identification problem, finding two or more TSAs belonging to the same physical object and an orbit fitting all the observations. This problem is well known for the determination of orbits of asteroids: we shall show how to apply the methods developed for preliminary orbit determination of heliocentric objects to geocentric objects. We shall focus on the definition of an admissible region for space debris, both in the case of optical observations and radar observations; then we shall outline a strategy to perform a full orbit determination.     

Chaotic scattering in the restricted three-body problem I. The Copenhagen problem

We consider the scattering motion of the planar restricted three-body problem with two equal masses on a circular orbit. Using the methods of chaotic scattering we present results on the structure of scattering functions. Their connection with primitive periodic orbits and the underlying chaotic saddle are studied. Numerical evidence is presented which suggests that in some intervals of the Jacobi integral the system is hyperbolic. The Smale horseshoe found there is built from a countable infinite number of primitive periodic orbits, where the parabolic orbits play a fundamental role.     

Computing with certainty individual members of families of periodic orbits of a given period

The accurate computation of families of periodic orbits is very important in the analysis of various celestial mechanics systems. The main difficulty for the computation of a family of periodic orbits of a given period is the determination within a given region of an individual member of this family which corresponds to a periodic orbit. To compute with certainty accurate individual members of a specific family we apply an efficient method using the Poincaré map on a surface of section of the considered problem. This method converges rapidly, within relatively large regions of the initial conditions. It is also independent of the local dynamics near periodic orbits which is especially useful in the case of conservative dynamical systems that possess many periodic orbits, often of the same period, close to each other in phase space. The only computable information required by this method is the signs of various function evaluations carried out during the integration of the equations of motion. This method can be applied to any system of celestial mechanics. In this contribution we apply it to the photogravitational problem.     

Orbit Determination Using Satellite—to—Satellite Tracking Data

1 INTRODUCTIONThe tracking arc-length should be increased in order to approve the accuracy in orbitdetermination of LEO (Low Earth Orbit) satellites. The local ground-based tracking networkdoes not provide sufficient orbit coverage for the user satellites. The most promising methodis to use high orbiting satellites, such as GPS and TDRS, as trackers to observe the usersatellites. For examPle, tWO geosynclironous satellites could cover more than 85% of the orbitof any given user sate…     

Study of lunar gravity assist orbits in the restricted four-body problem

In this paper, the lunar gravity assist (LGA) orbits starting from the Earth are investigated in the Sun–Earth–Moon–spacecraft restricted four-body problem (RFBP). First of all, the sphere of influence of the Earth–Moon system (SOIEM) is derived. Numerical calculation displays that inside the SOIEM, the effect of the Sun on the LGA orbits is quite small, but outside the SOIEM, the Sun perturbation can remarkably influence the trend of the LGA orbit. To analyze the effect of the Sun, the RFBP outside the SOIEM is approximately replaced by a planar circular restricted three-body problem, where, in the latter case, the Sun and the Earth–Moon barycenter act as primaries. The stable manifolds associated with the libration point orbit and their Poincaré sections on the SOIEM are applied to investigating the LGA orbit. According to our research, the patched LGA orbits on the Poincaré sections can efficiently distinguish the transit LGA orbits from the non-transit LGA orbits under the RFBP. The former orbits can pass through the region around libration point away from the SOIEM, but the latter orbits will bounce back to the SOIEM. Besides, the stable transit probability is defined and analyzed. According to the variant requirement of the space mission, the results obtained can help us select the LGA orbit and the launch window.