一个膺三阶辛积分器

在太阳系动力学中,辛积分器已成为研究哈密顿系统的长期定性演化的最佳工具.对于可积分离的哈密顿系统H=H0+∑i=1N∈iHi(∈≤1),构造了一个膺三阶辛积分器.它大约相当于Wisdom-Holman二阶辛积分器的一次校正或Forest-Ruth四阶辛算法的精度.此外,含力梯度的辛算法也适合处理哈密顿系统H=Ho(q,P)+∈H1(q),其精度好于原辛积分器,但不优越于相应膺高阶辛积分器.     

On the Structure of Symplectic Mappings. The Fast Lyapunov Indicator: a Very Sensitive Tool

It is already known (Froeschlé, Lega and Gonczi, 1997) that the Fast Lyapunov Indicator (FLI), that is the computation on a relatively short time of the largest Lyapunov indicator, allows to discriminate between ordered and weak chaotic motion. We have found that, under certain conditions, the FLI also discriminates between resonant and non-resonant orbits, not only for two-dimensional symplectic mappings but also for higher dimensional ones. Using this indicator, we present an example of the Arnold web detection for four and six-dimensional symplectic maps. We show that this method allows to detect the global transition of the system from an exponentially stable Nekhoroshevs like regime to the diffusive Chirikovs one.     

辛方法的校正公式

1996年Wisdom等提出了对辛方法进行校正的概念和实践，现在继续对辛校正进行详尽讨论和数值比较，尤其对哈密顿函数可分解为一个主要部分和多个次要部分的一般情形，用Lie级数推导任意阶的各种辛算法的一次和二次辛校正公式并对一些算法给出具体的辛校正公式。又以日、木、土三体问题为模型进行数值实验，结果表明一次辛校正能提高精度，改善数值稳定性。计算效率也比较高，因而值得推荐使用，辛方法通常用大步长数值积分，这时二次辛校正并没有显著提高结果的精度，却大大增加了计算时间，不应予以推荐。     

Three-body dynamics: intermittent chaos with strange attractor

We have studied the structure of chaos in three-body dynamics using the concept of intermittency, implying that violent states of a system alternate in time with quasi-regular states producing together a non-stationary and evolving pattern of unpredictable behaviour. Computer simulations are produced to demonstrate explicitly sporadic short violent bursts in quasi-regular hierarchical states of the systems. This is seen both in orbits and in the long time series generated by the system. The time series prove to be similar in shape to what is observed in various physical experiments with laboratory chaotic systems when they reveal the so-called type-III intermittency. The new effective methods of time series analysis enable us to discover a strange attractor with a fractal dimension slightly above 2. This shows that three-body dynamics has the same intrinsic qualitative structure and quantitative measure of chaos as the widely known chaotic system, the Lorenz attractor.     

On the Calculation of Lyapunov Indicators with Post-stabilization in a Weyl Field

We present details of a work aiming at the overestimation of Lyapunov exponents defined by the geodesic deviation equations in the previous work. The geodesic deviation vector with post-stabilization is used to compute the fast Lyapunov indicator, considered to be a very sensitive tool for discrimination between ordered or weakly chaotic motions. We make a detailed study of the dynamics in the superposed Weyl field between a black hole and shell of octopoles by using the fast Lyapunov indicator with the Poincare surface of section. In particular, we examine the effects on the dynamics around the fixed points, of varying one of the three parameters (specific energy E, specific angular momentum L and octopolar moment O), while keeping the other two fixed, and identify the intervals of the varying parameter where the motion is regular or chaotic.     

倾角函数的四精度计算试验

倾角函数是天体力学分析理论中一种常用的函数.当把摄动方程展开成时间和根数的形式时需要用到.历史上提出了很多经典的倾角函数递推算法,并在双精度平台下开发了Fortran程序.进行了1次四精度计算倾角函数的试验,结果表明:L平面递推方法的四精度计算精度可达10-22,计算速度比双精度Jacobi方法快6倍.     

An efficient computational algorithm for the initial value problem associated with the universal Y functions

