层流动能转捩模型改进方法研究 陈灿平.docx

中国工程热物理学会热机气动热力学学术会议论文编号：122097层流动能转捩模型改进方法研究陈灿平1,2，赵庆军1,3，孙小磊1中国科学院工程热物理研究所，北京 100190；2.中国科学院研究生院，北京 100190；3.中国科学院轻型动力重点实验室，北京 100190）（Tel：010Email:chencanping10@mails.gucas.ac.cn）摘要：作为一种完全从流动现象建模且基于局部变量的模型，层流动能转捩模型表现出了转捩过程中流向脉动逐渐发展成三维湍流脉动的过程而具有更加真实的物理意义。目前层流动能转捩模型对转捩点、转捩区长度预测的准确性还有待提高，模型中也没有考虑逆压梯度、表面粗糙度、来流湍流度等因素对转捩的影响。本文利用流线曲率因子和间歇因子对模型中小尺度湍流粘性系数和湍流粘性系数表达式进行改进。借助于ERCOFTAC T3系列平板、Zierke-Deutsch压气机叶栅和RAE2822跨音速翼型的实验数据对改进后的模型进行了验证。结果表明采用的湍流粘性表达式提高了Bypass转捩转捩点和转捩区长度预测的精度，同时引入流线曲率因子可以提高模型模拟分离流转捩的预测能力。关键词计算流动力学；层流动能转捩模型；Bypass转捩；分离流转捩0.前言转捩过程是流动状态从层流到湍流转变的过程，它对流动的摩擦阻力、热交换、流动分离以及边界层的发展都有影响。计算流体力学中，基于RANS方法和湍流模型已经能够很好的计算湍流流动。但是目前还没有成熟的模型能够模拟层流到湍流的转捩过程，因此转捩过程准确模拟仍是计算流体力学一个基础问题。然而在现代涡轮部件中，涡轮叶片表面边界层大部分处于层流或转捩状态。因此对叶片表面转捩流动的准确模拟是准确预测和分析叶片表面压力分布、摩擦阻力、流动损失规律的前提，进一步理解转捩流动的机制和特点对于现代涡轮设计也有很重要的指导意义。基于RANS方程的边界层转捩模型的发展远远落后于湍流模型的发展。目前计算转捩流动的主要模型有：（1）低雷诺数湍流模型；（2）经验公式的间歇因子转捩模型；（3）间歇因子输运模型；（4）层流动能转捩模型。在湍流模型的发展过程中，湍流模型的低雷诺数修正是利用粘性耗散方法来修正模型中湍流生成过多的问题，后来人们发现低雷诺数湍流模型具有一定模拟转捩的能力，但是其并未考虑转捩过程的物理机制，一般认为低雷诺数湍流型模拟转捩流动的能力只是一个巧合[1]。基于经验公式的间歇因子模型是通过经验关系式得到流向不同位置的间歇因子，利用间歇因子乘以湍流模型计算得到的湍流粘性将转捩前流动和湍流流动联系起来，从而模拟转捩的过程。由于利用经验关系式可以考虑影响转捩的很多因素，如逆压梯度、表面粗糙度、来流湍流度等，而且这些依据实验结果得到的关系式在某种程度上也表现了转捩流动的内在机理，是一种很适合工程应用的转捩模型。但是由于多数的关系式是基于边界层的积分厚度和边界层外部的流动信息，这给求解RANS计算方法带来不便，并且对于分离流动边界层外缘不易定义。一般经验公式的积分厚度都是二维概念，因此模型难以用于复杂流动。利用间歇因子建立间歇因子的输运方程来模拟转捩流动就是间歇因子输运模型。相比于经验公式的间歇因子转捩模型，间歇因子的输运方程不仅能模拟间歇因子流向发展过程，而且可以模拟横跨边界层的转捩流动。为了解决间歇因子模型的积分厚度不便获得的问题，Langtry和Menter[2]提出了基于流场局部信息的转捩模型，但是该模型本质上还是基于经验关联式，而不是基于转捩过程的内在物理规律。层流动能转捩模型通过在湍流模型中加入层流动能输运方程模拟转捩前脉动的发展来模拟转捩流动。1997年Mayle和Schulz[3]根据实验观察到转捩前区域存在流向脉动，且这种脉动在结构上和动力学上和湍流脉动有着显著的差别，建立了层流动能输运方程来模拟流向脉动生成、发展并最终成为湍流脉动的过程，第一次建立起了层流动能转捩模型。此后很多研究人员[3, 4]对于转捩前流向脉动增长过程进行了研究并建模，但是其内在规律还没有被完全了解。Volino[5]和Bradshaw[6]认为流向脉动的增长与“splat mechanism”有关，Walters[7, 8]等人成功对“splat mechanism”建模并认为从层流到湍流的转捩过程是雷诺应力输运方程中压力应变项中快速部分从抑制到逐渐激发的过程[9]。综上可以看出，层流动能转捩模型最大优势在于完全从流动现象建模，考虑了转捩过程的内在规律。目前Walters等人已经建立了完全基于局部变量的转捩模型，L.Cutrone[10]等人对比六种不同转捩模型模拟转捩流动的能力，计算结果显示层流动能转捩模型比其他模型具有更好的模拟效果。Salari[11]等人在Walters[8]模