Non-Keplerian orbits for electric sails

An electric sail is capable of guaranteeing the fulfilment of a class of trajectories that would be otherwise unfeasible through conventional propulsion systems. In particular, the aim of this paper is to analyze the electric sail capabilities of generating a class of displaced non-Keplerian orbits, useful for the observation of the Sun’s polar regions. These orbits are characterized through their physical parameters (orbital period and solar distance) and the spacecraft propulsion capabilities. A comparison with a solar sail is made to highlight which of the two systems is more convenient for a given mission scenario. The optimal (minimum time) transfer trajectories towards the displaced orbits are found with an indirect approach.     

Three-dimensional time optimal double angular momentum reversal trajectory using solar sails

A new concept of three dimensional non-Keplerian trajectories with double angular momentum reversal is investigated with high performance solar sails. The main discussion of this paper is about such 3D solar inverse orbits with inner constraints. The problem is addressed in a time optimal control framework solved by an indirect method. Two typical solar inverse orbits have been achieved and presented in a 3D non-dimensional dynamic model in the Heliocentric Inertial Frame. Starting from the Earth orbit ecliptic plane, a sailcraft in the inverse orbit exhibits a butterfly shape trajectory. As such, the new orbits are symmetrical with respect to a plane which contains the Sun-perihelion line. The relation of the sail attitude angles between the two symmetrical parts of the orbits are used to reduce the simulation effort. The quasi-heliostationary property at its aphelia is demonstrated with variation of the orbital radius. Evolutions of the orbital velocity and optimal sail orientations are also outlined and discussed to benefit future design work. As is suited for space observation guaranteed by its butterfly shape, the inverse orbits are thoroughly studied in terms of the concerned parameters. The discussion of the parametric influence is ranked in order as perihelion distance r _E , required maximum position z _max, perihelion position z _f and the sail lightness number β. Suitable ranges of each parameter are adopted to illustrate the orbital variation trend. Through numerical simulations the features of such inverse orbits are further emphasized to provide an initial reference for future researchers.     

Sun-perturbed earth-to-moon transfers with low energy and moderate flight time

We construct a spacecraft transfer with low cost and moderate flight time from the Earth to the Moon. The motion of the spacecraft is modeled by the planar circular restricted three-body problem including a perturbation due to the solar gravitation. Our approach is to reduce computation of optimal transfers to a non-linear boundary value problem. Using a computer software called AUTO, we solve it and continue its solutions numerically to obtain the optimal transfers. Our result also shows that the use of the solar gravitation can further lower the transfer cost drastically.     

Planar heliocentric roto-translatory motion of a spacecraft with a solar sail of complex shape

A complete treatment of the general motion of rotation and translation of a solar-sail spacecraft is proposed for the non-flat sail of complex shape. The planar heliocentric roto-translatory motion is considered, orbit-rotational coupling in the problem of altitude and orbital sail motion is investigated for the two-folding sail formed by two unequal reflective rectangular plates oriented at a right angle. The problem of orbit-rotational coupling is essentially a planar one: both sail plates are orthogonal to the orbital plane. The possibility of the non-controlled interplanetary transfer with such two-folding sail at its passive radiational orientation is established analytically from point of view of orbit-rotational coupling. Optimal geometric proportions of this sail are found at minimum-time interplanetary transfers.     

Laboratory Experiment of Plasma Flow Around Magnetic Sail

To propel a spacecraft in the direction leaving the Sun, a magnetic sail (MagSail) blocks the hypersonic solar wind plasma flow by an artificial magnetic field. In order to simulate the interaction between the solar wind and the artificially deployed magnetic field produced around a magnetic sail spacecraft, a laboratory simulator was designed and constructed inside a space chamber. As a solar wind simulator, a high-power magnetoplasmadynamic arcjet is operated in a quasisteady mode of 0.8 ms duration. It can generate a simulated solar wind that is a high-speed (above 20 km/s), high-density (10¹⁸ m⁻³) hydrogen plasma plume of ∼0.7 m in diameter. A small coil (2 cm in diameter), which is to simulate a magnetic sail spacecraft and can obtain 1.9-T magnetic field strength at its center, was immersed inside the simulated solar wind. Using these devices, the formation of a magnetic cavity (∼8 cm in radius) was observed around the coil, which indicates successful simulation of the plasma flow of a MagSail in the laboratory.     

