一种新的代数多重网格法及其在CFD中的应用3∫S∫SρuudS=∫S.PDFVIP

交通科学 第 33 卷 　第 1 期 武 汉 理 工 大 学 学 报 ( 与工程版) Vol . 33 　No . 1 2009 年 2 月 J our nal of Wuhan U niver sit y of Technolo gy Feb . 2009 ( Tran sport ation Science Engineering) 一种新的代数多重网格法及其在 CFD 中的应用 明平剑 　张文平 　朱明刚 ( 哈尔滨工程大学动力与能源工程学院　哈尔滨 　15000 1) ( ) 摘要 :提出了一种面向对象的代数多重网格 algebraic multigrid , A M G 算法 , 以每一层网格作为 研究单元. 网格粗化过程中 ,形成各单元 , 同时记录其前后单元 ,形成双向链表. 粗化过程采用 Ruge 和 St üben 算法 ,光滑算子用 Gau ssSeidel 迭代. 由于 A M G 算法与网格信息无关 ,可以作为“即插即 用”型的线性方程组求解器. 对 CFD 计算过程耗时最多的压力修正方程作了研究 ,分别对二维后 台阶流动模型在不同网格划分情况进行了计算 ,代数多重网格方法与单重网格的不完全分解共轭 梯度法对比发现 ,前者具有明显的优势 ,随着网格数 目增加 ,优势表现更为明显. 最后与 A M G1r5 相比 ,开发的程序内存 占用较少 ,最高只有 A M G1r5 的 36 %. 关键词 :代数多重网格 ;计算流体力学 ;面向对象算法 ;后台阶流动 中图法分类号 :O351 DOI :10 . 3963/ j . i ssn .. 2009 . 0 1. 023 0 　引　　言 1 　计算流体力学控制方程 流体力学数值计算耗时最多的就是线性方程 流体流动和能量传递可以用一组偏微分方程 组的求解 , 线性方程组的求解有 J acobi , Gau ss 来描述 ,此方程组根据质量 、动量以及能量守恒定 Seidel , 以及 Krylov 子空间等迭代法. 以上算法都 律推导出来. 属于单层迭代方法 ,对于大规模实际问题 ,单层迭 质量方程 代法的收敛速度已经接近其极限水平 ,上升空间 ρ ( ) u d S = 0 1 [ 1] S 很小 . 收敛速度主要受低频分量的制约 , 为此 ∫ 196 1 年 Fedorenko 首次采用几何多重网格方法 　　动量方程 将细网格上的低频分量转化为粗网格上的高频分 ρ μ( T)