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.     

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.     

Effect of Microaccelerations on the Distribution of a Dopant in a Semiconductor Melt during Space Flight

The effect of residual microaccelerations on the distribution of a dopant in a semiconductor melt located in a heated closed cavity onboard an Earth-orbiting satellite is considered in the context of a model problem of thermal convection. The amplitude–frequency characteristics of the response of this distribution to the perturbing microaccelerations are obtained. It is demonstrated that the effect of low-frequency microaccelerations decreases when the frequency increases. A comparison is made of the macroscopic inhomogeneities of the dopant concentration due to the actual low-frequency (quasi-static) component of microaccelerations onboard different spacecraft: the orbital station Mir, the satellite Foton-11, a Space Shuttle orbiter, and the International Space Station. A substantial effect of the rotational motion of the spacecraft on the character of the time behavior of a macroscopic inhomogeneity is demonstrated.     

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.     

Quasistatic Microaccelerations onboard the Foton Satellite during Its Flight at High-Altitude Orbits

The level of quasistatic microaccelerations onboard the Foton-M satellite is predicted for its flights in two orbits: the planned orbit with the altitudes in perigee h = 262 km and in apogee h = 304 km and the orbit with h = 262 km and h = 350 km. The prediction is based on mathematical simulation of the satellite motion with respect to its center of mass under the action of gravitational and aerodynamic moments. The model is represented by the system of equations of the satellite rotational motion. Parameters of this system are chosen from the condition of coincidence of the motion of preceding Foton satellites (h 220 km and h 400 km) calculated using this model with the results of determination of actual rotational motion of the Foton-11 and Foton-12 satellites. With the help of the model thus calibrated, a calculation is made of the rotational motion of the Foton-M satellite and of the quasistatic microaccelerations onboard it. As is shown by the results of simulation, the use of the first and the second orbits will result in reductions of microaccelerations by 30% and 60%, respectively.     

Uncontrolled Attitude Motion of the Foton-12 Satellite and Quasi-Steady Microaccelerations onboard It

The results of determination of the uncontrolled attitude motion of the Foton-12 satellite (placed in orbit on September 9, 1999, terminated its flight on September 24, 1999) are presented. The determination was carried out by the onboard measurement data of the Earth's magnetic field strength vector. Intervals with a duration of several hours were selected from data covering almost the entire flight. On each such interval the data were processed simultaneously using the least squares method by integrating the satellite's equations of motion with respect to the center of mass. The initial conditions of motion and the parameters of the mathematical model employed were estimated in processing. The results obtained provided for a complete representation of the satellite's motion during the flight. This motion, beginning with a small angular velocity, gradually sped up. The growth of the component of the angular velocity with respect to the longitudinal axis of the satellite was particularly strong. During the first several days of the flight this component increased virtually after every passage through the orbit's perigee. As the satellite's angular velocity increased, its motion became more and more similar to the regular Euler precession of an axisymmetric rigid body. In the last several days of flight the satellite's angular velocity with respect to its longitudinal axis was about 1 deg/s and the projection of the angular velocity onto the plane perpendicular to this axis had a magnitude of approximately 0.15 deg/s. The deviation of the longitudinal axis from the normal to the orbit plane did not exceed 60°. The knowledge of the attitude motion of the satellite allowed us to determine the quasi-steady microacceleration component onboard it at the locations of the technological and scientific equipment.     

An analysis of low-frequency component in microacceleration measurements made onboard the <Emphasis Type="Italic">Foton M-2</Emphasis> satellite

The low-frequency component is investigated in the data of measurements performed onboard the Foton M-2 satellite with the three-component accelerometer TAS-3. Investigations consisted in comparison of this component with its calculated analog found from a reconstruction of the satellite’s attitude motion. The influence of the Earth’s magnetic field on the accelerometer readings is discovered by way of spectral analysis of the functions representing the results of determining the low-frequency microacceleration by two methods. After making correction for this influence, the results obtained by these two methods coincided within a root-mean-square error of less than 10^?6 m/s².     

The Influence of Arrangement of Growth Setups onboard a Spacecraft on Microgravity Conditions of Experiments: An Example of Floating Zone Melting of InSb:Te onboard the Foton-3 Satellite

