绕翼型可压缩湍流流动的数值模拟.pdfVIP

第十届全国计算流体力学会议论文集 绕翼型可压缩湍流流动的数值模拟 杨明智陈作斌张玉伦 中国空气动力研究与发展中心，四川绵阳(邮政编码621000) 摘要车文运用Baldwin—L哪x代数模型、高雷诺数k-e两方程模型和低雷诺数k-s两方程模型 对NACA0012翼型和RAE2822翼型进行了计算，井对计算结果和实验结果进行了比较，两者吻台 较好。 关键词翼型．湍流，k一￡两方程模型 0引言 对绕翼型的粘性可压流计算，除了网格分布和离散格式外．湍流模型的选取也是获得 满意结果的重要因素。目前，在粘性可压流的计算中，通常采用的是Baldwin-Lomax代数 模型(下简称为B—L模型)。该模型虽然形式简单，使用方便，但是该模型只适于平衡边 界层的流动情况，当边界层分离时，用B-L模型计算得到的摩擦系数分布与实验差别很 大。为解决这一问题，可以考虑采用具有更强模拟能力的k1两方程模型来对绕翼型的粘 性可压缩流动进行数值模拟。 但是，在实际计算过程中发现．由于k-z两方程模型中椿项的存在，使得k一￡湍流方程 的时间步长被严格限制，结果使k-e湍流方程数值解收敛显著地滞后于流场方程，妨碍了 方程的收敛特性。为此，本文对k—e方程中的源项和耗散项进行了隐式处理，从而提高了 时间步长．改善了其收敛特性。并在亚跨超软件平台基础上，运用高雷诺数k一￡两方程模 了计算．和实验结果“1进行了比较，得到了有意义的结果。 1湍流模型 1．1B-L模型 B—L模型一类常用的零方程模型，目前已在可压缩湍流流动的数值模拟中得到了较为广泛 的应用，其具体形式参见文献[2]。 1 2k-c两方程模型 k一￡两方程模型是一类工程上常用的两方程模型．其雷诺应力通过布辛涅斯克 (Boussinesq)假设与流动应变相关联，即： ～68— 一面刊妾+》+知妻一；秘 n，， 按照对固壁边界条件处理方法的不同，这类模型又可进一步分为高雷诺数k-E两方程模型和 低雷诺数k一￡两方程模型，下面分别进行介绍。 (1)高雷诺数k-E两方程模型(ChenPatel[2])。 一般曲线坐标系下的高雷诺数b8两方程模型方程可表达为： 硎蓖j既i i+言2磊+马 (1‘2) 其中： 矿=艄 弘¨1IpvJ，习k V‘鸹雾 ∽，， 旦f堡，。+雾纠 吼L斑h E?2两1； (14) 盟f鱼一 以L盘h 十雾纠 一 蹦 胆 2一j 皿=‘ (1．5) C占一七 丛如丛船 一 跳 呜P 2—3 坠鼠勃 ，一七 ∽s， 尸=[(等+期一j2瓦chk岛1鲁s：考 方程中，k为湍流动能-占为湍流动能的耗散率一∥为分子粘性系数．岸，为湍流粘性 系数。p为平均密度，u，v为平均速度的x，y分量。