Families of halo orbits in the elliptic restricted three-body problem for a solar sail with reflectivity control devices

Solar sail halo orbits designed in the Sun-Earth circular restricted three-body problem (CR3BP) provide inefficient reference orbits for station-keeping since the disturbance due to the eccentricity of the Earth’s orbit has to be compensated for. This paper presents a strategy to compute families of halo orbits around the collinear artificial equilibrium points in the Sun-Earth elliptic restricted three-body problem (ER3BP) for a solar sail with reflectivity control devices (RCDs). In this non-autonomous model, periodic halo orbits only exist when their periods are equal to integer multiples of one year. Here multi-revolution halo orbits with periods equal to integer multiples of one year are constructed in the CR3BP and then used as seeds to numerically continue the halo orbits in the ER3BP. The linear stability of the orbits is analyzed which shows that the in-plane motion is unstable while the out-of-plane motion is neutrally stable and a bifurcation is identified. Finally, station-keeping is performed which shows that a reference orbit designed in the ER3BP is significantly more efficient than that designed in the CR3BP, while the addition of RCDs improve station-keeping performance and robustness to uncertainty in the sail lightness number.     

Bound orbits of solar sails and general relativity

We study general relativistic effects on the bound orbits of solar sails. The combined effects of spacetime curvature and solar radiation pressure (SRP) lead to deviations from Kepler’s third law. Such kind of deviations also arise from frame dragging, the gravitational multipole moments of the sun, a net electric charge on the sun, and a positive cosmological constant. The SRP increases these deviations by several orders of magnitude, possibly rendering some of them detectable. We consider how the SRP modifies the perihelion shift of non-circular orbits, as well as the Lense-Thirring effect involving the precession of polar orbits. We investigate how the pitch angle for non-Keplerian orbits changes due to the partial absorption of light, general relativistic effects, and the oblateness of the sun. It is predicted that there is an analog of the Lense-Thirring effect for non-Keplerian orbits, in that the orbital plane precesses around the sun. We also consider the Poynting–Robertson effect and show that this effect can, in principle, be compensated for by an extremely small tilt of the solar sail.     

Resonant quasi-periodic near-rectilinear Halo orbits in the Elliptic-Circular Earth-Moon-Sun Problem

Motivated by the near-future re-exploration of the cislunar space, this paper investigates dynamical substitutes of the Earth-Moon’s resonant Near-Rectilinear Halo Orbits (NRHOs) under the Elliptic-Circular Restricted Four-Body Problem formulation of the Earth-Moon-Sun system. This model considers that the Earth and Moon move in elliptical orbits about each other and that a third body, the Sun, moves in a circular orbit about the Earth-Moon barycenter. By making use of this higher-fidelity dynamical model, we are able to incorporate the Sun’s influence and the Moon’s eccentricity, two of the most significant perturbations of the cislunar environment. As a result of these perturbations, resonant periodic NRHOs of the Earth-Moon Circular Restricted Three-Body Problem (CR3BP) are hereby replaced by two-dimensional quasi-periodic tori that better represent the dynamical evolution of satellites near the vicinity of the Moon. We present the steps and algorithms needed to compute these dynamical structures in the Elliptic-Circular model and subsequently assess their utility for spacecraft missions. We focus on the planned orbit for the NASA-led Lunar Gateway mission, a 9:2 synodic resonant L₂ southern NRHO, as well as on the 4:1 synodic and 4:1 sidereal resonances, due to the proximity to the nominal orbit and their advantageous dynamical properties. We verify that the dynamical equivalents of these orbits preserve key dynamical attributes such as eclipse avoidance and near-linear stability. Furthermore, we find that the higher dimensionality of quasi-periodic solutions offers interesting alternatives to mission designers in terms of phasing maneuvers and low-altitude scientific observations.     

Periodic attitude motions along planar orbits in the elliptic restricted three-body problem

