Reconstruction of spacecraft rotational motion using a Kalman filter

Quasi-static microaccelerations of four satellites of the Foton series (nos. 11, 12, M-2, M-3) were monitored as follows. First, according to measurements of onboard sensors obtained in a certain time interval, spacecraft rotational motion was reconstructed in this interval. Then, along the found motion, microacceleration at a given onboard point was calculated according to the known formula as a function of time. The motion was reconstructed by the least squares method using the solutions to the equations of satellite rotational motion. The time intervals in which these equations make reconstruction possible were from one to five orbital revolutions. This length is increased with the modulus of the satellite angular velocity. To get an idea on microaccelerations and satellite motion during an entire flight, the motion was reconstructed in several tens of such intervals. This paper proposes a method for motion reconstruction suitable for an interval of arbitrary length. The method is based on the Kalman filter. We preliminary describe a new version of the method for reconstructing uncontrolled satellite rotational motion from magnetic measurements using the least squares method, which is essentially used to construct the Kalman filter. The results of comparison of both methods are presented using the data obtained on a flight of the Foton M-3.     

Estimation of dynamic characteristics of the <Emphasis Type="Italic">International Space Station</Emphasis> from measurements of microaccelerations

We describe the results of determining the mass of the International Space Station using the data of MAMS accelerometer taken during correction of the station orbit on August 20, 2004. The correction was made by approach and attitude control engines (ACE) of the Progress transporting spacecraft. The engines were preliminary calibrated in an autonomous flight using the onboard device for measuring apparent velocity increment. The method of calibration is described and its results are presented. The error in station mass determination is about 1%. The same data of MAMS and similar data obtained during the orbit correction on August 26, 2004 were used for the analysis of high-frequency vibrations of the station mainframe caused by operation of the ACE of Progress. Natural frequencies of the ACE are determined. They lie in the frequency band 0.024–0.11 Hz. ACE operation is demonstrated to result in a substantial increase of microaccelerations onboard the station in the frequency range 0–1 Hz. The frequencies are indicated at which disturbances increase by more than an order of magnitude. The study described was carried out as a part of the Tensor technological experiment.     

Observations of corotating solar wind structures at radio sounding by signals of the <Emphasis Type="Italic">Rosetta</Emphasis> and <Emphasis Type="Italic">Mars Express</Emphasis> spacecraft

In the implementation of the space projects Rosetta and Mars Express, a large-scale series of experiments has been carried out on radio sounding circumsolar plasma by decimeter (S-band) and centimeter (X-band) signals of the Rosetta comet probe (from October 3 to October 31, 2010) and the Mars Express satellite of Mars (from December 25, 2010 to March 27, 2011). It was found that in the phase of ingress the spacecraft behind the Sun, the intensity of the frequency fluctuations increases in accordance with a power function whose argument is the solar offset distance of radio ray path, and when the spacecraft is removed from the Sun (the egress phase), frequency fluctuations are reduced. Periodic strong increases in the fluctuation level, exceeding by a factor of 3–12 the background values of this value determined by the regular radial dependences, are imposed on the regular dependences. It was found that increasing the fluctuations of radio waves alternates with the periodicity m × T or n × T, where m = 1/2, n = 1, аnd T is the synodic period of the Sun’s rotation (T ≈ 27 days). It was shown that the corotating structures associated with the interaction regions of different speed fluxes are formed in the area of solar wind acceleration and at distances of 6–20 solar radii already have a quasi-stationary character.     

On the possibility of long-term maintenance of geostationary satellite attitude without using the external information sensors and inertial sensors

The problem of maintaining spacecraft attitude during the loss of information from attitude sensors and inertial sensors is considered. The problem is solved on the basis of the force gyro stabilization principle with producing an angular momentum in the plane of orbit. The redundant mode of attitude maintenance is developed for spacecraft of the Yamal series. The results of testing the mode during the in-flight tests of the Yamal-200 spacecraft are presented.     

Study of microaccelerations onboard the <Emphasis Type="Italic">International Space Station</Emphasis> with the DAKON-M convection sensor

The results of experiments with the DAKON-M convection sensor onboard the Russian orbital segment of the International Space Station are described. A comparison of the sensor measurements with the results of calculation of the quasistatic microacceleration component at the point of installation is made. For this comparison we have used three measurement intervals of the experiments in 2009, during which spacecraft were docked with the station, undocked from it, and actuation of jet engines of the attitude control system took place. When calculating microacceleration, we use the measurement data of the low-frequency MAMS accelerometer, installed on the American segment, and the telemetry data on the ISS rotational motion. This information allowed one to convert the MAMS measurements to the point of installation of the DAKON-M convection sensor. A comparison of sensor measurements with calculated microaccelerations showed sufficiently accurate coincidence between the calculated and measured data.     

