Formulation of the restricted three-body problem with two-center problem orbits as intermediate orbits

The restricted three-body Hamiltonian is partitioned into a two-center type principal part and its accompanying perturbational part. The mathematical analysis, involving the Jacobian elliptic functions, is adapted for the case of figure-eight orbits winding around the two given mass points. For many such orbits the elliptic function modulusk is small and can serve as a small parameter.Fourier expansions in terms of a parameter related tot are obtained for the intermediate orbit functions which provide representations in terms of elementary functions.     

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.     

Laplace's problem in mathematical aspect

In this paper we consider the mathematical aspect of Laplace's problem pertaining to the occurrence probability of elliptical and hyperbolical orbits for comets. We show that, among other things, if we use arbitrary velocity distributions in Laplace's problem, we may obtain an arbitrarily small probability for elliptical orbits and a probability neighbouring around one for hyperbolic orbits, or conversely. Our theorems hold, even for the variable initial conditions. The aspect presented here permits us to describe easily studies of Laplace's problem.Dedicated to my teacher, Professor W. Orlicz, on the occasion of 80th birthday.     

Periodic orbits in the gravity field of a fixed homogeneous cube

In the current study, the existence of periodic orbits around a fixed homogeneous cube is investigated, and the results have powerful implications for examining periodic orbits around non-spherical celestial bodies. In the two different types of symmetry planes of the fixed cube, periodic orbits are obtained using the method of the Poincaré surface of section. While in general positions, periodic orbits are found by the homotopy method. The results show that periodic orbits exist extensively in symmetry planes of the fixed cube, and also exist near asymmetry planes that contain the regular Hex cross section. The stability of these periodic orbits is determined on the basis of the eigenvalues of the monodromy matrix. This paper proves that the homotopy method is effective to find periodic orbits in the gravity field of the cube, which provides a new thought of searching for periodic orbits around non-spherical celestial bodies. The investigation of orbits around the cube could be considered as the first step of the complicated cases, and helps to understand the dynamics of orbits around bodies with complicated shapes. The work is an extension of the previous research work about the dynamics of orbits around some simple shaped bodies, including a straight segment, a circular ring, an annulus disk, and simple planar plates.     

High-order Analytical Solutions of Hill’s Equations

The purpose of this paper is the semi-analytical computation of the bounded orbits of Hill’s equations, describing the relative motion of two particles in their Keplerian motion around a central body. We have considered the case in which one of the particles moves along a circular reference orbit. The solutions obtained are the generalisation of the usual periodic orbits obtained for the linearised equations and also of the third-order solution computed by D.L. Richardson and J.W. Mitchell (2003). With the algorithm presented, those orbits can be computed in a fast and efficient way up to a high-order in the in-plane and out-of-plane amplitudes.     

Zonal harmonics of azimuthally averaged potential of rotating inhomogeneous ellipsoids

The distant retrograde orbits (DROs) can serve as the parking orbits for a long-term cis-lunar space station. This paper gives a comprehensive study on the transfer problem from DROs to Earth orbits, including low Earth orbits (LEOs), medium Earth orbits (MEOs), and geosynchronous orbits (GSOs), in the bicircular restricted four-body problem (BR4BP) via optimizations within a large solution space. The planar transfer problem is firstly solved by grid search and optimization techniques, and two types of transfer orbits, direct ones and low-energy ones, are both constructed. Then, the nonplanar transfer problem to Earth orbits with inclinations between 0 and 90 degrees are solved via sequential optimizations based on the planar transfers. The transfer characteristics in the cases of different destination orbit inclinations are discussed for both the direct and the low-energy transfer orbits. The important role of the lunar gravity in the low-energy transfers is also discussed, which can overcome the increase of transfer cost caused by the high inclination of Earth orbits. The distinct features of different transfer scenarios, including multiple revolutions around the Earth and Moon, the exterior phase, and the lunar flyby, are discovered. The energy of transfer orbits is exploited to discuss the effects of close lunar flybys. The results will be helpful for the transfer design in future manned or unmanned return missions, and can also provide valuable information for selecting proper parking DROs for cis-lunar space stations.

     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.     

月球卫星轨道力学综述

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

Periodic solutions around the collinear langrangian points in the photo gravitational restricted three-body problem: Sun-Jupiter case

The evolution of the periodic orbits around the collinear equilibrium positions, belonging to the Strömgren families a, b and c, with the radiation pressure parameter of the more massive body is studied in the Sun-Jupiter system. These families are determined for a single value of the radiation pressure parameter and particularly when the radiation force of the more massive body is equal to one half of the gravitational attraction. Then the critical stability orbits of each family are transferred with the radiation parameter. The stability of each periodic solution is also studied.     

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.     

