Application of a modified generalized sail model on a solar sail with wrinkles

In this paper, we present an analysis of effect of wrinkles on the solar sail performance. We describe different analytical, semi-analytical and numerical approaches to the calculation of general large-scale curvature of a solar sail as well as parameters of so-called wrinkled domains, and introduce the impact of such wrinkles on the thrust and torque of the solar sail. Finally, we present a model of an optically-orthotropic surface for such non-ideal sail, providing a connection with the Generalized Sail Model, and other solar sail thrust models.     

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.     

Solar sailing in realistic mode,based on a locally-optimal control laws

Ballistic design of solar sailing missions in the solar system is composed of defining the design parameters, the control programs, and the trajectories that provide performance goals of a flight. The use of a solar sail spacecraft imposes specific restrictions on mission parameters that include the degradation limit on the flight duration, the maximum temperature of solar sail's surface, the minimum distance from the Sun, the maximum angular velocity of the spacecraft's rotation and others.Many authors considered the impact of these restrictions on the design of the mission separately, but they used a sophisticated method of finding the exact optimal motion control or applied the most straightforward laws of motion control. This paper uses local-optimal control laws at the complete mathematical models of motion and functioning of solar sail spacecraft to describe a technique of designing interplanetary missions. The described method avoids the need to obtain an accurate optimal solution to the control problem and does not cause significant computational difficulties.     

Contingency target assessment,trajectory design,and analysis for NASA’s NEA scout solar sail mission

The presented study examines contingency target selection and trajectory design for NASA’s Near-Earth Asteroid Scout mission under the assumption of a missed lunar gravity assist. Two previously considered asteroids are selected as potential targets for the given scenario based on favorable orbital characteristics for launch dates ranging from June 27, 2020 through July 26, 2020. Initially, a simplified circular restricted 3-body problem + ideal solar sail model is utilized to survey trajectory options for a month-long launch window. Selected solutions from this data set are then converged in an N-body ephemeris + non-ideal sail model. Results suggest that NEA Scout can still perform asteroid rendezvous mission under the missed lunar gravity assist scenario with new targets, 2019 GF1, 2018 PK21, and 2007 UN12, based on the target launch dates. Further target assessment is carried out for 165 days beyond the current June 27, 2020 launch date.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

Statistical characteristics of low latitude Pc4s: Influence of solar wind and IMF parameters

Pc4 signatures for the year 2013, extracted from geomagnetic north–south and east–west components of induction coil magnetometer (LEMI 30) from low latitude station Desalpar (DSP), operated by Institute of Seismological Research (ISR), India have been investigated vis-à-vis the prevalent interplanetary parameters (IMF) as well as the geomagnetic activity indices. A clear dominance of Pc4-5 (467 events) over Pc3 (17 events) is observed. Local time variation of Pc4 shows a peak in the noon sector in both X and Y components. Our investigations show that the dominant peak frequency is 10 mHz at low latitude region. Correlations with solar wind and IMF parameters illustrate highest occurrence of Pc4 for a solar wind speed of 300–400 km/s and average IMF B field of 3–6 nT. The amplitude of Pc4s at DSP shows an increase with increasing solar wind speed, plasma density, solar wind dynamic pressure and average B field which is also reflected in the trend of frequency variation of these pulsations. We report that IMF clock angle at low latitude does not have influence on Pc4 occurrence. Based on the characteristics of these events, detected in latitudinally distributed stations from low and mid-latitudes from northern and southern hemisphere, we infer that modes were compressional, which could be driven by K-H instability or solar wind dynamic pressure, as compressional modes can propagate to low latitude with little attenuation.     

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.     

Plasma polarization electric field derived from radio sounding of solar wind acceleration region with spacecraft signals

Presented is the analytical approximation of averaged solar wind velocity radial dependence in the solar wind acceleration region at heliolatitudes below 60° under low and moderate solar activity. This empirical approximation is based on the data of radio sounding of the solar corona with radio signals from various spacecraft. Deduced is an equation connecting the solar wind velocity radial dependence and the radial dependence of solar wind plasma polarization electric field intensity. This allows constructing a semi-empirical radial dependence of plasma polarization electric field corresponding to the empirical radial dependence of solar wind velocity. Main properties of the semi-empirical dependence, which is based on radio sounding data, are described.     

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.