A method for the rapid estimation of transfer costs for the removal of debris in low Earth orbit is proposed. Debris objects among a population with similar inclination values are considered. The proposed approximate analysis can provide estimations of actual Deltav between any debris object pair as a function of time; these estimations allow for the rapid evaluation of the costs of large sequences of targets to be removed. The effect of Earth's oblateness perturbation (J2) is exploited to reduce transfer costs. The debris removal problem of the 9th edition of the Global Trajectory Optimization Competition is used to evaluate the estimation accuracy; Deltav estimations of the transfers between objects pairs are verified by comparing them with the GTOC9 solution proposed by the winning team from JPL. The comparison of the results demonstrates the very good accuracy of the simple approximation. Key words: Space debris; approximation; trajectory optimization; J2 perturbation

This paper summarizes the solution approach and the numerical methods developed by the joint team Sapienza University of Rome and Politecnico di Torino (Team Sapienza-PoliTo) in the context of the 10th Global Trajectory Optimization Competition. The proposed method is based on a preliminary partition of the galaxy into several small zones of interest, where partial settlement trees are developed, in order to match a (theoretical) optimal star distribution. A multi-settler stochastic Beam Best-First Search, that exploits a guided multi-star multi-vessel transition logic, is proposed for solving a coverage problem, where the number of stars to capture and their distribution within a zone is assigned. The star-to-star transfers were then optimized through an indirect procedure. A number of refinements, involving settle time re-optimization, explosion, and pruning, were also investigated. The submitted 1013-star solution, as well as an enhanced 1200-point rework, are presented.