Global search for planar and three-dimensional periodic orbits near Europa

A global grid search is performed to find axi- and doubly-symmetric periodic orbits in the restricted three-body problem using the dimensioned parameters associated with the Jupiter-Europa system. Local differential correctors are applied to regions of the initial condition phase space that appear to be near solutions. A three-dimensional initial condition mesh with billions of nodes is evaluated, and over 600,000 periodic solutions are identified. Families of direct and retrograde solutions, both new and previously published, are identified and discussed. Stability is analyzed for each solution and general regions of stability are noted. Of the most promising results is the observation and characterization of a large class of stable yet highly-inclined direct orbits. Finally, all of the solutions and associated properties are archived. The resulting database is a practical reference for preliminary design of missions to Europa.