Clustering Approach To Identifying Low Lunar Frozen Orbits In A High-Fidelity Model

Low lunar frozen orbits are of continued interest in the astrodynamics community for trajectory design and space domain awareness. This paper presents a data-driven approach to analyzing a wide variety of numerically-generated lunar trajectories in a 100 × 100 lunar gravity model with the point mass gravity of the Earth and Sun. First, clustering is used to extract a summary of these trajectories: within each cluster, trajectories possess a geometrically similar evolution of perilune but varying drift and lifetimes. Within some clusters, trajectories with a bounded perilune evolution are also identified to produce candidates for lunar frozen orbits of distinct geometries.