The Tisserand graph (TG) is a graphical tool commonly employed in the preliminary design of gravity-assisted trajectories. The TG is a two-dimensional map showing essential orbital information regarding the Keplerian orbits resulting from the close passage by one or more massive bodies, given the magnitude of the hyperbolic excess speed ($v_{\infty}$) and the minimum allowed pericenter height for each passage. Contours of constant $v_{\infty}$ populate the TG. Intersections between contours allow to link consecutive flybys and build sequences of encounters en route to a selected destination. When the number of perturbing bodies is large and many $v_{\infty}$ levels are considered, the identification of all the possible sequences of encounters through the visual inspection of the TG becomes a laborious task. Besides, if the sequences are used as input for a numerical code for trajectory design and optimization, an automated examination of the TG is desirable. This contribution describes an automatic technique to explore the TG and find all the encounter paths. The technique is based on a tree search method, and the intersections between contours are found using the regula-falsi scheme. The method is validated through comparisons with solutions available in the open literature. Examples are given of application to interplanetary mission scenarios, including the coupling with a trajectory optimizer.

The icy moons are in the focus of the exploration plans of the leading space agencies because of the indications of water-based life and geological activity observed in a number of these objects. In particular, the presence of geyser-like jets of water near Enceladus' south pole has turned this moon of Saturn into a priority candidate to search for life and habitability features. This investigation proposes a set of trajectories between Halo orbits about Lagrangian points L1 and L2 in the Saturn-Enceladus Circular Restricted Three-Body Problem as science orbits for a future in situ mission at Enceladus. The design methodology is presented, followed by the analysis of the observational performance of the solutions. The conclusion is that the proposed orbits exhibit suitable features for their use in the scientific exploration of Enceladus, i.e., long transfer times, low altitudes, wide surface visibility windows and long times of overflight.

This contribution investigates the properties of a category of orbits around Enceladus. The motivation is the interest in the in situ exploration of this moon following the detection on behalf of Cassini of plumes of water and organic compounds close to its south pole. In a previous investigation, a set of heteroclinic transfers were designed between Halo orbits around the equilibrium points L1 and L2 of the circular restricted three-body problem with Saturn and Enceladus as primaries. The kinematical and geometrical characteristics of those trajectories makes them good candidates as science orbits for the extended observation of the surface of Enceladus: they are highly inclined, they approach the moon and they are maneuver free. However, the low heights above the surface and the strong perturbing effect of Saturn impose a more careful look at their dynamics, in particular regarding the influence of the polar flattening of the primaries. Therefore, those solutions are here reconsidered by employing a dynamical model that includes the effect of the oblateness of Saturn and Enceladus, individually and in combination. Substitutes of the Halo orbits around the equilibrium points L1 and L2 and their stable and unstable hyperbolic invariant manifolds are obtained in the perturbed models, and maneuver-free heteroclinic transfers are identified in the new framework. A systematic comparison with the corresponding solutions of the unperturbed problem shows that qualitative and quantitative features are not significantly altered when the oblateness of the primaries is taken into account, and that J2 of Saturn plays a larger role than the oblateness of Enceladus. From a mission perspective, the results confirm the scientific value of the solutions obtained in the classical circular restricted three-body problem and suggests that this simpler model can be used in a preliminary feasibility analysis.

Orbit insertion at Saturn requires a large impulsive manoeuver due to the velocity difference between the spacecraft and the planet. This paper presents a strategy to reduce dramatically the hyperbolic excess speed at Saturn by means of deep-space electric propulsion. The interplanetary trajectory includes a gravity assist at Jupiter, combined with low-thrust maneuvers. The thrust arc from Earth to Jupiter lowers the launch energy requirement, while an ad hoc steering law applied after the Jupiter flyby reduces the hyperbolic excess speed upon arrival at Saturn. This lowers the orbit insertion impulse to the point where capture is possible even with a gravity assist with Titan. The control-law algorithm, the benefits to the mass budget and the main technological aspects are presented and discussed. The simple steering law is compared with a trajectory optimizer to evaluate the quality of the results and possibilities for improvement.

The description of the long-term dynamics of highly elliptic orbits under third-body perturbations may require an expansion of the disturbing function in series of the semi-major axes ratio up to higher orders. To avoid dealing with long series in trigonometric functions, we refer the motion to the apsidal frame and efficiently remove the short-period effects of this expansion in vectorial form up to an arbitrary order. We then provide the variation equations of the two fundamental vectors of the Keplerian motion by analogous vectorial recurrences, which are free from singularities and take a compact form useful for the numerical propagation of the flow in mean elements.

The minimum orbital intersection distance is used as a measure to assess potential close approaches and collision risks between astronomical objects. Methods to calculate this quantity have been proposed in several previous publications. The most frequent case is that in which both objects have elliptical osculating orbits. When at least one of the two orbits has low eccentricity, the latter can be used as a small parameter in an asymptotic power series expansion. The resulting approximation can be exploited to speed up the computation with negligible cost in terms of accuracy. This contribution introduces two asymptotic procedures into the SDG-MOID method presented in a previous article, it discusses the results of performance tests and their comparisons with previous findings. The best approximate procedure yields a reduction of 40% in computing speed without degrading the accuracy of the determinations. This remarkable result suggests that large benefits can be obtained in applications involving massive distance computations, such as in the analysis of large databases, in Monte Carlo simulations for impact risk assessment or in the long-time monitoring of the minimum orbital intersection distance between two objects.

We present an outline of a new reduction of the Hipparcos astrometric data, the justifications of which are described in the accompanying paper. The emphasis is on those aspects of the data analysis that are fundamentally different from the ones used for the catalogue published in 1997. The new reduction uses a dynamical modelling of the satellite's attitude. It incorporates provisions for scan-phase discontinuities and hits, most of which have now been identified. Solutions for the final along-scan attitude (the reconstruction of the satellite's scan phase), the abscissa corrections and the instrument model, originally solved simultaneously in the great-circle solution, are now de-coupled. This is made possible by starting the solution iterations with the astrometric data from the published catalogue. The de-coupling removes instabilities that affected great-circle solutions for short data sets in the published data. The modelling-noise reduction implies smaller systematic errors, which is reflected in a reduction of the abscissa-error correlations by about a factor 40. Special care is taken to ensure that measurements from both fields of view contribute significantly to the along-scan attitude solution. This improves the overall connectivity of the data and rigidity of the reconstructed sky, which is of critical importance to the reliability of the astrometric data. The changes in the reduction process and the improved understanding of the dynamics of the satellite result in considerable formal-error reductions for stars brighter than 8th magnitude.