Asymptotic solution for the two-body problem with constant tangential thrust acceleration

An analytical solution of the two body problem perturbed by a constant tangential acceleration is derived with the aid of perturbation theory. The solution, which is valid for circular and elliptic orbits with generic eccentricity, describes the instantaneous time variation of all orbital elements. A comparison with high-accuracy numerical results shows that the analytical method can be effectively applied to multiple-revolution low-thrust orbit transfer around planets and in interplanetary space with negligible error.     

Analysis of periodic orbits in the Saturn-Titan system using the method of Poincare section surfaces

We explore the periodic orbits and the regions of quasi-periodic motion around both the primaries in the Saturn-Titan system in the framework of planar circular restricted three-body problem. The location, nature and size of periodic and quasi-periodic orbits are studied using the numerical technique of Poincare surface of sections. The maximum amplitude of oscillations about the periodic orbits is determined and is used as a parameter to measure the degree of stability in the phase space for such orbits. It is found that the orbits around Saturn remain around it and their stability increases with the increase in the value of Jacobi constant C. The orbits around Titan move towards it with the increase in C. At C=3.1, the pericenter and apocenter are 358.2 and 358.5 km, respectively. No periodic or quasi-periodic orbits could be found by the present method around the collinear Lagrangian point L ₁ (0.9569373834…).     

Generalized Averaging Principle and the Secular Evolution of Planet Crossing Orbits

Planet crossing orbits give rise to mathematical singularities that make it not possible to apply the classical averaging principle to study the qualitative evolution of Near Earth Asteroids (NEAs). Recently this principle has been generalized to deal with crossings in a mathematical model with the planets on circular coplanar orbits. More accuracy is needed to compute the averaged evolution of planet crossing orbits for different purposes: computing reliable crossing times for the averaged motion, writing more precise proper elements and frequencies for NEAs, etc. In this paper we present the generalization of the averaging principle using a model where the eccentricity and the inclination of the planets are taken into account.     

Accurate Mars Express orbits to improve the determination of the mass and ephemeris of the Martian moons

The determination of the ephemeris of the Martian moons has benefited from observations of their plane-of-sky positions derived from images taken by cameras onboard spacecraft orbiting Mars. Images obtained by the Super Resolution Camera (SRC) onboard Mars Express (MEX) have been used to derive moon positions relative to Mars on the basis of a fit of a complete dynamical model of their motion around Mars. Since, these positions are computed from the relative position of the spacecraft when the images are taken, those positions need to be known as accurately as possible. An accurate MEX orbit is obtained by fitting two years of tracking data of the Mars Express Radio Science (MaRS) experiment onboard MEX. The average accuracy of the orbits has been estimated to be around 20–25 m. From these orbits, we have re-derived the positions of Phobos and Deimos at the epoch of the SRC observations and compared them with the positions derived by using the MEX orbits provided by the ESOC navigation team. After fit of the orbital model of Phobos and Deimos, the gain in precision in the Phobos position is roughly 30 m, corresponding to the estimated gain of accuracy of the MEX orbits. A new solution of the GM of the Martian moons has also been obtained from the accurate MEX orbits, which is consistent with previous solutions and, for Phobos, is more precise than the solution from the Mars Global Surveyor (MGS) and Mars Odyssey (ODY) tracking data. It will be further improved with data from MEX-Phobos closer encounters (at a distance less than 300 km). This study also demonstrates the advantage of combining observations of the moon positions from a spacecraft and from the Earth to assess the real accuracy of the spacecraft orbit. In turn, the natural satellite ephemerides can be improved and participate to a better knowledge of the origin and evolution of the Martian moons.     

Using space manifold dynamics to deploy a small satellite constellation around the Moon

The possibility of communicating with the far side of the Moon is essential for keeping a continuous radio link with lunar orbiting spacecraft, as well as with manned or unmanned surface facilities in locations characterized by poor coverage from Earth. If the exploration and the exploitation of the Moon must be sustainable in the medium/long term, we need to develop the capability to realize and service such facilities at an affordable cost. Minimizing the spacecraft mass and the number of launches is a driving parameter to this end. The aim of this study is to show how Space Manifold Dynamics can be profitably applied in order to launch three small spacecraft onboard the same launch vehicle and send them to different orbits around the Moon with no significant difference in the Delta-V budgets. Internal manifold transfers are considered to minimize also the transfer time. The approach used is the following: we used the linearized solution of the equations of motion in the Circular Restricted Three Body Problem to determine a first–guess state vector associated with the Weak Stability Boundary regions, either around the collinear Lagrangian point L1 or around the Moon. The resulting vector is then used as initial state in a numerical backward-integration sequence that outputs a trajectory on a manifold. The dynamical model used in the numerical integration is four-body and non-circular, i.e. the perturbations of the Sun and the lunar orbital eccentricity are accounted for. The trajectory found in this way is used as the principal segment of the lunar transfer. After separation, with minor maneuvers each satellite is injected into different orbits that lead to ballistic capture around the Moon. Finally, one or more circularization maneuvers are needed in order to achieve the final circular orbits. The whole mission profile, from launch to insertion into the final lunar orbits, is modeled numerically with the commercial software STK.     

