Recent works demonstrated that the dynamics caused by the planetary oblateness coupled with the solar radiation pressure can be described through a model based on singly-averaged equations of motion. The coupled perturbations affect the evolution of the eccentricity, inclination and orientation of the orbit with respect to the Sun--Earth line. Resonant interactions lead to non-trivial orbital evolution that can be exploited in mission design. Moreover, the dynamics in the vicinity of each resonance can be analytically described by a resonant model that provides the location of the central and hyperbolic invariant manifolds which drive the phase space evolution. The classical tools of the dynamical systems theory can be applied to perform a preliminary mission analysis for practical applications. On this basis, in this work we provide a detailed derivation of the resonant dynamics, also in non-singular variables, and discuss its properties, by studying the main bifurcation phenomena associated to each resonance. Last, the analytical model will provide a simple analytical expression to obtain the area-to-mass ratio required for a satellite to deorbit from a given altitude in a feasible timescale.

The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its mission lifetime. Few software packages perform this analysis, and most of them employ time-consuming ray-tracing methodology, where particles are randomly sampled from relevant distributions. In addition, they tend not to consider the risk associated with the secondary debris clouds. The paper presents the development of a vulnerability assessment model, which relies on a fully statistical procedure: the debris fluxes are directly used combining them with the concept of the vulnerable zone, avoiding the random sampling the debris fluxes. A novel methodology is presented to predict damage to internal components. It models the interaction between the components and the secondary debris cloud through basic geometric operations, considering mutual shielding and shadowing between internal components. The methodologies are tested against state-of-the-art software for relevant test cases, comparing results on external structures and internal components.

The coupled attitude and orbit dynamics of solar sails is studied. The shape of the sail is a simplified quasi-rhombic-pyramid that provides the structure helio-stablility properties. After adimensionalisation, the system is put in the form of a fast-slow dynamical system where the different time scales are explicitely related to the physical parameters of the system. The orientation of the body frame with respect to the inertial orbit frame is a fast phase that can be averaged out. This gives rise to a simplified formulation that only consists of the orbit dynamics perturbed by a flat sail with fixed attitude perpendicular to the direction of the sunlight. The results are exemplified using numerical simulations.

Among the mitigation measures introduced to cope with the space debris issue there is the de-orbiting of decommissioned satellites. Guidelines for re-entering objects call for a ground casualty risk no higher than 0.0001. To comply with this requirement, satellites can be designed through a design-for-demise philosophy. Still, a spacecraft designed to demise has to survive the debris-populated space environment for many years. The demisability and the survivability of a satellite can both be influenced by a set of common design choices such as the material selection, the geometry definition, and the position of the components. Within this context, two models have been developed to analyse the demise and the survivability of satellites. Given the competing nature of the demisability and the survivability, a multi-objective optimisation framework was developed, with the aim to identify trade-off solutions for the preliminary design of satellites. As the problem is nonlinear and involves the combination of continuous and discrete variables, classical derivative based approaches are unsuited and a genetic algorithm was selected instead. The genetic algorithm uses the developed demisability and survivability criteria as the fitness functions of the multi-objective algorithm. The paper presents a test case, which considers the preliminary optimisation of tanks in terms of material, geometry, location, and number of tanks for a representative Earth observation mission. The configuration of the external structure of the spacecraft is fixed. Tanks were selected because they are sensitive to both design requirements: they represent critical components in the demise process and impact damage can cause the loss of the mission because of leaking and ruptures. The results present the possible trade off solutions, constituting the Pareto front obtained from the multi-objective optimisation.

The paper is concerned with examining the effects that design-for-demise solutions can have not only on the demisability of components, but also on their survivability that is their capability to withstand impacts from space debris. First two models are introduced. A demisability model to predict the behaviour of spacecraft components during the atmospheric re-entry and a survivability model to assess the vulnerability of spacecraft structures against space debris impacts. Two indices that evaluate the level of demisability and survivability are also proposed. The two models are then used to study the sensitivity of the demisability and of the survivability indices as a function of typical design-for-demise options. The demisability and the survivability can in fact be influenced by the same design parameters in a competing fashion that is while the demisability is improved, the survivability is worsened and vice versa. The analysis shows how the design-for-demise solutions influence the demisability and the survivability independently. In addition, the effect that a solution has simultaneously on the two criteria is assessed. Results shows which, among the design-for-demise parameters mostly influence the demisability and the survivability. For such design parameters maps are presented, describing their influence on the demisability and survivability indices. These maps represent a useful tool to quickly assess the level of demisability and survivability that can be expected from a component, when specific design parameters are changed.

This work introduces two Monte Carlo (MC)-based sampling methods, known as line sampling and subset simulation, to improve the performance of standard MC analyses in the context of asteroid impact risk assessment. Both techniques sample the initial uncertainty region in different ways, with the result of either providing a more accurate estimation of the impact probability or reducing the number of samples required during the simulation with respect to standard MC techniques. The two methods are firstly described and then applied to some test cases, providing evidence of the increased accuracy or the reduced computational burden with respect to a standard MC simulation. Finally, a sensitivity analysis is carried out to show how parameter setting affects the accuracy of the results and the numerical efficiency of the two methods.

