广义Benjamin-Bona-Mahoney方程的多辛算法 A Multi-Symplectic Algorithm for Generalized Benjamin-Bona-Mahoney(BBM)Equation with Stable Long-Time Numerical Behavior.pdfVIP

2008年12月 西北工业大学学报 Dee． 2008 ofNorthwestern V01．26No．6 第26卷第651 Journal PolytechnicalUniversity 胡伟鹏1，邓子辰1’2 f 1．西北工业大学力学与土木建筑学院．陕西西安710072 1 I 2．大连理工大学工业装备结构分析国家重点实验室，辽宁大连116023／ 摘 微分方程组，利用变分原理得到了多种守恒律，构造了一种等价于Preissmann格式的隐式多辛格 式。钟状孤波解的数值模拟结果表明该多辛格式具有较好的长时间数值稳定性。 关键词：教值模拟，广义Benjamin—Bona—Mahoney方程，多辛算法，钟状孤波解 中图分类号：0241．82 文献标识码：A 文章编号：1000—2758(2008)06—0689—04 近年来，哈密尔顿系统逐渐成为学术界研究的 式中，肘，K∈R舣。为反对称矩阵；S：刑一尺是光滑 热点，主要原因是哈密尔顿系统具有一系列的几何 哈密尔顿函数。上述方程组是一个具多辛结构的多 特性，其中最具吸引力的是其具有辛结构形式。为了 辛哈密尔顿方程组，是因为方程组(1)满足多辛守 将这一理论应用于哈密尔顿体系下的偏微分方程求 恒律 解，Bridges提出了多辛以及多辛积分方法等概 ac￡，-I-a蟑=0 (2) 念‘u，多辛是偏微分方程的几何特性，而多辛积分方 式中，叫和志是预辛形式 dzAMdz 法正是保持偏微分方程这种几何特性的有效途径。 ，dzAKdz · ∞2—r一定2—百一 基于这一理论成果，许多保守型的偏微分方程。如 多辛方程组(1)具有能量守恒律(ECL) Schr6dinger方程、KdV方程等等，都存在能够较好 础+屯／。=0 (3) 保持守恒律的多辛格式[1~3]。 针对非线性广义Benjamin—Bona—Mahoney方 程(广义BBM方程)，本文导出了其多辛偏微分方 zTKaz 。 程组及其多种守恒律，随后构造了一种与 2 多辛方程组(1)还具有一个动量守恒律(MCL) Preissmann格式等价的隐式多辛离散格式。文末给 忍^+49=0 (4) 出了相应的数值算例，并对算例的结果进行了必要 的分析。 式中，h：z—rM百一a。z，g：s(z)一Z—rM了a—z。 考虑广义BBM方程 1广义BBM方程组的多辛形式 a雄+站。a．u一7屯U=0 大量描述守恒系统的偏微分方程可以改写成多 出的BBM方程‘‘3的更一般形式，它是比KdV方程 辛方程组形式 更合适的数学物理模型。 M3,z+Ka,z=V：S(z)z∈R。 (1) 引入正则动量五可=乱，W=以“，a“=P和g=