计算流体CFD2004-10.pptVIP

第十章 湍流模型 湍流模型分类 雷诺平均后的方程（RANS） 涡粘性系数模型： 10.2 Spalart-Allmaras 模型 10.3 标准、RNG和带旋流修正κ-ε模型 标准κ-ε模型 RNG κ-ε模型 Realizable κ-ε模型 10.4 标准κ-ω和SST κ-ω模型 10.4.2 剪切压力传输（SST） k-ω模型 10.5 雷诺应力模型 10.6 大涡模拟（LES）模型 10.7 湍流近壁的处理方法 10.7.2 壁面函数 10.7.3 增强壁面处理（双层模型） 增强壁面函数： * 华中科技大学能源与动力工程学院 N-S方程的简化处理（因直接求解困难） 脉动的速度场是湍流的主要特征，这种脉动使得流体介质之间相互交换动量、能量和浓度而引起了这些量的脉动。由于这些脉动规模小、频率高，直接模拟实际的工程对计算机的要求很高。因此，可以采用时间平均法、总体平均法或者用去掉那些小变化量的方法修正方程，以减少计算机的耗时。然而，修正后的方程又包含了未知的附加变量，湍流模型就是用来决定这些变量。 事实上，没有一个湍流模型对于所有的问题是通用的。选择哪种湍流模型，要根据流动中所包含的物理现象、给定的实际情况、精度要求、现有计算条件以及时间限制而定，要想选用最优的计算模型，首先要清楚不同模型的适用范围和限制。 10.1 简介 湍流 直接模拟 大涡模拟(大尺度涡直接模拟、小尺度工程模拟)见10.6 雷诺平均的工程模拟，见10.1(10.2~10.5) 时间平均的雷诺方程与N-S方程相比，除平均量外多出脉动项： 雷诺平均的工程模拟 式中的雷诺应力项 和粘性扩散项 由湍流模式得到 雷诺应力 模型 零方程模型 Spalart-Allmaras一方程模型 等二方程模型(κ-ω) 涡粘性模型 RSM（雷诺应力方程模型） 6个雷诺应力项： 或 3个粘性扩散项： 其中双方程模型考虑了两个独立湍流量（κ-ε或κ-ω）的对流、扩散及其随时间的变化。能较真实地描述许多流动的主要物理过程。既有较好的通用性和精度，又不需太大的计算量。因而是近二十年来研究得最深入、应用得最广泛的模式之一。 求湍流粘性系数 所用的输运微分方程数，根据各个模型水准而不同，按其增加输运方程的个数可分类为零方程、“1”方程( Spalart-Allmaras )、“2”方程模型等等 “1”方程( Spalart-Allmaras )模型考虑了1个独立湍流量--湍流运动粘性系数 的对流、扩散及其随时间的变化。 壁面附近处的流动 湍流受壁面的影响很大，很明显，平均流动区域将由于壁面不光滑而受到影响。当然，湍流还受到壁面其他的一些影响。在离壁面很近的地方，粘性将抑制流体切线方向速度的变化，而且流体运动受壁面阻碍从而抑制了正向的脉动。但近壁面的外部区域，湍流动能受平均流速的影响而增大，湍流运动加剧。κ-ε模型，RSM模型，LES模型都仅适用于湍流核心区域（一般都远离壁面），应该考虑怎样使这些模型适用于壁面边界层处的流动。 在湍流模型中利用Boussinesq近似，中心问题是怎样计算涡粘度 该模型相对简单，可以计算运动漩涡（湍流）粘性传输方程，在这个方程中，不必要计算对应于本地边界层厚度。该模型主要是为涉及壁面边界流动的航空应用设计的，对于受到逆向压力梯度的边界层，采用该模型能得到较好的结果，在涡轮机械中的应用也逐渐受到青睐。 主要的不同点是： 计算湍流粘性的方法； 控制κ和ε湍流耗散的湍流Prandtl数； 在ε方程中湍流的生成和耗散项。 每个模型计算湍流粘性的方法和模型常数不一样。 标准、RNG和带旋流修正κ-ε模型这三种模型有相似的形式，有κ方程和ε方程；有 包括方程中有许多项的模拟是一样的 是由于平均速度梯度产生的湍流动 能项，三个模型中都一样； 其中： ，Mt马赫数,可压缩流体考虑； C1?、C2?、C3?、?k 、??是常数； 湍流动能κ方程是个精确方程，湍流耗散率ε方程是由经验公式导出的方程 Prt为湍流能量普朗特数，在标准模型和Realizable 模型中，缺省为0.85；RNG模型中Prt=1/ ? ，?由下式计算： 由于浮力产生的湍流动能项 （如果重力场和温度梯度同时存在）； 热力膨胀系数； 上式中： 其中： 解析推导得到不同于标准κ-ε模型的常数，在κ、ε方程中有附加项和函数 增加项 对于涉及到旋转、逆压梯度下的边界层、分离流、二次流及回流等问题，此模型比另外的模型得到更优的结果。 对于流动中包括大流线曲率、旋涡及旋转的现象时，此模型和RNG模型