Determination of <Emphasis Type="Italic">Foton M-2</Emphasis> satellite attitude motion by the data of microacceleration measurements

The results of reconstruction of uncontrolled attitude motion of the Foton M-2 satellite using measurements with the accelerometer TAS-3 are presented. The attitude motion of this satellite has been previously determined by the measurement data of the Earth’s magnetic field and the angular velocity. The TAS-3 data for this purpose are used for the first time. These data contain a well-pronounced additional component which made impossible their direct employment for the reconstruction of the attitude motion and whose origin was unknown several years ago. Later it has become known that the additional component is caused by the influence of the Earth’s magnetic field. The disclosure of this fact allowed us to take into account a necessary correction in processing of TAS-3 data and to use them for the reconstruction of the attitude motion of Foton M-2. Here, a modified method of processing TAS-3 data is described, as well as results of its testing and employing. The testing consisted in the direct comparison of the motion reconstructed by the new method with the motion constructed by the magnetic measurements. The new method allowed us to find the actual motion of Foton M-2 in the period June 9, 2005–June 14, 2005, when no magnetic measurements were carried out.     

Analysis of Low-Frequency Microaccelerations onboard the Foton-11Satellite

The method and the results of investigating the low-frequency component of microaccelerations onboard the Foton-11satellite are presented. The investigation was based on the processing of data of the angular velocity measurements made by the German system QSAM, as well as the data of measurements of microaccelerations performed by the QSAM system and by the French accelerometer BETA. The processing was carried out in the following manner. A low-frequency (frequencies less than 0.01 Hz) component was selected from the data of measurements of each component of the angular velocity vector or of the microacceleration, and an approximation was constructed of the obtained vector function by a similar function that was calculated along the solutions to the differential equations of motion of the satellite with respect to its center of mass. The construction was carried out by the least squares method. The initial conditions of the satellite motion, its aerodynamic parameters, and constant biases in the measurement data were used as fitting parameters. The time intervals on which the approximation was constructed were from one to five hours long. The processing of the measurements performed with three different instruments produced sufficiently close results. It turned out to be that the rotational motion of the satellite during nearly the entire flight was close to the regular Eulerian precession of the axially symmetric rigid body. The angular velocity of the satellite with respect to its longitudinal axis was about 1 deg/s, while the projection of the angular velocity onto the plane perpendicular to this axis had an absolute value of about 0.2 deg/s. The magnitude of the quasistatic component of microaccelerations in the locations of the accelerometers QSAM and BETA did not exceed 5 × 10^–5–10^–4m/s²for the considered motion of the satellite.     

An Analysis of the Low-Frequency Component in Measurements of Angular Velocity and Microacceleration onboard the Foton-12 Satellite

We compare the results of two methods used to determine the angular velocity of the Foton-12 satellite and the low-frequency component of microaccelerations onboard it. The first method is based on reconstruction of the satellite's rotational motion using the data of onboard measurements of the strength of the Earth's magnetic field. The motion (time dependence of the orientation parameters and angular velocity) was found from the condition of best approximation of the measurement data by the functions calculated along the solutions to equations of attitude motion of the satellite. The solutions found were used to calculate the quasistatic component of microaccelerations at certain points of the satellite, in particular, at the point of location of an accelerometer of the QSAM system. Filtration of the low-frequency component of the angular velocity and microacceleration from the data of measurements by a sensor of angular velocity and by the accelerometer of this system served as a second method. The filtration was made using the discrete Fourier series. A spectral analysis of the functions representing the results of determining the angular velocity and microacceleration by both methods is performed. Comparing the frequencies and amplitudes of the harmonic component of these functions allowed us to estimate the accuracy of measurements made by the QSAM system in the low-frequency range.     

Low-Frequency Microaccelerations onboard the Foton-11 Satellite

We analyze the microacceleration measurements carried out onboard the Foton-11 satellite with the three-component accelerometer BETA. The microaccelerations were recorded virtually throughout the entire orbital flight of the Foton-11 satellite. The data obtained were analyzed in the following way. First they were used to determine the actual rotational motion of the satellite for several arbitrarily selected time intervals 4 h long. This problem was solved by constructing the approximation of the microacceleretation low-frequency component (previously determined from the data) by its calculated analog computed along the solutions to differential equations of rotational motion of the satellite. The approximation was made by the least squares method. As a result, those mathematical model parameters and the solutions to equations of motion were found that gave the best consistency of the microacceleretation low-frequency component and its calculated analog. Then the spectral analysis of the low-frequency component and its calculated analog was made. It was shown that, although basic harmonics of these functions coincided sufficiently well, some harmonics of the low-frequency component failed to be interpreted in terms of the satellite's rotational motion.     

Processing the data of measurements of angular velocity and microaccelerations onboard the <Emphasis Type="Italic">Foton-12</Emphasis> satellite

