We revisit the concept of sphere of gravitational activity, to which we give both a geometrical and physical meaning. This study aims to refine this concept in a much broader context that could, for instance, be applied to exo-planetary problems (in a Galactic stellar disc-Star-Planets system) to define a first order "border" of a planetary system. The methods used in this paper rely on classical Celestial Mechanics and develop the equations of motion in the framework of the 3-body problem (e.g. Star-Planet-Satellite System). We start with the basic definition of planet's sphere of activity as the region of space in which it is feasible to assume a planet as the central body and the Sun as the perturbing body when computing perturbations of the satellite's motion. We then investigate the geometrical properties and physical meaning of the ratios of Solar accelerations (central and perturbing) and planetary accelerations (central and perturbing), and the boundaries they define. We clearly distinguish throughout the paper between the sphere of activity, the Chebotarev sphere (a particular case of the sphere of activity), Laplace sphere, and the Hill sphere. The last two are often wrongfully thought to be one and the same. Furthermore, taking a closer look and comparing the ratio of the star's accelerations (central/perturbing) to that of the planetary acceleration (central/perturbing) as a function of the planeto-centric distance, we have identified different dynamical regimes which are presented in the semi-analytical analysis.

Modeling phenomena from experimental data, always begin with a \emph{choice of hypothesis} on the observed dynamics such as \emph{determinism}, \emph{randomness}, \emph{derivability} etc. Depending on these choices, different behaviors can be observed. The natural question associated to the modeling problem is the following : \emph{"With a finite set of data concerning a phenomenon, can we recover its underlying nature ?} From this problem, we introduce in this paper the definition of \emph{multi-scale functions}, \emph{scale calculus} and \emph{scale dynamics} based on the \emph{time-scale calculus} (see \cite{bohn}). These definitions will be illustrated on the \emph{multi-scale Okamoto's functions}. The introduced formalism explains why there exists different continuous models associated to an equation with different \emph{scale regimes} whereas the equation is \emph{scale invariant}. A typical example of such an equation, is the \emph{Euler-Lagrange equation} and particularly the \emph{Newton's equation} which will be discussed. Notably, we obtain a \emph{non-linear diffusion equation} via the \emph{scale Newton's equation} and also the \emph{non-linear Schr\"odinger equation} via the \emph{scale Newton's equation}. Under special assumptions, we recover the classical \emph{diffusion} equation and the \emph{Schr\"odinger equation}.

Experimental data suggest that the Earth short time dynamics is related to stochastic fluctuation of its shape. As a first approach to this problem, we derive a toy-model for the motion of a rotating ellipsoid in the framework of stochastic differential equations. Precisely, we assume that the fluctuations of the geometric flattening can be modeled by an admissible class of diffusion processes respecting some invariance properties. This model allows us to determine an explicit drift component in the dynamical flattening and the second zonal harmonic whose origin comes from the stochastic term and is responsible for short term effects. Using appropriate numerical scheme, we perform numerical simulations showing the role of the stochastic perturbation on the short term dynamics. Our toy-model exhibits behaviors which look like the experimental one. This suggests to extend our strategy with a more elaborated model for the deterministic part.

We study the Sharma-Parthasarathy stochastic two-body problem introduced by N. Sharma and H. Parthasarathy in \cite{sharma}. In particular, we focus on the preservation of some fundamental features of the classical two-body problem like the Hamiltonian structure and first integrals in the stochastic case. Numerical simulations are performed which illustrate the dynamical behaviour of the osculating elements as the semi-major axis, the eccentricity and the pericenter. We also derive a stochastic version of Gauss's equations in the planar case.