高精度有限差分在壁面湍流直接数值模拟中的应用+研究.pdfVIP

中国．L=程热物理学会 馅热传质分会 学术会议论文 编号：06324,1 高精度有限差分在壁面湍流直接数值模 拟中的应用 马良栋，李增耀，陶文铨 (西安交通大学能源与动力j-程学院动力Ji程多棚流国家霞点实验窜，陕西西安710049) (电话：029mai 1～．itu．edu．crl) 1：wqtao@mai 摘要：本文采用高精度有限差分法(对流项采用五阶迎风差分，扩散项采用六阶巾心差分)对矩形管道 内充分发展湍流换热进行了数值模拟．湍流雷诺数Re，和普朗特数分别为400(Re。=6200)和0．71．计算 结果表明，同传统的：阶巾心差分格式相比，高精度有限差分可以在较少的网格下得到较好的结果． 通过对压力泊松方程采用不同离散格式的实施，其结果表明，雎力波松方程的散度算子和梯度算子采 用四阶中心差分或二阶中心差分对计算结果的影响不大，但压力波松方程采用二阶中心差分时，可以 节省大量的计算时间。 关键词：高稽度有限差分，湍流，换热，克接数位模拟， 1引言 直接数值模拟(DNS)作为研究|丁．段，与试验技术、理论分析以及其它的数值方法相 互补充。因此随着计算技术的提高，湍流商接数值模拟住湍流研究中越来越发挥着重要 的作用【IJ。而高精度差分格式在研究多尺度的流动结构和流动机理上具有重要的作用。 是湍流直接数值模拟研究的重要方法。Rai和MoinpI最早采用高精度的有限差分法对槽 道湍流进行了直接数值模拟，结果表明高精度的有限差分法的计算精度可以与谱方法相 媲美。之后，许多研究者均采用了这一方法11j1。本文在具有两个方向各项异性的壁面湍 流实施了高精度有限差分法，并且对压力波松方程的散度算子和梯度算子分别采用四阶 中心差分和二阶中心差分进行计算，结果表明，压力波松方程的离散方法对计算结果的 影响不大，特别在并行计算中，采用二阶中心差分是经济的，可以大大节省计算时问。 2物理数学模型 本文以矩形通道湍流流动和换热为例，其物理模型见图l所示。管道的高为Ⅳ，长 为￡，工为主流流动方向，Y和z为展向，管道左右肇面绝热，上下壁面等壁温且下表面 足高温．上表面足低温，设矩形通道内充满4i可压缩流体，在均匀的，F均压力梯度驱动 下流动，U、’，和w分别表示流向和展向的速度分量。流体不考虑重力的影响，则其控 制方程足连续性方程，不可压缩的N．S方程以及能最方程，见文献【6】所示。 本课题受丽家自然科学基会资助，共r11并行计算足在卜海超级汁算巾心的曙光4000A超级 计算机卜完成，特此表示感谢。 I 167 计算中蹙面采用无滑移边界条件，流向方 向速度和温度为周期性边界条件。采用分布算 法解决速度与压力的耦合。时lHj项采用显示：二 ， 阶的Adams．Bashforth格式。对流项采用氕阶 迎风偏置格式，扩散项采用六阶中心差分。 其具体的离散方法见文献[1，2]。这里为了比舶扎 较：．j阶巾心差分与高阶迎风偏置格式之l’HJ的特 性，表l给出了算例的物理参数。2-order表示 高温7； j二阶巾心差分格式，high-order表示高阶迎风偏 图1管道几何模型及嫩标 置格式，且压力泊松方程的散度算子和梯度算 子为四阶中心差分，high-order-2p表示高阶迎风偏置格式，但压力泊松方程的散度算子 和梯度算子为：二阶中-Ii,差分。 表l算例的物理参数 3计算结果比较 ‘图2给出了湍流脉动强度的分布，图2(a)表明，对于主流脉动强度甜。，，在16y+120 的范围内，本文的计算均小于大涡模拟的计算结果141，但是两种高阶格式的结果一致且 均好于二阶中心差分格式，密网格下的二阶中，11,格式略好于疏网格的结果。同Huser和 Bigingenl51计算的结果的差异，主要在于雷诺数的不同引起的差异。图2Co)和(c)表明， 对于’，。。和M■。，本文直接模拟