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.     

太阳帆飞行器轨道动力学分析

通过分析对轨道要素影响最大的加速度分量,使太阳帆总是位于光压力沿这个加速度方向的分量最大的方位。通过分析轨道要素调整时的相互影响关系,提出了同时修正多种轨道偏差的控制方案。对处于地球静止轨道上太阳帆飞行器的轨道调整进行了数值仿真。结果证明利用太阳光压力进行轨道调整是可行的,而且有利的太阳方位是进行快速有效的轨道调整的必要条件。     

On the station keeping of a solar sail in the elliptic Sun–Earth system

In this work we focus on the dynamics of a solar sail in the Sun–Earth Elliptic Restricted Three-Body Problem with solar radiation pressure. The considered situation is the motion of a sail close to the L₁ point, but displacing the equilibrium point with the sail so that it is possible to have continuous communication with the Earth. In previous works we derived a station keeping strategy for this situation but using the Circular RTBP as a model.     

Mission analysis and performance comparison for an Advanced Solar Photon Thruster

The so-called “compound solar sail”, also known as “Solar Photon Thruster” (SPT), is a design concept, for which the two basic functions of the solar sail, namely light collection and thrust direction, are uncoupled. In this paper, we introduce a novel SPT concept, termed the Advanced Solar Photon Thruster (ASPT), which does not suffer from the simplified assumptions that have been made for the analysis of compound solar sails in previous studies. After having presented the equations that describe the force on the ASPT and after having performed a detailed design analysis, the performance of the ASPT with respect to the conventional flat solar sail (FSS) is investigated for three interplanetary mission scenarios: an Earth–Venus rendezvous, where the solar sail has to spiral towards the Sun, an Earth–Mars rendezvous, where the solar sail has to spiral away from the Sun, and an Earth-NEA rendezvous (to near-Earth asteroid 1996FG3), where a large change in orbital eccentricity is required. The investigated solar sails have realistic near-term characteristic accelerations between 0.1 and 0.2 mm/s². Our results show that an SPT is not superior to the flat solar sail unless very idealistic assumptions are made.     

H-reversal trajectory: A new mission application for high-performance solar sails

The aim of this paper is to quantify the performance of a flat solar sail to perform a double angular momentum reversal maneuver and produce a new class of two-dimensional, non-Keplerian orbits in the ecliptic plane. For a given pair of orbital parameters, the orbital period and the perihelion distance, it is possible to find the minimum solar sail characteristic acceleration required to fulfil a double angular momentum reversal trajectory. This problem is addressed using an optimal formulation and is solved through an indirect approach. The new trajectories are symmetrical with respect to the sun-perihelion line and exhibit a bean-like shape. Two main difficulties must be properly taken into account. On one side the sail is required to perform a rapid reorientation maneuver when it approaches the perihelion. Suitable simulations have shown that such a maneuver is feasible. In the second place the new trajectories require the use of high performance solar sails. For example, assuming an orbital period equal to 5 years, the required solar sail characteristic acceleration is greater than 3.4 mm/s². Such a value, although beyond the currently available sail performance, is comparable to what is required by the original concept of H-reversal maneuvers introduced by Vulpetti in 1996.     

A simple control law for reducing the effective characteristic acceleration of a solar sail

The direction and magnitude of a solar sail acceleration are strongly related. For this reason, once the characteristic acceleration has been fixed, it is not possible to modulate the acceleration in a particular direction. In this work, a semi-analytical switching control law is derived, enabling a solar sail to emulate a smaller effective characteristic acceleration (without changes in geometry or optical properties); by periodically changing the pitch (cone) angle of the sail, in average over time, the acceleration produced by the sail matches exactly (in both direction and magnitude) that of a “smaller” sail. The range in which this is possible is determined, and the limitations on this range due to the size difference is computed. The method is validated on optimal Earth-Mars trajectories.     

Spinning solar sail orbit steering via spin rate control

