In the last decades there have been an increasing interest in improving the accuracy of spacecraft navigation and trajectory data. In the course of this plan some anomalies have been found that cannot, in principle, be explained in the context of the most accurate orbital models including all known effects from classical dynamics and general relativity. Of particular interest for its puzzling nature, and the lack of any accepted explanation for the moment, is the flyby anomaly discovered in some spacecraft flybys of the Earth over the course of twenty years. This anomaly manifest itself as the impossibility of matching the pre and post-encounter Doppler tracking and ranging data within a single orbit but, on the contrary, a difference of a few mm$/$s in the asymptotic velocities is required to perform the fitting. Nevertheless, no dedicated missions have been carried out to elucidate the origin of this phenomenon with the objective either of revising our understanding of gravity or to improve the accuracy of spacecraft Doppler tracking by revealing a conventional origin. With the occasion of the Juno mission arrival at Jupiter and the close flybys of this planet, that are currently been performed, we have developed an orbital model suited to the time window close to the perijove. This model shows that an anomalous acceleration of a few mm$/$s$^2$ is also present in this case. The chance for overlooked conventional or possible unconventional explanations is discussed.

The flyby anomaly is a persistent riddle in astrodynamics. Orbital analysis in several flybys of the Earth since the Galileo spacecraft flyby of the Earth in 1990 have shown that the asymptotic post-encounter velocity exhibits a difference with the initial velocity that cannot be attributed to conventional effects. To elucidate its origin, we have developed an orbital program for analyzing the trajectory of the spacecraft in the vicinity of the perigee, including both the Sun and the Moon's tidal perturbations and the geopotential zonal, tesseral and sectorial harmonics provided by the EGM96 model. The magnitude and direction of the anomalous acceleration acting upon the spacecraft can be estimated from the orbital determination program by comparing with the trajectories fitted to telemetry data as provided by the mission teams. This acceleration amounts to a fraction of a mm/s$^2$ and decays very fast with altitude. The possibility of some new physics of gravity in the altitude range for spacecraft flybys is discussed.

The flyby anomaly is the unexpected variation of the asymptotic post-encounter velocity of a spacecraft with respect to the pre-encounter velocity as it performs a slingshot manoeuvre. This effect has been detected in, at least, six flybys of the Earth but it has not appeared in other recent flybys. In order to find a pattern in these, apparently contradictory, data several phenomenological formulas have been proposed but all have failed to predict a new result in agreement with the observations. In this paper we use a multivariate dimensional analysis approach to propose a fitting of the data in terms of the local parameters at perigee, as it would occur if this anomaly comes from an unknown fifth force with latitude dependence. Under this assumption, we estimate the range of this force around 300 km.

It is a well-known fact that the gravitational effect of dark matter in galaxies is only noticeable when the orbital accelerations drop below $a_0 \simeq 2\times 10^{-8}$ cm s$^{-1}$ (Milgrom's Law). This peculiarity of the dynamic behaviour of galaxies was initially ascribed to a modification of Newtonian dynamics (MOND theory) and, consequently, it was used as an argument to criticize the dark matter hypothesis. In our model, warm dark matter is composed by collisionless Vlasov particles with a primordial typical velocity $\simeq 330$ km s$^{-1}$ and, consequently, they evaporated from galactic cores and reorganized in halos with a cusp at a finite distance from the galactic center (in contrast with Cold Dark Matter simulations which predict a cusp at the center of galaxies). This is confirmed by mean-field N-body simulations of the self-gravitating Vlasov dark matter particles in the potential well of the baryonic core. The rest mass of these particles, $\mu$, is determined from a kinetic theory of the early universe with a cosmological constant. We find that $\mu$ is in the range of a few keV. This result makes sterile neutrinos the best suited candidates for the main component of dark matter.