The preliminary mission design of spacecraft missions to asteroids often involves, in the early phases, the selection of candidate target asteroids. The final result of such an analysis is a list of asteroids, ranked with respect to the necessary propellant to be used, that the spacecraft could potentially reach. In this paper we investigate the sensitivity of the produced asteroids rank to the employed trajectory model in the specific case of a small low-thrust propelled spacecraft beginning its journey from the Sun-Earth $L_2$ Lagrangian point and heading to a rendezvous with some near-Earth asteroid. We consider five increasingly complex trajectory models: impulsive, Lambert, nuclear electric propulsion, nuclear electric propulsion including the Earth's gravity, solar electric propulsion including the Earth's gravity and we study the final correlation between the obtained target rankings. We find that the use of a low-thrust trajectory model is of great importance for target selection, since the use of chemical propulsion surrogates leads to favouring less attractive options 19\% of times, a percentage that drops to 8\% already using a simple nuclear electric propulsion model that neglects the Earth's gravity effects and thrust dependence on the solar distance. We also find that for the study case considered, a small interplanetary CubeSat named M-ARGO, the inclusion of the Earth's gravity in the considered dynamics does not affect the target selection significantly.

Finding the optimal ordering of k-subsets with respect to an objective function is known to be an extremely challenging problem. In this paper we introduce a new objective for this task, rooted in the problem of star identification on spacecrafts: subsets of detected spikes are to be generated in an ordering that minimizes time to detection of a valid star constellation. We carry out an extensive analysis of the combinatorial optimization problem, and propose multiple algorithmic solutions, offering different quality-complexity trade-offs. Three main approaches are investigated: exhaustive search (branch and prune), goal-driven (greedy scene elimination, minimally intersecting subsets), and stateless algorithms which implicitly seek to satisfy the problem's goals (pattern shifting, base unrank). In practical terms, these last algorithms are found to provide satisfactory approximations to the ideal performance levels, at small computational costs.

We consider the interplanetary trajectory design problem posed by the 8th edition of the Global Trajectory Optimization Competition and present the end-to-end strategy developed by the team ACT-ISAS (a collaboration between the European Space Agency's Advanced Concepts Team and JAXA's Institute of Space and Astronautical Science). The resulting interplanetary trajectory won 1st place in the competition, achieving a final mission value of $J=146.33$ [Mkm]. Several new algorithms were developed in this context but have an interest that go beyond the particular problem considered, thus, they are discussed in some detail. These include the Moon-targeting technique, allowing one to target a Moon encounter from a low Earth orbit; the 1-$k$ and 2-$k$ fly-by targeting techniques, enabling one to design resonant fly-bys while ensuring a targeted future formation plane% is acquired at some point after the manoeuvre ; the distributed low-thrust targeting technique, admitting one to control the spacecraft formation plane at 1,000,000 [km]; and the low-thrust optimization technique, permitting one to enforce the formation plane's orientations as path constraints.

Weierstrass elliptic and related functions have been recently shown to enable analytical explicit solutions to classical problems in astrodynamics. These include the constant radial acceleration problem, the Stark problem and the two-fixed center (or Euler's) problem. In this paper we review the basic technique that allows for these results and we discuss the limits and merits of the approach. Applications to interplanetary trajectory design are then discussed including low-thrust planetary fly-bys and the motion of an artificial satellite under the influence of an oblate primary including $J_2$ and $J_3$ harmonics.

We present for the first time an explicit, complete and closed-form solution to the three-dimensional problem of two fixed centres, based on Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our solution is exact, valid for all initial conditions and physical parameters of the system (including unbounded orbits and repulsive forces), and expressed via a unique set of formulae. Various properties of the three-dimensional problem of two fixed centres are investigated and analysed, with a particular emphasis on quasi-periodic and periodic orbits, regions of motion and equilibrium points.

While the constant radial acceleration problem is known to be integrable and has received some recent attention in an orbital mechanics context, a closed form explicit solution, relating the state variables to a time parameter, has eluded all researchers so far. It is here shown how such a solution exists and is elegantly expressed in terms of the Weierstrass elliptic and related functions. Previously known facts can be derived from the new explicit solution and new insights are revealed.

The orbital boundary value problem, also known as Lambert Problem, is revisited. Building upon Lancaster and Blanchard approach, new relations are revealed and a new variable representing all problem classes, under L-similarity, is used to express the time of flight equation. In the new variable, the time of flight curves have two oblique asymptotes and they mostly appear to be conveniently approximated by piecewise continuous lines. We use and invert such a simple approximation to provide an efficient initial guess to an Householder iterative method that is then able to converge, for the single revoltuion case, in only two iterations. The resulting algorithm is compared to Gooding's procedure revealing to be numerically as accurate, while having a smaller computational complexity.

We present a new general, complete closed-form solution of the Stark problem in terms of Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our analysis is exact and valid for all values of the external force field, and it is expressed via unique formul\ae{} valid for all initial conditions and parameters of the system. The simple form of the solution allows us to perform a thorough investigation of the properties of the dynamical system, including the identification of quasi-periodic and periodic orbits, the formulation of a simple analytical criterion to determine the boundness of the trajectory, and the characterisation of the equilibrium points.

We summarise the scientific and technological aspects of the SAGAS (Search for Anomalous Gravitation using Atomic Sensors) project, submitted to ESA in June 2007 in response to the Cosmic Vision 2015-2025 call for proposals. The proposed mission aims at flying highly sensitive atomic sensors (optical clock, cold atom accelerometer, optical link) on a Solar System escape trajectory in the 2020 to 2030 time-frame. SAGAS has numerous science objectives in fundamental physics and Solar System science, for example numerous tests of general relativity and the exploration of the Kuiper belt. The combination of highly sensitive atomic sensors and of the laser link well adapted for large distances will allow measurements with unprecedented accuracy and on scales never reached before. We present the proposed mission in some detail, with particular emphasis on the science goals and associated measurements.

The Doppler-tracking data of the Pioneer 10 and 11 spacecraft show an unmodelled constant acceleration in the direction of the inner Solar System. Serious efforts have been undertaken to find a conventional explanation for this effect, all without success at the time of writing. Hence the effect, commonly dubbed the Pioneer anomaly, is attracting considerable attention. Unfortunately, no other space mission has reached the long-term navigation accuracy to yield an independent test of the effect. To fill this gap we discuss strategies for an experimental verification of the anomaly via an upcoming space mission. Emphasis is put on two plausible scenarios: nondedicated concepts employing either a planetary exploration mission to the outer Solar System or a piggybacked micro-spacecraft to be launched from a mother spacecraft travelling to Saturn or Jupiter. The study analyses the impact of a Pioneer anomaly test on the system and trajectory design for these two paradigms. It is found that both paradigms are capable of verifying the Pioneer anomaly and determine its magnitude at 10% level. Moreover the concepts can discriminate between the most plausible classes of models of the anomaly, a central force, a blueshift of the radio signal and a drag-like force. The necessary adaptions of the system and mission design do not impair the planetary exploration goals of the missions.