湍流模型发展综述.docVIP

湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；Navier-Stokes方程组；J-K模型 Abstract：On the basis of introducing the problems of turbulence, this paper briefly analyzed four kinds of turbulence models and compared their ability of simulation in different situations. At last, the paper expounded the development direction of the turbulence model. Key words：Turbulence model; Navier-Stokes equations; k?ε，k?ω等，它们都需要求解两个输运方程。 二方程模型即 k ?ε模型以其结构简单，使用方便经济的优点，在工程上得到了极为广泛的应用。但是标准的 k ?ε模型也存在着一些缺陷。如该模型中的许多模拟项建立在各项同性的假设之上，使得该模型对于一些各项异性较强的流动（如分离流动）不能很好的进行模拟以及在粘性底层中存在数值困难等。Wicox提出k?ω模型克服了k?ε模型的许多弱点，特别是在近壁面区域中，k?ω模型的性能比k?ε模型有较大的提高；然而，k?ω对自由来流的ω值非常敏感。Menter通过混合函数将这两种模型结合起来，既保留了k?ω模型在近壁面区域、 k?ε模型在自由剪切层中各自的优良特性，又克服了k?ω模型对自由来流的ω 值的敏感性，从而设计出一种新的BSL模型。 三、湍流模型模拟能力的对比 B-L模型广泛地应用于CFD领域，能够较好地模拟附体流动，对于较小的局部分离流动，该模型也有一定的模拟能力。当采用B-L模型计算跨声速分离流动时，其激波总是较实际激波位置靠后。 J-K模型需要求解一个关于最大粘性剪切应力分布的常微分方程。由Johnson D A和King L S于1985年首先提出，经过不断地完善和发展,陆续发表了几个版本。J-K85模型能较好地计算二维翼型分离流动,但是该模型在模拟附体平衡流动时激波位置过于靠前，确定边界层厚度非常困难。将J-K模型应用于机翼分离流动，提出了J-K90A模型，其中湍流耗散项仅用于恢复区域。同年,为改进J-K85模型模拟附体平衡流动激波位置靠前的不足，Johnson D A和Coakley T J提出了J-K90J模型，该模型考虑到压缩性影响，采用双曲正切函数形式混合的内层湍流粘性系数(由平衡流和J-K85模型内层粘性系数构成)，外层湍流粘性系数与J-K90A模型相同，只是合成流场湍流粘性系数的函数由指数函数修改为双曲正切函数。 肖志祥采用B-L模型,J-K90A及J-K92模型分别数值模拟了跨声速三维机翼、超声速细长旋成体流场。结果表明B-L、J-K90A和J-K92模型数值模拟跨声速和超声速流动，对于附体流动，三种湍流模型都能得到与实验吻合良好的结果；对中等分离流动J-K90A和J-K92模型模拟的结果与实验更接近，对大攻角、大分离流动J-K92模型则体现出较其他两个模型更好的模拟能力来。 刘学炎采用一方程的S-A模型、二方程的标准k-ε模型、RNG k-ε模型、Realizable k-ε模型和二方程的标准k-ω模型，五种不同的湍流模式对二维平板绕流进行数值模拟，并与经验公式相比较。结果显示，采用五种不同的湍流模式，数值模拟平板的粘性流场,由于平板两端形成的奇点的影响,平板前后端压力突变引起了局部摩擦阻力发生变化,但是这部分摩擦阻力很小对总的摩擦阻力影响不大。而S-A模型的计算结果和经验公式结果符合较好。 李超等借助CFD软件，采用k-ε模型、RNG k-ε模型和k-ω模型三种湍流模型对导管螺旋桨进行不同空化数条件下的空化性能计算，通过与试验数据对比，分析不同湍流模型、不同空化数对计算结果的影响。分析结果表明：k-ω模型计算更加稳定，更加容易收敛，精度也相对较高，因而与k-ε模型和RNG k-ε模型相比更适于对导管螺旋桨的空化性能进行计算；在空化数或进速系数较小的条件下，空化相对较剧烈，计算的相对误差较大，需要采用更加精确的空化模型或通过进一步提高网格质量来提高计算精度；另外，与k-ε模型和RNG k-ε模型相比，k-ω模型对空化数较敏感，对不同空化数条件下计算的相对误差变化较大。 四、总结 目前,湍流模型理论虽在理论严谨性方面不断提高,也能处理一些复杂的流动问题,但由于湍流本身的复杂性,现阶段还没有找到对任何问题数值模拟都很