The orbit of a solar sail can be controlled by changing the attitude of the spacecraft. In this study, we consider the spinning solar power sail IKAROS (Interplanetary Kite-craft Accelerated by Radiation Of the Sun), which is managed by Japan Aerospace Exploration Agency (JAXA). The IKAROS attitude, i.e., the direction of its spin-axis, is nominally controlled by the rhumb-line control method. By utilizing the solar radiation torque, however, we are able to change the direction of the spin-axis by only controlling its spin rate. With this spin rate control, we can also control indirectly the solar sail’s trajectory. The main objective of this study is to construct the orbit control strategy of the solar sail via the spin-rate control method. We evaluate this strategy in terms of its propellant consumption compared to the rhumb-line control method. Finally, we present the actual flight attitude data of IKAROS and the change of its trajectory.     

Variation of ionospheric total electron content in Taiwan region of the equatorial anomaly from 1994 to 2003

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionospheric anomaly (EIA) region is studied by analyzing dual-frequency signals of the Global Position System (GPS) acquired from a chain of nine observational sites clustered around Taiwan (21.9–26.2°N, 118.4–112.6°E). In this study, we present results from a statistical study of seasonal and geomagnetic effects on the EIA during solar cycle 23: 1994–2003. It is found that TEC at equatorial anomaly crests yield their maximum values during the vernal and autumnal months and their minimum values during the summer (except 1998). Using monthly averaged I_c (magnitude of TEC at the northern anomaly crest), semi-annual variations is seen clearly with two maxima occurring in both spring and autumn. In addition, I_c is found to be greater in winter than in summer. Statistically monthly values of I_c were poorly correlated with the monthly Dst index (r = −0.22) but were well correlated with the solar emission F_10.7 index (r = 0.87) for the entire database for the period during 1994–2003. In contrast, monthly values of I_c were correlated better with Dst (r ? 0.72) than with F_10.7 (r ? 0.56) in every year during the low solar activity period (1994–1997). It suggests that the effect of solar activity on I_c is a longer term (years), whereas the effect of geomagnetic activity on I_c is a shorter term (months).     

利用绳系太阳帆减缓小行星自转的技术研究

提出一种利用太阳帆绳系系统逐渐减小小行星自转速率的方法。在系绳长度不变的情况下,利用太阳帆受到的太阳光压力使其始终保持与小行星同步,避免了由于小行星自转而引起的系绳缠绕问题。通过控制太阳帆使系绳始终拉紧,系绳中的拉力便可以持续提供一个与小行星自转方向相反的力矩,从而减小其自转速率。仿真结果表明,面积10⁶ m²的太阳帆,经过约86天可将小行星的自转消除,验证了该方法的有效性。     

Status of solar sail technology within NASA

In the early 2000s, NASA made substantial progress in the development of solar sail propulsion systems for use in robotic science and exploration of the solar system. Two different 20-m solar sail systems were produced. NASA has successfully completed functional vacuum testing in their Glenn Research Center’s Space Power Facility at Plum Brook Station, Ohio. The sails were designed and developed by Alliant Techsystems Space Systems and L’Garde, respectively. The sail systems consist of a central structure with four deployable booms that support each sail. These sail designs are robust enough for deployment in a one-atmosphere, one-gravity environment and are scalable to much larger solar sails – perhaps as large as 150 m on a side. Computation modeling and analytical simulations were performed in order to assess the scalability of the technology to the larger sizes that are required to implement the first generation of missions using solar sails. Furthermore, life and space environmental effects testing of sail and component materials was also conducted.     

Numerical techniques for generating and refining solar sail trajectories

Like all applications in trajectory design, the design of solar sail trajectories requires a transition from analytical models to numerically generated realizations of an orbit. In astrodynamics, three numerical strategies are often employed. Differential correctors (also known as shooting methods) are perhaps the most common techniques. Finite-difference methods and collocation schemes are also employed and are successful in generating trajectories with pseudo-continuous control histories. These three numerical techniques are employed here to generate periodic trajectories displaced below the Moon in a circular restricted three-body system. All these approaches reveal trajectory options within the design space for solar sail applications.     