A way to improve the accuracy of the three-body problem model is taking into account the eccentricity of primary attractors. Elliptic Restricted Three-Body Problem (ER3BP) is a model for studying spacecraft trajectory within the three-body problem such that the orbital eccentricity of primaries is reflected in it. As the principal cause of perturbation in the employed dynamical model, the primaries eccentricity changes the structure of orbits compared to the ideal Circular Restricted Three-Body Problem (CR3BP). It also changes the attitude behavior of a spacecraft revolving along periodic orbits in this regime. In this paper, the coupled orbit-attitude dynamics of a spacecraft in the ER3BP are exploited to find precise periodic solutions as the spacecraft is considered to be in planar orbits around Lagrangian points and Distant Retrograde Orbits (DRO). Periodic solutions are repetitious behaviors in which spacecraft whole dynamics are repeated periodically, these periodic behaviors are the main interest of this study because they are beneficial for future mission designs and allow delineation of the system’s governing dynamics. Previous studies laid the foundation for spacecraft stability analysis or studying pitch motion of spacecraft in the ER3BP regime. While in this paper, at first, initial guesses for correction algorithms were derived through verified search methods, then correction algorithms were used to refine calculated orbit-attitude periodic behaviors. Periodic orbits and full periodic solutions are portrayed and compared to previous studies and simpler models. Natural periodic solutions are valuable information eventuate in the longer functional lifetime of spacecraft. Since the problem assumption considered in this paper is much closer to real mission conditions, these results may be the means to use natural bounded motions in the actual operational environment.     

Solar sails for perturbation relief: Application to asteroids

A major cause of spacecraft orbital variation comes from natural perturbations, which, in close proximity of a body, are dominated by its non-spherical nature. For small bodies, such as asteroids, these effects can be considerable, given their uneven (and uncertain) mass distribution. Solar sail technology is proposed to reduce or eliminate the net secular effects of the irregular gravity field on the orbit. Initially, a sensitivity analysis will be carried out on the system which will show high sensitivity to changes in initial conditions. This presents a challenge for optimisation methods which require an initial guess of the solution. As such, the Genetic Algorithm (GA) is proposed as the preferred optimisation method as this requires no initial guess from the user. A multi-objective optimisation is performed which aims to achieve a periodic orbit whilst also minimising the effort required by the sail to do so. Given the system sensitivity, the control law for one orbit is not necessarily applicable for any subsequent orbit. Therefore, a new method of updating the control law for subsequent orbits is presented, based on linearisation and use of a Control Transition Matrix (CTM). The techniques will later find application in a multiple asteroid rendezvous mission with a solar sail as the primary propulsion system.     

Formation flying for magnetic sails around artificial equilibrium points

This paper investigates the problem of magnetic sail-based spacecraft formation control around the artificial equilibrium points (AEPs), which can eliminate the requirement of the propellant. The thrusts are achieved by utilizing the interaction between the solar wind and the artificial magnetosphere generated by superconducting current coil onboard. The circular restricted three-body problem (CRTBP) of magnetic sail is discussed including the allowed regions and linear stability of AEPs, the locations of collinear AEPs and the possibility of existence of periodic orbits around the collinear AEPs. Next, the dynamical models of magnetic sail formation around the collinear AEPs are established. A novel fast fixed-time nonsingular terminal sliding mode controller (FFNTSM) based on fixed-time disturbance observer (FTDO) is developed to account for external disturbances. Several numerical simulations are conducted to substantiate that spacecraft formation can be precisely controlled by the proposed propellantless propulsion method in the presence of external disturbances.     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

DRO计算及其在地月系中的摄动力研究

针对远距逆行轨道（DRO）的航天工程应用问题,研究了DRO的计算方法以及轨道特性,分析了DRO在实际力环境中的主要摄动因素,为DRO的精确建模和标称轨道设计奠定一定的理论基础。首先,利用仿真算例验证流函数法在计算DRO周期轨道族中的有效性。然后,利用该方法,通过改变雅可比常数,延拓计算DRO周期轨道族,获得不同共振比的DRO,仿真结果表明整数共振比的DRO在地月惯性坐标系中的轨迹是封闭的曲线,而共振比非整数的DRO则不封闭。最后,通过轨道外推分析影响DRO稳定性的主要摄动因素,仿真结果表明太阳引力和月球轨道偏心率是影响DRO稳定性的主要摄动因素。在动力学模型中,使用标准星历表示行星的运动状态,当积分时间多于10天时模型误差为km量级,因此在地月系这样大尺度的空间范围内,可以使用星历模型近似的分析DRO在真实力环境中的运动状态,为任务轨道设计提供依据。      

