孤立波方程的保结构算法.PDF

第 21 卷 第 5 期 计 　　算 　　物 　　理 Vol . 21 ,No . 5 2004 年 9 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Sep . , 2004 　 ( ) [文章编号] 　1001246X 2004 孤立波方程的保结构算法 1 ,2 2 2 王雨顺 , 　王 　斌 , 　季仲贞 ( 1 南京师范大学数学与计算机学院 , 江苏 南京 　210097 ; 2Lasg , 中国科学院大气物理所 , 北京 　100029) [摘 　要] 　讨论了孤立波方程的保结构差分算法 , 以一些经典的孤立波方程为例 , 如 KdV , sineGordon , KP 方 程 , 给出了它们的辛和多辛结构 , 说明辛和多辛算法的可适用性. 提出局部守恒算法和广义保结构算法的概念 , 它们是保结构算法的概念 自然推广. 还给出一种能系统构造局部守恒格式的复合方法. 数值例子说明 , 保结构 数值能很好模拟各种孤立波的演化. [ 关键词] 　孤立波方程 ; 保结构算法 ; 复合方法 ; 局部守恒格式 ; 数值实验 [ 中图分类号] 　O241 [文献标识码] 　A 0 　引言 随着计算机的飞速发展 ,数值方法越来越显示出其优越性和生命力. 但一些具体复杂问题 ,如分子化学 , 复杂流体计算 ,它们计算规模也越来越大 ,积分时间也越来越长 ,因此对数值计算方法的要求也越来越高. 对 这些特殊的问题 ,人们希望找到较好的数值方法. 冯康于 1984 年提出 ,能否尽可能多的保持原问题的本质特 [ 1] 性是评价一个数值方法好坏的标准. 基于这样的认识 ,冯康提出了保结构算法的概念 . 这里的保结构是指 [2 ] 保几何结构. 具体有三类 ,分别是哈密尔顿系统的辛算法 、无源系统的保体积算法 、接触系统的接触算 法[3 ] . 有限维 Hamilton 系统的几何结构和辛算法的理论发展很快 ,现在已日益完善. 冯康课题组非常漂亮的工 作之一是发展并完善了能系统构造任意阶精度的辛算法 ,该方法从离散辛变换的生成函数入手 ,通过 Hamil [4 ] tonJ acobi 方程对应得到原辛变换的一个逼近 ,即数值方法 . 1988 年 ,SanzSerna 等人发现满足 GaussLegend re 条件的 RungeKutta 方法是辛算法 ,这不仅可以合理解释这些数值算法的优越性 ,也为构造新的辛算法提 供了思路[5 ] . 由于 RungeKutta 方法在数值计算方法中的特殊地位 ,辛算法引起了计算界的极大兴趣 ,兴起一 个研究热潮 ,得到许多好的研究结果. 另外还有两类重要的构造方法 ,组合方法[6 ,7 ] 和变分方法. 辛算法的理 论分析工作也得到了系统地发展 ,主要有 :离散的 KAM 理论[8 ]