The results of reconstruction of uncontrolled rotational motion of the Foton-12 satellite using the measurement data of onboard sensors are presented. This problem has already been solved successfully several years ago. The satellite motion was reconstructed using the data of measuring the Earth’s magnetic field. The data of measuring the angular velocity and microaccelerations by the QSAM system were actually not used for this purpose, since these data include a clearly seen additional component whose origin was at that time unclear. This component prevented one from using these data directly for reconstruction of the angular motion. Later it became clear that the additional component was caused by the Earth’s magnetic field. Discovery of this fact allowed us to make necessary corrections when processing the QSAM system data and to use them for reconstruction of rotational motion of Foton-12. Below, a modified method of processing the QSAM system data is described together with the results of its application. The main result is obtained by comparing the motion reconstructed from measurements of angular velocity or acceleration with that found by way of processing the measurements of the Earth’s magnetic field. Their coincidence turned out to be rather accurate.     

Uncontrolled motion of the <Emphasis Type="Italic">Foton M-2</Emphasis> satellite and quasistatic microaccelerations on its board

The results of determining the uncontrolled rotational motion of the Foton M-2 satellite (in orbit from May 31 to June 16, 2005) are presented. The determination was made using the data of onboard measurements of the Earth’s magnetic field strength. Segments 270 min long (three orbits) were selected from these data covering the first two thirds of the flight. On each such segment the data were processed jointly by the least squares method using integration of the equations of attitude motion of the satellite. In processing, the initial conditions of motion and parameters of the used mathematical model were estimated. The thus obtained results gave a complete overview of the satellite motion. This motion, having started with a small angular velocity, gradually accelerated, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 09, 2005 (the last day of measurements) the angular velocity of the satellite relative to its lengthwise axis was about 1.1 deg/s, while the projection of the angular velocity onto a plane perpendicular to this axis had an absolute value of about 0.11 deg/s. Deviations of the lengthwise axis from a normal to the orbit plane did not exceed 60°. Based on the results of determination of the rotational motion of the satellite, calculations of quasi-static microaccelerations on its board are performed.     

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.     

Observation of Terrestrial gamma-ray flashes in the RELEC space experiment on the <Emphasis Type="Italic">Vernov</Emphasis> satellite

The RELEС scientific payload of the Vernov satellite launched on July 8, 2014 includes the DRGE spectrometer of gamma-rays and electrons. This instrument comprises a set of scintillator phoswich-detectors, including four identical X-ray and gamma-ray detector with an energy range of 10 kev to 3 MeV with a total area of ~500 cm² directed to the atmosphere, as well as an electron spectrometer containing three mutually orthogonal detector units with a geometric factor of ~2 cm² sr. The aim of a space experiment with the DRGE instrument is the study of fast phenomena, in particular Terrestrial gamma-ray flashes (TGF) and magnetospheric electron precipitation. In this regard, the instrument provides the transmission of both monitoring data with a time resolution of 1 s, and data in the event-by-event mode, with a recording of the time of detection of each gamma quantum or electron to an accuracy of ~15 μs. This makes it possible to not only conduct a detailed analysis of the variability in the gamma-ray range, but also compare the time profiles with the results of measurements with other RELEC instruments (the detector of optical and ultraviolet flares, radio-frequency and low-frequency analyzers of electromagnetic field parameters), as well as with the data of ground-based facility for thunderstorm activity. This paper presents the first catalog of Terrestrial gamma-ray flashes. The criterion for selecting flashes required in order to detect no less than 5 hard quanta in 1 ms by at least two independent detectors. The TGFs included in the catalog have a typical duration of ~400 μs, during which 10–40 gamma-ray quanta were detected. The time profiles, spectral parameters, and geographic position, as well as a result of a comparison with the output data of other Vernov instruments, are presented for each of candidates. The candidate for Terrestrial gamma-ray flashes detected in the near-polar region over Antarctica is discussed.     

Determination of characteristic frequencies of elastic oscillations of the <Emphasis Type="Italic">International Space Station</Emphasis> construction

The results of investigating free oscillations of the International Space Station construction appearing during spacecraft docking and undocking are described. The study is carried out using the measurement data of the low-frequency MAMS accelerometer. Several intervals of measurements performed in 2005 and 2006 were chosen to be studied. For chosen intervals, only the data intervals corresponding to the process of free attenuation of the oscillations construction elements were analyzed. Characteristic frequencies of elastic oscillations of the station construction and attenuation coefficients corresponding to them are found. The comparative analysis of the results obtained for various docking ports (nodes) is carried out. The described study is performed as a part of the technical experiment “The ISS Environment” carried out onboard the station in accordance with the Russian program of scientific and engineering experiments.     

Rotational motion of <Emphasis Type="Italic">Foton M-4</Emphasis>

