Three methods for the computation of the probability of collision between two space objects are presented. These methods are based on the high order Taylor expansion of the time of closest approach (TCA) and distance of closest approach (DCA) of the two orbiting objects with respect to their initial conditions. The identification of close approaches is first addressed using the nominal objects states. When a close approach is identified, the dependence of the TCA and DCA on the uncertainties in the initial states is efficiently computed with differential algebra (DA) techniques. In the first method the collision probability is estimated via fast DA-based Monte Carlo simulation, in which, for each pair of virtual objects, the DCA is obtained via the fast evaluation of its Taylor expansion. The second and the third methods are the DA version of Line Sampling and Subset Simulation algorithms, respectively. These are introduced to further improve the efficiency and accuracy of Monte Carlo collision probability computation, in particular for cases of very low collision probabilities. The performances of the methods are assessed on orbital conjunctions occurring in different orbital regimes and dynamical models. The probabilities obtained and the associated computational times are compared against standard (i.e. not DA-based) version of the algorithms and analytical methods. The dependence of the collision probability on the initial orbital state covariance is investigated as well.

A high order method to quickly assess the effect that uncertainties produce on orbital conjunctions through a numerical high-fidelity propagator is presented. In particular, the dependency of time and distance of closest approach to initial uncertainties on position and velocity of both objects involved in a conjunction is studied. The approach relies on a numerical integration based on differential algebraic techniques and a high-order algorithm that expands the time and distance of closest approach in Taylor series with respect to relevant uncertainties. The modeled perturbations are atmospheric drag, using NRLMSISE-00 air density model, solar radiation pressure with shadow, third body perturbation using JPL's DE405 ephemeris, and EGM2008 gravity model. The polynomial approximation of the final position is used as an input to compute analytically the expansion of time and distance of closest approach. As a result, the analysis of a close encounter can be performed through fast, multiple evaluations of Taylor polynomials. Test cases with objects ranging from LEO to GEO regimes are considered to assess the performances and the accuracy of the proposed method.

No abstract

Abstract-As the amount of space debris orbiting the earth is continuously increasing, it is becoming essential to monitor and predict the debris' trajectories in order to avoid collision that could threaten space missions, i.e. operative satellites or manned spacecraft. An innovative bistatic radar system, that couples multi-beaming and ranging techniques, is proposed in this work as an instrument for the observation of Earth-orbiting objects. The main characteristics of the sensor are described. In addition, tailored algorithms have been implemented to perform orbit determination using the peculiar outputs provided by the sensor. The results of numerical simulations are illustrated to assess its performances for space debris detection and orbit determination.

No abstract