期刊文献+

太阳系动力学中几种辛和类辛方法的比较 被引量:2

Comparison of Several Symplectic and Quasi-Symplectic Methods in Solar System Dynamics
在线阅读 下载PDF
导出
摘要 辛方法是目前研究太阳系天体轨道长期演化的最佳数值方法 .它能保持哈密顿系统的主要性质而不会引入人为的耗散 ,因而可以用大步长进行长时期的数值积分 .国际文献中用于太阳系动力学定性研究的辛方法主要是Wisdom和Holman的SYA方法 ,它是SYP方法的一种近似 .它们都采用雅可比坐标系 .在质心坐标系中建立的辛方法有传统的将动能和势能分离的SYS方法 .提出一种在质心坐标系中构造的类辛方法SYQ ,并对这 4种方法进行了较为详尽的效率分析和误差比较 ,结论是SYA和SYQ是比较好的方法 。 Symplectic methods are so far the best numerical methods for qualitative exploration in solar system dynamics. They maintain the symplectic structure and key properties of Hamiltonian systems and do not bring in any artificial dissipation, making possible long_term numerical integrations with a large step size. The symplectic method that has been widely adopted in references on qualitative studies of solar system dynamics is the method worked out by Wisdom and Holman (abbreviated as SYA hereafter). It is built in the Jacobian coordinate system and takes an approximation of the hamiltonian. The Wisdom and Holman's method for an exact hamiltonian is abbreviated as SYP. Actually a symplectic integrator can be built in the barycentric coordinate system (abbreviated as SYS), which separates the hamiltonian into two parts, the potential energy and the kinetic energy. Here we propose a quasi_symplectic method SYQ in the barycentric coordinate system. An extensive comparative study of these four types of methods is given, especially on their computation effciency and error accumulation. This research draws the following conclusions. Considering that symplectic integrators are mainly used in exploring the qualitative evolution of dynamical systems and a high precision is not required, SYS should not be recommended in solar system dynamics for its low effciency. During a 10\+8 years integration, SYP methods cause almost the same errors on the positions of the planets but they take almost 40% more computing time. We thus believe that SYP can not compete with SYA or SYQ, but it is hard to tell SYA or SYQ is better. Our research has also shown that resonances play a role in keeping the orbit configuration of a planetary system during long_term numerical integrations.
出处 《南京大学学报(自然科学版)》 CAS CSCD 北大核心 2002年第4期462-470,共9页 Journal of Nanjing University(Natural Science)
基金 国家自然科学基金 (10 1730 0 7) 紫金山天文台小行星基金会
关键词 类辛方法 太阳系动力学 数值方法 辛方法 SYA方法 SYQ方法 天体轨道 质心坐标系 solar system dynamics, numerical methods, symplectic methods
  • 相关文献

参考文献11

  • 1Ruth R D. A canonical integration technique,IEEK Trans Nucl Sci, 1983, NS- 30 : 2669-2671.
  • 2Feng K. Proc 1984 Beijing symp on diff geometry and diff Eq. Beijing: Science Press, 1985, 42.
  • 3Wisdom J, Holman M J. Sympletic maps for the n-body problem. Astron J, 1991, 102: 1528-1538.
  • 4Sussman G J, Wisdom J. Numerical evidence that the motion of Pluto is chaotic. Science, 1988, 241:433-437.
  • 5Quinlan G D, Tremine S. Symmetric multistep methods for the numerical integration of planetary orbits. Astron J, 1990, 100:1694-1700.
  • 6Huang Tian-Yi, Wang Chang-Bin, Zhao Zhang-Yin. The relation between symmetric and symplectic methods.Chin Astron Astrophys, 1998,22:101-112(in Chinese).
  • 7Cohen C J, Hubbard E C. Libration of the dose approaches of Pluto to Neptune. Astron J, 1965, 70:10-13.
  • 8Williams J G,Benson G S. Resonances in the Neptune-Pluto system. Astron J , 1971, 76:167-176.
  • 9Kinoshita H, Nakai H. Long-term behavior of the motion of Pluto over 5.5 billion years. Earth Moon Planet,1996, 72:165-173.
  • 10Wan X-S, Huang T-Y, Innanen K. The 1:1 super resonance in Pluto's motion. Astron J, 2001, 121 : 1155-1162.

同被引文献107

  • 1朱耀仲,吴斌,彭碧波.月球平均运动和地球自转速率长期变化的潮汐耗散[J].测绘学报,2000,29(Z1):1-4.
  • 2杨远玲,聂清香,吴晓梅,徐顺福.N体问题的几种数值算法比较[J].计算物理,2006,23(5):599-603. 被引量:7
  • 3E Fehlberg. Classical Fifth - , Sixth -, Seventh -and Eighth Order Runge - Kutta Formulas with Step Size Control[ R]. NASA TR R2287,1968.
  • 4J Stoer, R Bulirsch. Introduction to Numerical Analysis[ M]. New York: Springer-Verlag, 1980.
  • 5G D Quinlan, S Tremaine. Symmetric muhistep methods for the numerical integration of planetary orbits[ J]. Astron J, 1990,100 : 1694 - 1700.
  • 6R D Ruth. A canonical integration technique[ J]. IEEE Trans. Nucl. Sci., 1983, 30:2669-2671.
  • 7Feng K, M Quin. The symplectic methods for the computation of Hamiltonian equations[ J]. Lecture Notes in Math, 1987,1297 : 1 -37.
  • 8Jack Wisdom, Matthew Holman. Symplectic maps for the n - body problem[J]. Astron. J. , 1991,102:1528 - 1538.
  • 9J E Chambers. A hybrid symplectic integrator that permits close encounters between massive bodies [ J]. Month. Not. Roy. Astron. Soc. , 1999,304:793 - 799.
  • 10D Gauchez, M Fouchard, J Souchay. Inter- comparisons between numerical integrations and analytical theory for the solar system. Modern Celestial Mechanics: From Theory to Applications[M]. Kluwer Academic Publishers, 2002.

引证文献2

二级引证文献1

相关作者

内容加载中请稍等...

相关机构

内容加载中请稍等...

相关主题

内容加载中请稍等...

浏览历史

内容加载中请稍等...
;
使用帮助 返回顶部