KdV方程的二次B样条有限元孤立子模拟.pdf

第15卷第4捌 济南大学学报(自然科学版) Vd．15．Nu4 2001年J2月 JOURNALOFJINAN Dec．．200I UNIVERSII’Y(Sdence) 文章编号：1671—3559(2001}04—0311—04 KdV方程的二次B样条有限元孤立子模拟 曹爱增1，蔡卫东1，王兵2 (1济南太学信息I删学院，2济南^学土术建筑学院．山东济南250002) 摘要：应用Bubnov Galerkln万法致二次n样杀孝限元方 法，『旱到一种数值计算解KdV方程的方法。对其产生的五 到两个三角矩阵。而后再用追赶法解这个矩阵方 对角矩哇万程用教值线代数的Ib,little三角静解方法求解． 程，然后编写详细过程及解决这个问题的程序。最 并对二1种谘武的线性稳定性进行了分析研宽，结名孤立子 后再刘格式进行稳定性分析。 模型编写了}方算洼的汁算机程序，从而得到了给定韧边值条 件1、姥KdV方程数值模拟站果一 1问题的描述 关键词：KdV方程；I卜样条；有限元方法；孤立于 KdV方程 中图分类号：11})={OI 文献标识码：A L，f十eUU。+一r上，。=fJ，“≤j≤b(J) bldV方稃：U+【j【，，1U。，=【1是由Ko一在这里∥，e是正参数，下标l{，t代表对r，，的微 Vries第一次提出的…、用米描述潜水分，边界条件将从下面的条件中选择： 】_teweg和de 波的无损耗传播。KdV方程演不r非线性扩散波 r，(a，￡)=U(b，，)=0． 的一个主要规律。它的非线性项和扩散项保持平 Uf(日，，)=L，，(b，f)=0 (2) 衡，使稳定的孤立了波能够形成，这些孤波即孤立 再加上初值条件： 于。它的明显特征是在相互碰撞后能保持原来的大 U(a，0)=r(z)(3) 小、形状和速度，因此，在许多考虑简单非线性和简 根据给定的方程和初边值条件，我们可以选择比较 单扩敞影响的研究巾，它是一个有用的逼近。实际 好的方法来进行计算，进而得到较为满意的结果。 上，KdV方程的建立可描述一大批重要的物理现 2数值方法的选择、描述和特点 象，而其中的两个典型例子是浅水波和等离子区的 离子声波。 有上面给出的方程和条件，我们应用Galerkln KdV方程的局部问题的解析解、用反散射方 方法结合权函数v(。)，并且分步积分得 r b 法”』求解KdV方程及KdV方程的相似解【21都已实 o d 现。在此前的一些论文中，也有用三次样条建立了 在被积嫡数中二阶空间导数的出现要求插值函数和 KdV方程的一种三次样条有限元解L3,4j，而且得到 它们的一阶导数在整个区域内必须连续．而二次B r比其他方法有更小误差的解。那种选择是由于样 样条函数正好满足这个要求，因此我们选择二次B 条具有的特殊性质和知识制定的，但他对初等连续