FLUENT常用的湍流模型及壁面函数处理

FLUENTFLUENT 常用的湍流模型及壁面函数处理常用的湍流模型及壁面函数处理 本文内容摘自 精通 CFD 工程仿真与案例实战 实际上也是帮助文档的 翻译 英文好的可直接参阅帮助文档 FLUENT 中的湍流模型很多 有单方程模型 双方程模型 雷诺应力模型 转捩模型等等 这里只针对最常用的模型 1 湍流模型描述 模型模型描述描述 Spalart AllmarasSpalart Allmaras 单方程模型 直接解出修正过的湍流粘性 用于有界壁面 流动的航空领域 需要较好的近壁面网格 尤其是绕流过 程 该模型也可用于粗网格 StandardStandard k ek e 双方程模型 是默认的 k e 模型 系数由经验公式给出 只对高 Re 的湍流有效 包含粘性热 浮力 压缩性等选项 RNGRNG k ek e 标准 k e 模型的变形 方程和系数来自解析解 在 e 方程 中改善了模拟高应变流动的能力 用来预测中等强度的旋 游和低雷诺数流动 RealizableRealizable k ek e 标准 k e 模型的变形 用数学约束改善模型的性能 能用 于预测中等强度的旋流 StandardStandard k wk w 两个输运方程求解 k 与 w 对于有界壁面和低雷诺数流动性 能较好 尤其是绕流问题 包含转捩 自由剪切和压缩性 选项 SSTSST k wk w 标准 k w 模型的变形 使用混合函数将标准 k e 模型与 k w 模型结合起来 包含了转捩和剪切选项 ReynoldsReynolds StressStress 直接使用输运方程来解出雷诺应力 避免了其它模型的粘 性假设 模拟强旋流相比其它模型有明显优势 2 湍流模型的选择 模型模型用法用法 SpalartSpalart 计算量小 对一定复杂的边界层问题有较好的效果 AllmarasAllmaras 计算结果没有被 广泛的测试 缺少子模型 典型的应用场合为航空领域的绕流模拟 StandardStandard k ek e 应用多 计算量适中 有较多数据积累和比较高的精度 对于曲率较大和压力梯度较强等复杂流动模拟效果欠佳 一般工程计算都使用此模型 其收敛性和计算精度能满足一般 的工程计算要求 但模拟旋流和绕流时有缺陷 RNGRNG k ek e 能模拟射流撞击 分离流 二次流和旋流等 中等复杂流动 受到涡旋粘性同性假设限制 除强旋流过程无法精确预测外 其它流动都可以使用此模型 RealizableRealizable k ek e 和 RNG 基本一致 还可以更好的模拟圆形射流 受到涡旋粘性同性假设限制 除强旋流过程无法精确预测外 其它流动都可以使用此模型 StandStand k wk w 对于壁面边界层 自由剪切流 低雷诺数流动性能较好 适合 于存在逆压力梯度时的边界层流动 分离与转捩 SSTSST k wk w 基本与标准 k w 模型相同 由于对壁距离依赖性较强 因此不 太适合于自由剪切流 ReynoldsReynolds StressStress 是最复杂的 RANS 模型 避免了同性的涡粘性假设 占用较多的 CPU 时间和内存 收敛较难 对于复杂 3D 流动适用 如弯曲管 道 旋转 旋流燃烧 旋风 尤其是强旋流运动 有两种方法处理近壁面区域 一种方法 不求解粘性影响内部区域 粘性 子层及过渡层 使用一种称之为 wall function 的半经验方法去计算壁面 与充分发展湍流区域之间的粘性影响区域 采用壁面函数法 省去了为壁面的 存在而修改湍流模型 另一种方法 修改湍流模型以使其能够求解近壁粘性影响区域 包括粘性 子层 此处使用的方法即近壁模型 近壁模型不需要使用壁面函数 如一些 低雷诺数模型 K W 湍流模型是一种典型的近壁湍流模型 所有壁面函数 除 scalable 壁面函数外 的最主要缺点在于 沿壁面法向 细化网格时 会导致使数值结果恶化 当 y 小于 15 时 将会在壁面剪切力及 热传递方面逐渐导致产生无界错误 然而这是若干年前的工业标准 如今 ANSYS FLUENT 采取了措施提供了更高级的壁面格式 以允许网格细化而不产生 结果恶化 这些 y 无关的格式是默认的基于 w 方程的湍流模型 对于基于 epsilon 方程的模型 增强壁面函数 EWT 提供了相同的功能 这一选项同样 是 SA 模型所默认的 该选项允许用户使其模型与近壁面 y 求解无关 实际 上是这样的 K W 方程是低雷诺数模型 采用网格求解的方式计算近壁面粘性 区域 所以加密网格降低 y 值不会导致结果恶化 k e 方程是高雷诺数模型 其要求第一层网格位于湍流充分发展区域 而此时若加密网格导致第一层网格 处于粘性子层内 则会造成计算结果恶化 这时候可以使用增强壁面函数以避 免这类问题 SA 模型默认使用增强壁面函数 只有当所有的边界层求解都达到要求了才可能获得高质量的壁面边界层数 值计算结果 这一要求比单纯的几个 Y 值达到要求更重要 覆盖边界层的最小 网格数量在 10 层左右 最好能达到 20 层 还有一点需要注意的是 提高边界 层求解常常可以取得稳健的数值计算结果 因为只需要细化壁面法向方向网格 与增加精度向伴随的是计算开销的增加 对于非结构网格 建议划分 10 20 层 棱柱层网格以提高壁面边界层的预测精度 棱柱层厚度应当被设计为保证有 15 层或更多网格节点 这可以在获得计算结果后 通过查看边界层中心的最大湍 流粘度 该值提供了边界层的厚度 最大值的两倍位置即边界层的边 棱柱 层大于边界层厚度是必要的 否则棱柱层会限制边界层的增长 一些建议 1 对于 epsilon 方程 使用 enhanced 壁面函数 2 若壁 面函数有助于 epsilon 方程 则可以使用 scalable 壁面函数 3 对于基于 w 方程的模型 使用默认的增强壁面函数 4 SA 模型 使用增强壁面处理 以上内容翻译自 Fluent 理论文档 P121 1 标准壁面函数 ANSYS FLUENT 中的标准壁面函数是基于 launder 与 spalding 的工作 在 工业上有广泛的应用 在标准壁面函数中 用得比较多的变量为 y y 的下限为 15 低于此值 将会导致结果精度恶化 标准壁面函数用于以下模型 k epsilon 模型与 reynolds stress 模型 这两个模型均为高雷诺数模型 2 Scalable 壁面函数 该壁面函数是 14 0 新加的 以前的版本中没有的 也是 CFX 软件中默认的 湍流壁面函数 该壁面函数能避免在 y 11 时 该壁面函数的表现与标准壁 面函数一致 scalable 壁面函数的目的在于联合使用标准壁面方法以强迫使用对数律 该功能是通过使用限制器 y max y y limit 来实现的 其中 y limit 11 225 3 非平衡壁面函数 非平衡壁面函数的特点 1 用于平均速度的 launder 及 spalding 的对 数律对于压力梯度效应敏感 2 采用双层概念以计算临壁面单元的湍流动能 对于平均温度及组分质量分数则与标准壁面函数处理方式相同 非平衡壁面函数考虑了压力梯度效应 因此对于涉及到分离 再附着 及 撞击等平均速度与压力梯度相关且变化迅速的复杂流动问题 推荐使用些壁面 函数 但是非平衡壁面函数不适合于低雷诺流动问题 非平衡壁面函数适用于高雷诺流动问题 适用于以下湍流模型 1 K epsilon 模型 2 reynolds stress transport 模型 4 壁面函数方法的局限 对于大多数壁面边界流动问题 标准壁面函数能给出合理的预测 非平衡 壁面函数考虑了压力梯度效应 扩展了标准壁面函数的功能 但是一些流动问 题不适合使用壁面函数 否则可能导致不合理的解 如以下一些情况 1 低雷诺数流动或近壁面效应 例如小缝出流 高粘性低速流动问题 2 通过壁面的大量沸腾 3 大的压力梯度导致的边界层分离 4 强体力 如旋转圆盘附近的流动 浮力驱动流动 5 近壁区域高度三维流动 如 ekman 螺旋流动 高度歪斜的 3D 边界层 若模型中出现了以上的情况 则必须使用近壁模型 ANSYS FLUENT 中提供 了增强壁面处理以应对这些情况 这一方法能够用于 K epsilon 模型及 RSM 模 型 5 增强壁面处理 EWT 不依赖于壁面法则 对于复杂流动尤其是低雷诺数流动问题很适合 该方 法要求近壁面网格很密 y 接近于 1 对于 epsilon 方程的近壁面处理结合了双层模型 若近壁面网格足够密以 致于可以求解粘性子层时 通常第一节