Turbulent fluctuations of plasma and magnetic field parameters in the magnetosheath and the low-latitude boundary layer formation: Multisatellite observations on March 2, 1996

The results of simultaneous analysis of plasma and magnetic field characteristics measured on the INTERBALL/Tail Probe, WIND and Geotail satellites on March 2, 1996, are presented. During these observations the INTERBALL/Tail Probe crossed the low-latitude boundary layer, and the WIND and Geotail satellites measured the solar wind’s and magnetosheath’s parameters, respectively. The plasma and magnetic field characteristics in these regions have been compared. The data of the Corall, Electron, and MIF instruments on the INTERBALL/Tail Probe satellite are analyzed. Fluctuations of the magnetic field components and plasma velocity in the solar wind and magnetosheath, measured onboard the WIND and Geotail satellites, are compared. The causes resulting in appearance of plasma jet flows in the low-latitude boundary layer are analyzed. The amplitude of magnetic field fluctuations in the magnetosheath for a studied magnetosphere boundary crossing is shown to exceed the magnetic field value below the magnetopause near the cusp. The possibility of local violation of pressure balance on the magnetopause is discussed, as well as penetration of magnetosheath plasma into the magnetosphere, as a result of magnetic field and plasma flux fluctuations in the magnetosheath.     

Sсheme of rendezvous mission to lunar orbital station by spacecraft launched from Earth

In recent years, great experience has been accumulated in manned flight astronautics for rendezvous in near-Earth orbit. During flights of Apollo spacecraft with crews that landed on the surface of the Moon, the problem of docking a landing module launched from the Moon’s surface with the Apollo spacecraft’s command module in a circumlunar orbit was successfully solved. A return to the Moon declared by leading space agencies requires a scheme for rendezvous of a spacecraft launched from an earth-based cosmodromee with a lunar orbital station. This paper considers some ballistic schemes making it possible to solve this problem with minimum fuel expenditures.     

Fluxes of Quasi-trapped Electrons with Energies >0.08 MeV in the Near-Earth Space on Drift Shells <Emphasis Type="Italic">L</Emphasis> < 2

We study the characteristics of fluxes of electrons with energy >80 keV in the near-Earth space regions corresponding to the drift shells L = 1.7, 1.4, and 1.1 observed during the entire period of the GRIF experiment onboard the Spectr module of the Mir orbital station from October 1995 to June 1997. The obtained geographic maps of the distribution of electron fluxes at the height of the station flight (400 km) and, also, the estimates of the spectra indicate that the South-Atlantic Anomaly provides for a mechanism of stable replenishment for shells with L < 1.5. The mechanism of stable replenishment of shells with L < 1.5 may be due to the scattering, in the residual atmosphere, of electrons from the inner radiation belt precipitating into the region of the South-Atlantic Anomaly.     

Study of notches in the Earth’s plasmasphere based on data of the <Emphasis Type="Italic">MAGION-5</Emphasis> satellite

Depleted narrow (localized in longitude) regions (field tubes) in the plasmasphere, recently discovered in He⁺ radiation measurements on the IMAGE spacecraft, were first directly observed by the Magion-5 satellite. The low-density regions (notches) occupy <～ 10–30° in longitude and extend from L ～ 2–3 to the plasmasphere boundary in neighboring plasmasphere regions with larger densities. The Magion-5 data give evidence that in the low-density regions temperature is enhanced as compared to the neighboring denser plasmasphere regions. Formation of notches in the plasmasphere is, apparently, associated with AE intensification during weak magnetic storms, while strong magnetic storms usually result in the overall reduction of plasmasphere dimensions. However, even a strong magnetic storm on April 6–7, 2000 (max K _p = 9-and min D _st ～ ?290 nT), but accompanied by an isolated AE impulse, resulted in a density decrease only in the longitudinally limited post-midnight sector of the plasmasphere.     

Electron flux variations at altitudes of 600–800 km in the second half of 2014. preliminary results of an experiment using RELEC equipment onboard the satellite <Emphasis Type="Italic">VERNOV</Emphasis>

The results of measurements of fluxes and spectra carried out using the RELEC (relativistic electrons) equipment onboard the VERNOV satellite in the second half of 2014 are presented. The VERNOV satellite was launched on July 8, 2014 in a sun-synchronous orbit with an altitude from 640 to 830 km and an inclination of 98.4°. Scientific information from the satellite was first received on July 20, 2014. The comparative analysis of electron fluxes using data from RELEC and using experimental data on the electron detection by satellites Elektro-L (positioned at a geostationary orbit) and Meteor-M no. 2 (positioned at a circular polar orbit at an altitude of about 800 km as the VERNOV satellite) will make it possible to study the spatial distribution pattern of energetic electrons in near-Earth space in more detail.     