The aim of this work is to provide an analytical model to characterize the equilibrium points and the phase space associated with the singly-averaged dynamics caused by the planetary oblateness coupled with the solar radiation pressure perturbations. A two-dimensional differential system is derived by considering the classical theory, supported by the existence of an integral of motion comprising semi-major axis, eccentricity and inclination. Under the single resonance hypothesis, the analytical expressions for the equilibrium points in the eccentricity-resonant angle space are provided, together with the corresponding linear stability. The Hamiltonian formulation is also given. The model is applied considering, as example, the Earth as major oblate body, and a simple tool to visualize the structure of the phase space is presented. Finally, some considerations on the possible use and development of the proposed model are drawn.

In this paper the planar orbit and attitude dynamics of an uncontrolled spacecraft is studied, taking on-board a deorbiting device. Solar and drag sails with the same shape are considered and separately studied. In both cases, these devices are assumed to have a simplified pyramidal shape that endows the spacecraft with helio and drag stable properties. The translational dynamics is assumed to be planar and hence the rotational dynamics occurs only around one of the principal axes of the spacecraft. Stable or slowly-varying attitudes are studied, subject to disturbances due to the Earth oblateness effect and gravity gradient torques, and either solar radiation pressure or atmospheric drag torque and acceleration. The results are analysed with respect to the aperture of the sail and the center of mass - center of pressure offset.

This paper introduces a new analytical method for the determination of the coverage area modeling the Earth as an oblate ellipsoid of rotation. Starting from the knowledge of the satellite's position vector and the direction of the navigation antenna line of sight, the surface generated by the intersection of the oblate ellipsoid and the assumed conical field of view is decomposed in many ellipses, obtained by cutting the Earth's surface with every plane containing the navigation antenna line of sight. The geometrical parameters of each ellipse can be derived analytically together with the points intersection of the conical field of view with the ellipse itself by assuming a proper value of the half-aperture angle or the minimum elevation angle from which the satellite can be considered visible from the Earth's surface. The method can be applied for different types of pointing (geocentric, geodetic and generic) according to the mission requirements. Finally, numerical simulations compare the classical spherical approach with the new ellipsoidal method in the determination of the coverage area, and also show the dependence of the coverage errors on some relevant orbital parameters.

Numerical integration of orbit trajectories for a large number of initial conditions and for long time spans is computationally expensive. Semi-analytical methods were developed to reduce the computational burden. An elegant and widely used method of semi-analytically integrating trajectories of objects subject to atmospheric drag was proposed by King-Hele (KH). However, the analytical KH contraction method relies on the assumption that the atmosphere density decays strictly exponentially with altitude. If the actual density profile does not satisfy the assumption of a fixed scale height, as is the case for Earth's atmosphere, the KH method introduces potentially large errors for non-circular orbit configurations. In this work, the KH method is extended to account for such errors by using a newly introduced atmosphere model derivative. By superimposing exponentially decaying partial atmospheres, the superimposed KH method can be applied accurately while considering more complex density profiles. The KH method is further refined by deriving higher order terms during the series expansion. A variable boundary condition to choose the appropriate eccentricity regime, based on the series truncation errors, is introduced. The accuracy of the extended analytical contraction method is shown to be comparable to numerical Gauss-Legendre quadrature. Propagation using the proposed method compares well against non-averaged integration of the dynamics, while the computational load remains very low.

In this work the orbital dynamics of Earth satellites about the geosynchronous altitude are explored, with primary goal to assess current mitigation guidelines as well as to discuss the future exploitation of the region. A thorough dynamical mapping was conducted in a high-definition grid of orbital elements, enabled by a fast and accurate semi-analytical propagator, which considers all the relevant perturbations. The results are presented in appropriately selected stability maps to highlight the underlying mechanisms and their interplay, that can lead to stable graveyard orbits or fast re-entry pathways. The natural separation of the long-term evolution between equatorial and inclined satellites is discussed in terms of post-mission disposal strategies. Moreover, we confirm the existence of an effective cleansing mechanism for inclined geosynchronous satellites and discuss its implications in terms of current guidelines as well as alternative mission designs that could lead to a sustainable use of the geosynchronous orbital region.

Euclid is a space-based optical/near-infrared survey mission of the European Space Agency (ESA) to investigate the nature of dark energy, dark matter and gravity by observing the geometry of the Universe and on the formation of structures over cosmological timescales. Euclid will use two probes of the signature of dark matter and energy: Weak gravitational Lensing, which requires the measurement of the shape and photometric redshifts of distant galaxies, and Galaxy Clustering, based on the measurement of the 3-dimensional distribution of galaxies through their spectroscopic redshifts. The mission is scheduled for launch in 2020 and is designed for 6 years of nominal survey operations. The Euclid Spacecraft is composed of a Service Module and a Payload Module. The Service Module comprises all the conventional spacecraft subsystems, the instruments warm electronics units, the sun shield and the solar arrays. In particular the Service Module provides the extremely challenging pointing accuracy required by the scientific objectives. The Payload Module consists of a 1.2 m three-mirror Korsch type telescope and of two instruments, the visible imager and the near-infrared spectro-photometer, both covering a large common field-of-view enabling to survey more than 35% of the entire sky. All sensor data are downlinked using K-band transmission and processed by a dedicated ground segment for science data processing. The Euclid data and catalogues will be made available to the public at the ESA Science Data Centre.