求解Kdv方程的一种显式模平方守恒格式.pdfVIP

第29卷 第2期 海 南 大 学 学 报 自然 科 学 版 Vo1．29N02 2011年 6月 NATURAL SCIENCEJOURNALOFHAINANUNIVERSITY Jun．201l 文章编号：1004—1729(2011)02—0120一O8 求解 Kdv方程的一种显式模平方守恒格式 胡勋锋 ，孙建强，李昊辰 (海南大学 信息科学技术学院，海南 海 口，570228) 摘 要：利用Magnus方法求解Kdv方程．Kdv方程具有模平方守恒特性，首先用适当差分格式对其进行模平 方守恒空间离散，转化成模平方守恒的常微分方程，再用 Magnus方法求解．数值结果表明，Magnus方法能保 Kdv方程模平方守恒特性． 关键词：Kdv方程；Magnus方法 ；模平方守恒 中图分类号：0241．82 文献标志码：A 1895年，Korteweg和Devries导出了描述表面远场特点的非线性波动方程，即著名的Kdv方程 +ozuu +62Z／船 = 0，t0，∈[a，b]， (1) 其中，“(0，)=u()，“(￡，a)=u(t，b)，C和 是2个实数．1965年，Zabusky与 Kruskal运用数值方法对 Kdv方程在周期性边界条件下进行数值求解时，发现满足Kdv方程的一系列孤子波相互碰撞后仍保持波 形和速度不变，并定义为孤立子的碰撞 ，从而迎来了研究孤立子的高潮 ，尤其对 Kdv方程 的特性、特征进 行了深人的研究，并取得 了突破性的进展．如Kdv方程有孤立子解 ，有无穷多个守恒律，有化为线性方程 的解析解法，完全可积等 ．同时求解Kdv方程的各种各样的数值方法也被提出，包括有限差分法和 Fourier展开方法 ．这些方法充分考虑了Kdv方程的守恒特性，使构造 的格式尽量保持Kdv方程的某个 守恒特性，然而构造的差分守恒格式常常是隐式且低阶的． 冯康院士与其保结构算法组成员提出的保辛结构的辛几何算法 ，以及在辛几何算法基础上提出的多 辛算法，能长时间精确模拟微分方程的变化，且能近似保方程的能量和动量守恒等特性．这些算法不仅能 使微分方程的离散格式具有高阶精度和稳定性，还能尽量保持微分方程的某些守恒特性，如能量守恒、动 量守恒、辛守恒、模守恒等，所有这些保持微分方程内部结构特性的数值方法称为保结构算法，近几年来， 在保结构算法领域 内产生了保微分方程 内部结构特性的李群算法，是一种指数积分算法，具有保方程的 模守恒特性、稳定性好等优点，通过把一个李群方程转化成相应的李代数方程进行数值求解 ，再代回得到 李群方程的数值解，常见的李群方法有 RKMK方法、Magnus方法和Fer方法等 ’． 本文用李群算法中的Magnus方法求解 Kdv方程．首先介绍保模平方守恒的Magnus方法，再对Kdv 方程进行离散转化成模平方守恒的常微分方程，最后利用高阶显式保结构的Magnus方法和相应阶的显 式Runge．Kutta方法求解离散的Kdv方程，并对数值结果进行分析． 1 Kdv方程的守恒特性 一 般对于一个偏微分方程系统(PDE)，运动常数的存在也等价于守恒律的存在．这里的守恒律是指 对下述形式的方程 +Jx：0， (2) 收稿 日期：2010—10—03 基金项 目：国家 自然科学基金 ；海南省 自然科学基金 (110002)；海南大学科研启动基金 (kyqd1053)；海南大学青年基金(qnjj1022)资助 作者简介：胡勋锋 (1987一)，男，湖北黄冈人 ，海南大学信息科学技术学院应用数学系2009级硕士研究生． 通信作者：孙建强，教授，qiangqiangi6@gmail．eom 第2期 胡勋锋等：求解 Kdv方程的一种显式模平方守恒格式 12l 使得当ll一