The actual controlled rotational motion of the Foton M-4 satellite is reconstructed for the mode of single-axis solar orientation. The reconstruction was carried out using data of onboard measurements of vectors of angular velocity and the strength of the Earth’s magnetic field. The reconstruction method is based on the reconstruction of the kinematic equations of the rotational motion of a solid body. According to the method, measurement data of both types collected at a certain time interval are processed together. Measurements of the angular velocity are interpolated by piecewise-linear functions, which are substituted in kinematic differential equations for a quaternion that defines the transition from the satellite instrument coordinate system to the inertial coordinate system. The obtained equations represent the kinematic model of the satellite rotational motion. A solution of these equations that approximates the actual motion is derived from the condition of the best (in the sense of the least squares method) match between the measurement data of the strength vector of the Earth’s magnetic field and its calculated values. The described method makes it possible to reconstruct the actual rotational satellite motion using one solution of kinematic equations over time intervals longer than 10 h. The found reconstructions have been used to calculate the residual microaccelerations.     

Attitude control system of the first Russian nanosatellite <Emphasis Type="Italic">TNS-0 no. 1</Emphasis>

The results of designing the attitude control system of the first Russian nanosatellite TNS-0 no. 1 providing orientation of its longitudinal axis along the local geomagnetic field induction vector are presented. The system consists of a permanent magnet and two sets of hysteresis rods. The magnetic and geometric parameters of the magnet and rods are calculated. The influence of the permanent magnet field on the hysteresis rods and mutual influence of the rods in the case of compact satellite packaging are analyzed. Examples of calculations of transient processes and steady-state angular satellite motion are presented.     

Design of satellite constellations with continuous coverage on elliptic orbits of <Emphasis Type="Italic">Molniya</Emphasis> type

New methods of choosing the structures of satellite constellations (SC) on elliptical orbits of the Molniya type are presented. The methods, using critical inclination and putting the orbit apogee in the Earth’s hemisphere with an area of continuous coverage, are based on geometrical analysis of two-dimensional representation of the coverage conditions and SC motion in the space of inertial longitude of the orbit ascending node and time. The coverage conditions are represented in the form of a certain region. Dynamics of all satellites in this space is represented by uniform motion along a straight line approximately parallel to the ordinate axis, while the satellite system forms a grid. The problem of choosing a minimal (as far as the number of satellites is concerned) SC configuration can be formulated as a search for the most sparse grid. The contemporary advanced methods of computational geometry serve as an algorithmic basis for the problem solution. Design of SC for continuous coverage of latitude belts with the use of kinematically regular systems is considered. A method of analyzing single-track systems for continuous coverage of arbitrary geographic regions is described, which makes a region at any time instant observable by at least one satellite of the system. As an example, SC on elliptical orbits are considered with periods of ～4, 12, and 24 hours.     

Monitoring of the time and spatial distribution of neutron-flux spectral density outside the Russian segment of the <Emphasis Type="Italic">International Space Station</Emphasis> based on data from the BTN-Neutron space experiment

Results of measurements of neutron-flux spectral density in the vicinity of the International Space Station (ISS) based on BTN-Neutron space experimental data acquired in 2007–2014 have been presented in this paper. It has been shown that, during the flight of the ISS over different regions of the Earth’s surface, neutron flux in the energy range of 0.4 eV–15 MeV varies from 0.1 n/sm²/s in equatorial regions to 50 n/sm²/s in the South Atlantic anomaly region. The measurements were used to estimate the contribution of the neutron component to the overall exposure dose rate. The total contribution of fast neutrons is about 0.1–0.4 μ Zv/h above the equator area and more than 50 μ Zv/h above the South Atlantic anomaly region. A data analysis of BTN-Neutron data also showed that the time profile of neutron flux has long-periodic variations. It was found that, under the influence of Galactic cosmic rays (GCRs), modulation during 24th solar cycle neutron flux changed almost twofold (above high latitude regions). Maximum values of neutron flux were observed in January 2010 and minimum values were observed in January 2014.     

Measurement of the Quasistatic Component of Microaccelerations Using a Convection Sensor onboard a Satellite

The problem of the interpretation of measurements made by means of a convection sensor is considered. The sensor is a cubic chamber filled by a viscous fluid (gas). Fixed and unequal temperatures are maintained on two opposite sides of the cube; the other sides are perfect heat conductors. Two differential thermocouples are placed inside the chamber to measure the temperature difference at two pairs of fixed points. The sensor is mounted aboard the Earth's satellite. Mathematical models of various degrees of complexity are proposed which describe processes of heat and mass transfer under the action of a quasistatic component of microaccelerations. The results of mathematical simulation of the data of sensor thermocouples presenting a response to the real quasistatic component of microaccelerations which took place aboard the Mirstation are given. It is shown that under usual conditions of an orbital mission the sensor presents a linear low-frequency filter. By combining the data of several identical sensors, tightly arranged and oriented in a certain way, it is possible to measure low-frequency components of the angular acceleration of the satellite and linear microaccelerations at the point of the sensor position.     