In the present paper, an efficient iterative method of arbitrary integer order of convergent ≥2 based on the homotopy continuation techniques for the solution of the initial value problem of space dynamics using the universal Y functions is presented. The method is of dynamic nature in the sense that, ongoing from one iterative scheme to the subsequent one, only additional instruction is needed. Most importantly, the method does not need any prior knowledge of the initial guess. This is a property which avoids the critical situations between divergent to very slow convergent solutions that may exist in other numerical methods which depend on initial guesses. A computational package for digital implementation of the method is given, together with numerical applications for elliptic, hyperbolic, and parabolic orbits. The accuracy of the results for all orbits is O(10^–16). (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mathematical Model of Unit Vector Method (MMUVM) and Convergence of Its Iteration Method in Simplified Form (PUVM1)

The unit vector method (PUVM1) by which the observed data of artificial satellites are used to determine the perturbed initial orbits has been widely applied. In order to further perfect and improve this method, first, a mathematical model (MMUVM) corresponding to this method is constructed on the basis of measuring errors. In essence, the MMUVM is a nonlinear optimization problem. In the light of the MMUVM, simulated and observed data of multiple circuits are taken to form specific objective functions, which are then processed by means of a direct search method and with of a tri-diagonal quadratic interpolation model for solving optimization problems. The calculated results show that the optimization model MMUVM is right and reasonable and the adopted direct method is practical and effective. Secondly, the relation between PUVM1 and MMUVM is further clarified: PUVM1 is essentially a simplified form of the MMUVM, and the primary reason is found from mathematical principles, why it is that PUVM1 can only be applied to short-arc data within one circuit and not to long-arc or multi-circuit data. Finally, a preliminary theoretical analysis of the convergence of the iteration algorithm of PUVM1 is carried out and instancesof numerical verification given. It is pointed out that the iterative scheme of PUVM1 is conditionally convergent. This means that sometimes, even though the quasi-normalization equation is reasonable, the iteration diverges.     

Long-term stability of planetary orbits between Jupiter and Saturn

We extend our two previous studies on the existence of stable orbits in the Solar System by examining the domain between Jupiter and Saturn. We place (1) a massless object, (2) a Moon-mass object, (3) a Mars-mass object, (4) an Earth-mass object, and (5) a Uranus-mass object in the said region. Note that these objects are considered separately in the framework of our simulations. Our goal is to explore the orbital stability of those objects. We employ the Lie-integration method, which is fast and well established, allowing us to solve the respective differential equations for the $N$-body system. Hence, we consider the celestial bodies spanning from Jupiter to Neptune, including the aforementioned test object, the main focus for our model simulations. The integrations indicate that in some models the test objects placed in the region between Jupiter and Saturn reside in that region for more than 600 Myr. Between 5 and 10 au, mean-motion resonances (MMRs) take place acting upon the test objects akin to simulations of Paper I and II. Our models indicate relatively small differences for the long-term stability of the five test objects notwithstanding their vastly different masses. Generally, it is found that between $a_{\text{ini}}=7.04$ and 7.13 au the orbits become unstable mostly within 5 million years and further out, that is, up to $a_{\text{ini}}=7.29$ au, the duration of stability lengthens to up to hundreds of millions of years.     

倾角函数递推稳定性分析

历史上曾经提出了较多的倾角函数递推算法,但有一些已经被证明在高阶是不稳定的.通过对递推方向上倾角函数的数量级分析,可以判断倾角函数递推的稳定性.对于常用的3项递推,只有Mk(l)递推是稳定的,其他递推均是不稳定的.但是对于多项递推比较复杂,还需深入分析.     

Formation of the Kuiper Belt by Long Time-Scale Migration of Jovian Planets

The orbital migration of Jovian planets is believed to have played an important role in shaping the Kuiper Belt. We investigate the effects of the long time-scale (2×107 yr) migration of Jovian planets on the orbital evolution of massless test particles that are initially located beyond 28 AU. Because of the slowness of the migration, Neptune's mean motion resonances capture test particles very efficiently. Taking into account the stochastic behavior during the planetary migration and for proper parameter values, the resulting concentration of objects in the 3:2 resonance is prominent, while very few objects enter the 2:1 resonance, thus matching the observed Kuiper Belt objects very well. We also find that such a long time-scale migration is favorable for exciting the inclinations of the test particles, because it makes the secular resonance possible to operate during the migration. Our analyses show that the v8 secular resonance excites the eccentricities of some test particles, so decreasing their perihelion distances, leading to close encounters with Neptune, which can then pump the inclinations up to 20℃.     

Research on the Method of Orbit Determination Based on the Self-adaptive Stable Least p-norm Estimate

The traditional least square estimation (LSE) method for orbit determination will not be optimal if the error of observational data does not obey the Gaussian distribution. In order to solve this problem, the least p-norm (Lp) estimation method is presented in this paper to deal with the non-Gaussian distribution cases. We show that a suitable selection of parameter p may guarantee a reasonable orbit determination result. The character of Lp estimation is analyzed. It is shown that the traditional Lp estimation method is not a robust method. And a stable Lp estimating based on data depth weighting is put forward to deal with the model error and outlier. In the orbit determination process, the outlier of observational data and coarse model error can be quantitatively described by their weights. The farther is the data from the data center, the smaller is the value of data depth and the smaller is the weighted value accordingly. The result of the new Lp method is stabler than that of the traditional Lp estimation and the breakdown point could be up to 1/2. In addition, the orbit parameter is adaptively estimated by residual analysis and matrix estimation method, and the estimation efficiency is enhanced. Finally, by taking the Space-based Space Surveillance System as an example and performing simulation experiments, we show that if there are system error or abnormal value in the observational data or system error in satellite dynamical model and space-based observation platform, LSE will not be optimal even though the observational data obeys the Gaussian distribution, and the orbit determination precision by the self-adaptive robust Lp estimation method is much better than that by the traditional LSE method.     

A Study on Orbit Determination by Using the Data of Crossovers from Satellite-borne Altimeter

Starting with the status of the developments of oceanic altimetry satellites, the significance of orbit determination by using altimeter data is introduced. Then the error correction model of the altimeter data and the calculational method of the data of crossovers are analyzed. The modification quantities of the errors and the adopted model concerned in the files of the altimeter data from the aircraft JASON-1 are also introduced in detail. Finally, through the calculations of the simulated data and the data of actual measurements, the highest accuracies of the orbit determination reached by solely using the altimeter data and the data of crossovers are analyzed, respectively. And this work provides a valuable reference to the practical applications in future.     

倾角函数及其导数的定积分计算方法

给出了一种倾角函数及其导数的定积分计算方法,表达式十分简单,其计算精度:倾角函数可达10^-15,导数可达10^-13,可与Gooding方法相媲美.该方法的稳定性和适用倾角范围均较好,可供倾角函数的最高阶数L_max≤50时使用.     

热海王星系统HD 106315的近$2:1$ 平运动共振捕获与轨道演化

通过结合理论分析和数值模拟方法,可以对热海王星系统HD 106315轨道迁移中的近2:1平运动共振捕获机制以及潮汐作用下的演化过程进行研究.在轨道迁移阶段,初始轨道半长径、初始偏心率以及行星c的偏心率衰减系数K会对系统轨道构型产生影响.数值模拟结果显示当初始轨道半长径分别为ab～0.4 au、ac～0.8 au,偏心率eb和ec均小于0.03时, HD 106315b和HD 106315c在中央恒星的引力作用以及原行星盘粘滞作用下向内迁移, 65000 yr左右两颗行星均可迁移至当前观测位置附近并形成近2:1平运动共振捕获.此外,中央恒星的潮汐效应也可能会对行星系统共振构型产生影响,理论分析表明当行星潮汐耗散系数Q=100时,潮汐效应造成的轨道半长径衰减使系统轨道周期比发生的变化可能是系统脱离共振构型的原因.数值模拟结果显示, HD 106315系统内两颗行星Q103时,来自中央恒星的潮汐效应并不会使行星系统产生明显的偏心率和轨道半长径衰减,不足以使HD 106315行星系统在剩余寿命内脱离2:1平运动共振轨道构型.     

数据深度加权半参数估计法在方差分量估计中的应用研究

针对天地基联合测控中多源异质数据存在多耦合性系统模型误差和粗差的情况,提出了基于数据深度加权半参数估计的方差分量估计法.利用半参数估计法补偿不确定性模型误差,同时构建数据深度权抑制粗差的影响,使半参数估计不受粗差的影响,为方差分量估计提供合理的权矩阵,使得方差分量估计能对存在不确定性模型误差和粗差的结构定权.进行了天地基对空间目标的联合测控仿真试验.结果表明,改进的方差分量估计法能有效补偿不确定性模型误差并抑制粗差的影响,确定更加合理的权值.     

Analysis of Land Water Storage in Southwest China Based on GRACE Data

The inversion of the variation in the land water storage in Southwest China is carried out by taking advantage of the data obtained by the earth gravity satellite GRACE (Gravity Recovery and Climate Experiment) for 64 months from January 2005 to April 2010. The result shows that by selecting an appropriate Gauss radius (R = 600 km) and taking the average gravitational field of the adopted data as the back-ground gravitational field, the land water storage in Southwest China inverted on the basis of the GRACE data reflects the drought in Southwest China at the beginning of 2010 very well.     

The total number of giant planets in debris discs with central clearings

Infrared spectra from the Spitzer Space Telescope ( SSC ) of many debris discs are well fit with a single blackbody temperature which suggest clearings within the disc. We assume that clearings are caused by orbital instability in multiple planet systems with similar configurations to our own. These planets remove dust-generating planetesimal belts as well as dust generated by the outer disc that is scattered or drifts into the clearing. From numerical integrations, we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10⁸ yr old debris disc with a dust disc edge at a radius of 50 au hosted by an A star must contain approximately five Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disc systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However, an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.