对数哈密顿方法及其应用

来源 :南昌大学 | 被引量 : 0次 | 上传用户:dozen
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
天体力学数值方法作为天体力学的重要领域之一在辛算法的提出后得到长足发展,辛算法保持哈密顿系统辛结构且计算过程中系统没有能量和角动量的长期误差累积。辛算法适用于哈密顿系统的长期定性演化研究同时也具有数值精度不高、显辛算法要求固定步长的不足。通常积分计算天体紧密交汇问题或大偏心率轨道运动都需缩短步长来克服天体受引力过大而剧增的加速度,直接变步长将丢失辛算法保持辛结构的优势,考虑时间变换的思路,原时间变量取变步长而新的时间变量仍为固定步长,则既能调节步长又能保持辛算法固有优势。本文的主要内容为构造针对不同哈密顿系统的对数哈密顿算法及论证其在具有更高的数值精度和保证获得有效的混沌判别结果方面的优势。  针对不同的哈密顿系统结构构造不同形式的时间变换辛算法。对于可分解为分别只含状态量广义动量和广义坐标的动能部分和势能部分的哈密顿函数,可构造取时间变换函数为形式不同但等价的两个函数得到显式对数哈密顿方法,其中时间变换作用于哈密顿函数,本文构造了由三个二阶蛙跳算子构成的显式对数哈密顿Yoshida四阶方法。对于动能部分具有广义动量和广义坐标的交叉项而势能部分仅含位置变量的系统构造显隐式混合对数哈密顿方法,对于动能部分应用隐式中点法。而对于更一般的系统则构造隐式对数哈密顿方法。隐式方法具有更广泛的应用但也由于算法构造中包括迭代需耗费更多的计算机时间降低计算效率。  本文详细论证了显式对数哈密顿方法在应用于牛顿圆型限制性三体问题及相对论圆型限制性三体问题时较于非时间变换辛算法更具数值精度优势。且在前一系统的精度优势独立于轨道偏心率的变化。对于后一系统这一现象未能发生但数值精度也明显优越于常规辛算法。特别对于高偏心率轨道,非时间变换算法得到的虚假的混沌判别指标,如Lyapunov指标和快速Lyapunov指数(FLI)。而通过对数哈密顿方法则可获得可靠地定性分析结果,彻底地解决后牛顿圆型限制性三体问题的高偏心率轨道Lyapunov指数的过度估计和FLI快速增大的问题。在得到论证后本文应用对数哈密顿方法讨论了动力学参数两主天体间距离的变化对动力学系统有序和混沌转化的影响。本文通过数值模拟验证了对数哈密顿方法具有更高的数值精度及可得到可靠的定性研究成果的优势。适用于定性研究和定量计算高偏心率问题,为天体力学研究开拓了新思路。在实际的天体紧密交汇处的动力学演化提供反映动力学实质的积分工具。
其他文献
自从上世纪初发现宇宙线以来,宇宙线的产生、加速机制问题一直是一个尚未解决的粒子天体物理的基本问题。宇宙线膝区(1014~1016eV)对这一百年难题十分敏感,对膝区的成因做出的正
M87作为第一个被发现的非blazar(耀变类星体)的AGN(活动星系核)类TeV(1012eV)源,它有着快速的TeV耀变,以及非常平的TeV光谱。类似blazars,M87作为一个射电星系,其核区的辐射虽然
大质量恒星在星系的演化中扮演着重要的角色。它们是星系中重子物质和UV辐射的主要来源。然而,对于大质量恒星的形成和演化现在还知之甚少。首先对其起始条件不清楚;其次,对
多年以来宇宙起源一直是一个热门话题,特别是在宇宙学奇点问题上。本文给出了一个可以避免宇宙奇点的模型,即反弹暴涨模型,并用最新的观测数据对模型进行限制。该模型描述了宇宙
如果奇异夸克物质为宇宙中最稳定的基态这一假设成立,那么完全由夸克物质组成的致密天体——奇异夸克星,将成为致密星旗下的另一个有别于中子星的分支。通过对r-模不稳定性的研
描述哈密顿系统的混沌运动需要依赖于可靠的数值方法和混沌识别方法。本文主要工作在于数值方法的开发与利用。   我们在Chin等提出的四阶力梯度辛算法基础上构造了两类显
据最近的研究表明:一个低阶的后牛顿近似拉格朗日函数动力学的方程式可以通过泰勒级数展开到无穷高阶的后牛顿欧拉拉格朗日方程式,并且所展开的无穷高阶的后牛顿欧拉拉格朗日方
在不引入暗物质的情况下,MOND可以成功地解释众多星系动力学现象。但是MOND只是一个唯象理论,为了给出一个完备的理论,后续发展了一些.非相对论与相对论理论。   宇宙大尺