Artificial halo orbits for low-thrust propulsion spacecraft

We consider periodic halo orbits about artificial equilibrium points (AEP) near to the Lagrange points L ₁ and L ₂ in the circular restricted three body problem, where the third body is a low-thrust propulsion spacecraft in the Sun–Earth system. Although such halo orbits about artificial equilibrium points can be generated using a solar sail, there are points inside L ₁ and beyond L ₂ where a solar sail cannot be placed, so low-thrust, such as solar electric propulsion, is the only option to generate artificial halo orbits around points inaccessible to a solar sail. Analytical and numerical halo orbits for such low-thrust propulsion systems are obtained by using the Lindstedt Poincaré and differential corrector method respectively. Both the period and minimum amplitude of halo orbits about artificial equilibrium points inside L ₁ decreases with an increase in low-thrust acceleration. The halo orbits about artificial equilibrium points beyond L ₂ in contrast show an increase in period with an increase in low-thrust acceleration. However, the minimum amplitude first increases and then decreases after the thrust acceleration exceeds 0.415 mm/s². Using a continuation method, we also find stable artificial halo orbits which can be sustained for long integration times and require a reasonably small low-thrust acceleration 0.0593 mm/s².     

Solar sail time-optimal interplanetary transfer trajectory design

The fuel consumption associated with some interplanetary transfer trajectories using chemical propulsion is not affordable. A solar sail is a method of propulsion that does not consume fuel. Transfer time is one of the most pressing problems of solar sail transfer trajectory design. This paper investigates the time-optimal interplanetary transfer trajectories to a circular orbit of given inclination and radius. The optimal control law is derived from the principle of maximization. An indirect method is used...     

Periodic and quasi-periodic motions of a solar sail close to SL 1 in the Earth–Sun system

Solar sails are a proposed form of spacecraft propulsion using large membrane mirrors to propel a satellite taking advantage of the solar radiation pressure. To model the dynamics of a solar sail we have considered the Earth–Sun Restricted Three Body Problem including the Solar radiation pressure (RTBPS). This model has a 2D surface of equilibrium points parametrised by the two angles that define the sail orientation. In this paper we study the non-linear dynamics close to an equilibrium point, with special interest in the bounded motion. We focus on the region of equilibria close to SL ₁, a collinear equilibrium point that lies between the Earth and the Sun when the sail is perpendicular to the Sun–sail direction. For different fixed sail orientations we find families of planar, vertical and Halo-type orbits. We have also computed the centre manifold around different equilibria and used it to describe the quasi-periodic motion around them. We also show how the geometry of the phase space varies with the sail orientation. These kind of studies can be very useful for future mission applications.     

Analytical results for solar sail optimal missions with modulated radial thrust

The aim of this paper is to analyze the optimal trajectories of a spacecraft subjected to a modulated radial thrust, whose magnitude is inversely proportional to the square of the distance from the primary body. This case is representative of a Sun-facing solar sail with a passive attitude control system. In this study the sailcraft is assumed to perform a finite number of reorientation maneuvers to set the propelling acceleration to zero and generate suitable coasting arcs along the trajectory. Accordingly, the resulting generalized orbit is a sequence of either propelled or ballistic conic arcs, whose main characteristics (in terms of semimajor axis, eccentricity, and perihelion radius) can be calculated in closed form. As a result, the sailcraft optimal performance can be studied using an analytical approach. In particular, some compact relationships are drawn and discussed that allow one to find the optimal sailcraft characteristics required to reach a prescribed final orbit.     

Controlled Solar Sail Transfers into Near-Sun Regions Combined with Planetary Gravity-Assist Flybys

Results of numerical simulations of 'local-optimal' (or 'instantaneously optimal') trajectories of a space probe with a flat solar sail which moves from the circular Earth orbit to near-Sun regions are presented. We examine planar (ecliptic) solar sail transfer with gravity-assist flybys of Earth, Venus and Mercury. Several complex control modes of the sail tilt orientation angle for near-Sun orbits and for some 'falling onto the Sun' trajectories are investigated. The numerical simulations are used to examine the flight duration of some sail missions and to investigate the evolution of osculating elliptical orbits.     

New Challenges to Trajectory Design by the Use of Electric Propulsion and Other New Means of Wandering in the Solar System

