Pressing Issues of the Day in Studying Deep Space

Two problems in studying deep space are discussed that are, in the author's opinion, the most important. The first is soil sampling from the smaller bodies of the Solar System, such as the Martian satellite Phobos and asteroids of groups C and S of the Main Asteroid Belt. This soil (so-called primordial substance) can help to elucidate some problems of the Solar System's formation; in particular, to construct a reliable model of the internal structure of the Earth. The second problem is to reveal all sufficiently large asteroids penetrating inside the Earth's orbit and to catalog those asteroids that are hazardous from the viewpoint of collision with the Earth. To this end, it is suggested to launch five or six Earth-orbiting spacecraft with telescopes capable of recording objects down to a brightness of 22– 25 ^m . It is pointed out that both problems can be solved in the near future using comparatively cheap standardized space vehicles launched into near-Earth orbits by the Soyuz carrier rocket and boosted further by electro-jet engines of small thrust.     

Regularization of the Perturbed Spatial Restricted Three-Body Problem by <Emphasis Type="Italic">L</Emphasis>-Transformations

Equations of motion for the perturbed circular restricted three-body problem have been regularized in canonical variables in a moving coordinate system. Two different L-matrices of the fourth order are used in the regularization. Conditions for generalized symplecticity of the constructed transform have been checked. In the unperturbed case, the regular equations have a polynomial structure. The regular equations have been numerically integrated using the Runge–Kutta–Fehlberg method. The results of numerical experiments are given for the Earth–Moon system parameters taking into account the perturbation of the Sun for different L-matrices.     

地月平动点拟周期轨道探测器自主导航方法研究

地月L2点的拟周期轨道可以用于实现与月球背面的持续通信,具有重要的科学研究价值和广阔的应用前景。针对地月L2点探测器所处的弱稳定拟周期轨道,论证了基于日地月信息的自主导航方法的可行性,并进行了深入分析。首先,推导了会合坐标系下带有星历的精确导航动力学方程;其次,针对弱稳定轨道不同于近地强稳定轨道的特性,在基于日地月方位信息导航的基础上,提出了三种敏感器组合方案。使用迭代最小二乘方法给出导航仿真结果,并结合非线性可辨识性理论对这三种情况下历元状态的可辨识性及可辨识度进行分析。最后,仿真结果表明,使用日地月敏感器以及对地多普勒雷达可以满足历元状态的可辨识性、导航的收敛性以及系统经济性的要求。     

从月球逃逸探测小行星的发射机会搜索

从月球逃逸探测小行星的发射机会搜索因需考虑日、地、月引力的影响而使问题变得复杂。针对该多体系统的发射机会搜索问题,提出了一种分层渐近的搜索方法。该方法首先通过分析地月系质心与小行星的几何关系,搜索从地月系质心到小行星的发射机会,进而以地月运动为研究对象,推导出了从月球轨道切向逃逸机会的判别条件,并基于此判别条件及等高线图法对逃逸机会进行了搜索。同时,为提高所得发射机会在多体模型下的轨道修正收敛性,给出了基于月心逃逸轨道参数为终端约束的日-地与日-地-月动力学模型的轨道渐近修正方法。最后,以近地小行星(3908)Nyx和(190491)2000 FJ20为例,搜索其从月球逃逸的发射机会,仿真计算表明了该方法的有效性。     

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.     

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.     

嫦娥五号平动点拓展任务轨道方案研究

中国嫦娥五号探测器成功实现月球采样返回任务,为最大限度利用任务资源,研究了利用嫦娥五号轨道器的平动点拓展任务轨道方案,设计了平动点轨道及其转移轨道.首先,给出了任务轨道设计的轨道动力学模型,包括圆型限制性三体问题模型和精确力模型.其次,基于嫦娥二号和嫦娥5T1平动点拓展任务设计经验,介绍了平动点轨道直接转移与入轨等轨道...     

Calculation and Study of Limited Orbits around the <Emphasis Type="Italic">L</Emphasis><Subscript>2</Subscript> Libration Point of the Sun–Earth System

