三维泊松方程的高精度多重网格解法.pdfVIP

应 　用 　数 　学 MA T H EMA TICA A PPL ICA TA 2006 ,19 (2) :3 13～3 18 三维泊松方程的高精度多重网格解法 1 1 2 葛永斌 , 田振夫 ,马红磊 ( 1. 宁夏大学应用数学与力学研究所 ,宁夏 银川 75002 1 ;2 . 航天医学工程研究所 ,北京 100094) 摘要 :利用对称网格点泰勒展开式中各阶导数项明显的对称性 ,得到了数值求解三维 泊松方程的四阶和六阶精度的紧致差分格式 ,其推导过程简便直接. 为了克服传统迭 代法在求解高维问题时计算量大 、收敛速度慢的缺陷 ,采用了多重网格加速技术 ,设 计了相应的多重网格算法 ,求解了三维泊松方程的 Dirichlet 边值问题. 数值实验结 果表明 ,本文所提出的高精度紧致格式达到了期望的精度并且多重网格方法的加速 效果是非常显著的. 关键词 :泊松方程 ;有限差分法 ;高阶紧致格式 ;多重网格 中图分类号 :O24 1. 82 　　AMS( 2000) 主题分类 :65M06 文献标识码 :A 　　文章编号 :100 19847 (2006) 0203 1306 对于二维泊松方程 , 已经存在不少高精度紧致格式[ 1 ,2 ] 及其多重网格算法[3 ,4 ] . 但是对于 三维泊松方程的高精度紧致差分方法 , 由于格式的构造比二维具有更高的代数复杂性要求 ,所 以研究报道并不多. 具代表性的是 Spotz 和 Carey[ 5 ] 基于中心差分余项修正的方法和 Gup t a 和 Kouat chou [ 6 ] 采用 Mat hmatica 数学软件包所得到的 7 点、15 点、19 点、2 1 点和 27 点格式. 本 文首先从泰勒展式入手 ,通过引进对称和的方法 ,消去中间的导数项 ,从而较容易的得到了数 ( ) 值求解方程 1 的两种精度分别为四阶和六阶的紧致差分格式. 然后 ,利用了多重网格加速技 术 ,采用多重网格的 V 循环算法对三维泊松方程的Dirichlet 边值问题进行了数值求解和方法 验证. 本文考虑如下三维泊松方程的Dirichlet 边值问题 : Δ ( ) ( ) ( ) u = f x , y , z , x , y , z ∈D , 1 　 ( ) ( ) ( ) ( ) u x , y , z = g x , y , z , x , y , z ∈D , 2 2 2 2 其中Δ= 2 + 2 + 2 为L ap lace 算子 , D 为求解区域 , D 为区域边界. 式中 u ( x , y , z) 是 x y z 待求的函数 , f ( x , y , z) 和 g ( x , y , z) 可以视为已知函数. 为了保证文中所得格式的完整性 ,不 失一般性 ,假定函数 u ( x , y , z) 和 f ( x , y , z) 均为充分光滑的函数. 收稿 日期 :2005053 1 ( ) 基金项 目: 国家 自然科学基金资助项 目 ;教育部“高等学校优秀青年教师教学科研奖励计 划”资助项 目;宁夏高等学校科学技术研究项 目 作者简介 :葛永斌 ,男 ,汉 ,宁夏人 ,副教授 ,博士 ,主要从事偏微分方程数值解和计算流体力学的研究. 3 14