The design of spacecraft trajectories is a crucial part of a space mission design. Often the mission goal is tightly related to the spacecraft trajectory. A geostationary orbit is indeed mandatory for a stationary equatorial position. Visiting a solar system planet implies that a proper trajectory is used to bring the spacecraft from Earth to the vicinity of the planet. The first planetary missions were based on conventional trajectories obtained with chemical engine rockets. The manoeuvres could be considered 'impulsive' and clear limitations to the possible missions were set by the energy required to reach certain orbits. The gravity-assist trajectories opened a new way of wandering through the solar system, by exploiting the gravitational field of some planets. The advent of other propulsion techniques, as electric or ion propulsion and solar sail, opened a new dimension to the planetary trajectory, while at the same time posing new challenges. These 'low thrust' propulsion techniques cannot be considered 'impulsive' anymore and require for their study mathematical techniques which are substantially different from before. The optimisation of such trajectories is also a new field of flight dynamics, which involves complex treatments especially in multi-revolution cases as in a lunar transfer trajectory. One advantage of these trajectories is that they allow to explore regions of space where different bodies gravitationally compete with each other. We can exploit therefore these gravitational perturbations to save fuel or reduce time of flight. The SMART-1 spacecraft, first European mission to the Moon, will test for the first time all these techniques. The paper is a summary report on various activities conducted by the project team in these areas.     

Formation flying solar-sail gravity tractors in displaced orbit for towing near-Earth asteroids

Several methods of asteroid deflection have been proposed in literature and the gravitational tractor is a new method using gravitational coupling for near-Earth object orbit modification. One weak point of gravitational tractor is that the deflection capability is limited by the mass and propellant of the spacecraft. To enhance the deflection capability, formation flying solar sail gravitational tractor is proposed and its deflection capability is compared with that of a single solar sail gravitational tractor. The results show that the orbital deflection can be greatly increased by increasing the number of the sails. The formation flying solar sail gravitational tractor requires several sails to evolve on a small displaced orbit above the asteroid. Therefore, a proper control should be applied to guarantee that the gravitational tractor is stable and free of collisions. Two control strategies are investigated in this paper: a loose formation flying realized by a simple controller with only thrust modulation and a tight formation realized by the sliding-mode controller and equilibrium shaping method. The merits of the loose and tight formations are the simplicity and robustness of their controllers, respectively.     

Periodic orbits of solar sail equipped with reflectance control device in Earth–Moon system

In this paper, families of Lyapunov and halo orbits are presented with a solar sail equipped with a reflectance control device in the Earth–Moon system. System dynamical model is established considering solar sail acceleration, and four solar sail steering laws and two initial Sun-sail configurations are introduced. The initial natural periodic orbits with suitable periods are firstly identified. Subsequently, families of solar sail Lyapunov and halo orbits around the \(L_{1}\) and \(L_{2}\) points are designed with fixed solar sail characteristic acceleration and varying reflectivity rate and pitching angle by the combination of the modified differential correction method and continuation approach. The linear stabilities of solar sail periodic orbits are investigated, and a nonlinear sliding model controller is designed for station keeping. In addition, orbit transfer between the same family of solar sail orbits is investigated preliminarily to showcase reflectance control device solar sail maneuver capability.     

Control of chaos in the vicinity of the Earth‐Moon L5 Lagrangian point to keep a spacecraft in orbit

The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic idea of the method of Ott, Grebogi, and York (OGY method) to stabilize the motion of a spacecraft around the triangular Lagrange point L₅ of the Earth‐Moon system. We have determined the escape rate of the trajectories in the general three‐ and four‐body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail the phase space around the L₅ point of the Earth‐Moon system. Using the phase space portrait our next goal was to apply a modified OGY method to keep a spacecraft close to the vicinity of L₅. We modified the equation of motions with the addition of a time dependent force to the motion of the spacecraft. In our orbit‐keeping procedure there are three free parameters: (i) the magnitude of the thrust, (ii) the start time, and (iii) the length of the control. Based on our numerical experiments we were able to determine possible values for these parameters and successfully apply a control phase to a spacecraft to keep it on orbit around L₅. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Evolution of the \mathcal {L}_1 halo family in the radial solar sail circular restricted three-body problem

