In this paper, we exploit the fact that the dynamics of homogeneous and isotropic Friedmann-Lemaitre universes is a special case of generalized Lotka-Volterra system where the competitive species are the barotropic fluids filling the Universe. Without coupling between those fluids, Lotka-Volterra formulation offers a pedagogical and simple way to interpret usual Friedmann-Lemaitre cosmological dynamics. A natural and physical coupling between cosmological fluids is proposed which preserve the structure of the dynamical equations. Using the standard tools of Lotka-Volterra dynamics, we obtain the general Lyapunov function of the system when one of the fluids is coupled to dark energy. This provides in a rigorous form a generic asymptotic behavior for cosmic expansion in presence of coupled species, beyond the standard de Sitter, Einstein-de Sitter and Milne cosmologies. Finally, we conjecture that chaos can appear for at least four interacting fluids.

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.

In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby called {\em global symplectic integrator}. In particular, the proposed method allows us to recover the correct orbits character with very large integration time steps, small energy losses and short CPU times. To illustrate the numerical performances of the global symplectic integrator we will apply it to two well-known and widely studied problems: the H\'enon-Heiles model and the restricted three-body problem.

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.