A procedure has been proposed for calculating limited orbits around the L₂ libration points of the Sun–Earth system. The motion of a spacecraft in the vicinity of the libration point has been considered a superposition of three components, i.e., decreasing (stable), increasing (unstable), and limited. The proposed procedure makes it possible to correct the state vector of the spacecraft so as to neutralize the unstable component of the motion. Using this procedure, the calculation of orbits around various types of libration points has been carried out and the dependence on the orbit type on the initial conditions has been studied.     

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.     

Overcoming Genesis mission design challenges

Genesis is the fifth mission selected as part of NASA's Discovery Program. The objective of Genesis is to collect solar wind samples for a period of approximately 2 years while in a halo orbit about the Sun–Earth colinear libration point, L1, located between the Sun and Earth. At the end of this period, the spacecraft follows a free-return trajectory with the samples delivered to a specific recovery point on the Earth for subsequent analysis. This type of sample return has never been attempted before and presents a formidable challenge, particularly with regard to planning and execution of propulsive maneuvers. Moreover, since the original inception, additional challenges have arisen as a result of emerging spacecraft design concerns and operational constraints. This paper will describe how these challenges have been met to date in the context of the better-faster-cheaper paradigm. [This paper addresses an earlier mission design, as of May 2000.]     

Synchronized SETI-the case for "opposition"

If the signals being sought in search for extraterrestrial intelligence (SETI) programs exist but are brief (for example, they are produced intermittently to conserve energy), then it is essential to know when these signals will arrive at the Earth. Different types of transmitter/receiver synchronization schemes are possible, which vary in the relative amount of effort required by the transmitter and the receiver. The case is made for a scheme that is extremely simple for the receiver: Make observations of a target when it is at maximum angular distance from the Sun (i.e., "opposition"). This strategy requires that the transmitter has accurate knowledge of the distance and proper motion of the Sun and the orbit of the Earth. It is anticipated that within the next 10-20 years it will be possible to detect directly nearby extrasolar planets of approximately terrestrial mass. Since extraterrestrial transmitters are expected to have significantly more advanced technology, it is not unreasonable to expect that they would be able to detect the presence of the Earth and measure its orbit at even greater distances. This strategy is simple to implement, and opposition is also typically the time when observations are easiest to make. Limited opposition surveys contained in a number of all-sky surveys have already been performed. However, full-sky opposition surveys are best suited to detectors with very large fields of view.     

On the Precession of Saturn

The precession of Saturn under the effect of the gravity of the Sun, Jupiter and planet’s satellites has been investigated. Saturn is considered to be an axisymmetric (A = B) solid body close to the dynamically spherical one. The orbits of Saturn and Jupiter are considered to be Keplerian ellipses in the inertial coordinate system. It has been shown that the entire set of small parameters of the problem can be reduced to two independent parameters. The averaged Hamiltonian function of the problem and the integrals of evolutionary equations are obtained disregarding the effect of satellites. Using the small parameter method, the expressions for the precession frequency and the nutation angle of the planet’s axis of rotation caused by the gravity of the Sun and Jupiter are obtained. Considering the planet with satellites as a whole preceding around the normal to the unmovable plane of Saturn’s orbit, the satellites effect on the Saturn rotation is taken into account via the corrections in the formula for the undisturbed precession frequency. The satellites are shown to have no effect on the nutation angle (in the framework of the accepted model), and the disturbances from Jupiter to make the main contribution to the nutation angle evolution. The effect of Jupiter on the nutation angle and the precession period is described with regard to the attraction of satellites.     

Parametric study of the diversion of a near Earth object on an Earth intersecting trajectory

