一类非线性伪抛物型方程的D格式.pdfVIP

第25卷第3期 阜阳师范学院学报(自然科学版) VoI．25．No．3 Joumalof Teacher8 Science) Jun．200Z 2007年9月 Fuyang colle鹳(NatⅢal 一类非线性伪抛物型方程的D格式 舒阿秀 (安庆师范学院教学与计算科学学院，安徽安庆246011) 摘要：针对一类非线性伪抛物型方程，构造了D格式，并给出了截断误差阶估计．该方法通过求解两个线性代数方 程得到原问题的解，避免了非线性迭代运算，提高了计算效率． 关键词；伪抛物型方程}D格式，截断误差阶估计 中圈分类号。0241．82 文献标识码。A 文章编号；1004·4329(2007)02—0011一02 值问题(O．1)的解存在且唯一． 0 引言 根据引理1．1，可对问题进行数值求解．将矩形 伪抛物型方程是一类含时间和空间的混合偏导 区域[o，1]×[o，1]进行等距剖分。选取正整数L和 数的高阶偏微分方程，描述了传热、传质及反应扩散 Ⅳ，并令^=1／L为空间方向上的网格步长，奄= 等多种物理现象[¨．考虑如下伪抛物型方程的初边 1／Ⅳ为时间方向上的网格步长，记 值问题t 而=访，f=0，1，…，LI r工缸=F(z，f，“)，z，t∈J×， “2疗五，疗=0，1，…，Ⅳ． ≮∞o，f)=Oz，￡∈ar×， (O．1) 设哌是定义在网格点(％，厶)，m=o，1，．．·，L，露= 【“(z，o)=厂(z)z∈J o，1，．．·，Ⅳ处的离散函数，“=p(‰)，“=g(‰)， 其中L弘一(户o)甜蹦)，+(po)‰)；一(go)+1)地 一口o)“，J一[o，1]，‘，=[o，1]．方程中的系数函 如下，分为两步： 数、右端项及初值满足如下条件； (H1)户(z)∈C1(J)，口(z)∈C(J)且户(z)O， 口o)≥0， (H2)厂o)∈C(j)，且初值满足连接性条件 (1．1) 厂(0)=0f 其中 (H。)F关于各变量存在连续的一阶偏导数． 本文将在文献[2]，[3]的基础上进一步研究问允(户-丸(岔吒))2壶(户一+考(嵋嚣1一嵋计1一 题(o．1)的高效差分格式．由于伪抛物型方程含有 吆+1+吒)一户一吾(畦一+1一嵋冀1一屹+秒：一1))， 时间与空间变量的混合偏导项，高效差分格式的建 立存在一定的难度．本文将针对初边值问题(o．1)∞姒芷≯))-赤(户州嵋#1一∥1+ 建立Douglas交替方向隐格式(即D格式)，提高了 嵋+1一眩)一p一丢(嵋计1一以才1+吒一嵋一1))， 时间方向上的精度，从而在数值求解过程中避免了 譬咤=华， 非线性迭代，减少了计算量． (I)以(夕。以(4屹))+丸(户。以(蝰掣))一 1 D格式 引理1．1乜】 若条件(H。)一(H。)成立，则初边 收稿日期：2006．12—23 作者简介：舒阿秀(1977一)，女，讲师，华东师范大学在职硕士生．研究方向：偏微分方程数值解． 万方数据万方数据 12 阜阳师范学院学报(自然科学版) 第24卷 (如+1)雠．g