The results of experiments on the growth of InSb:Te by the floating zone melting in the Zona-4 setup during a flight of the Foton-3 satellite are discussed in comparison with the data on the microgravity situation typical for satellites of this type. When analyzing inhomogeneities in the crystals obtained, we reveal the frequencies corresponding to periodical variations of the impurity channel position (the facet effect), to interleaving of packs in which the growth layers are grouped, and to location of subgroups of the growth layers inside the packs. These frequencies are close to those discovered during measurements of low-frequency (quasistationary) microaccelerations (g) onboard the spacecraft Foton-8, Foton-10, Foton-11, and Foton-12. Calculated values of g at the place where the Zona-4 setup is installed confirm the possibility of impact of such g on the heat and mass transfer in the melt.     

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.     

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.     

Determination of the Attitude Motion of a Satellite Based on a Kinematical Model of Motion and Using Data of Measurements of Its Angular Velocity and the Earth’s Magnetic Field Strength

An integral statistical procedure of determination of the attitude motion of a satellite using the data of onboard measurements of angular velocity vectors and the strength of the Earth’s magnetic field (EMF) is suggested. The procedure uses only the equations of kinematics of a solid body and is applicable to determining both controlled and uncontrollable motions of a satellite at any external mechanical moments acting upon it. When applying this procedure, the data of measurements of both types, accumulated during a certain interval of time, are processed jointly. The data of measuring the angular velocity are smoothed by discrete Fourier series, and these series are substituted into kinematical Poisson equations for elements of the matrix of transition from a satellite-fixed coordinate system to the orbital coordinate system. The equations thus obtained represent a kinematical model of the satellite motion. The solution to these equations (which approximate the actual motion of a satellite) is found from the condition of the best (in the sense of the least squares method) fit of the data of measuring the EMF strength vector to its calculated values. The results of testing the suggested procedure using the data of measurements of the angular velocity vectors onboard the Foton-12 satellite and measurements of EMF strengths are presented.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 295–305.Original Russian Text Copyright © 2005 by Abrashkin, Volkov, Voronov, Egorov, Kazakova, Pankratov, Sazonov, Semkin.     

The use of refined model of aerodynamic moment in the problem of reconstruction of the <Emphasis Type="Italic">Foton</Emphasis> satellite rotational motion

A new mathematical model of the uncontrolled rotational motion of the Foton satellite is presented. The model is based on the Euler dynamic equations of rigid body motion and takes into account the action upon the satellite of four external mechanical moments: gravitational, restoring aerodynamic, moment with constant components in the satellite-fixed coordinate system, and moment arising due to interaction of the Earth’s magnetic field with the satellite’s proper magnetic moment. To calculate the aerodynamic moment a special geometrical model of the outer satellite shell is used. Detailed form of the formulas giving above-mentioned moments in the equations of satellite motion is agreed with the form of the considered motion. Model testing is performed by determining with its help the rotational motion of the Foton M-2 satellite (it was in orbit from May 31, 2005 to June 16, 2005) using the data of the onboard measurements of the Earth’s magnetic field strength. The use of the new model has led to a relatively small improvement in the accuracy of the motion determination, but allowed us to obtain physically real estimates of some parameters.     

Determination of attitude motion of the <Emphasis Type="Italic">Foton M-3</Emphasis> satellite according to the data of onboard measurements of the Earth’s magnetic field

The results of reconstruction of rotational motion of the Foton M-3 satellite during its uncontrolled flight in September 2007 are presented. The reconstruction was performed by processing the data of onboard measurements of the Earth’s magnetic field obtained by the DIMAC instruments. The measurements were carried out continuously throughout the flight, but the processing technique dealt with the data portions covering time intervals of a few orbital revolutions. The data obtained on each such interval were processed jointly by the least squares method with using integration of the equations of satellite motion relative to its center of mass. When processing, the initial conditions of motion and the used mathematical model’s parameters were estimated. The results of processing 16 data sets gave us complete information about the satellite motion. This motion, which began at a low angular velocity, had gradually accelerated and in five days became close to the regular Euler precession of an axisymmetric solid body. At the end of uncontrolled flight the angular velocity of the satellite relative to its lengthwise axis was 0.5 deg/s; the angular velocity projection onto the plane perpendicular to this axis had a magnitude of about 0.18 deg/s.     