We present a detailed investigation of the dramatic changes that occur in the \(\mathcal {L}_1\) halo family when radiation pressure is introduced into the Sun–Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. in Int J Bifurc Chaos 1:493–520, 1991) the families can be fully mapped out as the parameter \(\beta \) is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at \(\beta \approx 0.0387\) acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at \(\beta \approx 0.289\) , resulting in a third new family. The linear stability of the families changes rapidly at low values of \(\beta \) , with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter \(\beta \) and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of \(\beta \) there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case.     

Linking Remote-Sensing and <Emphasis Type="Italic">In Situ</Emphasis> Observations of Coronal Mass Ejections Using STEREO

The twin STEREO spacecraft have been observing the Sun since 2006. Even though STEREO has only been active during solar minimum conditions so far, an important number of coronal mass ejections (CMEs) and their interplanetary counterparts (ICMEs) have been observed. Many of the ICMEs can be linked back to the corresponding CMEs on the Sun through the combination of remote-sensing and in situ observations. This paper aims to answer the question whether a CME observed by a coronagraph will be detected in situ by a spacecraft in a specific location in the heliosphere. We use a flux-rope-like model fit to the STEREO SECCHI/COR2 data to obtain the direction of CME propagation and its geometrical configuration in three dimensions. Based on model parameters, we then calculate their angular widths and determine whether they should have been detected by STEREO-A, STEREO-B, Wind or ACE. We compare the results with corresponding in situ observations of ICMEs. We find that predictions of ICME detections on the base of COR2 data generally match well the actual in situ observations.     

Optimization of the interplanetary trajectories of spacecraft with a solar electric propulsion power plant of minimal power

The problem of optimizing the interplanetary trajectories of a spacecraft (SC) with a solar electric propulsion system (SEPS) is examined. The problem of investigating the permissible power minimum of the solar electric propulsion power plant required for a successful flight is studied. Permissible ranges of thrust and exhaust velocity are analyzed for the given range of flight time and final mass of the spacecraft. The optimization is performed according to Portnyagin’s maximum principle, and the continuation method is used for reducing the boundary problem of maximal principle to the Cauchy problem and to study the solution/ parameters dependence. Such a combination results in the robust algorithm that reduces the problem of trajectory optimization to the numerical integration of differential equations by the continuation method.     

Solar Sail Halo Orbits at the Sun–Earth Artificial L 1 Point

Halo orbits for solar sails at artificial Sun–Earth L₁ points are investigated by a third order approximate solution. Two families of halo orbits are explored as defined by the sail attitude. Case I: the sail normal is directed along the Sun-sail line. Case II: the sail normal is directed along the Sun–Earth line. In both cases the minimum amplitude of a halo orbit increases as the lightness number of the solar sail increases. The effect of the z-direction amplitude on x- or y-direction amplitude is also investigated and the results show that the effect is relatively small. In case I, the orbit period increases as the sail lightness number increases, while in case II, as the lightness number increases, the orbit period increases first and then decreases after the lightness number exceeds ~0.01.     

On Motions Near the Lagrange Equilibrium Point L 4 in the Case of Routh's Critical Mass Ratio

We deal with the planar restricted circular problem of three bodies. We study trajectories in a small neighborhood of the Lagrange equilibrium point L ₄ when mass ratio is close to Routh's value. In particular, we show that the case of proper degeneracy takes place and for most initial conditions trajectories are conditionally-periodic. We obtain an approximate representation of families of periodic solution emanating from the equilibrium point L ₄. We also show that in the case of instability of L ₄ the trajectories starting in a vicinity of L ₄ remain in a finite domain forever. We give an upper bound of this domain. To carry out our investigation, we analyze the dynamics of a general Hamiltonian system with two degrees of freedom in the case of 1 : 1 resonance in detail.This revised version was published online in October 2005 with corrections to the Cover Date.     

Tracking a Ulysses High-latitude ICME Event Back to Its Solar Origins

High-latitude interplanetary mass ejections (ICMEs) observed beyond 1 AU are not studied very often. They are useful for improving our understanding of the 3D heliosphere. As there are only few such events registered by the Ulysses spacecraft, the task of detecting their solar counterparts is a challenge, especially during high solar activity periods, because there are dozens coronal mass ejections (CMEs) registered by SOHO that might be chosen as candidates. We analyzed a high-latitude ICME registered by the Ulysses spacecraft on 18 January 2002. Our investigation focused on the correlation between various plasma parameters that allow the identification to be made of the ICME and its components such as the forward shock, the magnetic cloud and the reverse shock.