The dynamics of solar sails with a non-point source of radiation pressure

The form of the solar radiation pressure on a heliocentric orbiting solar sail is obtained for a finite angular sized and limb darkened solar disk by the use of the radiation pressure tensor. It is found that the usual inverse square variation of the solar radiation pressure is modified by the finite angular size, and to a lesser extent by the solar limb darkening. The actual magnitude of the modification is in itself small, except at close heliocentric distances. However, its existence has implications for the dynamical stability of solar sails both in parked and circular orbital configurations and for the accuracy of trajectory calculations, particularly for sails in the inner solar system.     

Long-term evaluation of three-dimensional heliocentric solar sail trajectories with arbitrary fixed sail setting

The solar radiation effects upon the orbital behaviour of an arbitrarily shaped spacecraft (or a solar sail in particular) in a general fixed orientation with respect to the local coordinate frame are investigated. Through introduction of a quasi-angle in the osculating plane, the motion of the orbital plane becomes uncoupled from the in-plane perturbations. Exact solutions in the form of conic sections and logarithmic spirals can readily be formulated for certain specific initial conditions. An effective out-of-plane spiral transfer trajectory is obtained by reversing the force component normal to the orbital plane at specified positions in the orbit. By choosing the appropriate control angles for the sail orientation, any point in space can be reached eventually. In the case of general initial conditions, the long-term orbital behaviour is assessed asymptotically by means of the two-variable expansion procedure. An implicit expression for the eccentricity is derived and explicit results are established by an iteration scheme. The other orbital elements can be expressed in terms of the eccentricity and their asymptotic series for near-circular initial orbits are also obtained. While equations for the higher-order contributions as well as the periodic parts of their solutions can be formulated readily, their secular terms are determined only for a circular initial orbit.     

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.     

Logarithmic spiral trajectories generated by Solar sails

Analytic solutions to continuous thrust-propelled trajectories are available in a few cases only. An interesting case is offered by the logarithmic spiral, that is, a trajectory characterized by a constant flight path angle and a fixed thrust vector direction in an orbital reference frame. The logarithmic spiral is important from a practical point of view, because it may be passively maintained by a Solar sail-based spacecraft. The aim of this paper is to provide a systematic study concerning the possibility of inserting a Solar sail-based spacecraft into a heliocentric logarithmic spiral trajectory without using any impulsive maneuver. The required conditions to be met by the sail in terms of attitude angle, propulsive performance, parking orbit characteristics, and initial position are thoroughly investigated. The closed-form variations of the osculating orbital parameters are analyzed, and the obtained analytical results are used for investigating the phasing maneuver of a Solar sail along an elliptic heliocentric orbit. In this mission scenario, the phasing orbit is composed of two symmetric logarithmic spiral trajectories connected with a coasting arc.     

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.     

Utilization of an H-reversal trajectory of a solar sail for asteroid deflection

Near Earth Asteroids have a possibility of impacting the Earth and always represent a threat. This paper proposes a way of changing the orbit of the asteroid to avoid an impact. A solar sail evolving in an H-reversal trajectory is utilized for asteroid deflection. Firstly, the dynamics of the solar sail and the characteristics of the H-reversal trajectory are analyzed. Then, the attitude of the solar sail is optimized to guide the sail to impact the target asteroid along an H-reversal trajectory. The impact...     

Asteroid body-fixed hovering using nonideal solar sails

The problem of body-fixed hovering over an asteroid using a compact form of nonideal solar sails with a controllable area is investigated.Nonlinear dynamic equations describing the hovering problem are constructed for a spherically symmetric asteroid.Numerical solutions of the feasible region for body-fixed hovering are obtained.Different sail models,including the cases of ideal,optical,parametric and solar photon thrust,on the feasible region is studied through numerical simulations.The influence of the asteroid spinning rate and the sail area-to-mass ratio on the feasible region is discussed.The required orientations for the sail and their corresponding variable lightness numbers are given for different hovering radii to identify the feasible region of the body-fixed hovering.An attractive scenario for a mission is introduced to take advantage of solar sail hovering.     

A continuum model for the orbit evolution of self-propelled ‘smart dust’ swarms

A continuity equation is developed to model the evolution of a swarm of self-propelled ‘smart dust’ devices in heliocentric orbit driven by solar radiation pressure. These devices are assumed to be MEMs-scale (micro-electromechanical systems) with a large area-to-mass ratio. For large numbers of devices it will be assumed that a continuum approximation can be used to model their orbit evolution. The families of closed-form solutions to the resulting swarm continuity equation then represent the evolution of the number density of devices as a function of both position and time from a set of initial data. Forcing terms are also considered which model swarm sources and sinks (device deposition and device failure). The closed-form solutions presented for the swarm number density provide insights into the behaviour of swarms of self-propelled ‘smart dust’ devices an can form the basis of more complex mission design methodologies.     

Automatic Determination of the Conic Coronal Mass Ejection Model Parameters