Experiments with the DAKON-M convection sensor

We have described an express technique for processing the results of experiments with a DAKON-M convection sensor on board the Service Module of the International Space Station (ISS) in 2011. The technique uses a certain rule to compare the sensor measurements with the calculated data on the quasistatic component of microacceleration at the point of installation. The sensor experiments have been conducted during shuttle docking and undocking, when low-frequency microaccelerations on the ISS were significant. The microaccelerations have been calculated using measurement data of the MAMS low-frequency accelerometer installed in the Lab module and the telemetry data on the ISS rotational motion. This has made it possible to convert the MAMS measurement data to the DAKON-M convection sensor installation point. A comparison of sensor readings with calculated microaccelerations has revealed fairly good agreement between them.     

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.     

Estimating the accuracy of the technique of reconstructing the rotational motion of a satellite based on the measurements of its angular velocity and the magnetic field of the Earth

The paper has studied the accuracy of the technique that allows the rotational motion of the Earth artificial satellites (AES) to be reconstructed based on the data of onboard measurements of angular velocity vectors and the strength of the Earth magnetic field (EMF). The technique is based on kinematic equations of the rotational motion of a rigid body. Both types of measurement data collected over some time interval have been processed jointly. The angular velocity measurements have been approximated using convenient formulas, which are substituted into the kinematic differential equations for the quaternion that specifies the transition from the body-fixed coordinate system of a satellite to the inertial coordinate system. Thus obtained equations represent a kinematic model of the rotational motion of a satellite. The solution of these equations, which approximate real motion, has been found by the least-square method from the condition of best fitting between the data of measurements of the EMF strength vector and its calculated values. The accuracy of the technique has been estimated by processing the data obtained from the board of the service module of the International Space Station (ISS). The reconstruction of station motion using the aforementioned technique has been compared with the telemetry data on the actual motion of the station. The technique has allowed us to reconstruct the station motion in the orbital orientation mode with a maximum error less than 0.6° and the turns with a maximal error of less than 1.2°.     

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.     

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.     

Determination of parameters of attitude motion of a small communication satellite using the data of measurements of current of solar panels

A communication satellite (small spacecraft) injected into a geosynchronous orbit is considered. Flywheel engines are used to control the rotational spacecraft motion. The spacecraft after the emergency situation has passed into a state of uncontrolled rotation. In this case, no direct telemetric information about parameters of its rotational motion was accessible. As a result, the problem arose to determine the rotational satellite motion according to the available indirect information: current taken from the solar panels. Telemetric measurements of solar panel current obtained on the time interval of a few hours were simultaneously processed by the least squares method integrating the equations of rotational satellite motion. We present the results of processing 10 intervals of the measurement data allowing one to determine the real rotational spacecraft motion and to estimate the total angular momentum of flywheel engines.     

Study of the mode of angular velocity damping for a spacecraft at non-standard situation

Non-standard situation on a spacecraft (Earth’s satellite) is considered, when there are no measurements of the spacecraft’s angular velocity component relative to one of its body axes. Angular velocity measurements are used in controlling spacecraft’s attitude motion by means of flywheels. The arising problem is to study the operation of standard control algorithms in the absence of some necessary measurements. In this work this problem is solved for the algorithm ensuring the damping of spacecraft’s angular velocity. Such a damping is shown to be possible not for all initial conditions of motion. In the general case one of two possible final modes is realized, each described by stable steady-state solutions of the equations of motion. In one of them, the spacecraft’s angular velocity component relative to the axis, for which the measurements are absent, is nonzero. The estimates of the regions of attraction are obtained for these steady-state solutions by numerical calculations. A simple technique is suggested that allows one to eliminate the initial conditions of the angular velocity damping mode from the attraction region of an undesirable solution. Several realizations of this mode that have taken place are reconstructed. This reconstruction was carried out using approximations of telemetry values of the angular velocity components and the total angular momentum of flywheels, obtained at the non-standard situation, by solutions of the equations of spacecraft’s rotational motion.