对流项紧致与非紧致中心差分格式的对比研讨.pdfVIP

中国工程热物理学会 传热传质 学术会议论文 编号：083025 对流项紧致与非紧致中心差分格式的 对比研究 1 1 1* 2 3 4 王艺 张红娜 宇波 窦丹 魏进家 李凤臣 ( 1. 中国石油大学石油天然气工程学院，城市油气输配技术北京市重点实验室，北京 102249; 2. 中国石油天然气勘探开发公司，北京 100034; 3. 西安交通大学能源与动力工程学院，陕西 西安 710049; 4. 哈尔滨工业大学能源科学与工程学院，黑龙江 哈尔滨 150001) *通信作者：电话：010 E-mail:yubobox@cup.edu.cn 摘要：对于低雷诺数二维方腔顶盖驱动流动，分别采用均分网格下具有四阶精度和非均分网格下具有 三阶精度的对流项紧致与非紧致格式进行对比研究。研究表明对流项非紧致中心差分格式的稳定性高 于同阶截差精度的紧致中心差分格式。稀疏非均分网格下非紧致格式的准确度明显高于紧致格式，在 密网格下两格式均能得到与基准解吻合很好的解。 关键词：中心差分，紧致格式，对流项，顶盖驱动流 0 前言 紧致格式在数值计算中应用于非定常多尺度流动[1,2]时表现出比普通差分格式(非紧 致格式)更好的特性，主要表现在对数值解中的高波分量的模拟能力、达到相同精度时 所用的网格点较少以及对流动细微结构的模拟。傅德薰等比较了几种常用的迎风格式和 迎风紧致格式的特性，发现在同阶精度的格式中，紧致格式对高波分量的模拟能力比非 紧致格式强[3] 。那么对于非多尺度流动(此时数值解含有的高波分量很少) ，紧致格式是 否仍优于非紧致格式是一个值得研究的问题。另外，相同截差等级的紧致格式和非紧致 格式是否具有相同的准确性以及数值稳定性，关系到数值计算中离散格式的选取，也是 值得研究的问题。由于对流项的差分格式在离散格式中尤为重要，本文以低雷诺数下二 维方腔顶盖驱动层流流动问题为例(属于非多尺度流动) ，分别考察两种相同截差等级的 对流项中心型差分格式：五点中心差分格式和三点紧致中心差分格式，对以上问题进行 初步探索。 1 无量纲控制方程组和数值计算方法 引入无量纲量：U u/ Utop ，V v/ Utop ，X x / H ，Y y / H ，τ t / (H / Utop ) ， P p / ρU 2 和Re ρU H / μ( U 为顶盖速度， 为方腔边长) ，得到二维方腔顶盖驱动 ( top ) top top H *基金项目：教育部科学技术研究重点项目(106034)、高等学校博士学科点专项科研基金(20050425006)、 国家自然科学基金50521604) 。 流动的无量纲控制方程组如下： ∂U +∂V 0 (1) ∂X ∂Y ⎛ 2 2 ⎞ ∂U ∂U ∂U ∂P 1 ∂ U ∂ U (2) +U +V =− + ⎜ 2 + 2 ⎟ ∂τ ∂X ∂Y ∂X Re ⎝∂X ∂Y ⎠ ∂V ∂V ∂V ∂P 1 ⎛∂2 V ∂2 V ⎞ +U +V =−