二维波动方程的高精度隐格式及其多重网格算法-厦门大学学报自然.PDF

第 42 卷 　第 6 期 ( ) Vol . 42 　No . 6 厦门大学学报 自然科学版 2003 年 11 月 Journal of Xiamen University (Natural Science) Nov. 2003 　 ( ) 文章编号 2003 二维波动方程的高精度隐格式及其 多 重 网 格 算 法 1 1 2 葛永斌 , 吴文权 , 田振夫 ( 1. 上海理工大学动力工程学院 ,上海 200093 ; 2 . 上海大学 ,上海市应用数学和力学研究所 ,上海 200072) 摘要 :提出了数值求解二维波动方程的两种高精度三层紧致隐格式 ,利用 Fourier 分析方法证明了格式均是无条件 稳定的. 并在此基础上提出了求解该问题的多重网格算法 ,从而克服了传统迭代法在求解隐格式时收敛速度慢的缺 陷 ,大大加快了迭代收敛速度 ,提高了求解效率. 数值实验结果验证了方法的精确性和可靠性. 关键词 :二维波动方程 ;高精度 ;紧致隐格式 ;多重网格法 中图分类号 :O 241. 82 文献标识码 :A 　　　　　　　　　　　　　　　 　　波动方程是一种随时间发展的双 曲型微分方 其中 D = { ( x , y) ∶0 x , y 1} , 程 ,对其进行数值求解有非常重要的理论和实际意 初始条件给定为 : ( ) φ( ) ( ) ( ) 义. 和其它的高维非定常问题的求解一样 ,显格式由 u x , y , 0 = x , y , 　　 x , y ∈D 2 ( ) 于受到非常苛刻的稳定性条件限制 ,计算中很少采 u x , y , 0 ψ( ) ( ) ( ) t = x , y , 　　x , y ∈D 3 用. 隐式格式尽管稳定好 ,但采用传统迭代法也有计 算量大 、收敛速度慢的缺陷. 而多重网格方法是从 边界条件给定为 : 70 年代后才开始迅速发展起来的一种行之有效的