Influence of optical parameters on a solar sail motion

Detailed dynamic modeling of a solar sail requires recording of solar radiation pressure influence. A photon-solar sail is determined by the thrust value and the direction. We define the solar sail’s reflectivity depending on the film materials, the sail design and temperature, the thickness of multiple layers, and degradation factor, with a reasonable degree of accuracy. Thus, this work is devoted to the identification of optical characteristics of thin multilayer films in space flight conditions, i.e. to finding its reflectance, absorbance, and transmittance. In particular, the paper asks whether the solar sail simulates by a mathematical model of the optical characteristics of a multilayer epitaxial thin film. The temperature change effect and optical properties of solar sail degradation are considered as well. Solar sail flight from Earth to Mercury is designed as a simulation of the flight change in optical parameters.     

The investigation of powerful solar flares characteristics by analysis of excited states of C and various neutrons capture lines

Accelerated energetic particles in solar flares produced nuclear γ-lines in interactions with ambient solar atmosphere. Analysis of intensity of ratios between various γ-lines allows us to make estimations of abundance of elements, parameters of surrounding media and other solar characteristics. In this article we discuss the flux ratio between two lines from excited states of ¹²C (f_15.11/f_4.44) and our results of preliminary calculation of intensity ratio between two neutron capture lines at ³He and ¹H (f_20.58/f_2.223). In particular we consider the opportunity to obtain n(³He)/n(¹H) ratio during solar flares and using high-energy gamma-emission studying, based on the satellite data. Possible interpretation of spectral features observed during the January 20, 2005 solar flare is discussed. Preliminary analysis of energy spectrum in the band of 2–21 MeV gives n(³He)/n(¹H) ∼ 8 × 10⁻⁴ for January 20, 2005 solar flare.     

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.     

Changes in sea-level pressure over South Korea associated with high-speed solar wind events

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: first, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed V, daily SLP difference ΔSLP, and daily log(BV²) using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 km s⁻¹. We find that the daily ΔSLP averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ΔSLP is +2.5 hPa. The duration of deviation is a few days. We also find that ΔSLP is well correlated with both the speed of solar wind and log(BV²). The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are r ? 0.81 with P > 99.9%, and r ? 0.84 with P > 99.9%, respectively. We conclude by briefly discussing future direction to pursue.     

Solar sail near the Sun: Point-like and extended models of radiation source

Some modifications of solar sail radiation pressure forces on a plate and on a sphere for use in the numerical simulation of ‘local-optimal’ (or ‘instantaneously optimal’) trajectories of a spacecraft with a solar sail are suggested. The force model development is chronologically reviewed, including its connection with solar sail surface reflective and thermal properties. The sail surface is considered as partly absorbing, partly reflective (specular and diffuse), partly transparent. Thermal balance is specified because the spacecraft moves from circular Earth orbit to near-Sun regions and thermal limitations on the sail film are taken into account. A spherical sail-balloon can be used in near-Sun regions for scientific research beginning with the solar-synchronous orbit and moving outward from the Sun. The Sun is considered not only as a point-like source of radiation but also as an extended source of radiation which is assumed to be consequently as a point-like source of radiation, a uniformly bright flat solar disc and uniformly bright solar sphere.     

Particle acceleration and gamma rays in solar flares: Recent observations and new modeling

