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.     

Evolution of Motion of a Binary Planet

A model of a binary planet, consisting of a material point of small mass and a deformable viscoelastic sphere, is suggested. The center of mass of the binary planet moves in the gravitational field of a central body in the plane, which contains planets forming the binary planet. A deformable spherical planet rotates around the axis orthogonal to the plane of planetary motion. Planet deformations are described by the linear theory of viscoelasticity. It is shown that with an appropriate approximation of the gravitational potential, there is a class of quasicircular orbits, when the eccentricities of an orbit of the center of mass of a binary planet and an orbit, describing mutual planet motion, are equal to zero. The further evolution of motion is investigated in this class of orbits with the use of the canonical Poincare–Andoyer variables. Corresponding averaged equations are found, and phase pictures are constructed; the stability of stationary solutions is investigated on the basis of equations in variations. For the Solar system planets with their satellites, forming binary planets, the application of the presented model allows us to conclude that satellites sooner or later will fall on the corresponding planets.     

Spatial approaches to moons from resonance relative to invariant manifolds

In this study, the final approach to a moon or other body from resonance is explored and compared to the invariant manifolds of unstable periodic orbits. It is shown that the stable manifolds of planar Lyapunov orbits can act as a guide for the periods or resonances that are required for the final approach in both the planar and spatial problems. Previously developed techniques for the planar problem are expanded for use with resonances and used for comparison with trajectories approaching a moon from these resonances. A similar technique is then used for exploring the relationship of invariant manifolds to approach trajectories in the spatial problem. It is shown that the invariant manifolds of unstable periodic orbits provide insight into the trajectory design, and they can be used as a guide to the more direct approach trajectories.     

Calculation of quasisynchronous orbits of a spacecraft around Phobos for solving the problem of landing on its surface

The problem of selecting quasi-synchronous orbits of a spacecraft around Phobos is considered. These quasi-synchronous orbits are far (with respect to the Hill’s sphere) quasi-satellite orbits with retrograde rotation in the restricted three body problem. The orbit should pass through a given point at a specified time instant. It should also possess a property of minimum distance from the Phobos surface at every passage above the region of planned landing. The equations of dynamics are represented in the form describing the orbit as a combination of motions in two drifting ellipses, inner and outer ellipses. The center of the outer ellipse is located on the inner ellipse. A formula is derived that relates averaged values of half-axes of the inner and outer ellipses. It is used for construction of the first approximation of numerically designed orbit, which makes it possible to simplify and speed up the computing process. The tables of initial conditions obtained as a result of calculations are presented.     

解析求解非线性相对运动的最优重构问题

椭圆轨道相对运动模型的线性化导致其在大尺度相对运动应用中精度不能满足需求。针对任意椭圆轨道上的大尺度航天器编队最优重构问题,提出一种基于椭圆轨道非线性相对运动模型的近似解析求解方法。首先通过变分法建立了非线性最优重构问题的数学模型;然后采用摄动法,以偏近点角为积分变量求得了不含特殊积分的解析开环最优控制,有效地避免了真近点角域下最优控制解所含有的特殊积分。仿真验证了所求最优控制的有效性和优越性,结果表明在相对运动尺度较大时,相比基于椭圆轨道线性化模型的最优控制,在燃耗保持相近的情况下,所求非线性控制有效地降低了重构误差。     

Stable spatial orbits around collinear libration points

A technique of generation of spatial periodic solutions to the restricted circular three-body problem from periodic orbits of the planar problem has been used for the families of orbits around collinear libration points L ₁ and L ₂. Developing the families obtained at the 1: 1 resonance, we have obtained stable solutions both in the Earth-Moon system and in the Sun-Earth system. Of course, the term “around the libration point” is rather conventional; the obtained orbits become more similar to the orbits around the smaller attracting body. The further development of the family of orbits “around” the libration point L ₂ in the Sun-Earth system made it possible to find the orbits satisfying the new, much more rigorous constraints on cooling the spacecraft of the Millimetron project.     

