解抛物型方程的一类高精度显式差分格式.pdfVIP

第29卷第6期 泰 山学 院学 报 Vo1．29 NO．6 2007年 11月 JOURNAL OF TAISHAN UNIVERSITY Nov． 2Oo7 解 抛 物 型 方 程 的一 类 高 精 度 显 式 差 分 格 式 杨 萍，周生田，郭希秀 (中国石油大学(华东)数学与计算科学学院，山东 东营 257061) [摘 要】 利用待定系数法对一雏抛物型方程构造了一类高精度的三层七点显式差分格式，格式的截断 误差达到 O(r +h‘)，稳定性条件是0r~_4／5．当r取特定值0．1335或0．5118时，格式的截断误差可提高到 O(1-4+h。)． [关键词】 一维抛物型方程；待定系数法；截断误差；稳定性条件；显式差分格式 [中图分类号】 O175．2 [文献标识码】 A [文章编号】 1672—259o(2o07)o6—0012—04 0 引言 在研究粘性流动，扩散方程(热传导方程或气体扩散现象)等耗散性传播问题时，常遇到抛物型方 程．在一维情形，其模型为如下初边值问题： ： ≥ 。 ( ，0)： ( )，0 (1) (0，f)： (L，f)：0，t 0 对于模型(1)的差分格式在文献中n 。 已经列出了很多，对于三层七点格式，金承日It]给出了一个 精度为D(r +h )的显式差分格式，其稳定性条件为0r1／2．张天德 对其进行了改进，得到了一个 精度为0(r +h )的显式格式，但其稳定性条件只有．马明书等构造了精度为的绝对稳定差分格式，它 是一个三层九点隐格式，计算量较大，为此本文用待定系数法在点数确定的前提下，确定出了最高阶的 差分格式，又在格式含有一个参数的情况下，确定出了稳定性条件最好的参数r=0．1335或0．5118．该 格式具有精度高，点数少，计算量小的优点，可用追赶法进行求解． 1 差分格式的构造 设时间步长为r，空间步长为h：L／M(M为正整数)，对区域[0， ]X[0， ]作矩形剖分，取局部节 点集为( )(如下图所示)． {( 一-，t +-)，( ，t +-)，( +。，t +。)，( ．。，t ( ，t )，( +。，t )，( ，t 一。)，其中： ：jh，t ：盯，并 令[ ]，= ( ，f )， ，是[ 的近似值． * * * * * * * (口) [收稿扫期]2007—10—25 [作者简介]杨 萍(1982一)，女，山东新泰人，中国石油大学数学与计算科学学院在读硕士 第6期 杨 萍等：解抛物型方程的一类高精度显式差分格式 当问题(1)的解充分光滑时，有如下关系式成立： ： n ，m’ 刀 F贝跫鳅‘为非负整数． (2) ，m ， △ ： ： 鲁 I nj+可O．F2 h4 ·，+ ·， r ) △ ，： ： · 一鲁警· + · 一 · ： 华 ： ·，+ · +鲁(·+r) · +鲁(· ) ·，+署(· + 4，) I +署(1+l0H 20r2)0出7~u ah6(1+l5r+60r+30r3) o(r “。) △ ，： 华 ： 鲁(·+r) ·，一鲁(· ) ·，+署(· + 4， )警I，一鲁(1+lOr+2