Determination of attitude motion of the Foton M-2 satellite using the data of onboard measurements of its angular velocity

We present the resutls of a prompt determination of the uncontrolled attitude motion of the Foton M-2 satellite, which was in orbit from May 31 to June 16, 2005. The data of onboard measurements of the angular velocity vector were used for this determination. The measurement sessions were carried out once a day, each lasting 83 min. Upon terminating a session, the data were transmitted to the ground to be processed using the least squares method and integrating the equations of motion of the satellite with respect to its center of mass. As a result of processing, the initial conditions of motion during a session were estimated, as well as parameters of the mathematical model used. The satellite’s actual motion is determined for 12 such sessions. The results obtained in flight completely described the satellite’s motion. This motion, having begun with a small angular velocity, gradually became faster, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 14, 2005 the angular velocity of the satellite with respect to its longitudinal axis was approximately 1.3 degrees per second, and the angular velocity projection onto a plane perpendicular to this axis had a magnitude of about 0.11 degrees per second. The results obtained are consistent with more precise results obtained later by processing the data on the Earth’s magnetic field measured on the same satellite, and they complement the latter in determination of the motion in the concluding segment of the flight, when no magnetic measurements were performed.     

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.     

Topological model of the solar event including a flare and coronal mass ejection on October 19, 2001

Based on the analysis of a strong solar flare X1.6/2B on October 19, 2001 in the active region 9661, accompanied by a coronal mass ejection (CME) of the halo type, a topological model of development of this solar event is suggested. The model considers a unified process of development of CME and a chromosphere flare. According to the model, this process has a common source of energy supply: the turbulent current layer lying between the arcade of flare loops and the surface of CME going away. The structures on the ends of flare bands (SEFB) represent in this model chromosphere feet of the system of large-scale coronal magnetic arches at the initial stage of the dynamic processes whose evolution results in CME. Peripheral structures (PS) of the flare (elongated double bright emission strips beyond the limits of the active region) are interpreted as chromosphere bases of magnetic field lines that form an external shell (braid) of the CME at the late stage of the flare.     

Features of coronal mass ejections of minimum size

The coronal mass ejections (CME) with small angular dimensions (d ≤ 10°) have the simplest form, much simpler than large CME. This fact simplifies the problem of analyzing the CME structure and studying their origin. On the basis of the analysis of the LASCO C2 (SOHO) data, we show in this paper that the motion of a CME having small dimensions proceeds within a magnetic tube (a ray with increased brightness) of the streamer belt and leads to an “explosion-like” increase in the angular dimensions (rapid expansion) of the tube. A hypothesis is put forward that a small CME represents a “plasmoid” (a plasma bunch bounded in space, with its own magnetic field) thrown into the base of the magnetic tube and moving along it away from the Sun.     

Modes of uncontrolled rotational motion of the <Emphasis Type="Italic">Progress M-29M</Emphasis> spacecraft

We have reconstructed the uncontrolled rotational motion of the Progress M-29M transport cargo spacecraft in the single-axis solar orientation mode (the so-called sunward spin) and in the mode of the gravitational orientation of a rotating satellite. The modes were implemented on April 3–7, 2016 as a part of preparation for experiments with the DAKON convection sensor onboard the Progress spacecraft. The reconstruction was performed by integral statistical techniques using the measurements of the spacecraft’s angular velocity and electric current from its solar arrays. The measurement data obtained in a certain time interval have been jointly processed using the least-squares method by integrating the equations of the spacecraft’s motion relative to the center of mass. As a result of processing, the initial conditions of motion and parameters of the mathematical model have been estimated. The motion in the sunward spin mode is the rotation of the spacecraft with an angular velocity of 2.2 deg/s about the normal to the plane of solar arrays; the normal is oriented toward the Sun or forms a small angle with this direction. The duration of the mode is several orbit passes. The reconstruction has been performed over time intervals of up to 1 h. As a result, the actual rotational motion of the spacecraft relative to the Earth–Sun direction was obtained. In the gravitational orientation mode, the spacecraft was rotated about its longitudinal axis with an angular velocity of 0.1–0.2 deg/s; the longitudinal axis executed small oscillated relative to the local vertical. The reconstruction of motion relative to the orbital coordinate system was performed in time intervals of up to 7 h using only the angularvelocity measurements. The measurements of the electric current from solar arrays were used for verification.     

Coronal mass ejection (CME) activity of low mass M stars as an important factor for the habitability of terrestrial exoplanets. I. CME impact on expected magnetospheres of Earth-like exoplanets in close-in habitable zones

