四次B样条Galerkin有限元法求解KdV方程.pdfVIP

维普资讯 西 北 师 范 大 学 学 报 (自然科学版) 第 44卷 2008年第 4期 66 JournalofNorthwestNormalUniversity(NaturalScience) VoL44 2008 NO．4 四次B样条 Galerkin有限元法求解 KdV方程 刘兴霞，孙建安 ，豆福全 ，刘万海 (西北师范大学 物理与 电子工程学 院，甘肃 兰州 730070) 摘 要 ：采用有 限元方法进行空间离散 ，提 出了解一维非线性 KdV方程 的四次 B样簪 Galerkin方法．通过两个数值 算例来体现这种算法的精确度 ，对该方法得到 的数值解与精确解 以及二次 B样条 Galerkin有限元解进行 比较 ，结果表 明所求得 的数值解与精确解符合得很好． 关键词：KdV方程；四次B样条Galerkin有限元方法；数值解 中图分类号 ：O175．29 文献标识码 ：A 文章编号：1001—988X(2008)04～0066—05 GalerkinquarticB—splinefiniteelementmethod forSOlvingKdV equation LIU Xing—xia，SUN Jian—an，DOU Fu—quan，LIU Wan-hal (CollegeofPhysicsandElectronicEngineering，NorthwestNormalUniversity，Lanzhou730070，Gansu，China) Abstract：Numericalsolutionsfortheone—dimensionalnonlinearKorteweg—deVriesequationisobtainedby usingthemethodofdiscretizationinspaceandfiniteelementapproximation Todemonstratetheefficiency ofthe method tWO testproblemsare considered． The numericalsolutionsoftheKdV equation are comparedwithboththeexactsolutionsandGalerkinquadraticB—splinefiniteelementsolutions．The numericalsolutionsarefoundtobeingoodagreementwith theexactsolutions． Keywords：Korteweg～deVriesequation；GalerkinquarticB—splinefiniteelementmethod；numerical sO】Ihi．tjOn KdV 方程在浅水波、转动大气 中的 Rossby 1 方法简述 波、等离子体 中的离子声波 、低温非线性晶格声子 波包的热激发和非谐晶格的振动等问题中有着广泛 考虑如下 KdV方程的初边值问题 的应用．文献[1]、[2]研究了 KdV方程 的解析 t … 一 0 解．文献[3]用热平衡积分法对 KdV方程进行 了 (＆≤ z≤ b，t 0)， (1) 求解．文献[4]～[6]构造 了解 KdV方程 的有 限 U(z，0)一 U。(z) (＆≤ 32≤ 6)， (2) 差分法．文献[7]～[9]建立 了解 KdV方程 的有 U (＆，)一 U (6，)一 0 ( 0) (3) 限元法 ，分别采用了立方 B样条形 函数 ，二次 B样