Characterization of the three-dimensional structure of solar transients using incomplete plane of sky data is a difficult problem whose solutions have potential for societal benefit in terms of space weather applications. In this paper transients are characterized in three dimensions by means of conic coronal mass ejection (CME) approximation. A novel method for the automatic determination of cone model parameters from observed halo CMEs is introduced. The method uses both standard image processing techniques to extract the CME mass from white-light coronagraph images and a novel inversion routine providing the final cone parameters. A bootstrap technique is used to provide model parameter distributions. When combined with heliospheric modeling, the cone model parameter distributions will provide direct means for ensemble predictions of transient propagation in the heliosphere. An initial validation of the automatic method is carried by comparison to manually determined cone model parameters. It is shown using 14 halo CME events that there is reasonable agreement, especially between the heliocentric locations of the cones derived with the two methods. It is argued that both the heliocentric locations and the opening half-angles of the automatically determined cones may be more realistic than those obtained from the manual analysis.     

Lissajous motion near Lagrangian point $$ L_{2} $$ in radial solar sail

An attempt was made to study the dynamics close to the collinear libration point \( L_{2} \) of the radial solar sail circular-restricted three-body problem (RSCRTBP) in the Sun–Jupiter System, where the third massless body is a solar sail. We analyse the qausi-periodic (Lissajous solutions) orbits about the libration point \( L_{2} \). The Lindstedt–Poincaré approximation for the qausi-periodic orbits was used for numerical simulations. We utilized linear quadratic regulator (LQR) to stabilize the full nonlinear model, and linear state-feedback controller was designed to stabilize the trajectory.     

Artificial equilibrium points for a generalized sail in the circular restricted three-body problem

This paper introduces a new approach to the study of artificial equilibrium points in the circular restricted three-body problem for propulsion systems with continuous and purely radial thrust. The propulsion system is described by means of a general mathematical model that encompasses the behavior of different systems like a solar sail, a magnetic sail and an electric sail. The proposed model is based on the choice of a coefficient related to the propulsion type and a performance parameter that quantifies the system technological complexity. The propulsion system is therefore referred to as generalized sail. The existence of artificial equilibrium points for a generalized sail is investigated. It is shown that three different families of equilibrium points exist, and their characteristic locus is described geometrically by varying the value of the performance parameter. The linear stability of the artificial points is also discussed.     

Interfacing MHD Single Fluid and Kinetic Exospheric Solar Wind Models and Comparing Their Energetics

An exospheric kinetic solar wind model is interfaced with an observation-driven single-fluid magnetohydrodynamic (MHD) model. Initially, a photospheric magnetogram serves as observational input in the fluid approach to extrapolate the heliospheric magnetic field. Then semi-empirical coronal models are used for estimating the plasma characteristics up to a heliocentric distance of 0.1 AU. From there on, a full MHD model that computes the three-dimensional time-dependent evolution of the solar wind macroscopic variables up to the orbit of Earth is used. After interfacing the density and velocity at the inner MHD boundary, we compare our results with those of a kinetic exospheric solar wind model based on the assumption of Maxwell and Kappa velocity distribution functions for protons and electrons, respectively, as well as with in situ observations at 1 AU. This provides insight into more physically detailed processes, such as coronal heating and solar wind acceleration, which naturally arise from including suprathermal electrons in the model. We are interested in the profile of the solar wind speed and density at 1 AU, in characterizing the slow and fast source regions of the wind, and in comparing MHD with exospheric models in similar conditions. We calculate the energetics of both models from low to high heliocentric distances.     

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.     

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.     

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.     

Planetary close encounters: geometry of approach and post-encounter orbital parameters

Öpik's assumptions on the geometry of particle trajectories leading to and through planetary close encounters are used to compute the distribution of changes in heliocentric orbital elements that result from such encounters for a range of initial heliocentric orbits. Behaviour at encounter is assumed to follow two-body (particle—planet) gravitational scattering, while before and after encounter particle motion is only governed by the force of the Sun. Derivation of these distributions allows precise analysis of the probability of various outcomes in terms of the physical characteristics of the bodies involved. For example, they allow an explanation and prediction of the asymmetry of the extreme energy perturbations for different initial orbits. The formulae derived here may be applied to problems including the original accumulation of planets and satellites, and the continuing evolution of populations of small bodies, such as asteroids and comets.     

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.     

Interplanetary disturbances in the solar wind produced by density,temperature, or velocity pulses at 0.08 AU

Time-dependent solutions of a one-fluid model of the interplanetary medium are investigated. This set of unsteady hydrodynamic equations has been written in conservation form in order to apply the Lax-Wendroff method for the solution of this problem. The initial disturbance is specified by a pulse at 0.08 AU (astronomical units). Physically, this pulse can be interpreted as having been caused by a solar flare, surge, or any other solar disturbance. The equilibrium condition is determined to be the steady solution of the governing equations and represents the quiet solar wind. The results are presented in terms of density, temperature, and velocity profiles of the interplanetary gas flow at heliocentric distances up to 6 AU at several times. Also, the trajectories of disturbances for various initial pulses are shown. Finally, we have used some June 1972 interplanetary observational data to compare with these theoretical calculations. On the basis of these results, the effects of solar disturbances on the interplanetary environment (such as the generation of large non-linear wave trains in the shocks' wakes) can be inferred.     

Energy Gain of Positive Ions in Solar Polar Coronal Holes

We have investigated heating of solar polar coronal holes and acceleration of fast solar wind by means of lower hybrid (LH) waves. A three-fluid Maxwell model comprising electrons, protons, and α-particles is employed at around two solar radii heliocentric distance, where wave dissipation starts to be dominated by collisionless processes. We suggest specific wavenumber ranges corresponding to LH as well as stochastic instabilities and find that these instabilities may bring about a significant energy gain in positive ions.