Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

On the periodic orbits and the integrability of the regularized Hill lunar problem

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune-Triton system

The problem of a spacecraft orbiting the Neptune-Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune's satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a x I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton's orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L-2, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincare section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method. Published by Elsevier Ltd. on behalf of COSPAR.

Distant stable direct orbits around the Moon

In the present work we explore regions of distant direct stable orbits around the Moon. First, the location and size of apparently stable regions are searched for numerically, adopting the approach of temporary capture time presented in Vieira Neto & Winter (2001). The study is made in the framework of the planar, circular, restricted three-body problem, Earth-Moon-particle. Regions of the initial condition space whose trajectories are apparently stable are determined. The criterion adopted was that the trajectories do not escape from the Moon during an integration period of 10(4) days. Using Poincare surface of sections the reason for the existence of the two stable regions found is studied. The stability of such regions proved to be due to two families of simple periodic orbits, h1 and h2, and the associated quasi-periodic orbits that oscillate around them. The robustness of the stability of the larger region, h2, is tested with the inclusion of the solar perturbation. The size of the region decreases, but it is still significant in size and can be useful in spacecraft missions.

The Liapunov exponent as a tool for exploring phase space

We have used the Liapunov exponent to explore the phase space of a dynamical system. Considering the planar, circular restricted three-body problem for a mass ratio mu = 10(-3) (close to the Jupiter/Sun case), we have integrated similar to 16,000 starting conditions for orbits started interior to that of the perturber and we have estimated the maximum Liapunov characteristic exponent for each starting condition. Despite the fact that the integrations, in general, are for only a few thousand orbital periods of the secondary, a comparative analysis of the Liapunov exponents for various values of the 'cut-off' gives a good overview of the structure of the phase space. It provides information about the diffusion rates of the various chaotic regions, the location of the regular regions associated with primary resonances and even details such as the location of secondary resonances that produce chaotic regions inside the regular regions of primary resonances.

