Motion of a space station. I

In this first part of the work we develop the equations of motion of a triaxial space station in orbit around the oblate Earth. A first order solution of the problem is presented and the method of complete integration of the system is outlined up to second order of approximation. The zero order part of the Hamiltonian includes both the kinetic and potential energy (Earth's Newtonian attraction) of the station, while the motion in the vicinity of a specific configuration is assumed.The solution leads to deviations on the attitude introduced by the oblateness of the Earth. Such attitude is an exact solution of the equations of the station when its center of mass moves in an elliptic Keplerian orbit.The explicit expressions of the complete solution, discussion of other possible effects on the motion and numerical comparisons will be presented in the second part of the work.     

Three-dimensional subsatellite motion under air drag and oblateness perturbations

The three-dimensional relative motion of a subsatellite with respect to a reference station in an elliptical orbit is studied. A general theory based on the variation of the relative elements, i.e. the instantaneous differences between the orbital parameters of the subsatellite and those of the station, is formulated in order to incorporate arbitrary perturbing forces acting on both satellites. The loss of precision inherent in the subtraction of almost identical quantities is avoided by the consistent use of difference variables. In the absence of perturbations exact analytical representations can be obtained for the relative state parameters. The influences of air drag and Earth's oblateness on the relative motion trajectories are investigated and illustrated graphically for a number of cases.     

An exact solution to determination of an open orbit

We present an exact solution of the equations for orbit determination of a two body system in a hyperbolic or parabolic motion. In solving this problem, we extend the method employed by Asada, Akasaka and Kasai (AAK) for a binary system in an elliptic orbit. The solutions applicable to each of elliptic, hyperbolic and parabolic orbits are obtained by the new approach, and they are all expressed in an explicit form, remarkably, only in terms of elementary functions. We show also that the solutions for an open orbit are recovered by making a suitable transformation of the AAK solution for an elliptic case.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.     

估算有中子星作伴星的黑洞质量

在本文中，我们讨论了由中子星和黑洞构成的密近双星系统的引力效应。首先，我们清楚地解释了引力辐射存在引起的轨道圆化，并估算了圆化的时间尺度。求解了围绕史瓦西黑洞作圆轨道运动的矢量进动方程，得到了一个精确解。由进动周期和中心黑洞的质量的关系，我们得到了估算中子星伴随的黑洞质量的方法。随后我们讨论了以中子星脉冲信号随自旋进动变化来寻找黑洞的方法。     

East-west stationkeeping requirements of nearly synchronous satellites due to Earth's triaxiality and luni-solar effects

A closed form solution, for longitude and semimajor axis deviations in the neighborhood of a prespecified station, is obtained for nearly synchronous satellites. The model use includes the important terms in Earth's zonal and tesseral harmonics as well as the luni-solar perturbations. The initial semimajor axis for two-maneuver east-west stationkeeping is then deduced. Due to the luni-solar effects, it is found that the initial semimajor axis deviation from synchronous orbit value is highly dependent on the initial position of the satellite relative to the Moon and the Sun. Verifications of the results by means of numerical integrations are also included.     

Orbit determination of space debris: admissible regions

The main problem in the orbit determination of the space debris population orbiting our planet is identifying which separate sets of data belong to the same physical object. The observations of a given object during a passage above an observing station are collectively called a Too Short Arc (TSA): data from a TSA cannot allow for a complete determination of an orbit. Therefore, we have to solve first the identification problem, finding two or more TSAs belonging to the same physical object and an orbit fitting all the observations. This problem is well known for the determination of orbits of asteroids: we shall show how to apply the methods developed for preliminary orbit determination of heliocentric objects to geocentric objects. We shall focus on the definition of an admissible region for space debris, both in the case of optical observations and radar observations; then we shall outline a strategy to perform a full orbit determination.     

Instabilities in the Sun–Jupiter–Asteroid three body problem

We consider dynamics of a Sun–Jupiter–Asteroid system, and, under some simplifying assumptions, show the existence of instabilities in the motions of an asteroid. In particular, we show that an asteroid whose initial orbit is far from the orbit of Mars can be gradually perturbed into one that crosses Mars’ orbit. Properly formulated, the motion of the asteroid can be described as a Hamiltonian system with two degrees of freedom, with the dynamics restricted to a “large” open region of the phase space reduced to an exact area preserving map. Instabilities arise in regions where the map has no invariant curves. The method of MacKay and Percival is used to explicitly rule out the existence of these curves, and results of Mather abstractly guarantee the existence of diffusing orbits. We emphasize that finding such diffusing orbits numerically is quite difficult, and is outside the scope of this paper.     

A Research on the Orbit Determination Method by Means of Sparse Data of Electronic Fences

Based on the Lambert equation and knowledge of space geometry a method of orbit determination is given using the sparse observational data provided by the space monitoring electronic fence device. Our simulated experiment of a large number of targets shows that the initial orbit determined by this method can be improved and can converge to a final accuracy better than 100 m, so proving that the method can be applied to the orbit determination of an overwhelming majority of space targets with the observed data of the electronic fence. Finally, the effect of the latitude of the observing station on the application of the method is discussed.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

Dispersion of ensembles of non-interacting particles

The dynamics of an ensemble of particles emanating from a common point with a distribution of velocities is modeled as a continuum of particles described by a phase space distribution function. A general solution for the distribution function and the associated spatial density function is obtained for a general dynamical system. The special cases of linear dynamical systems and slow dispersion from a circular orbit are treated in detail. A transcendental equation is derived, the roots of which determine the time since initial dispersion from knowledge of the spatial density function at later times.     

The eccentric case of a fast-rotating,gravity-gradient-perturbed satellite attitude solution

A closed-form first-order perturbation solution for the attitude evolution of a triaxial space object in an elliptical orbit is presented. The solution, derived using the Lie–Deprit method, takes into account gravity-gradient torque and is facilitated by an assumption of fast rotation of the object. The formulation builds on the earlier implementation of Lara and Ferrer, which assumes a circular orbit. The previously presented work—which assumes spin about an object’s axis of maximum inertia—is further extended by the explicit presentation of the transformations required to apply the solution to an object spinning about its axis of minimum inertia. Additionally, several numerical analyses are presented to more completely assess the utility of the solution. These studies (1) validate the elliptical solution, (2) assess the impact of varying the small parameter of the perturbation procedure, (3) analyze the assumption of fast rotation, and (4) apply the solution to the common and important scenario of a tumbling rocket body.     

Constructive-analytical solution of the problem of the secular evolution of polar satellite orbits

The well-known twice-averaged Hill problem is considered by taking into account the oblateness of the central body. This problem has several integrable cases that have been studied qualitatively by many scientists, beginning with M.L. Lidov and Y. Kozai. However, no rigorous analytical solution can be obtained in these cases due to the complexity of the integrals. This paper is devoted to studying the case where the equatorial plane of the central body coincides with the plane of its orbital motion relative to the perturbing body, while the satellite itself moves in a polar orbit. A more detailed qualitative study is performed, and an approximate constructive-analytical solution of the evolution system in the form of explicit time dependences of the eccentricity and pericenter argument of the satellite orbit is proposed. The methodical accuracy for the polar orbits of lunar satellites has been estimated by comparison with the numerical solution of the system.     

高轨空间碎片的观测模拟

对高轨空间碎片进行了观测模拟。在观测模拟中，分别选取南京和云南两个观测站；采用2种跟踪方式：跟踪恒星和跟踪卫星；对于同一地理经度上的同步卫星进行了观测；提供了发现新目标的方法；并给出了新目标的初轨。这为我国空间碎片的数据库的建立作了有益的尝试，对于卫星和航天器的飞行安全具有特别重要的意义。     

基于转发式测轨系统的GEO卫星机动检测初探

通过对利用C波段转发式测轨网观测"鑫诺一号"卫星得到的原始资料进行初步分析,探讨卫星轨道机动时卫星在测站与卫星连线方向的距离、速度和加速度的变化规律,并提出一种通过对原始观测数据的拟合比对,实现GEO卫星(地球同步卫星)轨道机动检测的方法。分析比较表明,用该方法检测到的卫星机动开始和结束时刻与星载推力器实际喷火开始和结束时刻之间的差别小于3 min,该方法对处于机动期间的轨道确定有一定的积极意义。     

Application of Hamiltonian structure-preserving control to formation flying on a J 2-perturbed mean circular orbit

The bounded quasi-periodic relative trajectories are investigated in this paper for on-orbit surveillance, inspection or repair, which requires rapid changes in formation configuration for full three-dimensional imaging and unpredictable evolutions of relative trajectories for non-allied spacecraft. A linearized differential equation for modeling J ₂ perturbed relative dynamics is derived without any simplified treatment of full short-period effects. The equation serves as a nominal reference model for stationkeeping controller to generate the quasi-periodic trajectories near the equilibrium, i.e., the location of the chief. The developed model exhibits good numerical accuracy and is applicable to an elliptic orbit with small eccentricity inheriting from the osculating conversion of orbital elements. A Hamiltonian structure-preserving controller is derived for the three-dimensional time-periodic system that models the J ₂-perturbed relative dynamics on a mean circular orbit. The equilibrium of the system has time-varying topological types and no fixed-dimensional unstable/stable/center manifolds, which are quite different from the two-dimensional time-independent system with a permanent pair of hyperbolic eigenvalues and fixed-dimensions of unstable/stable/ center manifolds. The unstable and stable manifolds are employed to change the hyperbolic equilibrium to elliptic one with the poles assigned on the imaginary axis. The detailed investigations are conducted on the critical controller gain for Floquet stability and the optimal gain for the fuel cost, respectively. Any initial relative position and velocity leads to a bounded trajectory around the controlled elliptic equilibrium. The numerical simulation indicates that the controller effectively stabilizes motions relative to the perturbed elliptic orbit with small eccentricity and unperturbed elliptic orbit with arbitrary eccentricity. The developed controller stabilizes the quasi-periodic relative trajectories involved in six foundational motions with different frequencies generated by the eigenvectors of the Floquet multipliers, rather than to track a reference relative configuration. Only the relative positions are employed for the feedback without the information from the direct measurement or the filter estimation of relative velocity. So the current controller has potential applications in formation flying for its less computation overload for on-board computer, less constraint on the measurements, and easily-achievable quasi-periodic relative trajectories.     

The computation of relative motion with increased precision

Encke's method as modified by Potter to increase the accuracy of orbit computations of gravitationally interacting bodies is applied to the problem of relative motion of non-interacting space vehicles. This technique is then combined with a simple transformation of the independent variable to arrive at a system of equations from which the relative motion may be determined with increased precision.     

A fast tour design method using non-tangent v-infinity leveraging transfer

The announced missions to the Saturn and Jupiter systems renewed the space community interest in simple design methods for gravity assist tours at planetary moons. A key element in such trajectories are the V-Infinity Leveraging Transfers (VILT) which link simple impulsive maneuvers with two consecutive gravity assists at the same moon. VILTs typically include a tangent impulsive maneuver close to an apse location, yielding to a desired change in the excess velocity relative to the moon. In this paper we study the VILT solution space and derive a linear approximation which greatly simplifies the computation of the transfers, and is amenable to broad global searches. Using this approximation, Tisserand graphs, and heuristic optimization procedure we introduce a fast design method for multiple-VILT tours. We use this method to design a trajectory from a highly eccentric orbit around Saturn to a 200-km science orbit at Enceladus. The trajectory is then recomputed removing the linear approximation, showing a Δv change of <4%. The trajectory is 2.7 years long and comprises 52 gravity assists at Titan, Rhea, Dione, Tethys, and Enceladus, and several deterministic maneuvers. Total Δv is only 445 m/s, including the Enceladus orbit insertion, almost 10 times better then the 3.9 km/s of the Enceladus orbit insertion from the Titan–Enceladus Hohmann transfer. The new method and demonstrated results enable a new class of missions that tour and ultimately orbit small mass moons. Such missions were previously considered infeasible due to flight time and Δv constraints.