Families of periodic solutions in three-dimensional restricted three-body problem

Planar orbits of three-dimensional restricted circular three-body problem are considered as a special case of three-dimensional orbits, and the second-order monodromy matrices M (in coordinate z and velocity v _z ) are calculated for them. Semi-trace s of matrix M determines vertical stability of an orbit. If |s| ≤ 1, then transformation of the subspace (z, v _z ) in the neighborhood of solution for the period is reduced to deformation and a rotation through angle φ, cosφ = s. If the angle ? can be rationally expressed through 2π,φ = 2π·p/q, where p and q are integer, then a planar orbit generates the families of three-dimensional periodic solutions that have a period larger by a factor of q (second kind Poincareé periodic solutions). Directions of continuation in the subspace (z, v _z ) are determined by matrix M. If |s| < 1, we have two new families, while only one exists at resonances 1: 1 (s = 1) and 2: 1 (s = ?1). In the course of motion along the family of three-dimensional periodic solutions, a transition is possible from one family of planar solutions to another one, sometimes previously unknown family of planar solutions.     

Peculiarities of long-wave radio bursts from solar flares preceding strong geomagnetic storms

Radio bursts in the frequency range of 100–1500 kHz, recorded in 1997–2000 on the INTERBALL-1 satellite during the solar flares preceding the strong geomagnetic storms with D _st < ?100 nT, are analyzed in this paper. The observed long-wave III-type radio bursts of solar origin at frequencies of 1460 and 780 kHz were characterized by large values of the flux S _f = 10^?15 ?10^?17 W/m² Hz and duration longer than 10 min. The rapid frequency drift of a modulated radio burst continued up to a frequency of 250 kHz, which testified that the exciting agent (a beam of energetic electrons) propagated from the Sun to the Earth. All such flares were characterized by the appearance of halo coronal mass ejections, observed by the LASCO/SOHO, and by the presence of a southward B_z-component of the IMF, measured on the ACE and WIND spacecraft. In addition, shortly after radio bursts, the INTERBALL-1 satellite has recorded the fluxes of energetic electrons with E > 40 keV.     

Observations of electromagnetic emissions inside the Earth’s plasmasphere from the <Emphasis Type="Italic">Interball-1</Emphasis> satellite

The spectrum analyzer AKR-X onboard the Interball-1 satellite at the beginning (August–October 1995) and at the end (August–October 2000) of satellite operation in perigees of its orbital motion recorded and analyzed electromagnetic emissions of the inner regions of the Earth’s plasmasphere in the frequency band 100–1500 kHz at distances of 1.1–1.8 R _E. The observations have shown that the electromagnetic modes (the Z and LO modes escaping the magnetosphere) which are formed at the altitudes 600–4000 km are associated with the subauroral nonthermal continuum and with the recently discovered kilometric continuum. There are noticeable differences in the spectral character of these emissions during the minimum (1996) and maximum (2000) solar activity, when, as a rule, the LO mode escaping the plasmaphere and the continua are not present.     

Forecasting the velocity of quasi-stationary solar wind and the intensity of geomagnetic disturbances produced by it

A brief review is given of contemporary approaches to solving the problem of medium-term forecast of the velocity of quasi-stationary solar wind (SW) and of the intensity of geomagnetic disturbances caused by it. At the present time, two promising models of calculating the velocity of quasi-stationary SW at the Earth’s orbit are realized. One model is the semi-empirical model of Wang-Sheeley-Arge (WSA) which allows one to calculate the dependence V(t) of SW velocity at the Earth’s orbit using measured values of the photospheric magnetic field. This model is based on calculation of the local divergence f _S of magnetic field lines. The second model is semi-empirical model by Eselevich-Fainshtein-Rudenko (EFR). It is based on calculation in a potential approximation of the area of foot points on the solar surface of open magnetic tubes (sources of fast quasistationary SW). The new Bd-technology is used in these calculations, allowing one to calculate instantaneous distributions of the magnetic field above the entire visible surface of the Sun. Using predicted V(t) profiles, one can in EFR model calculate also the intensity of geomagnetic disturbances caused by quasi-stationary SW. This intensity is expressed through the K _p index. In this paper the EFR model is discussed in detail. Some examples of epignosis and real forecast of V(t) and K _p (t) are discussed. A comparison of the results of applying these two models for the SW velocity forecasting is presented.     

Study of singular surfaces of systems of unsupported gyrodynes by the parameter continuation method

