1、湍流模型,Introductory FLUENT Training,湍流是什么?,非定常，无规律 (无周期) 运动，输运量 (质量, 动量, 组分) 在时间和空间中波动 湍流漩涡. 增强的混合（物质，动量 能量，等等）效果 流动属性和速度呈现随机变化 统计平均结果 湍流模型 包括一个大范围的湍流漩涡尺寸 (比例频谱). 大涡的尺寸和速率与平均流动在一个量级 大涡流动从平均流动中得到能量 能量从大涡向小涡转移 在最小尺度的涡中，湍流能量随着粘性耗散转移为内能,流动是否为湍流,外部流动,内部流动,自然对流,沿着表面,沿着障碍物,where,where,其它因素比如自由流动湍流,，表面条件，扰动等，
2、在低雷诺数下可能导致转变为紊流,is the Rayleigh number,is the Prandtl number,湍流结构,Energy Cascade Richardson (1922),计算方法总览,雷诺德平均NS模型(RANS) 解总体均值（或者时间均值）纳维斯托克斯方程 在RANS方法中，所有湍流尺度都进行模拟 在工业流动计算中使用得最为广泛 大涡模拟 (LES) 解算空间平均 N-S 方程，大涡直接求解, 比网格尺度小的涡通过模型得到 计算消耗小于DNS，但是对于大多数的实际应用来说占用计算资源还是太大了 直接数值模拟 (DNS) 理论上来说，所有的紊流流动能够由数值解出所有
3、的N-S方程来模拟 解出尺寸频谱，不需要任何模型 花费太高! 对工程流动不实用 ，目前 DNS 在 Fluent中不可用。 现在没有一种简单而实用的湍流模型能够可靠的预测出具有充分精度的所有湍流流动,在FLUENT中可用的湍流模型,基于RANS的模型,1方程模型Spalart-Allmaras 2方程模型标准 kRNG krealizable k标准 kSST k 雷诺德应力模型 分离涡模拟 大涡模拟,RANS 模拟 时间平均,将NS方程中的瞬时变量分解成平均量和脉动量:Reynolds-averaged 动量方程如下Reynolds 应力是由附加的平均过程引起的，因此为了封闭控制方程组，必须
4、对Reynolds应力建模,波动项,时均项,Example: 完全发展 湍流管流 速度分布,瞬时项,(Reynolds 应力张量),方程封闭,RANS 模型能够用下列方法封闭 (1) 涡粘模型 (通过 Boussinesq 假设)Boussinesq假设 Reynolds 应力 通过使用涡流粘性（湍流粘性）T模拟， 对简单湍流剪切流来说假设是合理的，例如 边界层、 圆形射流、 混合层、 管流 等等。（S-A， k ）(2) 雷诺德应力模型 (通过雷诺应力输运方程) RSM 对复杂的 3D湍流流动更有效，但是模型更加复杂， 计算强度更大， 比涡粘模型更难收敛,基于量纲分析， T 能够由 湍流时间
5、尺度 (或速度尺度) 和空间尺度来决定 湍流动能 L2/T2 湍流耗散率 L2/T3 比耗散率 1/T 每种湍流模型用不同的方法计算 T Spalart-Allmaras 解模拟湍流粘性的输运方程 标准 k, RNG k, Realizable k 解关于 k 和 的输运方程. 标准 k, SST k 解关于 k 和 的输运方程.,计算湍流粘性,Spalart-Allmaras 模型,Spalart-Allmaras 是一种低耗的求解关于改进的涡粘输运方程的RANS 模型 主要用于空气动力学/涡轮机， 比如机翼上的超音速/跨音速流动, 边界层流动 等等对于有壁面边界空气动力学流动应用较好 在有
6、逆压梯度的情况下给出了较好的结果 在涡轮机应用中很广泛 相对较新的模型 还没有应用于各种复杂的工程流动 对流动尺度变换较大的流动不太合适（平板射流，自由剪切流）,k 湍流模型,标准 k (SKE) 模型 在工程应用中使用最为广泛的湍流模型 稳定而且相对精确 包括可压缩性、 浮力、 燃烧等子模型 局限性 方程包括一个 不能在壁面上计算的项， 因此 必须使用壁面函数 在流动有强分离、大压力梯度情况下结果不太准确 RNG k模型 k 方程中 的常数通过renormalization group 定理得到 包括以下子模型 解决低雷诺数下的differential viscosity模型 由解析方法得到
7、的 Prandtl / Schmidt数的代数公式 旋流修正 对更复杂的剪切流来说比SKE 表现更好，比如剪切流、旋涡和分离流,k 湍流模型,Realizable k (RKE) 模型 realizable 意味着这个模型满足在雷诺应力上的特定数学约束, 与物理湍流流动一致. 法向应力为正 关于 Reynolds 剪切应力的Schwarz不等式 : 耗散率更能体现能量在谱空间的传输 优点： 对平面射流和圆形射流的散布率预测得更加精确. 对包括旋转、逆压梯度下的边界层、 分离, 循环流动提供较好性能三种模型区别：计算湍流粘性方法不同；控制湍流扩散的Pr数不同；耗散项的形式不同,k 湍流模型,k
8、湍流模型得到广泛特点： 模型方程不包括在壁面上没有定义的项，例如不需要壁面函数可以在壁面积分 对于有压力梯度的大范围边界层流动是精确稳定的 FLUENT 提供k 模型下的两个子模型 标准k (SKW) 模型 在航天和涡轮机械领域得到最广泛的应用 几个k子模型选项：压缩效果，转錑，剪切流修正. 剪切应力输运k (SSTKW) 模型(Menter, 1994) SST k 模型使用混合函数从壁面附近的标准k 模型逐渐过渡到边界层的外部的高雷诺数k模型. 包含修正的湍流粘性公式来解决湍流剪应力引起的输运效果,雷诺应力模型 (RSM),RSM 是最复合物理现象的模型: 各向异性，输运中的雷诺应力可以直
9、接计算出来 RSM 对控制方程需要更多的建模（其中压应力是最关键和有难度的参数之一） RSM 比2方程模型需要时间长且较难收敛 适合有大弯曲流线、漩涡和转动的3维流动,Turbulent diffusion,Stress production,Rotation production,Pressure Strain,Dissipation,Modeling required for these terms,大涡模拟 (LES),过滤NS方程中的湍流涡频谱： 通过网格尺寸筛选 比网格尺寸小的涡被忽略，用subgrid scale (SGS) 建模 较大尺度涡用数值方法直接求解NS方程,修正 N-S
10、 方程,N-S 方程,Subgrid Scale,Resolved Scale,Instantaneous component,(Subgrid scale Turbulent stress),大涡模拟,大涡模拟 (LES) LES 非常成功的应用于 RANS 模型不能满足要求的高端应用 对N-S方程在物理空间进行过滤，大涡直接求解，小涡各向同性模拟 方法 亚网格尺度(SGS) 湍流模型 Smagorinsky-Lilly 模型 Wall-Adapting Local Eddy-Viscosity(WALE) 壁面适应局部涡粘模型 Dynamic Smagorinsky-Lilly 模型 Dy
11、namic Kinetic Energy Transport 动能传输 分离涡 (DES) 模型 LES在FLUENT中对所有燃烧模型适用 有基本统计学工具：对求解值进行时均分析，内置快速傅立叶变换 (FFT) 在运行 LES之前, 参考帮助中对 LES方法的指导 (包括网格建议，亚网格模型， 数值方法， 边界条件等),分离涡流模拟 (DES),产生原因 对于高雷诺数壁面边界流动， LES 在解近壁面区域时显得比较耗费时间 在近壁面区域使用RANS 可以降低对网格的要求 基于Spalart-Allmaras turbulence 模型的RANS/LES 混合模型 :一方程SGS 湍流模型 在平
12、衡状态下，简化为代数模型 在高雷诺数的外部空气动力流动方面，DES是LES 的有效替代,RANS 湍流模型描述,RANS 湍流模型用法,壁面和近壁面处理原则,壁面对湍流有明显影响,where y is the normal distance from the wall,近壁区域分为粘性底层，过渡区和完全湍流层。 处理方法：半经验公式（壁面函数）以及改进湍流模型,Wall shear stress,壁面边界条件,k 系列模型和 RSM 模型在近壁面区域不可用，而Spalart-Allmaras 和 k 模型对所有区域都有效 (假设网格足够好) 壁面函数法 标准壁面函数法利用对数校正法提供了必需的
13、壁面边界条件（对于平衡湍流边界层） 非平衡壁面函数法用来改善高压力梯度、分离、再附、滞止等情况下的结果 对能量和组分方程采用同样的方法 优点： 壁面函数允许在近壁面区域上使用相对较粗的网格 增强壁面处理选项 把混合边界模型和两层边界模型结合起来. 对低雷诺数流动或者复杂近壁面现象很适合 湍流模型在内层上得到修正.,标准壁面函数法,标准壁面函数 动量边界条件遵循Launder-Spaulding 的关于壁面的法则对于能量和组分使用相同方法 附加公式考虑到 k, , .和 当流动偏离预先假定的条件时，会产生误差 例如，变化大的 p 或接近壁面的非平衡流动，高度蒸腾和大的体积力，低雷诺数和高速三维流
14、动等,where,几种壁面处理方法比较,第一个网格点的布置,对于标准或者非平衡的壁面函数法，每个壁面相邻的单元体中心必须位于对数层（loglaw layer）中 对于增强的壁面处理 (EWT), 每个与壁面相邻的单元体中心应该位于粘性亚层上在生成网格之前怎样估计壁面相邻的单元体大小:表面摩擦系数可以从经验公式中估算出来 :在建立好流动模型之后使用后处理工具（XY图或者等值线图）来仔细检查近壁面网格布置,平板：,管道：,近壁面建模：推荐策略,对于大多数高雷诺数情况使用标准的或者非平衡的壁面函数(Re 106) 在分离、再附着或者射流流动中使用非平衡壁面函数 考虑使用加强壁面处理（EWT）的情况：
15、 特征雷诺数很低或者贴体特征需要解出来 大部分壁面区域上y+变化明显. 使用大小合适的网格，避免将近壁面网格放置在过渡区中 (5 y+ 30).,入口和出口边界条件,根据所选湍流模型，需要在入口或者出口给定湍流边界条件 。 直接或者间接的定义湍流参数的四种方式：直接输入 k, , 这是唯一允许定义分布的方式 通过帮助文件查看详细设置 湍流强度 和 长度尺度 长度尺度与包含的大涡尺度相关. 对于边界层流动 l 0.4 99 对于网格下游的流动 l opening size 湍流强度 和 水力直径 和内部流动 (管流) 相匹配 湍流强度 和 湍流粘性比 对外部流动 1 mt/m 10 湍流强度依赖
16、于上游条件,湍流模型选项,近壁面处理,无粘,层流, 或者湍流,附加选项,Boundary Conditions,Define,湍流模型的 GUI,Viscous,Define,Models,示例1 流过平面湍流,使用四种不同的湍流模型可以模拟湍流通过一个钝平面 8,700 四边形网格, 在前沿和再附着位置加密网格 非平衡边界层处理,N. Djilali and I. S. Gartshore (1991), “Turbulent Flow Around a Bluff Rectangular Plate, Part I: Experimental Investigation,” JFE, Vo
17、l. 113, pp. 5159.,示例1 流过平面湍流,Inlet,Outlet,Wall,Wall,Symmetry,示例1 流过平面湍流,RNG k,Standard k,Reynolds Stress,Realizable k,湍流动能云图 (m2/s2),Experimentally observed reattachment point is at x / D = 4.7,预测分离区:,示例1 流过平面湍流,Skin Friction Coefficient Cf 1000,比较而言，RKE模型对分离区的预测较为准确。,Distance Along Plate, x / D,示例2
18、 旋风分离器湍流,40,000-六面体网格单元 使用高阶迎风格式 使用SKE, RNG, RKE 和RSM 模型以及标准壁面函数 强旋流 (Wmax = 1.8 Uin),示例2 旋风分离器湍流,0.41 m处切向速度剖面图,Iso-Contours of Instantaneous Vorticity Magnitude,Time-averaged streamwise velocity along the wake centerline,CL spectrum,示例3 通过正方形棱柱的流动 (LES),(ReH = 22,000),Streamwise mean velocity along the wake centerline,Streamwise normal stress along the wake centerline,示例3 通过正方形棱柱的流动 (LES),小结: 湍流模型指南,湍流模型选择 物理流动 可用计算机资源 工程需要 准确性 时间 近壁面处理 建模程序 计算特征值 Re 和确定流动是否为湍流 在生成网格之前估计贴体网格中心 y+ 从 SKE (standard k-)开始计算 ，如果需要的话，之后改成 RKE, RNG, SKW, SST 。 强旋流，三维，旋转流动下使用RSM 模型 大多数情况下使用壁面函数功能,
