热流计算中数值格式的网格敏感性的研究.pdfVIP

第十四届全国计算流体力学会议论文集 热流计算中数值格式的网格敏感性研究 张亮，姜贵庆 (中国航天空气动力技术研究院，北京100074) 摘 要：采用数值求解雷诺平均N_s方程方法对钝锥和钝双锥两种实验模型进行了计算，比较了在 不同壁面网格尺度下不同数值格式在热流计算中的表现，认为在热流计算中，不同数值格式对网格尺度的 敏感程度是不同的．线性波分辨率高的格式受网格尺度影响较小，反之则受网格尺度影响较小。据此，本 文提出在CFD热流计算中应将网格尺度与格式物理特性综合考虑以提高热流计算的精度． 关键词：数值格式；热流；网格敏感性；线性波 0 引 言 高超声速气动热的数值模拟是CFD研究中的难点问题。热流是由流体粘性起主导作用的物理现 象，在计算中受数值格式、网格分布、收敛过程甚至热流的后处理计算等各方面因素的影响，多种 影响因素之间的交错决定了热流计算的复杂性。因此迄今为止，热流cFD计算中还存在一些没有解 决的问题。 在国内外，研究者对影响热流的各种因素进行了讨论。在网格尺度方面，Klopfer和Yeeul指出， 热流计算结果与物面附近网格尺寸有密切关系，并且对于钝头体给出了通过减小网格法向间距改善 热流计算的方法；李君哲比1的研究指出：影响热流CFD计算的是沿物面的法向网格，而与切向网格 分布基本无关，并结合来流参数，给出了网格划分标准。 在计算方法方面，贺国宏嘲讨论了计算格式对表面热流计算的影响，认为壁面边界附近采用二 阶中心型格式更有利于热流的计算。谢锦睿H1分析了TVD格式中熵修正函数的耗散本质，提出了改 进的各向异性熵修正函数以提高热流计算的精度。 多年研究表明，虽然人们在气动热的数值模拟方面有了相当大的进步，但仍旧存在诸多问题， 数值研究多侧重于在相同数值格式下网格尺度的确定，或在相同网格尺度下数值格式的选择，对于 不同格式对网格敏感性方面的研究尚不充分。本文采用数值求解雷诺平均N-S方程方法对钝锥和钝 双锥两种实验模型进行了数值模拟，比较了在不同网格尺度下，不同空间离散格式在热流模拟中的 表现，研究了不同数值格式对于网格尺度的敏感程度，将数值格式的线性波分辨率与网格敏感性建 立关系，从而更好地指导热流计算。 1 数学模型及算法 1．1控制方程 在笛卡儿坐标系下，采用层流三维非定常可压缩雷诺平均N-S方程： 丝+望+塑+塑：盟4。-塑+盟 a 苏 砂 色 叙 砂 昆 ·__。一…’。。。’一‘-·-‘一 其中，豆为守恒变量；户，召，豆为三个方向的无粘矢通量；丘，色，鼠b--+方向的粘 性矢通量，具体形式可参阅文献[5] 1．2数值方法及定解条件 控制方程采用有限体积法进行离散。为使格式达到二阶精度，选用minmod限制器通过MUSCL方 法进行插值，无粘通量选择以下四种经典格式进行离散，即， 345 第十四月々目”算楫i#力学☆Ⅸ*女{ (I)vanLeer格式 (2)Roe的FDS格式 (3)AUSM格式 (4)AUSM up格式 粘性通量采用二阶中心格式离散，时间方向则采用隐式LuSGS方法推进。 边界条件的处理}，物面处采用无滑移、等壁温、零压力梯度边界条件；八u边界赋来流参数 出口边界则采用外插处理。 2结果分析和讨论 对于钝锥和钝烈锥两种实验模型，在相同网格条件下研究不同数值格式的热流计算表现。同时． 对丁每种数值格式，选择系列网格进行格式的网格敏感性研究。 2 1数值格式对热流计算的影响 影响热流计算的因素众多，为晟大程度地排除其他因素对于格式敛应研究的影响．这里选用有 详尽实验数据的简单层流算例””进行计算。在两组算例中，壁面法向第层网格高度均取为1x10■。 2 1 I钝锥绕流 696唧，头部曲率半径Rn=2794mm，半 钝锥对称面及其壁面网格如图I所示，模型长度￡：488