Collisions of stars by oscillating orbits of the second kind

Ordinary estimations of the number of star collisions in our galaxy—by simple kinematic considerations—lead to a very small number of such collisions: about one or even less every millions of years. However star collisions can occur through the following indirect way which has a much higher probability. (a) Binary stars are very common in our galaxy, about 30–50% of the stars. (b) If two binary stars meet a triple system can be formed by an ordinary exchange type motion. (c) A triple system is generally decomposed into the “inner orbit” (i.e. the relative orbit of the two nearest stars) and the “outer orbit” (i.e. the relative orbit of the third star with respect to the center of mass of the two nearest stars). The major axes of these two orbits have generally small perturbations and it is the same for the eccentricity of the outer orbit. On the contrary, if the relative inclination of the two orbits is large, the perturbations of the eccentricity of the inner orbit are important and can even in some cases lead to an eccentricity equal to one, that is to a collision of the two stars of the inner orbit.Such orbits can be called “oscillating orbits of the second kind”, indeed the first oscillating orbits—conceived by Khilmi and described for the first time in an example by Sitnikov—have unbounded mutual distances r_ij, but the system always come back to small sizes, it has an infinite number of very large expansions followed by strong contractions and, in the three-body case, an upper bound of lim inf (r_1.2 + r_1.3 + r_2.3) can be given in terms of the three masses and the integrals of motion. For the oscillating orbits of the second kind the mutual distances r_ij are bounded, but the velocities are unbounded (i.e. lim inf r_ij = 0 for at least one r_ij) and the system goes to a collision if the bodies have non-zero radius even small. The analytical study of the oscillating orbits of the second kind is a part of the general analytical study of the three-body problem, a part which must be valid for large eccentricities and large inclinations. The use of Delaunay's variables and of a Von Zeipel transformation lead to a first order integrable approximation, valid for any eccentricities and any inclinations, and giving the following results: (a) The oscillating orbits of the second kind occur when the angular momentum of the outer orbit has a modulus sufficiently close to the modulus of the total angular momentum of the three-body system. Hence these orbits occur for inclinations in the vicinity of 90°. (b) The oscillating orbits represent a set of positive measure of phase space and the first order study allows to give a rough estimation of the probability of collisions—even for stars of infinitely small radius. This probability, for given initial major axes and eccentricities and for isotropic arbitrary initial orientations, is generally of the order of m₃$\text{R}\text{M}$ (m₃ being the mass of the outer star, M the total mass and R the ratio of the period of the inner orbit to that of outer orbit).One question remains to be solved: how many collisions of stars are due to that phenomenon? That question is difficult because the probability of formation of a triple system by a random meeting of two binaries is very uneasy to estimate. However it seems that, compared to the usual evaluations based on pure kinematic considerations without gravitational effects, the number of collisions must be multiplied by a factor between one thousand and one million.     

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.     

Collision Probabilities for Keplerian Orbits

We present a new derivation of the probability of collisions between spherical satellites occupying Keplerian orbits. The equations follow from the central concept of the instantaneous collision rate, an expression that describes the occurrence of collisions by using a Dirac -function. The derivation proceeds by showing how this instantaneous collision rate can be averaged over orbital mean anomaly angles and, additionally, over orbital precession angles to generate expressions appropriate for intermediate and long time scales. Collision rates averaged over mean anomalies tend to be non-zero during relatively brief collision seasons, when the peak collision probability may exceed the long-term average by several orders of magnitude. Derived precession-angle averages have a functional form similar but not identical to the collision probability expression derived using the spatial density approach of Kessler (Icarus, 48: 39–48, 1981), and the two methods have been found to yield numerical results to within 1% for all cases examined.     

Dynamics of a plasma ring rotating in the magnetic field of a central body: Magneto-gravitational waves

