求解高维扩散方程的一种恒稳定的半显式差分格式.pdfVIP

第24卷第 12期 商 丘 师 范 学 院学 报 Vol_24 N0．12 2008年 12月 J0URNAL0FSHANGQIUTEACHERSC0LLEGE Dec． 2008 求解高维扩散方程的一种恒稳定的半显式差分格式 田巧娴，杨国锋，葛永斌 (宁夏大学 应用数学与力学研究所，宁夏 银川 750021) 摘 要：提出了数值求解二维和三维扩散方程的一种半显式差分格式，其截断误差为D(。+ )，并且是无条 件稳定的．数值算例验证 了方法的精确性和可靠性． 关键词：高维扩散方程；差分方法；半显式；无条件稳定 中图分类号：0241．82 文献标识码 ：A 文章编号：1672—3600(2008)12—0015—05 Anunconditj0nallystablesemi—explicitdiff_erenceschemef0rthehigh—dimensionaldimlsionequation TIANQia0一xian，YANGGuo—feng，GEY0ng-bin (InstituteofAppliedMathematicsandMechanics，NingxiaUniversity，Yinchuan750021，China) Abstract：A semi—explicitdifflerencescheme isprop0sedf0rnumerically solvingtwo— andthree—dimensional diffusi0nequation．Themethodisunc0nditi仰allystableandthetmncatione丌0rsareD( + )．Thenumerical experimentsprovetheaccuI．acyandreliabilityofpresentmeth0d． Keyw0rds：highdimensionaldiflusi0nequation；din-erencemeth0d；semi—exp1icitscheme；unc0nditi0nallystable O 引 言 对扩散方程数值解法的研究一直是科研工作者关注的热点问题之一，特别是对高维情形．其数值方法一 般可分为显式方法和隐式方法两种．隐式方法一般稳定性好、求解精度高，但不适合并行计算或计算量较大． 为此，葛永斌等人提出了基于高精度隐格式的多重网格算法 ¨0j，加快了传统迭代法的收敛速度，提高了计 算效率．与此同时，显式方法计算简单、且利于并行计算，但一般稳定性较差．为此，很多学者进行改进，提出 了具有较高稳定性的各种显式差分格式 【卜 ，较古典显格式，稳定性范围有了较大的改善，但一般仍是条件 稳定的．对所提出的格式 ，其网格比A=J『／ 均不能超过 l／2．最近冯青华和王文洽在Evans和 Abudul— lah 所提出的组显差分格式的基础上，提出了二维扩散方程 的一类交替分组方法_】。。，格式是绝对稳定的， 并且适合并行计算． 本文从经典的crank—Nicolson(C—N)格式出发，对其进行改造，得到一种求解二维和三维扩散方程的 半显式差分格式．格式可以显式计算，具有并行性，并且是绝对稳定的．文中给出数值算例的计算结果，验证 了本文的结论． 1 半显式差分格式 考虑二维扩散方程的初边值问题 收稿 日期 ：2008—05—14 基金项 目：国家 自然科学基金资助项 目(1O502026 作者简介：田巧娴 (1985一)，女，宁夏盐池人 ，宁夏大学硕士研究生 ，主要从事偏微分方程数值解法和计算流体力学的研究 16 商丘师范学 院学报 2008年 ：。(( ++ ))++／(，)，，，¨)