To date, NASA's “Near Earth Object Program” has discovered over 5500 comets and asteroids on trajectories that bring them within “the neighborhood” of Earth's orbit. Nearly 1000 of these objects are classified as “potentially hazardous,” passing within 0.05 astronomical units of Earth's orbit. Discovery rates of such threatening bodies increase each year. Given this multitude of threats, in addition to evidence that the planet has absorbed many impacts over its history, it is reasonable to assume that another object will strike the Earth at some point in the future. Consequently, researchers have studied and proposed several mitigation techniques for such an occurrence. This study seeks to determine how effectively the attachment of a tether and ballast mass would divert the trajectory of such threatening objects. Specifically, the study analyzes the effects over time of such a system on objects of varying orbital semimajor axis and eccentricity, using various tether lengths and ballast masses. It was determined that the technique is most effective for NEOs with high eccentricity and small semimajor axis, and that system performance increases as tether length and ballast mass increase.     

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.     

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.     

地-月系平动点及Halo轨道的应用研究

地-月系统的平动点L1点及L2点的Halo轨道在探月工程中有重要的应用价值，可分别用于地月连续通信覆盖和月球背面的探测。由于在地-月系统中太阳的引力不可忽略，特别是在长时间作用以后，其动力学行为与摄动力较小的日-地系统有明显的不同。本文分析了如何利用太阳引力进入地-月系统的L1点及L2点的Halo轨道、以及由Halo轨道进入近月轨道的问题，两者综合起来构成了一条完整的地月低能转移轨道。研究结果对探月轨道设计有一定的参考价值。     

Application of periodic solutions of three-dimensional three-body problem to designing the orbit of a space telescope

A method of constructing three-dimensional orbits with a necessary evolution in the system the Sun — (Earth + Moon) is described. The orbit (promising from the viewpoint of solving formulated research problems) of the Millimetron spacecraft is suggested. Feasibility of such an orbit is demonstrated, as well as a possibility to observe with its help the majority of objects on the celestial sphere and to transmit the data to the Earth.     

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.     

Halo orbits in the vicinity of the L 2 libration point in the Sun-Earth system

Methods are proposed for constructing the orbits of spacecraft remaining for long periods of time in the vicinity of the L ₂ libration point in the Sun-Earth system (so-called halo orbits), and the trajectories of uncontrolled flights from low near-Earth orbits to halo orbits. Halo orbits and flight trajectories are constructed in two stages: A suitable solution to a circular restricted three-body problem is first constructed and then transformed into the solution for a restricted four-body problem in view of the real motions of the Sun, Earth, and Moon. For a halo orbit, its prototype in the first stage is a combination of a periodic Lyapunov solution in the vicinity of the L ₂ point and lying in the plane of large-body motion, with the solution for the linear second-order system describing small deviations of the spacecraft from this plane along the periodic solution. The desired orbit is found as the solution to the three-body problem best approximating the prototype in the mean square. The constructed orbit serves as a similar prototype in the second stage. In both stages, the approximating solution is constructed by continuation along a parameter that is the length of the approximation interval. Flight trajectories are constructed in a similar manner. The prototype orbit in the first stage is a combination of a solution lying in the plane of large-body motion and a solution for a linear second-order system describing small deviations of the spacecraft from this plane. The planar solution begins near the Earth and over time tends toward the Lyapunov solution existing in the vicinity of the L ₂ point. The initial conditions of both prototypes and the approximating solutions correspond to the spacecraft’s departure from a low near-Earth orbit at a given distance, perigee, and inclination.     

Magnetic trapping in orbit

A toroidal solenoid in orbit can represent a magnetic trap for bodies carrying dipole moment. We find the dependence of the trap dimensions on the electrical parameters and the orbit height. There is no trap if the solenoid radius exceeds a certain value. Various trap dimensions are calculated for a non-spinning space station in Earth or Sun orbit. The trap dimension is nearly independent of the solenoid orientation so that the solenoid may be fastened to the station rigidly. If the solenoid radius is in the range of some meters, the trap dimension is one order of magnitude greater for a distant Earth orbit and two orders of magnitude for an orbit around the Sun at distances of the order of one AU. The calculations are performed for a non-superconducting solenoid. The trap can be used for simplifying the docking on of maintenance vehicles or building modules, as safety mechanism for crewmembers working in outer space, and for the collection of paramagnetic or ferromagnetic dust stopped by mechanical means.