Nanosatellite <Emphasis Type="Italic">REFLECTOR</Emphasis>: Choice of parameters of the attitude control system

We consider methodological and applied problems of choosing parameters and developing a passive gravity-gradient attitude control system for the REFLECTOR nanosatellite. The attitude control system consists of appropriately distributed satellite mass and additional masses forming the required tensor of inertia, as well as of a set of hysteresis rods manufactured from a magnetically soft material. The combination of these elements allows one, independent of the satellite’s initial angular motion, to maintain its three-axial orientation in the orbital coordinate system with preset performance in precision and fast response. The satellite’s dynamics is investigated using asymptotic and numerical methods. The problems of arrangement of the additional masses and hysteresis rods are discussed, as well as their interaction with magnetizable elements.     

静电加速度计对重力卫星平台的要求分析

静电加速度计是重力卫星的核心有效载荷,其非保守力测量精度是地球重力场反演任务成功执行的有效保证。文章分析了静电加速度计的测量信号以及对重力卫星平台的指标要求,重点指出满足这些要求应解决的关键技术,最后建议我国开展重力卫星研究。     

Uncontrolled rotational motion of the <Emphasis Type="Italic">Aist</Emphasis> small spacecraft prototype

The results of reconstructing the uncontrolled rotational motion of the Aist small spacecraft prototype during its flight in early 2014 have been presented. The reconstruction was carried out by processing data from onboard measurements of the Earth’s magnetic field. The processing procedure used portions of data covering intervals of time with durations ranging from a few dozen minutes to three hours. Data obtained in each such interval were processed jointly by the least-squares method by integrating the equations of the satellite motion relative to the center of mass. The initial conditions of the motion and the parameters of the used mathematical model during processing have been estimated. The results of processing for several data intervals have provided a fairly complete picture of the satellite motion. This was the weakly disturbed Euler–Poinsot motion.     

Single-impulse transfers from orbits of artificial satellites to orbits around <Emphasis Type="Italic">L</Emphasis><Subscript>1</Subscript> or <Emphasis Type="Italic">L</Emphasis><Subscript>2</Subscript> libration points

In the context of the restricted circular three-body problem a method for constructing families of periodic orbits is described. Each orbit contains a segment of transfer from artificial satellite orbit of a smaller body to an orbit around L ₁ or L ₂ points of the Sun-Earth and Earth-Moon systems, a segment of multiple flyby of this libration point, and a segment of return to the artificial satellite orbit. Dependences of velocities at the pericenter on the pericenter radius are given.     

Modelling and simulation of the space mission MICROSCOPE

MICROSCOPE is a French space mission for testing the weak equivalence principle (WEP). The mission goal is the determination of the Eötvös parameter $\eta$ with an accuracy of 10^?15. The French space agency CNES is responsible for the satellite which is developed and produced within the Myriade series. The satellite's payload T-SAGE (Twin Space Accelerometer for Gravitation Experimentation) is developed and built by the French institute ONERA. It consists of two high-precision capacitive differential accelerometers. One accelerometer is used as reference sensor with two test masses of platinum, the science sensor contains a platinum and a titanium proof mass. The detection of the test mass movement and their control is done via a complex electrode system. As a member of the MICROSCOPE performance team, the German department ZARM will be involved in the data analysis of the MICROSCOPE mission. For this purpose, mission simulations and the preparation of the mission data evaluation in close cooperation with the French partners CNES, ONERA and OCA are realised. The development status of the simulation tool which will represent the complex spacecraft dynamics and all error sources in order to design and test data reduction procedures is presented and some features are discussed in detail.     

Longitudinal dependence of the H<Superscript>+</Superscript> concentration distribution in the plasmasphere according to <Emphasis Type="Italic">INTERBALL-1</Emphasis> satellite data

The paper has presented a study of the dependence of the H⁺ ions concentration in the plasmasphere on geographic longitude. A vast database of measurements of the cold plasma density by the Alpha-3 instrument on board the INTERBALL-1 satellite has been used for the study. Based on these measurements, a dependence of the H⁺ ions concentration in the filled magnetic flux tube in the plasmasphere in the equatorial plane under quiet geomagnetic conditions has been obtained as a function of geographic longitude. Studies have been performed for two seasons, summer and winter. It has been shown that, during the summer in the near-midnight sector, the minimum in the H⁺ concentration falls within geographic longitudes of 270°–315°. The ratio of the concentration of H⁺ ions at various longitudes could reach a factor of three. During the winter, in the near-noon sector, the maximum of the H⁺ ions concentration falls within longitudes of 180°–225°, whereas the concentration ratio could reach a factor of 2.2.