一种新的三阶非力梯度辛积分算法  

A New Solution to the Third-Order Non-Gradient Symplectic Integration Algorithm

作　　者：汪文帅[1,2] 李小凡[2] 

机构地区：[1]宁夏大学数学计算机学院,宁夏银川750021 [2]中国科学院地质与地球物理研究所,北京100029

出　　处：《武汉大学学报（理学版）》2012年第3期221-228,共8页Journal of Wuhan University:Natural Science Edition

基　　金：国家自然科学基金项目(41174047;40874024);国家重点基础研究发展计划(973计划)(2007CB209603)资助项目

摘　　要：基于分部的Runge-Kutta离散形式,引入了相位误差最小的思想,给出了一种新的三级三阶非力梯度辛积分算法,并通过数值试验与经典的Ruth、McLachlan&Atela以及Iwatsu的三级三阶非力梯度辛算法从稳定性、长时程、保结构性等方面进行了对比.结果显示新推导的三级三阶非力梯度辛算法稳定性较好、长时程运算误差小,表明该算法具有好的保结构性和较强的长时程跟踪能力.进一步通过数值试验与力梯度辛算法比较,也显示出该算法的有效性和具有较高的精度.In this paper, a new solution to the three-stage third-order non-gradient symplectic integration algorithm is given based on the minimum phase error principle and the partitioned Runge-Kutta form. Several other nongradient methods, such as Ruth＇s, McLachlan＆Atela＇s and Iwatsu＇s symplectic integration algorithm of three-stage third-order, are employed to compare the performance with the present algorithm by numerical experiments. The numerical results show that the present algorithm is good at stability, and much superior to the other methods in the features of long-time computing ability. The further numerical experiments is used to compare the algorithm with above non-gradient symplectic methods and some force-gradient symplectic methods, the results also show that the algorithm is effective with high accuracy. These appealing features of the algorithm would make it effective to have the ability of structure-preserving and long-time tracing.

关 键 词：辛算法 力梯度 稳定性 相位误差 

分 类 号：O343[理学—固体力学]