Divergence Evaluation Criteria for Lunar Departure Trajectories Under Bi-Circular Restricted Four-Body Problem

This study focuses on the nonlinear departure dynamics of spacecraft from the Near Rectilinear Halo Orbit (NRHO) to the outer regions of Selenocentric Space. By carefully selecting the combination of orbital parameters and the order of the evaluation process, it becomes possible to precisely identify the divergence moment and to reliably classify the subsequent dynamical space. An empirical divergence detection algorithm is proposed by integrating multiple parameters derived from multi-body dynamical models, including gravitational potentials and related quantities. In an applied analysis using this method, it is found that the majority of perturbed trajectories diverge into the outer Earth–Moon Vicinity, while transfers into the inner Earth–Moon Vicinity are relatively limited. Furthermore, transfers to Heliocentric Space are found to be dependent not on the magnitude of the initial perturbation but on the geometric configuration of the Sun, Earth, and Moon during the transfer phase. The investigation of the Sun’s initial phase reveals a rotationally symmetric structure in the perturbation distribution within the Sun–Earth–Moon system, as well as localized conditions under which the destination space varies significantly depending on the initial state. Identifying the divergence moment allows for comparative evaluation of the spacecraft’s nonlinear dynamical state, providing valuable insights for the development of safe and efficient transfer strategies from selenocentric orbits, including those originating from the NRHO.