四阶抛物型方程一个两层的高精度隐式格式.docxVIP

2003 年 11 月Journal of Quanzhou Normal University (Natural Science)Nov. 2003四阶抛物型方程一个两层的高精度隐式格式陈世平(泉州师范学院 数学系 ,福建 泉州 362000)摘 要 :对四阶抛物型方程构造一族新的含参数隐式差分格式 . 适当选取参数时 ,可得到 2003 年 11 月 Journal of Quanzhou Normal University (Natural Science) Nov. 2003 四阶抛物型方程一个两层的高精度隐式格式 陈 世 平 (泉州师范学院 数学系 ,福建 泉州 362000) 摘 要 :对四阶抛物型方程构造一族新的含参数隐式差分格式 . 适当选取参数时 ,可得到一个 高精度恒稳格式 ,其截断误差达到 O [ ( △t) 2 + ( △x) 8 ] ,数值例子表明该格式是有效的 . ① 关键词 :四阶抛物型方程 ;绝对稳定 ;隐式差分格式 中图分类号 :O241. 82 文献标识码 :A 文章编号 :1009 - 8224 (2003) 06 - 006 - 03 对四阶抛物型方 程 笔 者 曾 构 造 一 族 含 参 数 隐 式 差 分 格 式[1 ] , 其 精 度 达 到 O [ ( △t) 2 + ( △x) 6 , 本文对同类问题构造一族含参数的两层高精度隐式格式 , 其截断误差达到 O [ ( △t) 2 + ( △x) 8 , 且数值例子表明该格式是有效的 . 1 差分格式的构造 考虑下列四阶抛物型方程周期初值问题 54 u 5 u = 0 , ( x ∈ R , t 0) ; + 5 x4 5 t ( 1) u ( x , 0) = f ( x) , ( x ∈ R) ; u ( x + L , t) = u ( x , t) , ( x ∈ R , t ≥0) . 设时间步长为 △t , 空间步长为 △x , 网域由点集 ( xj , t n) ( j = 0 , 1 , 组成 , 其中 , xj = j △x , t n = n △t , 又设 r = △t /( △x) 4 为网格比 . 对方程 (1) 建立如下含参数差分格式 , J ; n = 0 , 1 , 2 , ) 所 un +1 n n +1 n un +1 n un +1 n un +1 n - uj u j +1 - uj +1 j - 1 - uj - 1 j +2 - uj +2 j - 2 - uj - 2 j s1 + s2 [ + ] + s3 [ + ] + △t △t △t △t △t un +1 un +1 δ4 n +1 δ4 n n n j +3 - uj +3 j - 3 - uj - 3 x u j x uj ] + θ + (1 - θ) ( △x) 4 (2) s4 [ + = 0 , △t △t ( △x) 4 其中 ,δ4 un +1 un +1 n +1 n +1 n +1 n +1 4 n n n n n n 4 j - 1 + u j - 2 ,δx uj uj +2 - = j +2 - 4 u j +1 + 6 u j - = 4 uj +1 + 6 uj - 4 uj - 1 + uj - 2 . x j un j 表示在节点 ( j △x , n △t) 上的值 , 而 sk ( k = 1 , 2 , 3 , 4) 与θ都是待定参数 , 通过适当选 取这些参数 , 可以使差分方程 ( 2) 逼近微分方程 ( 1) 具有尽可能高阶的离散误差. 收稿日期 :2003 - 05 - 12 作者简介 :陈世平 (1965 - ① ) ,男 ,福建惠安人 ,副教授 ,主要从事计算数学研究 . 第 6 期陈世平 :四阶抛物型方程一个两层的高精度隐式格式7 第 6 期 陈世平 :四阶抛物型方程一个两层的高精度隐式格式 7 泉州师范学院学报 (自然科学)2003 年 11 月8由 98542co sθ -·2co s2θ +·2co s2θ =(9956 + 5328co sθ -+151 泉州师范学院学报 (自然科学) 2003 年 11 月 8 由 9854 2679 51 5 1 ·2co sθ - ·2co s2θ + ·2co s2θ = (9956 + 5328co sθ - + 15120 15120 15120 15120 15120 1 204co s2θ + 40co