Improving magnetosphere in situ observations using solar sails

Past and current magnetosphere missions employ conventional spacecraft formations for in situ observations of the geomagnetic tail. Conventional spacecraft flying in inertially fixed Keplerian orbits are only aligned with the geomagnetic tail once per year, since the geomagnetic tail is always aligned with the Earth-Sun line, and therefore, rotates annually. Solar sails are able to artificially create sun-synchronous orbits such that the orbit apse line remains aligned with the geomagnetic tail line throughout the entire year. This continuous presence in the geomagnetic tail can significantly increase the science phase for magnetosphere missions. In this paper, the problem of solar sail formation design is explored using nonlinear programming to design optimal two-craft, triangle, and tetrahedron solar sail formations, in terms of formation quality and formation stability. The designed formations are directly compared to the formations used in NASA’s Magnetospheric Multi-Scale mission.     

Influence of curved thin-film device on deformation of a solar sail

A spinning solar sail IKAROS’s membrane is estimated to unexpectedly deform into an inverted pyramid shape due to thin-film devices with curvature, such as thin-film solar cells and steering devices on the membrane. It is important to investigate the deformation caused by the curved thin-film devices and predict the sail shape because the out-of-plane deformation greatly affects solar radiation pressure (SRP) and SRP torque. The purpose of this paper is to clarify the relationship between the global shape and orientation and position of curved thin-film devices and to evaluate SRP torque on the global shape using finite element analysis. The global shape is evaluated based on the out-of-plane displacement and the SRP torque. When the curved thin-film devices make the membrane shrink in the circumferential, diagonal, and radial direction, the sail deforms into a pyramid shape, an inverted pyramid one, and a saddle one, respectively. The saddle shape is more desirable for solar sails than the inverted pyramid shape and the pyramid one from the viewpoint of shape stability to SRP and control of SRP torque in the normal direction of the sail (windmill torque). The position of the thin-film device tends to increase the absolute value of windmill torque when it is biased circumferentially from the petal central axis. The suggested design principles for the arrangement of thin-film devices is that the curved thin-film devices should be directed so that the sail shrinks in the radial direction in order to deform the sail into a saddle shape with high shape stability, and the position of the thin-film devices should be biased in the circumferential direction paying attention to the absolute value of windmill torque to determine the direction of windmill torque.     

A torus-shaped solar sail accelerated via thermal desorption of coating

A torus-shaped sail consists of a reflective membrane attached to an inflatable torus-shaped rim. The sail’s deployment from its stowed configuration is initiated by introducing inflation pressure into the toroidal rim with an attached circular flat membrane coated by heat-sensitive materials that undergo thermal desorption (TD) from a solid to a gas phase. Our study of the deployment and acceleration of the sail is split into three steps: at a particular heliocentric distance a torus-shaped sail is deployed by a gas inflated into the toroidal rim and the membrane is kept flat by the pressure of the gas; under heating by solar radiation, the membrane coat undergoes TD and the sail is accelerated via TD of coating and solar radiation pressure (SRP); when TD ends, the sail utilizes thrust only from SRP. We study the stability of the torus-shaped sail and deflection and vibration of the flat membrane due to the acceleration by TD and SRP.     

地月L2点周期轨道的月球背面覆盖分析

