四阶抛物方程一类变步长本性并行格式的稳定性分析.pdfVIP

维普资讯 应 用 数 学 MATHEMATICA APPLICATA 2008，2l(3)：485～489 四阶抛物方程一类变步长本性 并行格式的稳定性分析 郭阁阳 ，刘播 (1．天津工程师范学院数学系，天津 300222；2．吉林大学数学学院，吉林 长春 130012) 摘要 ：本文给 出了四阶抛物方程一类变步长 的本性并行差分格式 ，并用发展 的能量法 给出了该格式的绝对稳定性证明 ，得到 了绝对稳定性估计． 关键词 ：四阶抛物方程；变步长；本性并行 ；绝对稳定 中图分类号：0246 AMS(2000)主题分类 ：65M06；65M12；65Y05 文献标识码：A 文章编号 ：1001—9847(2008)03—0485—05 1．引言 近年来 ，抛物方程的并行算法已经得到广泛的发展 【1。。．对于扩散方程 ，Evans和 Abdul— lab首先提出交替分组方法 (AGE)[¨．随后 ，张宝琳设计 了交替分段显一隐式法 (ASE—I)r引．袁 光伟等在文献[3~63中把 AGE方法应用到三阶色散方程 ，并用能量法证明了无条件稳定性 ， 得到稳定性估计．王文恰等又把 AGE方法应用到色散方程和变系数对流一扩散方程 ]，并得 到了令人满意的结果．此外，文献[10]给出并行求解 四阶抛物方程的两种并行算法：AGE和 ASE—I算法 ，该算法具有较高的收敛阶． 如果使用某种格式 ，可以编制能直接在并行机上 自然实施的并行计算程序 ，我们就称这样 的格式为本性并行格式．本文给出了四阶抛物方程一个变步长本性并行的差分格式 ，该格式包 含了多个变步长并行算法．最后用发展 的能量法给出了该变步长并行差分格式的绝对稳定性 证明． 考虑如下四阶抛物方程 挈十 ：0，∈[o，z]，￡∈[o，明， (1) ‘ 口工： 初值和边值条件为 “(，O)一 “。()， ∈[O，z]， (2) “(0，￡)一 “(O，￡)一 “(2，￡)一 “(z，￡)一0，t∈[O，T]， (3) 其中 ()是给定的函数． 2．变步长并行交替差分格式的构造 * 收稿 日期 ：2007-09-20 基金项 目：天津市高校科技发展基金资助项 目 作者简介 ：郭阁阳，男 ，汉 ．吉林人 ，讲师 ．研究方 向：并行数值算法． 维普资讯 应 用 数 学 首先做网格削分．用一组平行线z— = (=0，1，2，…，J)，t— t”一t-+r坩，n一0， 1，2，…，N，将区域[0，z]×[0，T]划分成矩形网格，其中h—l／J和r分别是空间和时间网格 步长，-，和N是 自然数．数值求解周期初边值问题 (1)～ (3)，只需计算网格点(z，t”)(= 1， 2，…，-，一 1；一 1，2，…，N)上的近似解 “7． 引进算予 = I__ “；=一 ±～ ． 并记 = ／h， 一 ￡A；“ +(1_一￡) “，则一般的交替差分格式可表示为 + △i(抛 抖 )一 0，