NIMASTEP is a dedicated numerical software developed by us, which allows one to integrate the osculating motion (using cartesian coordinates) in a Newtonian approach of an object considered as a point-mass orbiting a homogeneous central body that rotates with a constant rate around its axis of smallest inertia. The code can be applied to objects such as particles, artificial or natural satellites or space debris. The central body can be either any terrestrial planet of the solar system, any dwarf-planet, or even an asteroid. In addition, very many perturbations can be taken into account, such as the combined third-body attraction of the Sun, the Moon, or the planets, the direct solar radiation pressure (with the central body shadow), the non-homogeneous gravitational field caused by the non-sphericity of the central body, and even some thrust forces. The simulations were performed using different integration algorithms. Two additional tools were integrated in the software package; the indicator of chaos MEGNO and the frequency analysis NAFF. NIMASTEP is designed in a flexible modular style and allows one to (de)select very many options without compromising the performance. It also allows one to easily add other possibilities of use. The code has been validated through several tests such as comparisons with numerical integrations made with other softwares or with semi-analytical and analytical studies. The various possibilities of NIMASTEP are described and explained and some tests of astrophysical interest are presented. At present, the code is proprietary but it will be released for use by the community in the near future. Information for contacting its authors and (in the near future) for obtaining the software are available on the web site http://www.fundp.ac.be/en/research/projects/page_view/10278201/

We study the dynamics of the 3-D three-body problem of a small body moving under the attractions of a star and a giant planet which orbits the star on a much wider and elliptic orbit. In particular, we focus on the influence of an eccentric orbit of the outer perturber on the dynamics of a small highly inclined inner body. Our analytical study of the secular perturbations relies on the classical octupole hamiltonian expansion (third-order theory in the ratio of the semi-major axes), as third-order terms are needed to consider the secular variations of the outer perturber and potential secular resonances between the arguments of the pericenter and/or longitudes of the node of both bodies. Short-period averaging and node reduction (Laplace plane) reduce the problem to two degrees of freedom. The four-dimensional dynamics is analyzed through representative planes which identify the main equilibria of the problem. As in the circular problem (i.e. perturber on a circular orbit), the "Kozai-bifurcated" equilibria play a major role in the dynamics of an inner body on quasi-circular orbit: its eccentricity variations are very limited for mutual inclination between the orbital planes smaller than ~40^{\deg}, while they become large and chaotic for higher mutual inclination. Particular attention is also given to a region around 35^{\deg} of mutual inclination, detected numerically by Funk et al. (2011) and consisting of long-time stable and particularly low eccentric orbits of the small body. Using a 12th-order Hamiltonian expansion in eccentricities and inclinations, in particular its action-angle formulation obtained by Lie transforms in Libert & Henrard (2008), we show that this region presents an equality of two fundamental frequencies and can be regarded as a secular resonance. Our results also apply to binary star systems where a planet is revolving around one of the two stars.

In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth's gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth's shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.

The aim of Dawn mission is the acquisition of data from orbits around two bodies, (4)Vesta and (1)Ceres, the two most massive asteroids. Due to the low thrust propulsion, Dawn will slowly cross and transit through ground-track resonances, where the perturbations on Dawn orbit may be significant. In this context, to safety go the Dawn mission from the approach orbit to the lowest science orbit, it is essential to know the properties of the crossed resonances. This paper analytically investigates the properties of the major ground-track resonances (1:1, 1:2, 2:3 and 3:2) appearing for Vesta orbiters: location of the equilibria, aperture of the resonances and period at the stable equilibria. We develop a general method using an averaged Hamiltonian formulation with a spherical harmonic approximation of the gravity field. If the values of the gravity field coefficient change, our method stays correct and applicable. We also discuss the effect of one uncertainty on the C20 and C22 coefficients on the properties of the 1:1 resonance. These results are checked by numerical tests. We determine that the increase of the eccentricity appearing in the 2:3 resonance is due to the C22 and S22 coefficients. Our method can be easily adapted to missions similar to Dawn because, contrarily to the numerical results, the analytical formalism stays the same and is valid for a wide range of physical parameters of the asteroid (namely the shape and the mass) as well as for different spacecraft orbits. Finally we numerically study the probability of the capture in resonance 1:1. Our paper reproduces, explains and supplements the results of Tricarico and Sykes (2010).

We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lie in the equatorial plane (Sun or Jupiter for example) using an Hamiltonian formalism. We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on average, to characterize their stability/unstability and to compute the periods of the equilibria. The proposed theory is general enough, to be applied to a wide range of probes around planet or natural planetary satellites. The BepiColombo mission is used to motivate our analysis and to provide specific numerical data to check our analytical results. Finally, we also bring to the light that the coefficient J_2 is able to protect against the increasing of the eccentricity due to the Kozai-Lidov effect.

In this paper we provide an extensive analysis of the global dynamics of high-area-to-mass ratios geosynchronous (GEO) space debris, applying a recent technique developed by Cincotta et al. (2000), Mean Exponential Growth factor of Nearby Orbits (MEGNO), which provides an efficient tool to investigate both regular and chaotic components of the phase space. We compute a stability atlas, for a large set of near-geosynchronous space debris by numerically computing the MEGNO indicator, to provide an accurate understanding of the location of stable and unstable orbits as well as the timescale of their exponential divergence in case of chaotic motion. The results improve the analysis presented in Breiter et al. (2005a) notably by considering the particular case of high-area-to-mass ratios space debris. The results indicate that chaotic orbits region can be highly relevant, especially for very high area-to-mass ratios. Then, we provide some numerical investigations and an analytical theory which lead to a detailed understanding of the resonance structures appearing in the phase space. These analyses bring to the fore a relevant class of secondary resonances on both sides of the well-known pendulum-like pattern of geostationary space debris, leading to complex dynamics of such objects.