附着和分离流动中几种湍流模型的性能.pdfVIP

中困。r程热物理学会 学术会议论文 热机气动热力学 编I∽1 062095 附着和分离流动中几种湍流模型的性能 陈钱林智荣袁新 清lf；大学热能工程系，热科学与动力工程教育部重点实验室，北京。100084 Tel：01 0E-mail：chenqian04@mails．tSinghua．edu．crl 摘要 在三种附者和分离流动中数值研究厂SA模型、SST模型、k一缈模型、q一彩模型的预测性 能。数值计算程序采用新型Lu隐式格式和改良型MUSCL TVD格式。计算值．。j经验公式值以及实验 值的比较表明，这几种湍流模型能较准确地顶测、I￡板边界层流动的近壁参数、后向台阶流动的分离． “附参数和透’r-转了iFmint’栅流动的跨声速特征。 关键词 湍流模型 甲板边界层流动后向台阶流动透、r-转了平面叶栅流动模型比较 1 引言 随着叶轮机械计算流体力学的趋势向多级计算和1F定常计算以及更细致地考虑二 次流路流动的计算的转移?，湍流和转捩模型在这些越来越精细的计算中发挥着越米越 重要的作用‘1·31。 湍流和转捩模型的研究降13】使得目前使用不太多的计算网格(相对于人涡模拟和直 接数值模拟的计算网格)就能较精确地模拟旋转影响、曲率影响、分离、再附、激波． 边界层相互作m等较复杂的流动现象。将各种形式的湍流模型应用于大量实际』=程中内 流利外流的研究U4-2H表明，没有一种湍流模型能适用丁：所有实际的流动，但各种湍流模 艰有着不同的适用范同，而且除了计算精确性，在应用中还需要考虑不同的湍流模型有 不同的计算稳定性和计算收敛性。 本文将在三种流动中考察几种得到J“泛应用的湍流模型应用丁本研究组的二维叶 轮机械计算流体力学程序时的预测性能。湍流模型包括Spalart-Allmaras模型【22．23】(以 ．卜．简称SA模型)、Shear-stress transport模型【2425】(以卜．简称SST模型)、k一力模型12t 25l、q—co模型【26’”】，对其中的SA模型还采用了文献【28】中的另一种公式，另外，本文 对其中的SST模型还进行了封仞J常数的改进。三种流动包括平板边界层流动【29】和后向 台阶流动p训，以及透平转子平面叶栅流动【311。程序中的数值算法采用新型LU隐式格式 lj¨J和改良型MUSCLTVD格式p引。 2基本方程和湍流模型 对可压缩流动完整形式的Navier-Stokes方程忽略质苗力和辐射热后进行Favre平均 pl，得剑控制流动的基本方程。限于篇幅，本文省略这些基本方程，详细公式见文献[5】。 基丁．相同原因，本文也省略以卜．湍流模型的方程。 SA模型 SA模拟依据经验、带纲分析、伽利略不变性、对分子粘性的选择性依赖得到涡粘 性系数的输运方程。(模型构建的详细过程和公式见文献[22】)。 本文除了采心SA模型的常用公式外，还采用了文献【28】的中计算S的新公式。 SST模型 SST模型是一种分区模型：在边界层的不同层和白由剪切层中采用不同的二方程模 型；相对丁传统的涡粘性系数的公式，SST模型还对其进行了修正以改进该模型在逆压 力梯度流动中的预测性能(模型构建的详细过程和公式见文献[25】)。 k一∞模犁 302 使SST模型中的E恒等于1，将封闭常数吒，的值还原为0．5，并使涡粘性系数的 公式还原为坼：一k，则得到后’一∞模型。 co q～∞模型 q一∞模型的详细公式见文献【32】。， 3数值算法和实施 本研究组原有的数值计算程序中已包含本文所考察的g一缈模型，该湍流模型由两 个单独的模块来执行，并与核心土模块之间有数据传递，采用的数值算法为新型Lu隐 式格式【32】和改良型MUSCLTVD格式【33】。由下本文所考察的另几种湍流模型(SA模型、 SST模型、k一03模型)和q一彩模型一样都属于．涡粘性模型，其中的输运方程与流动 基本方程具有相同的特征，所以当在数值计算程序中加入新的模块以包含这几种模型 时，可采用与文献[32】相同的数值算法。 4算例描述和网格 平板边界层流动算例与文献【29】的计算给出的流动条件相同，便于与该文献采用显 式代数应力模型计算得到的结果进行比较。三维网格节点数为65×97×5，其在x．Y平 面内也与文献[29】的二维网格相同，其中的流向网格均匀分布，法向网格第二层的节点 距壁面的距离为1x10-6m(Y+约为0．2)。本文研究过程中曾在平板前缘附近区域将流 向网格加密，并将其计算结果与流向均匀的网格的计算结果比较，二者的著别可忽略。 后向台阶流动算例与文献【30】的实验之一给出的流动条件相同。计算域进口取在台 阶上游15倍台阶高度。文献[341计算了文献【35】中的后向台阶流动，指出后向台阶流动 以往的人量计算由丁二计算域出口距离台阶不够远，导