The model problem of the dynamics of a planar plasma ring rotating in the dipole magnetic field of a central body is considered. A finite-dimensional mathematical model of the system is synthesized by the Boubnov-Galerkin method. The class of solutions corresponding to magneto-gravitational waves associated with deformations of the ring boundaries is investigated.     

从地球到日-火系平动点轨道的转移

平动点轨道特殊的空间位置及动力学特征,使其在深空探测中具有重要的应用。以日-火系平动点轨道(Lissajous与Halo轨道)任务为目标,结合平动点轨道的不变流形理论,研究了小推力转移问题。首先给出了圆型限制性三体动力学模型下平动点附近不变流形(稳定和不稳定流形)高阶分析解以及相应的计算实例。接着以流形分析解为基础,建立了初始小推力轨道优化模型,并利用改进的协作进化算法求解初始小推力轨道。最后将初始轨道离散,采用多点打靶法将最优控制问题转化为参数优化问题,并用序列二次规划方法(SQP)求解。仿真结果证明轨道设计方法的有效性。     

Biobjective planning of an active debris removal mission

The growth of the orbital debris population has been a concern to the international space community for several years. Recent studies have shown that the debris environment in Low Earth Orbit (LEO, defined as the region up to 2000 km altitude) has reached a point where the debris population will continue to increase even if all future launches are suspended. As the orbits of these objects often overlap the trajectories of satellites, debris create a potential collision risk. However, several studies show that about 5 objects per year should be removed in order to keep the future LEO environment stable. In this article, we propose a biobjective time dependent traveling salesman problem (BiTDTSP) model for the problem of optimally removing debris and use a branch and bound approach to deal with it.     

Mapping orbits with low station keeping costs for constellations of satellites based on the integral over the time of the perturbing forces

The present paper is concerned with the search for orbits that have potential to require low fuel consumption for station-keeping maneuvers for constellations of satellites. The method used to study this problem is based on the integral over the time of the undesired perturbing forces. This integral measures the change of velocity caused by the perturbation forces acting on the satellite, so mapping orbits that are less perturbed, which generates good candidates for orbits that requires low fuel consumption for station-keeping maneuvers. The integral over the time depends only on the orbit of the spacecraft and the dynamical system considered. The type of engine and the control technique applied to the spacecraft are not considered to search for those orbits. It can be a good strategy to be applied for a first mapping of orbits. For this search, it is analyzed the integral of orbits with different values of the Keplerian elements in order to find the best ones with respect to this criterion. The perturbations considered are the ones caused by the third body, which includes the Sun and the Moon, and the J₂ term of the geopotential. The results presented here show numerical simulations to obtain the integral of those perturbing forces for different orbits. The GPS and the Molniya constellations are used as examples for those calculations.     

Qualitative Properties of Orbits for Minimum-Fuel Impulsive Transfers between Circular Coplanar Orbits with a Given Launch Time

The coplanar problem of minimizing propellant consumption in impulsive transfer between circular boundary orbits is investigated. The launch time and the initial configuration of objects on the boundary orbits are specified arbitrarily. The qualitative properties of optimal two-impulse trajectories and their optimality in the class of multi-impulse transfers are studied.     

Modifying the atlas of low lunar orbits using inert tethers

For long enough tethers, the coupling of the attitude and orbital dynamics may show non-negligible effects in the orbital motion of a tethered satellite about a central body. In the case of fast rotating tethers the attitude remains constant, on average, up to second order effects. Besides, for a tether rotating in a plane parallel to the equatorial plane of the central body, the attitude–orbit coupling effect is formally equal to the perturbation of the Keplerian motion produced by the oblateness of the central body and, therefore, may have a stabilizing effect in the orbital dynamics. In the case of a tethered satellite in a low lunar orbit, it is demonstrated that feasible tether lengths can help in modifying the actual map of lunar frozen orbits.     

地月高能共振循环轨道的快速计算及最优选择

