Bounded lunar relative orbits

Finding relative orbits for satellite formations flying around the Earth, which are long-term bounded under various perturbations, has been a vibrant field of study. However, much less attention has been given to detecting such orbits for circumlunar formation flying missions. As opposed to low-Earth orbits, in circumlunar missions the third-body effect is large, and the magnitude of the C ₂₂ sectorial harmonic has the same order as the J ₂ zonal harmonic. This renders the analysis of bounded relative orbits more challenging. In this paper, we detect a new family of long-term bounded lunar relative orbits, which can be used for lunar formation flying missions. The main idea is to find an analytical solution for the mean relative distance among the formation satellites, and use it to derive formation geometries that are resilient to the gravitational and third-body perturbations. Two invariant mean-distance conditions and one bounded mean-distance conditions are derived. Numerical simulations indicate that the newly-derived invariance conditions yield long-term bounded relative motion.