New Definition of Discovery for Solar System Objects

We propose a new formal definition of discovery for a Solar System object. It is based on an objective and mathematically rigorous algorithm to assess when a set of observations is enough to constitute a discovery. When this definition is satisfied, in almost all cases the orbit is defined well enough to establish the nature of the object discovered (Main Belt vs. Near Earth Asteroid, Trans-Neptunian vs. long period comet). The frequency of occurrence of exceptions is estimated by a set of numerical experiments. The availability of a non-subjective definition of discovery allows some rules to be adopted for the assignment of discovery credit with a minimum risk of dispute. Such rules should be fair, encourage good practice by the observers and acknowledge the contribution of the orbit computers providing the identifications and the orbits, as well as the one of all the contributing observers.     

Effective Hamiltonian for the D'Alembert Planetary Model Near a Spin/Orbit Resonance

The D'Alembert model for the spin/orbit problem in celestial mechanics is considered. Using a Hamiltonian formalism, it is shown that in a small neighborhood of a p:q spin/orbit resonance with (p,q) different from (1,1) and (2,1) the 'effective' D'Alembert Hamiltonian is a completely integrable system with phase space foliated by maximal invariant curves; instead, in a small neighborhood of a p:q spin/orbit resonance with (p,q) equal to (1,1) or (2,1) the 'effective' D'Alembert Hamiltonian has a phase portrait similar to that of the standard pendulum (elliptic and hyperbolic equilibria, separatrices, invariant curves of different homotopy). A fast averaging with respect to the 'mean anomaly' is also performed (by means of Nekhoroshev techniques) showing that, up to exponentially small terms, the resonant D'Alembert Hamiltonian is described by a two-degrees-of-freedom, properly degenerate Hamiltonian having the lowest order terms corresponding to the 'effective' Hamiltonian mentioned above.     

Numerical Exploration of Chermnykh’s Problem

We study the equilibrium points and the zero-velocity curves of Chermnykh’s problem when the angular velocity ω varies continuously and the value of the mass parameter is fixed. The planar symmetric simple-periodic orbits are determined numerically and they are presented for three values of the parameter ω. The stability of the periodic orbits of all the families is computed. Particularly, we explore the network of the families when the angular velocity has the critical value ω = 2√2 at which the triangular equilibria disappear by coalescing with the collinear equilibrium point L₁. The analytic determination of the initial conditions of the family which emanate from the Lagrangian libration point L₁ in this case, is given. Non-periodic orbits, as points on a surface of section, providing an outlook of the stability regions, chaotic and escape motions as well as multiple-periodic orbits, are also computed. Non-linear stability zones of the triangular Lagrangian points are computed numerically for the Earth–Moon and Sun–Jupiter mass distribution when the angular velocity varies.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

月球重力场对绕月低轨卫星的长期影响

从解析形式出发,利用月球重力场模型JGL165P1,分析了月球重力场(带谐项)对绕月低轨卫星的长期影响。为了减少计算误差,保证计算精度,在分析解中使用循环公式来计算倾角函数。结果指出对于一个高度为100km的极月轨道卫星,冻结轨道存在的可能性不大,但是当轨道倾角在i=90°附近或者高度再高一些,则有可能存在冻结轨道;对于100km高的初始圆轨道,卫星在无控的情况下半年内将会坠落到月球表面,如果高度增加到200km,则不进行轨道控制也不会坠落到月面上。利用仿真软件GEODYN解算出来的结果证实了上述结论。     

地月系L_2点中继星激光测距成功率的估算

为了模拟位于地月系L2点的中继星"鹊桥"与月球的位置关系,进而估算中继星激光测距的成功率,按照轨道周期约为14天的要求对中继星所在的晕轨道进行计算,建立了一个综合考虑望远镜抖动、大气抖动和预报轨道横向偏离的模型。从数值上给出了一条轨道周期为14. 78天,X方向(地月连线方向)振幅为12 493 km,Y方向为34 596 km,Z方向(垂直于地月轨道平面方向)为11 916 km的周期轨道。由于晕轨道的最小振幅远大于月球遮挡的临界振幅4 000 km,因此月球对中继星不存在遮挡问题。基于建立的测距成功率模型,根据昆明站(国际编号:7820)的激光测距系统对运行在该轨道上的中继星进行测距成功率分析,结果表明:测距成功率随着中继星横向轨道标准差的增大呈快速降低的趋势。对于中继星到测站的平均距离而言,当中继星没有横向偏离时,探测器产生的光电子数为0. 151,成功率为14. 07%;横向偏离2 km时,光电子数降为0. 035,成功率降为3. 46%。对比最近距离与最远距离的情况,无横向偏离的情况下,探测器产生的光电子数从0. 174降为0. 139,成功率从16. 01%降为13. 02%。该计算结果可为云南天文台1. 2 m望远镜实现中继星激光测距提供参考。     

Automated design of gravity-assist trajectories to Mars and the outer planets

In this paper, a new approach to planetary mission design is described which automates the search for gravity-assist trajectories. This method finds all conic solutions given a range of launch dates, a range of launch energies and a set of target planets. The new design tool is applied to the problems of finding multiple encounter trajectories to the outer planets and Venus gravity-assist trajectories to Mars. The last four-planet grand tour opportunity (until the year 2153) is identified. It requires an Earth launch in 1996 and encounters Jupiter, Uranus, Neptune, and Pluto. Venus gravity-assist trajectories to Mars for the 30 year period 1995–2024 are examined. It is shown that in many cases these trajectories require less launch energy to reach Mars than direct ballistic trajectories.Assistant Professor, School of Aeronautics and AstronauticsGraduate Student, School of Aeronautics and Astronautics     

初轨的稳健估计算法

讨论单测角资料(数据中包含异常值)初轨的稳健估计算法．该算法既可用于提供可靠的初轨,又可用作资料预处理,剔除异常值．稳健估计算法的崩溃点为16％至25％,可剔除6σ以上的异常值．     

推荐一组对长、短积分间隔皆优的线性多步积分公式

推荐了一组对短、中、长积分间隔以及带耗散力或较大偏心率等多种类型卫星轨道数值积分皆可获取较高计算精度的线性多步积分公式（MTM）,供读者选择使用。