针对现有高能共振循环轨道计算方法存在计算量大、有可能改变轨道共振特性和不能构造共振比大于2.3的地月循环轨道等缺点,本文提出了一种地月圆型限制性三体问题下高能共振循环轨道的快速计算方法。首先根据轨道在月球附近的组成弧段对高能共振循环轨道进行分类;然后根据轨道类型构建二体开普勒椭圆轨道;再进一步计算圆型限制性三体问题下的地月高能共振循环轨道;最后根据能量、稳定性、时间周期、近地点高度和近月点高度对所计算出的地月高能共振循环轨道进行最优选择。仿真结果表明,本文所提出的方法简单有效,能够计算出共振比为5:2的地月高能共振循环轨道。     

On dynamics of a plasma ring rotating in the magnetic field of a central body: Magneto-gyroscopic waves. Problems of stability and quantization

Based on a mathematical model described in [1], some new aspects of the dynamics of a thin planar plasma ring rotating in the magnetic field of a central body are considered. The dipole field is considered assuming that the dipole has a small eccentricity, and the dipole axis is inclined at a small angle to the central body’s axis of rotation. Emphasis is placed on the problem of stability of the ring’s stationary rotation. Unlike [1], the disturbed motion is considered which has a character of eddy magneto-gyroscopic waves. The original mathematical model is reduced to a system of finite-difference equations whose asymptotic analytical solution is obtained. It is demonstrated that some “elite” rings characterized by integral quantum numbers are long-living, while “lethal” or unstable rings (antirings) are associated with half-integer quantum numbers. As a result, an evolutionally rife rotating ring of magnetized plasma turns out to be stratified into a large number of narrow elite rings separated by gaps whose positions correspond to antirings. The regions of possible existence of elite rings in near-central body space are considered. Quantum numbers determining elite eigenvalues of the mean sector velocity (normalized in a certain manner) of a ring coincide with the quantum numbers appearing in the solution to the Schrödinger equation for a hydrogen atom. Perturbations of elite orbits corresponding to these quantum numbers satisfy the de Brogli quantum-mechanical condition. This is one more illustration of the isomorphism of quantization in microcosm and macrocosm.     

The Investigation of Possible Approaches of Cataloged Space Objects to Manned Spacecraft

Initial orbital parameter errors are used to examine the miss distance between a spacecraft and an ensemble of tracked objects by a Monte Carlo-type analysis. The radial separation between orbits is evaluated and a keep-out zone is determined, which reduces the risk of collision to an acceptable level.An operational prediction methodology is suggested based on a catalog database, which identifies potentially hazardous approaches and computes the probability of collision for selected spacecraft. An example for the Mir Space Station is presented, which estimates the collision probability and the cross-sectional flux of cataloged objects for the time frame of interest. The results appear to be in good agreement with those of other space debris models.     

Particle swarm optimization applied to impulsive orbital transfers

The particle swarm optimization (PSO) technique is a population-based stochastic method developed in recent years and successfully applied in several fields of research. It mimics the unpredictable motion of bird flocks while searching for food, with the intent of determining the optimal values of the unknown parameters of the problem under consideration. At the end of the process, the best particle (i.e. the best solution with reference to the objective function) is expected to contain the globally optimal values of the unknown parameters. The central idea underlying the method is contained in the formula for velocity updating. This formula includes three terms with stochastic weights. This research applies the particle swarm optimization algorithm to the problem of optimizing impulsive orbital transfers. More specifically, the following problems are considered and solved with the PSO algorithm: (i) determination of the globally optimal two- and three-impulse transfer trajectories between two coplanar circular orbits; (ii) determination of the optimal transfer between two coplanar, elliptic orbits with arbitrary orientation; (iii) determination of the optimal two-impulse transfer between two circular, non-coplanar orbits; (iv) determination of the globally optimal two-impulse transfer between two non-coplanar elliptic orbits. Despite its intuitiveness and simplicity, the particle swarm optimization method proves to be capable of effectively solving the orbital transfer problems of interest with great numerical accuracy.