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

2007年8月 安庆师范学院学报(矗然科学版)A岵栅 第”卷第3期 Joum刚ofAn咖g 一类非线性伪抛物型方程的D—R格式 舒阿秀t2 《1．华东师范大学数学系，上海200062：2．安庆师范学院数学与计算科学学院，安徽安庆246133) 摘要：针对一类非线性伪抛物型方程．构造了D-R格式，并给出了截断误差阶估计。该方法通过求解两个 线性代数方程得到原问腰的解．避免了非线性迭代运算，提高了计算效率。 关键词：伪抛物型方程；D-R格式；截断误差阶估计 中囤分类号：0241．82文献标识码：A 文章编号：l007_4260(2007)03—000l—02 O引言 伪抛物型方程是一类含时间和空间的混合偏导数的高阶偏微分方程，描述了传热、传质及反应扩 散等多种物理现象”l。考虑如下伪抛物型方程的初边值问题： f工¨=F(x，f，“)，，，fE』xJ {“(J，f)=ox，fea，×， … ‘’ ，∈， l口(J，0)=，(x) 接性条件，(O)=O；(H3)，关于各变量存在连续的一阶偏导数。由于伪抛物型方程含有时问与空问变量 的混合偏导项，高敛差分格式的建立存在一定的难度。本文将针对初边值问题(1)建立D—d∞一Rach— f矾交替方向隐格式(即D—R格式)，在数值求解过程中避免非线性迭代，减少了计算量。 1 D—R格式 引理1Ⅲ若条件(HI)一(H3)成立．则初边值问题(1)的解存在且唯一。 由引理1．将矩形区域【0，l】×【o，l】进行等距剖分，选取正整数L和，v，并令h=l儿为空间方向上的 0。．o，m∞，l，…，￡；n=0，1，…“处的离散函数妒御缸J，g捌0_)，构造D—R格式如下，分为两步： (2) (I)t(p。疋(△：v二))+疋(p，t口“)一(“+1)A：v：一‰v：=，(‰，，。，v：) 疋(JD，以v?“) =古(巳+∥。1)_巳一÷(∥《1))△．v：=华， (Ⅱ)t(p．正(△，嵋))+疋(pm瓦口“)一0，+1)矗．《一“喀“=，(‘，‘。穆“) 其中： t(p。t(△，v二))．去。。+∥：一《“一屹一v二)-p，一∥“一V葛一v二+v二心 哦(p卅正碟+1)=嘉(巳{(啦!一口+1卜p，。；(。“一V≥。))’色嵋=!：；笠 收稿日期：200“12—30 作者简介．舒阿秀“977一)，女，安徽旌德人，安庆师范学院数学与计算科学学院讲师，华东师范大学在职硕士生。 万方数据 ．2． 安庆师范学院学报(自然科学版)20∞年 初边值条件为：“：=丘， m=o，l，…，上； “：=略=o，H=o，l，…，Ⅳ． 格式第(I)步(2)可化简为：n11嚣1+6110“+cl蜗12q嵋+l+届v三+^《-1+F(靠，‘“，噱) (4) 2(古+扣巳一；， 其中：铲(古+寺p。+；，铲一(嘉+扫。。+；¨。一；)一妻(”1)q 口 = 上脯 ■0： 属=一嘉(巳+；+巳一；)一÷(g_+1)圯，^2古p，一÷， 格式第(Ⅱ)步(3)可化伺为： 1 1 1 1 其中：铲奇p。+÷，6z一斋(巳+；¨。一÷)一言(g*“)，铲方p。一÷，a舢 岛一奔(匕+；+4’一{(‰+1)，扎2c2，