任意阶显式精细积分多步法的常用形式及其高阶次数值计算

计算物理 ›› 2004, Vol. 21 ›› Issue (3): 333-338.

任意阶显式精细积分多步法的常用形式及其高阶次数值计算

闫海青, 唐晨, 张皞, 刘铭, 张桂敏

天津大学应用物理系, 天津 300072

收稿日期:2003-02-08 修回日期:2003-08-05 出版日期:2004-05-25 发布日期:2004-05-25
作者简介:闫海青(1962-),男,辽宁,副教授,硕士,主要从事凝聚态物理与计算物理方面的研究工作.
基金资助:
上海交通大学振动、冲击、噪声国家重点实验室基金(VSN-2003-03)资助项目

Common Formulae for Free-order Explicit Multistep Method of Precise Time Integration and the Higher Order Numerical Simulation

YAN Hai-qing, TANG Chen, ZHANG Hao, LIU Ming, ZHANG Gui-min

Department of Applied Physics, Tianjin University, Tianjin 300072, China

Received:2003-02-08 Revised:2003-08-05 Online:2004-05-25 Published:2004-05-25

摘要/Abstract

摘要： 基于任意阶显式精细积分多步法的一般公式,给出其几种常用形式,并实现了高阶次数值计算,将新算法应用于射线方程和双原子系统经典轨迹数值计算中.数值计算结果表明任意阶显式精细积分多步法是一种高精度、高效率、稳定性较好的方法,并且可方便地进行高阶次的运算.

关键词: 任意阶显式精细积分多步法, 高阶次数值计算, 经典轨迹, 稳定性分析

Abstract: Common formulae for the free-order explicit multistep method of precise time integration are proposed. When the higher order explicit algorithms of precise time integration are applied for calculating the ray equation and the classical trajectories of diatomic system, the effect is admirable. The numerical results reveal that the pressent method is higher accurated and efficient, capable of keeping computational stability for long time simulation, and suitable for higher order numerical computation.

Key words: free-order multi-step method of precise integration, higher order numerical computation, classical trajectory stability analysis

中图分类号:

O414

闫海青, 唐晨, 张皞, 刘铭, 张桂敏. 任意阶显式精细积分多步法的常用形式及其高阶次数值计算[J]. 计算物理, 2004, 21(3): 333-338.

YAN Hai-qing, TANG Chen, ZHANG Hao, LIU Ming, ZHANG Gui-min. Common Formulae for Free-order Explicit Multistep Method of Precise Time Integration and the Higher Order Numerical Simulation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 21(3): 333-338.

[1]	庞军刚, 宋琳, 唐娜, 杨雪滢, 李晓霖, 席忠红, 石玉仁. 磁化尘埃等离子体中的冲击波及其横向扰动稳定性[J]. 计算物理, 2019, 36(6): 682-690.
[2]	沈智军, 胡立军, 闫伟. 二维浅水波方程的数值激波不稳定性[J]. 计算物理, 2012, 29(1): 25-35.
[3]	郭静, 刘世兴, 徐天赋, 刘学深, 丁培柱. 激光场中一维和三维氢分子离子模型经典动力学行为的比较[J]. 计算物理, 2008, 25(4): 470-476.
[4]	赵慧勇, 乐嘉陵. 双时间步方法的应用分析[J]. 计算物理, 2008, 25(3): 253-258.
[5]	匙玉华, 刘学深, 丁培柱. 双原子分子在啁啾激光作用下的经典解离[J]. 计算物理, 2006, 23(4): 489-493.
[6]	刘世兴, 王怀民, 祁月盈, 刘学深, 丁培柱. 强激光场中CO分子经典轨迹的辛算法[J]. 计算物理, 2005, 22(4): 325-328.
[7]	刘嘉, 雷麦芳, 姚文秀, 王发民. 高超声速边界层流动转捩研究[J]. 计算物理, 2004, 21(1): 61-67.
[8]	李桂琴, 周茜, 刘淑琴. 应用Murrell-Sorbie双原子势计算H₂和O₂的经典轨迹[J]. 计算物理, 1999, 16(4): 358-361.
[9]	王强, 傅德薰, 马延文. 可压缩平面混合层稳定性数值计算[J]. 计算物理, 1997, 14(S1): 413-416.
[10]	石爱民, 母英魁, 丁培柱. N₂双原子系统经典轨迹的辛算法计算[J]. 计算物理, 1997, 14(S1): 433-434.