We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first establish a baseline for evaluating various techniques for close-encounter prediction and collision-probability calculation (Hoots et al. 1984, Gronchi 2005, JeongAhn and Malhotra 2015) by carrying out brute-force numerical simulations and using a sequence of filters to greatly reduce the computational expense of the algorithm. Next, we estimate conjunction events in the orbital environment following the anticipated deployments of the OneWeb LEO and SpaceX Starlink mega-constellations. As a final step, we investigate Minimum Space Occupancy (MiSO) orbits (Bombardelli et al. 2018), a generalization of the well-known frozen orbits that account for the perturbed-Keplerian dynamics of the Earth-Moon-Sun-satellite system. We evaluate the ability of MiSO configurations of the proposed mega-constellations, as suggested by Bombardelli et al. 2018, to reduce the risk of endogenous (intra-constellation) collisions. The results indicate that the adoption of the MiSO orbital configuration can significantly reduce risk with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.

An approximate analytical solution for the two body problem perturbed by a radial, low acceleration is obtained, using a regularized formulation of the orbital motion and the method of multiple scales. Formulating the dynamics with the Dromo special perturbation method allows us to separate the two characteristic periods of the problem in a clear and physically significative way, namely the orbital period and a period depending on the magnitude of the perturbing acceleration. This second period becomes very large compared to the orbital one for low thrust cases, allowing us to develop an accurate approximate analytical solution through the method of multiple scales. Compared to a regular expansion, the multiple scales solution retains the qualitative contributions of both characteristic periods and has a longer validity range in time. Looking at previous solutions for this problem, our approach has the advantage of not requiring the evaluation of special functions or an initially circular orbit. Furthermore, the simple expression reached for the long period provides additional insight on the problem. Finally, the behavior of the asymptotic solution is assessed through several test cases, finding a good agreement with high-precision numerical solutions. The results presented not only advance in the study of the two body problem with constant radial thrust, but confirm the utility of the method of multiple scales for tackling problems in orbital mechanics.

This paper is concerned with the comparison of semi-analytical and non-averaged propagation methods for Earth satellite orbits. We analyse the total integration error for semi-analytical methods and propose a novel decomposition into dynamical, model truncation, short-periodic, and numerical error components. The first three are attributable to distinct approximations required by the method of averaging, which fundamentally limit the attainable accuracy. In contrast, numerical error, the only component present in non-averaged methods, can be significantly mitigated by employing adaptive numerical algorithms and regularized formulations of the equations of motion. We present a collection of non-averaged methods based on the integration of existing regularized formulations of the equations of motion through an adaptive solver. We implemented the collection in the orbit propagation code THALASSA, which we make publicly available, and we compared the non-averaged methods to the semi-analytical method implemented in the orbit propagation tool STELA through numerical tests involving long-term propagations (on the order of decades) of LEO, GTO, and high-altitude HEO orbits. For the test cases considered, regularized non-averaged methods were found to be up to two times slower than semi-analytical for the LEO orbit, to have comparable speed for the GTO, and to be ten times as fast for the HEO (for the same accuracy). We show for the first time that efficient implementations of non-averaged regularized formulations of the equations of motion, and especially of non-singular element methods, are attractive candidates for the long-term study of high-altitude and highly elliptical Earth satellite orbits.

A realistic mission scenario for the deflection of fictitious asteroid 2015PDC is investigated that makes use of the ion beam shepherd concept as primary deflection technique. The article deals with the design of a low thrust rendezvous trajectory to the asteroid, the estimation of the propagated covariance ellipsoid and the outcome of a slow-push deflection starting from three worst case scenarios (impacts in New Delhi, Dhaka and Teheran). Displacing the impact point towards very low populated areas, as opposed to full deflection, is found to be the simplest and most effective mitigation approach. Mission design, technical and political aspects are discussed.