A system of two-degree-of-freedom force gyroscopes (gyrodynes) is considered to be used for spacecraft attitude control. Possible values of its total angular momentum form some finite region P in the frame of reference rigidly connected with the spacecraft. Near the boundary of this region and singular surfaces located inside it the control of the angular momentum is complicated or impossible. The program angular momentum of the gyrodynes, realizing the law of variation of the spacecraft orientation, should lie inside P and outside some singular surfaces, and due to this fact the boundary and internal singular surfaces should be studied. This work is dedicated to the numerical construction of region P and its internal singular surfaces by the method of parameter continuation. Using the results by E.N. Tokar’ we formalize sufficient conditions which in some cases allow one to determine the type of the singular surface. As an example, a system of six gyrodynes is considered, for which the regions of variations of the intrinsic angular momentum and singular surfaces are constructed. The possibilities of the system are demonstrated for the case when one gyrodyne fails.     

Designing corrections for the trajectory of the <Emphasis Type="Italic">Spektr-R</Emphasis> spacecraft in the event of immersions into the Moon’s sphere of influence

The results of updating the parameters of motion of the Spektr-R spacecraft at the end of 2016 have shown that, in January 2018, with a probability close to unity, the condition that a spacecraft stay in the Earth’s shadow is violated; however, in May of the same year, the ballistic life of the spacecraft will be terminated. Thus, in 2017, the question arose of how to design the correction of flight of this spacecraft using its onboard propulsion system. The correction was designed with allowance for the fact that, for the first time since it was launched, the spacecraft in the course of several years, beginning with 2017, repeatedly approaches the Moon, deeply immersing into its sphere of influence. This paper presents the technologically and organizationally convenient, allowable versions of upcoming correction of the Spektr-R spacecraft trajectory and justifies the particular scheme of its implementation.     

Solar activity index for long-term ionospheric forecasts

Based on the comparison of solar activity indices (annual average values of the relative number of sunspots Rz₁₂ and solar radio emission flux at a wavelength of 10.7 cm F₁₂) with the ionospheric index of solar activity IG₁₂ for 1954–2013, we have found that the index F₁₂ is a more accurate (than Rz₁₂) indicator of solar activity for the long-term forecast of foF2 (the critical frequency of the F2-layer). This advantage of the F₁₂ index becomes especially significant after 2000 if the specific features of extreme ultraviolet radiation of the Sun are additionally taken into account in the minima of solar cycles, using an appropriate correction to F₁₂. Qualitative arguments are given in favor of the use of F₁₂ for the long-term forecast of both foF2 and other ionospheric parameters.     

A method of calculation of vertical cutoff rigidity in the geomagnetic field

On the basis of generalization of the results of extensive trajectory calculations for trial charged particles moving in the geomagnetic field the method of calculation of effective vertical cutoff rigidity, taking into account the values of K _p -index and local time, is developed. The IGRF and Tsyganenko-89 models are used for the geomagnetic field. A comparison of the results of model simulations with the experimental data on penetration of charged particles into near-Earth space is made, and penetration functions for typical spacecraft orbits are calculated.     

A Concept of Guidance for an Air-Launch Vehicle with Compensation of Initial Downrange and Launch Time Errors at Direct Insertion into a Rendezvous Point in the Orbit

A flying launcher (airplane carrier) can generate initial errors in position and time of launch. In order to compensate for these errors, one should have two control parameters in addition to those that provide for a spacecraft's injection into a preset orbit. We suggest the concept of controlling the trajectory of injection by choosing thrust values (within allowable regions of control) of second-stage engines or/and of a space booster of the Polyot carrier launcher. As an example, a rendezvous of the spacecraft at the end of its boost phase with the International Space Station (ISS) is considered. The methodology of the suggested approach can be extended to other mobile systems of launch to rendezvous orbits.     

A simplified technique for determining the rotational motion of a satellite based on the onboard measurements of the angular velocity and magnetic field of the Earth

The mathematical model, which allowed us to reconstruct the rotational motion of the Bion M-1 and Foton M-4 satellites by processing the measurements of onboard magnetometers and the angular velocity sensor, is sufficiently detailed and accurate. If we slightly lower the requirements for accuracy and transfer to a rougher model, i.e., we will not update the biases in measurements of the angular velocity component, then the measurement processing technique can be significantly simplified. The volume of calculations in minimizing the functional of the least-square technique is reduced; the most complicated part of calculations is performed using the standard procedure of computational linear algebra. This simplified technique is described below, and the examples of its application for reconstructing the rotational motion of the Foton M-4 satellite are presented. A noticeable distinction in the reconstructions of motion, constructed by simplified and more exact techniques, is revealed in processing the measurements over time intervals longer than 4 hours.