Periodic orbits in the Chermnykh-like restricted problem of oblate bodies with radiation

In this paper, the periodic orbits around triangular points in the range of linear stability of the restricted three body problem, when the smaller primary and the test particle have the shape of an oblate spheroid and the larger primary is a radiation emitter with the allowance for the gravitational potential from the belt, is studied. It is observed that the orbits around these points are elliptical and have long and short periodic orbits. The period, orientation, eccentricities, the semi-major and semi-minor axis of the elliptic orbits are found. The study includes some numerical examples in the case of the Sun-Earth and Sun-Jupiter systems.     

Oblateness effect of Saturn on periodic orbits in the Saturn-Titan restricted three-body problem

We explore the effect of oblateness of Saturn (more massive primary) on the periodic orbits and the regions of quasi-periodic motion around both the primaries in the Saturn-Titan system in the framework of planar circular restricted three-body problem. First order interior and exterior mean motion resonances are located. The effect of oblateness is studied on the location, nature and size of periodic and quasi-periodic orbits, using the numerical technique of Poincare surface of sections. Some of the periodic orbits change to quasi-periodic orbits due to the effect of oblateness and vice-versa. The stability of the orbits around Saturn, Titan and both varies with the inclusion of oblateness. The centers of the periodic orbits around Titan move towards Saturn, whereas those around Saturn move towards Titan. For the orbit around Titan at C=2.9992, x=0.959494, the apocenter becomes pericenter. By incorporating oblateness effect, the orbit around Titan at C=2.99345, x=0.924938 is captured by Saturn, remains in various trajectories around Saturn, and as time progresses it spirals away around both the primaries.     

A Hamiltonian approach to the planar optimization of mid-course corrections

Lawden’s primer vector theory gives a set of necessary conditions that characterize the optimality of a transfer orbit, defined accordingly to the possibility of adding mid-course corrections. In this paper a novel approach is proposed where, through a polar coordinates transformation, the primer vector components decouple. Furthermore, the case when transfer, departure and arrival orbits are coplanar is analyzed using a Hamiltonian approach. This procedure leads to approximate analytic solutions for the in-plane components of the primer vector. Moreover, the solution for the circular transfer case is proven to be the Hill’s solution. The novel procedure reduces the mathematical and computational complexity of the original case study. It is shown that the primer vector is independent of the semi-major axis of the transfer orbit. The case with a fixed transfer trajectory and variable initial and final thrust impulses is studied. The acquired related optimality maps are presented and analyzed and they express the likelihood of a set of trajectories to be optimal. Furthermore, it is presented which kind of requirements have to be fulfilled by a set of departure and arrival orbits to have the same profile of primer vector.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

"Quasi Satellite Orbits" to observe a possible small moon of Pallas

The purpose of this paper is to make a numerical search for natural orbits that can be used for a spacecraft to study a possible small moon of Pallas. There are many speculations about the existence of a small companion around this large asteroid, so finding and classifying orbits around this possible celestial body is an interesting problem in astrodynamics and that can be used for a spacecraft to observe this body. It is assumed that this moon has a radius that can vary from 0.125 to 1 km and that is located 750 or 500 km away from the center of Pallas. The idea is to show the effects of this parameter in the orbits around this moon. It means that the moon is much smaller than Pallas, so Keplerian orbits are not possible around it. To solve this problem, it is possible to find some special orbits that are called "Quasi Satellite Orbits" (QSO). They are orbits dominated by the gravity of Pallas, but that use the smaller perturbation from the moon to keep the spacecraft close to it. The present work searches for orbits that make the spacecraft to remain at given limits in its distance from the moon, like in the range from 3 to 50 km, the values used as an example in the present paper. This value is used because it is a good range to observe the body without getting to close to it, so reducing the risks of collisions. Each trajectory can be identified by the initial conditions of the spacecraft with respect to the moon, which means its initial position and velocity. The dynamics considers the restricted three-body problem and the influence of the solar radiation pressure, because some spacecraft may have higher values for the area-to-mass ratio, which gives a non-negligible effect in the trajectory of the spacecraft.