非均匀网格泊松方程高阶紧致离散及加速方法-西安交通大学CFD.PDFVIP

第 卷第 期 工 程 热 物 理 学 报 ， 年 月 ， 非均匀 网格泊松方程高阶紧致离散及加速方法 兆 辉 封 永亮 陶 文铨 西 安交通 大学 热流科 学与工程教育部重点 实验 室 ， 西 安 摘 要 本文将非均 勻 网格直接离 散的 高阶紧致格式从二维推广 到三维 结合 附加修 正多重 网格方法提高 了传统迭 代方法 ， ， 且 证 现 。 明： 法 以 的 程 的 压 的收敛 效率 并 验 了该格式在 不同边 界条件 的数值表 结果表 该方 可 有效 求解 方 中 力泊松方程 关键词 高阶紧致格式； 附近 修 正多重 网格方法 ； 压 力泊 松方程 中 图分类 ： 文献标识码： 文章 编 ： 号 号 ， 引 言 不可压缩 方程求解 中常常 需要求解 压力 泊 同 时 网格变换时还 有 可能产生数值误差。 因此， 本 松方程 过去二十年 中 ， 出现 了 大量求解压力 泊松 文将文献 中 的非均 匀 网格上直接构造 的二维 高 方程 的方法 ⑵ 。 特别 是在湍流 的直接模拟 和大涡模 阶 紧致格式推广 到三维 ， 并采用 了基 于 附加修正 的 拟 中 需要求解高阶离散精度 的压 力 泊松方程， 近年 循环 多重 网格算法加速迭代收敛 。 来提 出 的高阶紧致格式 （ 该格 式可 以在任意坐标方 向上的 非均 匀 网 格上 结合 多重 网格方法受到人们关注 。 例如， 采 实施 ， 依赖 于 网 格 的 间距， 计算精度至 少三 阶或 四 用 及 循环 的多重 网格算法求解 了 二维压力 阶 。 并且可 处理第一类边界 、 齐次边界 以 及 周期性 方程， 对 比五点 二阶 中心差分 ， 计算结果表 明， 无 边界 多种边界条件 ， 可应用 于 方程 的求解过程 论在 数值精度还是在计算效率上 都有 了非 常大 的提 中 。 基于 附近修 正的 多重 网格方法， 系数 自然生成 ， 高 很多文献研 究 了均 勻 网 格 下的 格式实施， 不需要特别 的构造 限定及 延拓算子， 因而实施 简单， 然而考虑 精度及 节省 资源 ，