摄动有限体积法重构近似高精度意义.pdfVIP

第 21 卷 第 2 期 计 　　算 　　物 　　理 Vol . 21 ,No . 2 2004 年 3 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Mar . , 2004 　 [文章编号] 　1001246X 2004 0座机电话号码 摄动有限体积法重构近似高精度的意义 高 　智 , 　向　华 , 　申义庆 中国科学院力学研究所 , 北京 　100080 [摘 　要] 　研讨有限体积 FV 方法重构近似高精度的作用问题. FV 方法中积分近似采用中点规则为二阶精 [ 1 ,2 ] 度时 , 重构近似高精度 精度高于二阶 的意义和作用是一个有争议的问题. 利用数值摄动技术 构造了标量 输运方程的积分近似为二阶精度 、重构近似为任意阶精度的迎风型和中心型摄动有限体积 PFV 格式. 迎风 PFV 格式无条件满足对流有界准则 CBC , 中心型 PFV 格式为正型格式 , 两者均不会产生数值振荡解. 利用 PFV 格式求解模型方程的数值结果表明: 与一阶迎风和二阶中心格式相比 , PFV 格式精度高 、对解的间断分辨率 高 、稳定性好 、雷诺数的适用范围大 , 数值地 “证实”重构近似高精度和 PFV 格式的实际意义和好处. [ 关键词] 　计算流体力学 ; 有限体积方法 ; 摄动有限体积方法 [ 中图分类号] 　O35 [文献标识码] 　A 0 　引言 在计算流体力学 CFD 、数值传热学 NHT 中 ,有限体积 FV 方法是最主要的实用方法 ,它在 CFDNHT 软件中得到了普遍的应用. 但一般认为它的精度不如有限差分方法和有限元方法的精度高. 这是因为 FV 方 法是二级近似 :积分近似和插值 或称重构 近似. 积分近似常取中点规则 ,具有二阶精度 ,精度高于二阶的积 分近似需要过多的基点 ,因此很难使用. 精度高于二阶的重构近似需要同时使用相邻和较远的控制单元 ,对 非结构网格 特别是三维情况 重构近似为高精度的格式既复杂又难构造. 何况当积分近似为二阶精度时 ,重 构近似的精度再高 、FV 格式的整体精度仍为二阶. 因此重构近似高精度的 FV 格式并不多见. 摄动有限体积 PFV 方法[ 1 ,2 ] 的提出 ,使我们有可能构造节点数与一阶迎风和二阶中心有限体积格式一致 , 即积分近似精 度为二阶但重构近似为任意阶精度的 PFV 格式. 由于 PFV 格式的整体精度理论上仍为二阶 ,此外还应当强 调 ,在 FV 方法中积分方程对任一小体积单元列出 , 因此使用更多基点明显偏离积分方程的数学规定 ,物理