期刊文献+
共找到7篇文章
< 1 >
每页显示 20 50 100
一个膺三阶辛积分器 被引量:4
1
作者 刘福窑 伍歆 陆本魁 《天文学报》 CSCD 北大核心 2004年第4期402-412,共11页
在太阳系动力学中,辛积分器已成为研究哈密顿系统的长期定性演化的最佳工具.对于可积分离的哈密顿系统H=H0+∑i=1N∈iHi(∈≤1),构造了一个膺三阶辛积分器.它大约相当于Wisdom-Holman二阶辛积分器的一次校正或Forest-Ruth四阶辛算法的精... 在太阳系动力学中,辛积分器已成为研究哈密顿系统的长期定性演化的最佳工具.对于可积分离的哈密顿系统H=H0+∑i=1N∈iHi(∈≤1),构造了一个膺三阶辛积分器.它大约相当于Wisdom-Holman二阶辛积分器的一次校正或Forest-Ruth四阶辛算法的精度.此外,含力梯度的辛算法也适合处理哈密顿系统H=Ho(q,P)+∈H1(q),其精度好于原辛积分器,但不优越于相应膺高阶辛积分器. 展开更多
关键词 辛算法 哈密顿系统 可积 四阶 二阶 高阶 梯度 积分器 校正
下载PDF
辛方法的校正公式 被引量:6
2
作者 伍歆 黄天衣 万晓生 《天文学报》 CSCD 北大核心 2002年第4期391-402,共12页
1996年Wisdom等提出了对辛方法进行校正的概念和实践.现在继续对辛校正进行详尽讨论和数值比较,尤其对哈密顿函数可分解为一个主要部分和多个次要部分的一般情形,用Lie级数推导任意阶的各种辛算法的一次和二次辛校正公式并对一些算法给... 1996年Wisdom等提出了对辛方法进行校正的概念和实践.现在继续对辛校正进行详尽讨论和数值比较,尤其对哈密顿函数可分解为一个主要部分和多个次要部分的一般情形,用Lie级数推导任意阶的各种辛算法的一次和二次辛校正公式并对一些算法给出具体的辛校正公式.又以日、木、土三体问题为模型进行数值实验,结果表明一次辛校正能提高精度,改善数值稳定性,计算效率也比较高,因而值得推荐使用.辛方法通常用大步长数值积分,这时二次辛校正并没有显著提高结果的精度,却大大增加了计算时间,不应予以推荐. 展开更多
关键词 天体力学 数值积分 辛方法 哈密顿函数
下载PDF
约束条件和数值积分 被引量:4
3
作者 伍歆 黄天衣 《天文学报》 CSCD 北大核心 2004年第3期310-319,共10页
自治的哈密顿系统存在约束条件,例如能量积分或广义相对论中的4速度大小为常数.它能否在数值积分过程中始终满足将直接影响数值稳定性.在牛顿力学中哈密顿系统的动能一般为椭圆型,直接运用约束条件对方程进行降阶存在开平方判断正负号... 自治的哈密顿系统存在约束条件,例如能量积分或广义相对论中的4速度大小为常数.它能否在数值积分过程中始终满足将直接影响数值稳定性.在牛顿力学中哈密顿系统的动能一般为椭圆型,直接运用约束条件对方程进行降阶存在开平方判断正负号的困难,导致应用高精度的经典数值积分器时能量存在耗散.然而相对论力学的度规为双曲型,利用约束条件有可能实行方程降阶.在时空具有一定对称性的情况下,能够找到整个时空的一个全局变换使变换后的度规的主对角线某一元素为零,于是从约束方程中不需开平方能够解出某一动量,顺利实现运动方程的降阶.相对论力学中另一个可以降阶的模型是Mixmaster宇宙模型.数值实验表明将经典算法用于降阶后的运动方程能够严格地满足约束,但不一定能保持辛结构. 展开更多
关键词 约束条件 数值积分 哈密顿系统 广义相对论 Mixmaster宇宙模型
下载PDF
两类辛算法的稳定性比较 被引量:1
4
作者 刘福窑 钱晓明 《南昌大学学报(理科版)》 CAS 北大核心 2009年第2期123-126,共4页
主要对一阶隐式Euler法、二阶隐式Euler中点公式分别嵌入到一阶显辛Euler法和二阶leapfrog显辛积分器所得算法M1和M2进行了稳定性比较。对于坐标与速度耦合的线性哈密顿系统,理论分析寻找稳定区比较困难,依靠数值方法定性地比较了他们... 主要对一阶隐式Euler法、二阶隐式Euler中点公式分别嵌入到一阶显辛Euler法和二阶leapfrog显辛积分器所得算法M1和M2进行了稳定性比较。对于坐标与速度耦合的线性哈密顿系统,理论分析寻找稳定区比较困难,依靠数值方法定性地比较了他们的稳定性。结果表明M1算法一般优于M2算法。 展开更多
关键词 天体力学 辛算法 稳定性
下载PDF
天文动力学方程数值积分中的一种有效变步法 被引量:2
5
作者 刘承志 崔斗星 《天文学报》 CSCD 北大核心 2002年第4期387-390,共4页
利用积分曲线的曲率控制步长的技巧,使天文动力学方程数值解法的精度和速度有较大提高.这种方法适用于天体精密定轨以及一些精度要求高的常微分方程初值问题的数值积分.
关键词 天文动力学方程 天体精密定轨 曲率 数值积分 天体力学
下载PDF
稳定化方法与后稳定化方法的比较
6
作者 朱军芳 伍歆 《天文学报》 CSCD 北大核心 2006年第3期291-300,共10页
对BaumgaLrte的稳定化和Chin的后稳定化进行了详尽讨论与数值比较.用经典数值方法并结合这两种稳定化方式都能提高数值精度和改善数值稳定性.在最佳稳定参数下稳定化精度一般不等价于后稳定化.两者精度优劣并无常定.考虑到Baumgart... 对BaumgaLrte的稳定化和Chin的后稳定化进行了详尽讨论与数值比较.用经典数值方法并结合这两种稳定化方式都能提高数值精度和改善数值稳定性.在最佳稳定参数下稳定化精度一般不等价于后稳定化.两者精度优劣并无常定.考虑到Baumgarte的稳定化使得数值积分的右函数更复杂和增加计算耗费,尤其是存在稳定参数最佳选取的麻烦,故推荐后稳定化投入实算.但值得注意的是用后稳定化与没有经过稳定化处理的经典积分器来比不宜扩大积分步长. 展开更多
关键词 天体力学 数值积分 稳定化 后稳定化
下载PDF
Optimized third-order force-gradient symplectic algorithms 被引量:3
7
作者 LI Rong & WU Xin School of Science,Nanchang University,Nanchang 330031,China 《Science China(Physics,Mechanics & Astronomy)》 SCIE EI CAS 2010年第9期1600-1609,共10页
With the natural splitting of a Hamiltonian system into kinetic energy and potential energy,we construct two new optimal thirdorder force-gradient symplectic algorithms in each of which the norm of fourth-order trunca... With the natural splitting of a Hamiltonian system into kinetic energy and potential energy,we construct two new optimal thirdorder force-gradient symplectic algorithms in each of which the norm of fourth-order truncation errors is minimized.They are both not explicitly superior to their no-optimal counterparts in the numerical stability and the topology structure-preserving,but they are in the accuracy of energy on classical problems and in one of the energy eigenvalues for one-dimensional time-independent Schrdinger equations.In particular,they are much better than the optimal third-order non-gradient symplectic method.They also have an advantage over the fourth-order non-gradient symplectic integrator. 展开更多
关键词 SYMPLECTIC INTEGRATORS SYMPLECTIC scheme-shooting METHOD celestial mechanics time-independent Schrdinger equation energy eigenvalues numerical stability BISECTION METHOD topological structure
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部