对流扩散方程的摄动有限体积(PFV)方法和讨论.pdfVIP

维普资讯 第 36 卷 第 1期 力 学 学 报 V01．36．NO．1 2004 年 1 月 ACTA MECHANICA SINICA Jan．，2004 对流扩散方程的摄动有限体积 (PFV)方法及讨论 高 智 柏 威 (中国科学院力学研究所高温气体动力学重点实验室，北京 100080) 摘要 在有限体积 (FV)方法的重构近似中，引入数值摄动处理，即把界面数值通量摄动展开成网格间距的幂 级数，并利用积分方程 自身的性质求出幂级数的系数，同时获得高精度迎风和中心型摄动有限体积 (PFV)格 式．对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型 PFV格式，这些PFV格式 的结构形式及使用基点数与一阶迎风格式完全一致，迎风 PFV格式满足对流有界准则；二阶和四阶中心 PFV 格式对网格Peclet数的任意值均为正型格式，比常用的二阶中心格式优越．用一维标量输运和方腔流动算例 说明PFV格式的优 良性能，并把 PFV方法与性质相近的摄动有限差分 (PFD)方法及相关的高精度方法作 了对 比分析． 关键词 计算流体力学，有限体积方法，数值摄动，摄动有限体积方法，对流扩散方程 引 言 心型 PFV格式，讨论 了PFV格式的性能，并 以数 值算例说明迎风和中iX,PFV格式的良好性能． 有限体积 (FV)方法能够方便地利用多种类型 1标量输运方程的摄动有限体积 (PFV)格式 的网格 (结构及非结构网格)，适于计算区域复杂的 情况，因此在计算流体力学 (CFD)和 CFD软件 中 摄动有限体积 (PFV)格式的导出以一阶迎风格 得到了普遍的应用．成为与有限差分 (FD)方法和有 式为出发点．标量输运方程的积分形式为 限元 (FE)等方法互为补充的有效数值方法． Ot + ~u rids= a (1) FV方法为二级近似法，对第 1级积分近似，通 常采用中点规则、精度为二阶．第 2级重构近似， 其中 为扩散量， ll为控制体 (CV)体积， 兰 则采用差分格式，最简单的是一阶迎风格式和二阶 ， 和 分别是 的CV平均值和CV中心值． 中心格式．一阶迎风FV格式满足对流有界准则 _l』’ 对二维和三维笛卡儿 (Cartesian)网格和二维非结构 不会产生振荡解但数值耗散太大 [1．引．二阶中心格式 (三角形)网格 (见图1和图2)．在 CV和相邻 CV中 在对流效应占优、网格 Peclet大时发生数值振荡． 心连线垂直于交界面的条件下，方程 (1)的半离散 高阶多基点格式需要求解更复杂的方程、对边界的 一 阶迎风格式为 处理更复杂，特别是对三维非结构网格，很难构造 多基点高精度格式．因此，正如文献 [2]评论的，基 Ot：∑M[([-2d．S) ，广 于精度、效率、简单和可靠的综合考虑，二阶中心 格式 比其他格式好．另一方面，根据摄动有限差分