二维抛物型方程的高精度分支稳定隐格式-安徽师范大学学报.PDFVIP

38 3 ( ) Vol． 38 No． 3 第 卷 期 安徽师范大学学报 自然科学版 2 0 15 5 Journal of Anhui Normal University (Natural Science) May ． 2 0 1 5 年 月 DOI :10 ． 14 182 / J． cnki． 100 1 － 2443． 20 15 ． 03． 004 二维抛物型方程的高精度分支稳定隐格式 1 2 3 ， ， 詹涌强 杨小辉 谭志明 (1． ， 510800 ;2 ． ， 510440 ;3 ． 华南理工大学广州学院计算机工程学院 广东广州 广东警官学院计算机系 广东广州 广东理工职业学院数学 ， 528458) 教研室 广东中山 : ． 摘 要 用待定系数法构造了求解二维抛物型方程的高精度分支稳定隐式差分格式 格式的截断 2 4 1 O (t + x )． r ， ． ， 误差达到 Δ Δ 证明了当 ≥ 时 差分格式是稳定的 通过数值试验 比较了差分格式的 6 ， ． 解和精确解的区别 说明了差分格式的有效性 : ; ; 关键词 二维抛物型方程 隐式差分格式 截断误差 中图分类号:O24 1． 82 文献标志码:A 文章编号:100 1 － 2443 (20 15)03 － 0233 － 06 引 言 、 、 ， ， 在渗流 扩散 热传导等领域中经常会遇到求解二维抛物型方程的初边值问题 在二维情形 其模型问 题为 2 2 u =  u +  u ，0 ＜ x ，y ＜ 1，t ＞ 0 2 2 t x y  u (x ，y ，0)= φ (x ，y )，0 ≤ x ，y ≤ 1 (1) u (0 ，y ，t)= f 1 (y ，t)，u (1，y ，t)= f 2 (y ，t)，0 ≤ y ≤ 1，0 ≤ t ≤ T  