半隐Euler法和隐中点法嵌入混合辛积分器的比较
Comparison of second-order mixed symplectic integrator between semi-implicit Euler method and implicit midpoint rule作者机构：南昌大学物理系南昌330031
出 版 物：《物理学报》 (Acta Physica Sinica)
年 卷 期：2011年第60卷第9期
页 面：104-115页
核心收录：
学科分类：07[理学] 070201[理学-理论物理] 0702[理学-物理学]
基 金：国家自然科学基金(批准号:10873007) 南昌大学创新团队项目资助的课题~~
摘 要：当Hamilton函数分解为可积和不可积两部分时,前者能用分析方法给出解析解,而后者可借助一阶半隐Euler法或二阶隐中点法等数值求解,将这种解析和数值解法组合能构造二阶混合辛积分器.理论分析表明Euler嵌入法的稳定区要小于中点嵌入法的.再分别以圆形限制性三体问题和相对论自旋致密双星后牛顿Hamilton构型为例,详细比较了两嵌入法的性能特点.二者的数值精度、稳定性及计算效率与Hamilton的分解方式和轨道类型有关.就圆形限制性三体问题而言,当Hamilton采用势能和含坐标与动量混合项在内的动能分解时,无论对有序还是混沌轨道两嵌入法的数值性能基本相同;然而,当将含坐标与动量混合项做一个独立部分而其余的为另一个部分去分解Hamilton时Euler嵌入法数值稳定性劣于中点嵌入法.针对自旋致密双星来说,Euler嵌入法比中点嵌入法也具有明显的计算效率优势;当两嵌入法均保持数值稳定时二者精度没有显著差别;可是,在混沌存在等情况下Euler嵌入法大多变得数值不稳定.综合表明中点嵌入法具有较好的数值稳定性,尤其值得推荐用来求解后牛顿Hamilton动力学问题.
维普期刊数据库
同方期刊数据库