地月L2点附近轨道具备独特的动力学和运动学特性,是月球背面探测任务的中继卫星首选布设位置。面向未来月球背面探测任务的中继通信需求,分析并研究了地月L2点周期轨道(halo轨道)对月球背面的覆盖。在圆型限制性三体问题模型下,研究并给出了halo轨道族延拓计算方法,基于延拓法设计了地月系大范围南北halo轨道族;给出了中继卫星的月球背面覆盖计算模型,定义了相应的时间覆盖因子;数值仿真了地月系南北halo轨道族的月球背面覆盖情况。研究结果表明:地月L2点周期轨道幅值和类型决定其对月面的覆盖性,幅值较小的轨道的月面整体覆盖性较好,幅值较大的轨道对月球南北极覆盖较好,南北族轨道分别有利于月球南北半球的覆盖。文章研究可为我国"嫦娥4号"月球背面探测任务的中继星轨道设计提供有益参考和借鉴。     

Sailing with solar and planetary radiation pressure

Literature on solar sailing has thus far mostly considered solar radiation pressure (SRP) as the only contribution to sail force. However, considering a sail in a planetary mission scenario, a new contribution can be added. Since the planet itself emits radiation, this generates a radial planetary radiation pressure (PRP) that is also exerted on the sail. Hence, this work studies the combined effects of both SRP and PRP on a sail for two case studies, i.e. Earth and Venus. In proximity of the Earth, the effect of PRP can be significant under specific conditions. Around Venus, instead, PRP is by far the dominating contribution. These combined effects have been studied for single- and double-sided reflective coating and including eclipse. Results show potential increase in the net acceleration and a change in the optimal attitude to maximise the acceleration in a given direction. Moreover, an increasing semi-major axis manoeuvre is shown with and without PRP, to quantify the difference on a real-case scenario.     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

Drag force on solar sails due to absorption of solar radiation

We consider a special relativistic effect, known as the Poynting–Robertson effect, on various types of trajectories of solar sails. Since this effect occurs at order v^?/c, where v^? is the transversal speed relative to the sun, it can dominate over other special relativistic effects, which occur at order v²/c². While solar radiation can be used to propel the solar sail, the absorbed portion of it also gives rise to a drag force in the transversal direction. For escape trajectories, this diminishes the cruising velocity, which can have a cumulative effect on the heliocentric distance. For a solar sail directly facing the sun in a bound orbit, the Poynting–Robertson effect decreases its orbital speed, thereby causing it to slowly spiral towards the sun. We also consider this effect for non-Keplerian orbits in which the solar sail is tilted in the azimuthal direction. While in principle the drag force could be counter-balanced by an extremely small tilt of the solar sail in the polar direction, periodic adjustments are more feasible.     

平动点周期轨道间小推力转移的Gauss伪谱法

针对三体问题共线平动点附近周期轨道间的小推力转移问题,构造了一种新的形状函数,在此基础上提出了一种基于Gauss伪谱法的优化设计方法。首先,建立小推力轨道转移动力学模型,参考初始轨道和目标轨道的类型,构造一种新的形状函数以近似小推力转移轨道。为满足不同的约束要求,提出了振幅和相位按多项式变化的假设,推导了小推力转移轨道的近似解析解;然后利用Gauss伪谱法将小推力轨道转移的最优控制问题转化为非线性规划问题,并对推导的近似解析解进行解算和处理,为Gauss伪谱法求解非线性规划问题提供较为有效的控制变量的初始猜测值;最后以地月系统L1点附近Halo轨道间的小推力转移问题为例进行了仿真分析。仿真结果表明,小推力转移轨道近似解析解具备有效性和普适性,使得Gauss伪谱法的迭代效率提高55%以上,同时也表明Gauss伪谱法可有效解决平动点周期轨道间的小推力转移轨道优化设计问题。     

Paths not taken – The Gossamer roadmap’s other options

