共查询到19条相似文献，搜索用时 140 毫秒
1.
《天文学进展》2021,(2)
辛算法作为研究哈密顿系统长期定性演化的最佳积分工具,自问世以来就受到了很大的关注。通过对哈密顿函数的截断误差分析,可以从不同角度构造出较高精度的辛算法,也可以通过引入正规化技术实现自动调整积分步长和改善数值稳定性。从辛算法的表现形式可以将它分为显式和隐式两种。当哈密顿系统能够分解为几个可积部分且每部分的解能用时间显函数来表示时,可以构造显式算法。显式算法有非力梯度显式辛算法、力梯度辛算法、辛校正、类高阶辛算法四种。当哈密顿系统变量不能分离时,适合应用隐式辛算法和扩充相空间对称算法求解。分别对这些算法的构造方法及其适用的物理模型进行归纳对比,分析了各种辛算法的优劣性和发展趋势,对如何选择辛算法高效高精度地解决实际问题提供了一定的理论和数值计算依据。 相似文献
2.
当史瓦西黑洞周围存在渐近均匀的外部磁场时, 描述带电粒子在史瓦西黑洞附近运动的哈密顿系统会变为不可积系统. 类似于这样的相对论哈密顿系统不存在有显式分析解的2部分分离形式, 给显式辛算法的构建和应用带来困难. 近一年以来的系列工作提出将相对论哈密顿系统分解为具有显式分析解的2个以上分离部分形式, 成功解决了许多相对论时空构建显式辛算法的难题. 最近的工作回答了哈密顿系统显式可积分离数目对长期数值积分精度有何影响、哪种显式辛算法有最佳长期数值性能这两个问题, 指出哈密顿有最小可积分离数目即3部分分裂解形式并且应用于优化的4阶分段龙格库塔显式辛算法可取得最好精度. 由此选择上述数值积分方法并利用庞加莱截面、最大李雅普诺夫指数和快速李雅普诺夫指标研究在磁化史瓦西黑洞附近运动的带电粒子轨道动力学. 结果显示: 针对某特定的粒子能量和角动量, 较小的外部磁场很难形成混沌轨道; 较大的正磁场参数容易使轨道产生混沌, 并且随着磁场的增大, 轨道的混沌程度也随之加强; 粒子能量适当变大也可以加剧混沌程度, 但负磁场参数和粒子角动量变大都会减弱混沌. 相似文献
3.
4.
约束条件和数值积分 总被引：3，自引：1，他引：2
自治的哈密顿系统存在约束条件，例如能量积分或广义相对论中的4速度大小为常数，它能否在数值积分过程中始终满足将直接影响数值稳定性．在牛顿力学中哈密顿系统的动能一般为椭圆型，直接运用约束条件对方程进行降阶存在开平方判断正负号的困难，导致应用高精度的经典数值积分器时能量存在耗散．然而相对论力学的度规为双曲型，利用约束条件有可能实行方程降阶．在时空具有一定对称性的情况下，能够找到整个时空的一个全局变换使变换后的度规的主对角线某一元素为零，于是从约束方程中不需开平方能够解出某一动量，顺利实现运动方程的降阶．相对论力学中另一个可以降阶的模型是Mixmaster宇宙模型．数值实验表明将经典算法用于降阶后的运动方程能够严格地满足约束，但不一定能保持辛结构。 相似文献
5.
刘林 《紫金山天文台台刊》1997,16(2):82-94
太阳系小天体的运动对应—哈密顿(Hamilton)系统,对其轨道演化的数值研究宜采用哈密顿算法(即辛算法)。本文将仔细讨论这一问题,并以主带小行星的运动为例,较系统地介绍几种辛算法对应的显式辛差分格式。 相似文献
6.
7.
辛算法在动力天文中的应用（Ⅲ） 总被引：3，自引：0，他引：3
文［１］和文［２］从哈密顿系统的整体结构保持一角度阐明了辛算法［３－６］的主要功能，本文将从定量的角度进一步表明辛算法的另一独特优点－可以控制天体运动沿迹误差的快速增长，并对可分离哈密顿系统的显式辛差分格式稍加改进，推广应用到一般动力系统，该系统含有小耗散项或小的不可分离项，计算结果表明，效果极佳，因此，辛算法与传统的数值解法相比，确有很多优点。 相似文献
8.
李焱 《中国天文和天体物理学报》1989,(4)
本文发展了一种解恒星线性非绝热非径向脉动问题的退耦化方法。这个方法把非绝热非径向脉动问题的六阶线性微分方程,分解为由一个代数方程联系起来的一个四阶线性微分方程和另一个二阶线性微分方程进行数值求解。这样的一个退耦处理,有利于克服以前在数值解这类问题时常常遇到的收敛域小和收敛速度慢等困难,并且为数值解方程时所采用的Henyey方法提供了一个自然和方便的初始猜测解。 相似文献
9.
本文说明对称法在一定意义下是一种积分线性保守动力系统的辛方法．并指出在ｈλ的左半复平面上存在有一个对称法的相对稳定区域．这里，ｈ是步长而λ是动力系统的特征根．对称法适用于不显含速度的非线性Ｈａｍｉｌｔｏｎ系统，但不适用于显含速度的系统． 相似文献
10.
考虑周期解的数值延拓问题并提出基于Broyden拟牛顿法来延拓周期解的一种有效算法,先后以布鲁塞尔振子、平面圆型限制性三体问题(Planar Circular Restricted Three-Body Problem, PCRTBP)的周期解为例进行了验证.这里的Broyden方法包含线性搜索、正交三角分解求线性方程组的步骤.对一般的周期解,周期性条件方程组中含有周期作为待延拓参数,可用周期来决定积分时长,将解代入周期性条件得到积分型的非线性方程组,利用Broyden方法迭代延拓直至初值收敛.根据两次垂直通过一个超平面的轨道是对称周期轨道的性质,可采用插值的方法求得再次抵达超平面的解分量,得到周期性条件方程组,再用Broyden方法求解.结合哈密顿系统的对称性和PCRTBP周期轨道的一些分类,对2/1、3/1的内共振周期解族进行了数值研究.最后,对算法和计算结果做了总结和讨论. 相似文献
11.
《Chinese Astronomy and Astrophysics》2007,31(2):172-186
In this paper, we analyze the linear stabilities of several symplectic integrators, such as the first-order implicit Euler scheme, the second-order implicit mid-point Euler difference scheme, the first-order explicit Euler scheme, the second-order explicit leapfrog scheme and some of their combinations. For a linear Hamiltonian system, we find the stable regions of each scheme by theoretical analysis and check them by numerical tests. When the Hamiltonian is real symmetric quadratic, a diagonalizing by a similar transformation is suggested so that the theoretical analysis of the linear stability of the numerical method would be simplified. A Hamiltonian may be separated into a main part and a perturbation, or it may be spontaneously separated into kinetic and potential energy parts, but the former separation generally is much more charming because it has a much larger maximum step size for the symplectic being stable, no matter this Hamiltonian is linear or nonlinear. 相似文献
12.
We construct an explicit reversible symplectic integrator for the planar 3-body problem with zero angular momentum. We start with a Hamiltonian of the planar 3-body problem that is globally regularised and fully symmetry reduced. This Hamiltonian is a sum of 10 polynomials each of which can be integrated exactly, and hence a symplectic integrator is constructed. The performance of the integrator is examined with three numerical examples: The figure eight, the Pythagorean orbit, and a periodic collision orbit. 相似文献
13.
The time-transformed leapfrog scheme of Mikkola Aarseth was specifically designed for a second-order differential equation with two individually separable forms of positions and velocities.It can have good numerical accuracy for extremely close two-body encounters in gravitating few-body systems with large mass ratios,but the non-time-transformed one does not work well.Following this idea,we develop a new explicit symplectic integrator with an adaptive time step that can be applied to a time-dependent Hamiltonian.Our method relies on a time step function having two distinct but equivalent forms and on the inclusion of two pairs of new canonical conjugate variables in the extended phase space.In addition,the Hamiltonian must be modified to be a new time-transformed Hamiltonian with three integrable parts.When this method is applied to the elliptic restricted three-body problem,its numerical precision is explicitly higher by several orders of magnitude than the nonadaptive one's,and its numerical stability is also better.In particular,it can eliminate the overestimation of Lyapunov exponents and suppress the spurious rapid growth of fast Lyapunov indicators for high-eccentricity orbits of a massless third body.The present technique will be useful for conservative systems including N-body problems in the Jacobian coordinates in the the field of solar system dynamics,and nonconservative systems such as a time-dependent barred galaxy model in a rotating coordinate system. 相似文献
14.
辛积分器中沿迹误差的一种补偿方法 总被引：2，自引：0，他引：2
辛积分器严格描述了一摄动Ｈａｍｉｌｔｏｎ系统的流，因而导致天体轨道的沿迹误差随时间呈线性增长趋势。本文利用这一特点，提出了一种对其沿迹误差进行估算的数值方法，从而达到了对数值结果进行沿迹误差补偿的目的，数值结果证实了此方法在较大积分步长和较长积分时间的数值计算中是有效的。 相似文献
15.
A recurrent method of solving the formal integrals of symplectic integrators is given. The special examples show that there are no long-term variations in all integrals of the Hamiltonian system in addition to the energy one when symplectic integrators are used in the numerical studies of the system. As an application of the formal integrals, the relation between them and the linear stability of symplectic integrators is discussed. 相似文献
16.
Jean-Marc Petit 《Celestial Mechanics and Dynamical Astronomy》1998,70(1):1-21
We investigate the numerical implementation of a symplectic integrator combined with a rotation (as in the case of an elongated
rotating primary). We show that a straightforward implementation of the rotation as a matrix multiplication destroys the conservative
property of the global integrator, due to roundoff errors. According to Blank et al. (1997), there exists a KAM-like theorem
for twist maps, where the angle of rotation is a function of the radius. This theorem proves the existence of invariant tori
which confine the orbit and prevent shifts in radius. We replace the rotation by a twist map or a combination of shears that
display the same kind of behaviour and show that we are able not only to recover the conservative properties of the rotation,
but also make it more efficient in term of computing time. Next we test the shear combination together with symplectic integrator
of order 2, 4, and 6 on a Keplerian orbit. The resulting integrator is conservative down to the roundoff errors. No linear
drift of the energy remains, only a divergence as the square root of the number of iterations is to be seen, as in a random
walk. We finally test the three symplectic integrators on a real case problem of the orbit of a satellite around an elongated
irregular fast rotating primary. We compare these integrators to the well-known general purpose, self-adaptative Bulirsch–Stoer
integrator. The sixth order symplectic integrator is more accurate and faster than the Bulirsch–Stoer integrator. The second-
and fourth- order integrators are faster, but of interest only when extreme speed is mandatory.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
17.
We present a generalisation of the Levi-Civita and Kustaanheimo-Stiefel regularisation. This allows the use of more general time rescalings. In particular, it is possible to find a regularisation which removes the singularity of the equations and preserves scaling invariance. In addition, these equations can, in certain cases, be integrated with explicit symplectic Runge-Kutta-Nyström methods. The combination of both techniques gives an explicit adaptive symplectic (EASY) integrator. We apply those methods to some perturbations of the Kepler problem and illustrate, by means of some numerical examples, when scaling invariant regularisations are more efficient that the LC/KS regularisation. 相似文献
18.
An explicit symplectic integrator is constructed for the problem of a rotating planetary satellite on a Keplerian orbit. The
spin vector is fixed perpendicularly to the orbital plane. The integrator is constructed according to the Wisdom-Holman approach:
the Hamiltonian is separated in two parts so that one of them is multiplied by a small parameter. The parameter depends on
the satellite’s shape or the eccentricity of its orbit. The leading part of the Hamiltonian for small eccentricity orbits
is similar to the simple pendulum and hence integrable; the perturbation does not depend on angular momentum which implies
a trivial ‘kick’ solution. In spite of the necessity to evaluate elliptic function at each step, the explicit symplectic integrator
proves to be quite efficient.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
19.
By adding force gradient operators to symmetric compositions, we build a set of explicit fourth-order force gradient symplectic algorithms, including those of Chin and coworkers, for a separable Hamiltonian system with quadratic kinetic energy T and potential energy V . They are extended to solve a gravitational n-body Hamiltonian system that can be split into a Keplerian part H 0 and a perturbation part H 1 in Jacobi coordinates. It is found that the accuracy of each gradient scheme is greatly superior to ... 相似文献