Xiaoyu Fu

We present the solution approach developed by the team “TheAntipodes” during the 12th edition of the Global Trajectory Optimization Competition (GTOC12). An overview of the approach is as follows: (1) generate asteroid subsets, (2) chain building with beam search, (3) convex low-thrust trajectory optimization, (4) manual refinement of rendezvous times, and (5) optimal solution set selection. The generation of asteroid subsets involves a heuristic process to find sets of asteroids that are likely to permit high-scoring asteroid chains. Asteroid sequences “chains” are built within each subset through a beam search based on Lambert transfers. Low-thrust trajectory optimization involves the use of sequential convex programming (SCP), where a specialized formulation finds the mass-optimal control for each ship’s trajectory within seconds. Once a feasible trajectory has been found, the rendezvous times are manually refined with the aid of the control profile from the optimal solution. Each ship’s individual solution is then placed into a pool where the feasible set that maximizes the final score is extracted using a genetic algorithm. Our final submitted solution placed fifth with a score of 15,489.

The dynamical environment around the asteroid (162173) Ryugu is analyzed in detail using a constant-density polyhedron model based on the measurements from the Hayabusa2 mission. Six exterior equilibrium points (EPs) are identified along the ridge line of Ryugu, and their topological classifications fall into two distinctive categories. The initial periodic orbit (PO) families are computed and analyzed, including distant retrograde/prograde orbit (DRO/DPO) families and fifteen PO families emanating from the exterior EPs. The fifteen PO families are further divided into three categories: seven converge to an EP, seven reach Ryugu’s surface, and one exhibits cyclic behavior during its progression. The existence of initial PO families converging to an EP is analyzed using the bifurcation of a degenerate EP. Connection between these families and similar ones in the circular restricted three-body problem (CRTBP) is also examined. Bifurcated PO families are identified and computed from the initial PO families, including ten families from the DROs, fifteen from the DPOs, and twenty-five associated with the EPs. The bifurcated families are separately analyzed and categorized in terms of their corresponding initial families. Connections established by the same bifurcation points between different bifurcated families are identified. A comparison is made for the dynamical environments of Ryugu and Bennu to evaluate the similarities and differences in the evolution of EPs and the progression of PO families in top-shaped asteroids.