Low mass M- and K-type stars are much more numerous in the solar neighborhood than solar-like G-type stars. Therefore, some of them may appear as interesting candidates for the target star lists of terrestrial exoplanet (i.e., planets with mass, radius, and internal parameters identical to Earth) search programs like Darwin (ESA) or the Terrestrial Planet Finder Coronagraph/Inferometer (NASA). The higher level of stellar activity of low mass M stars, as compared to solar-like G stars, as well as the closer orbital distances of their habitable zones (HZs), means that terrestrial-type exoplanets within HZs of these stars are more influenced by stellar activity than one would expect for a planet in an HZ of a solar-like star. Here we examine the influences of stellar coronal mass ejection (CME) activity on planetary environments and the role CMEs may play in the definition of habitability criterion for the terrestrial type exoplanets near M stars. We pay attention to the fact that exoplanets within HZs that are in close proximity to low mass M stars may become tidally locked, which, in turn, can result in relatively weak intrinsic planetary magnetic moments. Taking into account existing observational data and models that involve the Sun and related hypothetical parameters of extrasolar CMEs (density, velocity, size, and occurrence rate), we show that Earth-like exoplanets within close-in HZs should experience a continuous CME exposure over long periods of time. This fact, together with small magnetic moments of tidally locked exoplanets, may result in little or no magnetospheric protection of planetary atmospheres from a dense flow of CME plasma. Magnetospheric standoff distances of weakly magnetized Earth-like exoplanets at orbital distances 相似文献   

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.     

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.     

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.     

Statistical analysis of attitude motion of a light capsule entering the atmosphere

The laws of distribution of the space angle of attack are analyzed for a light descent capsule on the atmosphere’s conventional boundary and on the segment of motion in dense layers of the atmosphere until the moment of reaching the maximum velocity head. It is assumed that upon detachment from a base spacecraft the angular velocity components of the descent capsule represent independent random quantities distributed according to the normal law. The shape of the descent capsule is close to a body of revolution whose instantaneous aerodynamic characteristic has a sufficiently simple form (one stable and one unstable positions of static equilibrium). It is demonstrated that there is a possibility to approximate the distribution densities by well-known laws and by approximate functions constructed on the basis of simplified models. Evolution of the distribution laws at increasing mass-inertia asymmetry of the descent capsule is studied.     

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.     

A study of angular motion of the Chibis-M microsatellite with three-axis flywheel control

We study the controlled angular motion of the Chibis-M microsatellite. Executive elements are three pairs of flywheels, whose axes are mutually perpendicular. The task of the control system is realization of a required program motion and support of its asymptotic stability. In this paper, we synthesize a control algorithm and study the evolution of the angular momentum of flywheels on long time intervals. The attitude accuracy is estimated for the case when disturbances act upon the spacecraft.     

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.     

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.     

A specific case of stability of cylindrical precession of a satellite

In a central Newtonian gravitational field, the motion of a dynamically symmetrical satellite along an elliptical orbit of arbitrary eccentricity is considered. The particular motion of the satellite is known when its axis of symmetry is perpendicular to the orbit plane, and the satellite rotates about this axis with a constant angular velocity (cylindrical precession). A nonlinear analysis of stability of this motion has been performed under the assumption that the geometry of the satellite mass corresponds to a thin plate. At small values of orbit eccentricity e the analysis is analytical, while numerical analysis is used for arbitrary values of e.     

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.     

Experimental investigation of the modes of operation of uncontrolled attitude motion of the Progress spacecraft

Results of in-flight tests of three modes of uncontrolled attitude motion of the Progress spacecraft are described. These proposed modes of experiments related to microgravity are as follows: (1) triaxial gravitational orientation, (2) gravitational orientation of the rotating satellite, and (3) spin-up in the plane of the orbit around the axis of the maximum moment of inertia. The tests were carried out from May 24 to June 1, 2004 onboard the spacecraft Progress M1-11. The actual motion of this spacecraft with respect to its center of mass, in the above-mentioned modes, was determined by telemetric information about an electric current tapped off from solar batteries. The values of the current obtained during a time interval of several hours were processed jointly using the least squares method by integration of the equations of the spacecraft’s attitude motion. The processing resulted in estimation of the initial conditions of motion and of the parameters of mathematical models used. For the obtained motions the quasi-static component of microaccelerations was computed at a point onboard, where installation of experimental equipment is possible.     

皮卫星分离参数优化

为实现“浙大皮星二号”卫星无干涉分离,达到分离初始姿态要求,对皮卫星分离机构进行了设计、理论分析及优化。通过对分离系统动力学特性分析,得出了影响该分离系统的关键因素,同时对星箭分离机构运动系统进行了优化,并进行了试验验证。仿真分析及试验结果表明：该设计能够完全满足分离机构分离速度和角速度要求,星箭分离过程的仿真分析和试验校验了“浙大皮星二号”星箭分离机构可实现无干涉分离,皮卫星初始速度、角速度均满足所提出的各项初始分离姿态要求。     

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.     

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°.