MKdV方程的多辛Fourier拟谱方法.pdfVIP

2009年 6月 应用数学与计算数学学报 第 23卷 第 1期 June，2009 C0MM．0N APPL．MATH．ANDC0NIPUT V01．23 No．1 MKdV方程的多辛 Fourier拟谱方法 鞠瑞亮 马和平 张中强 摘要基于谱微分矩阵方法，给出MKdV方程的多辛Fourier拟谱格式及其相应多辛 离散守恒律，证明了它等价于通常的Fourier拟谱格式．数值结果表明，格式对于长时间 计算具有稳定性与高精度． 关键词 多辛，谱微分矩阵，Fourier拟谱 M ultisym plecticFourierPseud0spectralM ethods fortheM KdV Equation JuRuiliang MaHeping ZhangZhongqiang Abstract A multisymplecticFourierpseudospectralscheme，basedonthedifferenti— ationmatrixmethod，ispresentedfortheMKdV equation．Itisshownthattheschemeis equivMenttoatraditionalFourierpseudospectralscheme．Numericalexperimentsshow theschemeisstable，ofhigh accuracy，and especially suitableforlargetime behavior simulation． Keywordsmultisymplectic，spectraldifferentiationmatrix，Fourierpseudospectral 1 引 言 上世纪80年代，冯康先生 [】提出基于Hamiltion系统辛几何的辛算法．此后，辛算 法得到了很大的发展．由于辛算法较传统数值方法具有长时间的数值稳定性，因此， 越来越受到人们的重视．目前 已经推广到许多科学领域，如量子力学、天体力学、海 洋科学、分子动力学等，取得 了很好的数值结果． 作为辛算法的延伸，Bridges等 [2-5j提出了偏微分方程的多辛算法．许多方程，例 如KdV方程，Kadomtsev。Petviashvili方程，Zabolotskaya-Khokhlov方程，sine．Gordon 方程等均可表示为多辛形式，并且可构造多辛格式取得很好的计算效果． 多辛 Hamilton系统的一般形式为 MO￡Z+KOxz=VS(z) (X，t)∈R。， ∈R ， (1．1) 收稿 日期： 2008年 4月 21日． 国家 自然基金 (No资助项 目． 1．上海大学理学院，上海 200444；SciencesCollege，ShanghaiUniversity，Shanghai200444，China 82 应用数学与计算数学学报 23卷 其 中吖， 为 R (n≥3)上的反对称矩阵，S：R __}R是一个光滑函数 ．如果 (1．1) 满足多辛守恒律 at+‘-g-xi (1～_2) 其中 =~dzAMdz，=~dzAKdz为予辛格式，那么称 (1．1)为多辛Hamilton系统， 而能保持多辛Hamilton系统离散多辛守恒律的格式，称之为多辛格式． 到 目前为止，许多学者对多辛Hamilton系统 (1．1)构造了一系列的多辛格式，如 多辛Preissman格式，多辛有限体积法，多辛有限元方法等 [4】．当问题具有周期边界条 件时，可以用Fourier谱方法来构造多辛格式．基于微分矩阵理论，Bridges等人率先