Experiments on SMM, GAMMA, Yohkoh, GRANAT, Compton GRO, INTEGRAL, RHESSI and CORONAS-F satellites over the past three decades have provided copious data for fundamental research relating to particle acceleration, transport and energetics of flares and to the ambient abundance of the solar corona, chromosphere and photosphere. We summarize main results of solar gamma-astronomy (including some results of several joint Russian–Chinese projects) and try to appraise critically a real contribution of those results into modern understanding of solar flares, particle acceleration at the Sun and some properties of the solar atmosphere. Recent findings based on the RHESSI, INTEGRAL and CORONAS-F measurements (source locations, spectrum peculiarities, ³He abundance etc.) are especially discussed. Some unusual features of extreme solar events (e.g., 28 October 2003 and 20 January 2005) have been found in gamma-ray production and generation of relativistic particles (solar cosmic rays, or SCR). A number of different plausible assumptions are considered concerning the details of underlying physical processes during large flares: (1) existence of a steeper distribution of surrounding medium density as compared to a standard astrophysical model (HSRA) for the solar atmosphere; (2) enhanced content of the ³He isotope; (3) formation of magnetic trap with specific properties; (4) prevailing non-uniform (e.g., fan-like) velocity (angular) distributions of secondary neutrons, etc. It is emphasized that real progress in this field may be achieved only by combination of gamma-ray data in different energy ranges with multi-wave and energetic particle observations during the same event. We especially note several promising lines for the further studies: (1) resonant acceleration of the ³He ions in the corona; (2) timing of the flare evolution by gamma-ray fluxes in energy range above 90 MeV; (3) separation of gamma-ray fluxes from different sources at/near the Sun (e.g., different acceleration sources/episodes during the same flare, contribution of energetic particles accelerated by the CME-driven shocks etc.); (4) asymmetric magnetic geometry and new magnetic topology models of the near-limb flares; (5) modeling of self-consistent time scenario of the event.     

The dynamics and control of the CubeSail mission: A solar sailing demonstration

The CubeSail mission is a low-cost demonstration of the UltraSail solar sailing concept (, ,  and ), using two near-identical CubeSat satellites to deploy a 260 m-long, 20 m² reflecting film. The two satellites are launched as a unit, detumbled, and separated, with the film unwinding symmetrically from motorized reels. The conformity to the CubeSat specification allows for reduction in launch costs as a secondary payload and utilization of the University of Illinois-developed spacecraft bus. The CubeSail demonstration is the first in a series of increasingly-complex missions aimed at validating several spacecraft subsystems, including attitude determination and control, the separation release unit, reel-based film deployment, as well as the dynamical behavior of the sail and on-orbit solar propulsion. The presented work describes dynamical behavior and control methods used during three main phases of the mission. The three phases include initial detumbling and stabilization using magnetic torque actuators, gravity-gradient-based deployment of the film, and steady-state film deformations in low Earth orbit in the presence of external forces of solar radiation pressure, aerodynamic drag, and gravity-gradient.     

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.     

Charging of dust grains by anisotropic solar wind multi-component plasma

In this paper we study the charging process of small grain particles by anisotropic multi-component solar wind plasmas (electrons, protons and heavy ions), versus two-component (electron/proton) plasmas. We are focusing attention on the important characteristics of the charging process, namely the charging time, floating potential and current content as functions of plasma parameters such as He⁺⁺/H⁺ (α/p) number density and T_α/T_p temperature ratios of alpha particles to protons, as well as plasma streaming velocity v₀. Measured statistical properties of solar wind plasma parameters at 1 AU show considerable variations in α/p-temperature ratios from 1 to 10, in α/p-number density ratio from 0.01 to 0.35, as well as in values of streaming velocity v₀ from 200 km/s to 1000 km/s and more. Periods of these variations could last for several days each, leading to significant variability in the charging process, according to newly derived general analytical expressions. Numerical calculations performed for protons/alphas plasmas showed large disparity in the charging characteristics. For example, in anisotropic plasma, grain charging time varies up to 90% depending on α/p-particles temperature and number density ratios, whereas changes in floating potential are up to 40%. In contrast, in isotropic plasma, charging characteristic for grains do not change very much for the same plasma parameters variations, with charging time varying about 12% and floating potential only varying about 4%. It is also shown that in highly anisotropic plasma, with all ballistic electrons and ions, dust grains could not hold their charges, and characteristic discharged time is calculated. We note that the analysis is equally applicable to any sized body immersed in solar wind plasma.