Highly efficient low-thrust propulsion is increasingly applied beyond commercial use, also in mainstream and flagship science missions, in combination with gravity assist propulsion. Another recent development is the growth of small spacecraft solutions, not in size but in numbers and individual capabilities.Just over ten years ago, the DLR-ESTEC Gossamer Roadmap to Solar Sailing was set up to guide technology developments towards a propellant-less and highly efficient class of spacecraft for solar system exploration and applications missions: small spacecraft solar sails designed for carefree handling and equipped with carried application modules.Soon, in three dedicated Gossamer Roadmap Science Working Groups it initiated studies of missions uniquely feasible with solar sails such as Displaced L₁ (DL1) space weather advance warning and monitoring, Solar Polar Orbiter (SPO) delivery to very high inclination heliocentric orbit, and multiple Near-Earth Asteroid (NEA) rendezvous (MNR). Together, they demonstrate the capability of near-term solar sails to achieve at least in the inner solar system almost any kind of heliocentric orbit within 10 years, from the Earth-co-orbital to the extremely inclined, eccentric and even retrograde. Noted as part of the MNR study, sail-propelled head-on retrograde kinetic impactors (RKI) go to this extreme to achieve the highest possible specific kinetic energy for the deflection of hazardous asteroids.At DLR, the experience gained in the development of deployable membrane structures leading up to the successful ground deployment test of a (20 m)², i.e., 20 m by 20 m square solar sail at DLR Cologne in 1999 was revitalized and directed towards a 3-step small spacecraft development line from as-soon-as-possible sail deployment demonstration (Gossamer-1) via in-flight evaluation of sail attitude control actuators (Gossamer-2) to an envisaged proving-the-principle flight in the Earth-Moon system (Gossamer-3). First, it turned the concept of solar sail deployment on its head by introducing four separable Boom Sail Deployment Units (BSDU) to be discarded after deployment, enabling lightweight 3-axis stabilized sailcraft. By 2015, this effort culminated in the ground-qualified technology of the DLR Gossamer-1 deployment demonstrator Engineering Qualification Model (EQM). For mission types using separable payloads, such as SPO, MNR and RKI, design concepts can be derived from the BSDU characteristic of DLR Gossamer solar sail technology which share elements with the separation systems of asteroid nanolanders like MASCOT. These nano-spacecraft are an ideal match for solar sails in micro-spacecraft format whose launch configurations are compatible with ESPA and ASAP secondary payload platforms.Like any roadmap, this one contained much more than the planned route from departure to destination and the much shorter distance actually travelled. It is full of lanes, narrow and wide, detours and shortcuts, options and decision branches. Some became the path taken on which we previously reported. More were explored along the originally planned path or as new sidings in search of better options when circumstance changed and the project had to take another turn. But none were dead ends, they just faced the inevitable changes when roadmaps face realities and they were no longer part of the road ahead. To us, they were valuable lessons learned or options up our sleeves. But for future sailors they may be on their road ahead.     

Attitude stability criteria of axisymmetric solar sail

Passive attitude stability criteria of a solar sail whose membrane surface is axisymmetric are studied in this paper under a general SRP model. This paper proves that arbitrary attitude equilibrium position can be designed through adjusting the deviation between the pressure center and the mass center of the sail. The linearized method is applied to inspect analytically the stability of the equilibrium point from two different points of views. The results show that the attitude stability depends on the membrane surface shape and area. The results of simulation with full dynamic equations confirm that the two stability criteria are effective in judging the attitude stability for axisymmetric solar sail. Several possible applications of the study are also mentioned.     

地月系L2平动点轨道长期维持过程研究

针对地月系L2平动点附近两种常用轨道——拟halo轨道和Lissajous轨道,研究了轨道长期维持过程。采用穿越xz平面单圈控制的轨道维持策略,控制量通过单步预测法进行寻优。在全摄动动力学模型下对拟halo轨道和Lissajous轨道受控维持下的一年飞行过程分别进行了仿真和比较,结果显示Lissajous轨道的维持控制代价优于拟halo轨道,从轨道拓扑构型保持方面对该结果进行解释,并设计对比实验加以验证。研究结论对地月系L2平动点轨道选择、维持控制策略设计、长期飞行效果分析等方面具有工程参考意义。     

A scalable bus-based attitude control system for Solar Sails

Recently, there has been a renewed interest in Solar Sails as an alternative means of space propulsion. Many different attitude control systems have been designed for Solar Sails taking advantage of the centre-of-mass (CM)/centre-of-pressure (CP) offset while utilising the main sail structure to position the actuators. However, by attaching actuators to the main sail, these systems increase the risks involved in the deployment subsystem.