Gravity assist pump-down trajectories that leverage \(v_{\infty }\) through sequential planetary moon flybys are enabling scientific missions to previously intractable destinations such as Saturn’s Enceladus, Jupiter’s Europa, and Uranus’s Miranda. These missions typically require flybys with multiple moons in the system. With each small moon requiring on the order of ten flybys and dozens of revolutions to enable transfers to the next moon in the series, the number of trajectory paths to consider quickly expands. To combat this combinatorial increase of the search space, a new method is developed using a low-dimensioned directed graph of the pump-down tour space of a planetary moon. The graph and its discretized information on the moon’s resonant and non-resonant transfer space are used to efficiently evaluate same-body only pump-down transfer sequences for time of flight and maneuver costs. Promising same-body only transfer sequences are converted to trajectories in a patched conics model, the associated solutions are converted to integrated trajectories in a high-fidelity model, and a Pareto front is identified. The full process is automated and enables a thorough search of the design space, resulting in families of optimized high-fidelity same-body only pump-down flyby tour trajectories.