Resonance and chaos .1. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted three-body problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincare surfaces of section with a mass ratio of 10(-3) (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. on theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Irregular Satellites of Jupiter: a study of the capture direction

The irregular satellites of Jupiter are believed to be captured asteroids or planetesimals. In the present work is studied the direction of capture of these objects as a function of their orbital inclination. We performed numerical simulations of the restricted three-body problem, Sun-Jupiter-particle, taking into account the growth of Jupiter. The integration was made backward in time. Initially, the particles have orbits as satellites of Jupiter, which has its present mass. Then, the system evolved with Jupiter losing mass and the satellites escaping from the planet. The reverse of the escape direction corresponds to the capture direction. The results show that the Lagrangian points L1 and L2 mainly guide the direction of capture. Prograde satellites are captured through these two gates with very narrow amplitude angles. In the case of retrograde satellites, these two gates are wider. The capture region increases as the orbital inclination increases. In the case of planar retrograde satellites the directions of capture cover the whole 360 degrees around Jupiter. We also verified that prograde satellites are captured earlier in actual time than retrograde ones.

Controlling the Eccentricity of Polar Lunar Orbits with Low-Thrust Propulsion

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Resonance and chaos: II. Exterior resonances and asymmetric libration

The motion of a test particle in the vicinity of exterior resonances is examined in the context of the planar, circular, restricted three-body problem. The existence of asymmetric periodic orbits associated with the 1 : n resonances (where n = 2, 3, 4, 5) is confirmed; there is also evidence of asymmetric resonances associated with larger values of n. A detailed examination of the evolution of the family of orbits associated with the 1:2 resonance shows the sequence that leads to asymmetric libration. On the basis of numerical studies of the phase space it is concluded that the existence of asymmetric libration means that the region exterior to the perturbing mass is more chaotic than the interior region. The apparent absence of 'particles' in 1 : n resonances in the solar system may reflect this inherent bias.

Application of the GPS system for preliminary satellite orbit determination

In this paper, we discuss a method of preliminary orbit determination for an artificial satellite based on the navigation message of the GPS constellation. Orbital elements are considered as state variables and a simple dynamic model, based on the classic two-body problem, is used. The observations are formed by range and range and range-rate with respect to four visible GPS. A discrete Kalman filter with simulated data is used as filtering technique. The data are obtained through numerical propagation (Cowell's method), which considers special perturbations for the GPS satellite constellation and a user satellite. © 1997 COSPAR. Published by Elsevier Science Ltd.

The stability evolution of a family of simply periodic lunar orbits

The present work deals with a family of simply periodic orbits around the Moon in the rotating Earth Moon-particle system. Taking the framework of the planar, circular, restricted three-body problem, we follow the evolution of this family of periodic orbits using the numerical technique of Poincaré surface of section. The maximum amplitude of oscillation about the periodic orbits are determined and can be used as a parameter to measure the degree of stability in the phase space for such orbits. Despite the fact that the whole family of periodic orbits remain stable, there is a dichotomy in the quasi-periodic ones at the Jacobi constant Cj = 2.85. The quasi-periodic orbits with Cj < 2.85 oscillate around the periodic orbits in a different way from those with Cj > 2.85. © 1999 Elsevier Science Ltd. All rights reserved.

Interplanetary natural routes via Moon and Lagrangiam points L1 and L2

Swing-by techniques are extensively used in interplanetary missions to minimize fuel consumption and to raise payloads of spaceships. The effectiveness of this type of maneuver has been proven since the beginning of space exploration. According to this premise, we have explored the existence of a natural and direct links between low Earth orbits and the lunar sphere of influence, to obtain low-energy interplanetary trajectories through swing-bys with the Moon and the Earth. The existence of these links are related to a family of retrograde periodic orbits around the Lagrangian equilibrium point L1 predicted for the circular, planar, restricted three-body Earth-Moon-particle problem. The trajectories in these links are sensitive to small disturbances. This enables them to be conveniently diverted reducing so the cost of the swing-by maneuver. These maneuvers allow us a gain in energy sufficient for the trajectories to escape from the Earth-Moon system and to stabilize in heliocentric orbits between the Earth and Venus or Earth and Mars. On the other hand, still within the Earth sphere of influence, and taking advantage of the sensitivity of the trajectories, is possible to design other swing-bys with the Earth or Moon. This allows the trajectories to have larger reach, until they can reach the orbit of other planets as Venus and Mars.(3σ)Broucke, R.A., Periodic Orbits in the Restricted Three-Body Problem with Earth-Moon Masses, JPL Technical Report 32-1168, 1968.

Optimal low-cost transfer to L4 and L5 lagrangian points based on G-family of periodic orbits

An alternative transfer strategy to send spacecrafts to stable orbits around the Lagrangian equilibrium points L4 and L5 based in trajectories derived from the periodic orbits around LI is presented in this work. The trajectories derived, called Trajectories G, are described and studied in terms of the initial generation requirements and their energy variations relative to the Earth through the passage by the lunar sphere of influence. Missions for insertion of spacecrafts in elliptic orbits around L4 and L5 are analysed considering the Restricted Three-Body Problem Earth- Moon-particle and the results are discussed starting from the thrust, time of flight and energy variation relative to the Earth. Copyright© (2012) by the International Astronautical Federation.

Nebular gas drag and planetary accretion with eccentric high-mass planets

Aims.We investigate the dynamics of pebbles immersed in a gas disk interacting with a planet on an eccentric orbit. The model has a prescribed gap in the disk around the location of the planetary orbit, as is expected for a giant planet with a mass in the range of 0.1-1 Jupiter masses. The pebbles with sizes in the range of 1 cm to 3 m are placed in a ring outside of the giant planet orbit at distances between 10 and 30 planetary Hill radii. The process of the accumulation of pebbles closer to the gap edge, its possible implication for the planetary accretion, and the importance of the mass and the eccentricity of the planet in this process are the motivations behind the present contribution. Methods. We used the Bulirsch-Stoer numerical algorithm, which is computationally consistent for close approaches, to integrate the Newtonian equations of the planar (2D), elliptical restricted three-body problem. The angular velocity of the gas disk was determined by the appropriate balance between the gravity, centrifugal, and pressure forces, such that it is sub-Keplerian in regions with a negative radial pressure gradient and super-Keplerian where the radial pressure gradient is positive. Results. The results show that there are no trappings in the 1:1 resonance around the L 4 and L5 Lagrangian points for very low planetary eccentricities (e2 < 0.07). The trappings in exterior resonances, in the majority of cases, are because the angular velocity of the disk is super-Keplerian in the gap disk outside of the planetary orbit and because the inward drift is stopped. Furthermore, the semi-major axis location of such trappings depends on the gas pressure profile of the gap (depth) and is a = 1.2 for a planet of 1 MJ. A planet on an eccentric orbit interacts with the pebble layer formed by these resonances. Collisions occur and become important for planetary eccentricity near the present value of Jupiter (e 2 = 0.05). The maximum rate of the collisions onto a planet of 0.1 MJ occurs when the pebble size is 37.5 cm ≤ s < 75 cm; for a planet with the mass of Jupiter, it is15 cm ≤ s < 30 cm. The accretion stops when the pebble size is less than 2 cm and the gas drag dominates the motion. © 2013 ESO.

Stable regions around Pluto

In a previous work, GiuliattiWinter et al. found several stable regions for test particles in orbit around Pluto associated with families of periodic orbits obtained in the circular, restricted three-body problem. They have shown that a possible eccentricity of the Pluto-Charon binary slightly reduces but does not destroy any of these stable regions. In thiswork, we extended their results by analysing the cases with the orbital inclination (I) equal to zero and considering the argument of pericentre (w) equal to 90°, 180° and 270°. We explore the influence of the orbital inclination of the particles in these stable regions. In this case, the initial inclination varies from 10° to 170° in steps of 10°. We also present a sample of results for the longitude of the ascending node Ω = 90°, considering the cases I = 20°, 50°, 130° and 180°. Our results show that stable regions are present in all of the inclined cases, except when the initial inclination of the particles is equal to 110°. A sample of 3D trajectories of quasi-periodic orbits were found related to the periodic orbits obtained in the planar case by Giuliatti Winter et al. © 2013 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.