四次b样条galerkin有限元方法数值求解kuramoto-兰州交通大学.pdfVIP

第29 卷　第 6 期 兰 州 交 通 大 学 学 报 Vol.29 No.6 2010 年12 月 Journal of Lanzhou Jiaotong University Dec.2010 　　文章编号:1001-4373(2010)06-0172-06 四次 样条Galerkin 有限元方法数值求解 ＊ Kuramoto-Sivashinsky 方程 张 锋, 　孙建安, 　陶　娜, 　陈继宇, 　石玉仁 (西北师范大学 物理与电子工程学院, 甘肃 兰州　730070) 摘　要:采用四次 样条Galerkin 有限元近似和对时间离散的方法, 构造了求解(1+1)维 Kuramoto-Sivashinsky 方程的数值格式.通过4 个算例和相关文献的对比, 表明了该格式精度高, 适应性强. 关键词:Kuramoto-Sivashinsky 方程;四次 样条;Galerkin 方法;有限元方法;数值解. 中图分类号:O242.21　　　　　文献标志码:A 　　Kuramoto-Sivashinsky(简称K-S)方程是一类 方法用到工程分析中取到很好的效果[15-16] .但目前 重要的非线性演化方程.1975 年, Kuramoto 和 利用 样条Galerkin 有限元方法求解 K-S 方程的 Tsuzuki 在研究空间三维 elousov-Zhabotinski 型 相关文献较少.本文利用四次 样条有限元方法求 [ 1] 反应扩散系统的角-相位湍流时导出了此方程 . 解K-S 方程的数值解, 并与文献[6, 10]作了比较. 1977 年, Sivashinsky 在研究空间二维层焰面的微 [2] 1　Kuramoto-Sivashinsky(K-S)方程 热扩散不稳定性时也导出了此方程 .K-S 方程在 离子物理问题, 热传导及氧化反应的扩散、动力学、 K-S 方程的一般形式如下: 自由膜表面的流动等问题中有着广泛的应用, 很多 ut +uu x +αu xx +δu xx x +βu xxx x =0, x ∈ [a, b], 学者对该方程的精确解和数值解做了大量研究: t 0, (1) Hooper 和 Grimshaw 使用行波变换法得到了该方 其中:α、β、δ是常数.本文考虑该方程初边值问题的 [3] 程的行波解 ;Manickam 和 Moudgalya 用正交立 求解, 初值和边值条件为 方样条配置法与二阶分裂法的组合得到了该方程的 u(x , 0)=f (x ), x ∈ [a, b], (2) [4] 数值解 ;Fan 运用双曲正切函数法得到了该方程 u(a, t )=g1 (t ), u (