A class of symplectic integrators with adaptive time step for separable Hamiltonian systems

被引:106
作者
Preto, M [1 ]
Tremaine, S [1 ]
机构
[1] Princeton Univ Observ, Princeton, NJ 08544 USA
关键词
celestial mechanics; stellar dynamics; methods : numerical; solar system : general;
D O I
10.1086/301102
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
Symplectic integration algorithms are well suited for long-term integrations of Hamiltonian systems, because they preserve the geometric structure of the Hamiltonian flow. However, this desirable property is generally lost when adaptive time step control is added to a symplectic integrator. We describe an adaptive time step, symplectic integrator that can be used if the Hamiltonian is the sum of kinetic and potential energy components and the required time step depends only on the potential energy (e.g., test-particle integrations in fixed potentials). In particular, we describe an explicit, reversible, symplectic, leapfrog integrator for a test particle in a near-Keplerian potential; this integrator has a time step proportional to distance from the attracting mass and has the remarkable property of integrating orbits in an inverse-square force field with only "along-track" errors; i.e., the phase-space shape of a Keplerian orbit is reproduced exactly, but the orbital period is in error by O(N-2), where N is the number of steps per period.
引用
收藏
页码:2532 / 2541
页数:10
相关论文
共 22 条
[1]  
[Anonymous], ANN NUMER MATH
[2]  
Arnold V. I., 1984, NONLINEAR TURBULENT, V3, P1161
[3]   Variable step implementation of geometric integrators [J].
Calvo, MP ;
Lopez-Marcos, MA ;
Sanz-Serna, JM .
APPLIED NUMERICAL MATHEMATICS, 1998, 28 (01) :1-16
[4]   SYMPLECTIC INTEGRATION OF HAMILTONIAN-SYSTEMS [J].
CHANNELL, PJ ;
SCOVEL, C .
NONLINEARITY, 1990, 3 (02) :231-259
[5]  
Danby JMA., 1988, FUNDAMENTALS CELESTI
[6]   A multiple time step symplectic algorithm for integrating close encounters [J].
Duncan, MJ ;
Levison, HF ;
Lee, MH .
ASTRONOMICAL JOURNAL, 1998, 116 (04) :2067-2077
[7]   Linear multistep methods for integrating reversible differential equations [J].
Evans, NW ;
Tremaine, S .
ASTRONOMICAL JOURNAL, 1999, 118 (04) :1888-1899
[8]   The adaptive Verlet method [J].
Huang, WZ ;
Leimkuhler, B .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1997, 18 (01) :239-256
[9]   BUILDING A BETTER LEAPFROG [J].
HUT, P ;
MAKINO, J ;
MCMILLAN, S .
ASTROPHYSICAL JOURNAL, 1995, 443 (02) :L93-L96
[10]  
Hut P, 1997, ASTR SOC P, V123, P26