4.2 核心 RANS 模型详解

4.2 核心 RANS 模型详解

方法论框架立好后，本节打开工具箱，把工程里出场率最高的几个 RANS 成员逐一过堂：每个模型的两个方程在算什么、它的参数标定源自哪类流动、因此它天然擅长什么、必然在哪里翻车。读完你应当能像老法师一样，看一眼工况就报出候选模型与理由。

两方程家族：k 与 epsilon 或 omega 的分工

两方程模型用湍动能 k（脉动能量）与一个特征长度相关量封闭方程组。分歧在于第二个变量选谁，而选谁决定了模型的性格。

标准 k-epsilon。第二个变量是湍动能耗散率 epsilon。参数标定源自均匀湍流与自由剪切流，性格豪放：对自由剪切流、充分发展管道流表现稳健，计算便宜、皮实耐糙。软肋在近壁与逆压梯度：需要壁面函数配合，对分离起始点的预测偏晚，回流区偏小。它至今仍是通用外流初算与通风空调类工程的常青树，前提是你不需要精确的分离位置。

Realizable k-epsilon。数学上强制雷诺应力满足可实现性约束（某些应变组合下正应力不得为负），并给耗散率方程换了更合理的公式。性格比标准版聪明：平面射流与圆射流的扩散率不再需要两套参数，含强应变的流动（分离、二次流）明显更稳。Fluent 里它几乎取代了标准版成为 epsilon 系默认，除非做文献对标必须复现标准版行为。

k-omega 与 k-omega SST。第二个变量换成比耗散率 omega。omega 系在近壁低雷诺数区天然友好——不需要壁面函数也能算。但纯 k-omega 对来流湍流度敏感，自由流里 omega 的微小误差被放大成剪切层预测的大误差。Menter 的 SST 方案是工程智慧的样板：近壁用 k-omega 保精度，外区切换成 k-epsilon 的变体抗自由流污染，中间用混合函数平滑过渡，再加一个雷诺应力限制项压住逆压梯度下的过高湍流黏度。于是 SST 性格全面：分离起点预测准、逆压边界层可靠、近壁无需特殊处理。代价是分离区湍流黏度偏低、回流区略大，以及对来流湍流度设定仍然讲究。

Spalart-Allmaras。一方程模型，直接解一个修正的涡黏性输运方程。出身航空，对附体流动与轻度分离的机翼效率极高：方程少、数值稳、网格宽容。跨音速翼型、直升机旋翼、外流的初轮筛选，它都是性价比之王。软肋是拿它算复杂射流与大规模分离时底气不足——它没有为这些场景标定过。

图 4-2：RANS 模型地盘图

控制台会话：一次完整的模型设置与检查

> define/models/viscous/model k-omega sst
  壁面处理：Shear Stress Transport（自动启用低雷诺数近壁解析）
> define/models/viscous/turbulence-expert-parameters
  curvature-correction? yes        # 叶轮通道曲率修正
> define/boundary-conditions/velocity-inlet inlet-1
  turbulence specification: intensity and viscosity ratio
  intensity 4.5 percent            # 按风洞实测来流湍流度
  viscosity ratio 10
> report/definitions/add/surface-report/area-weighted-avg
  场变量：turbulent-viscosity-ratio  监测面：出口
> solve/initialize/hybrid-initialization
​

会话里两处最容易被忽视、也最值得写进验收单：来流湍流度不是装饰——SST 的转捩相关行为、混合段的掺混速率都对它敏感，取值要有出处（实测、规范或同类算例）；湍流黏度比是收敛健康的哨兵，出口监测值若在迭代中爬升到数千量级，往往是分离区失控的前兆。

案例：三种模型算同一个闸阀，结论为何打架

背景。某阀门厂评估新款闸阀流阻，三位工程师分别用标准 k-epsilon 加壁面函数、SST、RSM 算同一模型，压降结果相差百分之八，评审会现场各执一词。

操作。先对齐变量：三份计算用同一网格（SST 路线要求的低雷诺数边界层网格对壁面函数版本来说过密但可用）、同一离散格式、同一收敛判据。再逐项排查差异来源：阀门全开时流道存在轻度分离与再附，正是各模型分歧的敏感区。提取分离点位置与再附长度做横向对比，epsilon 版分离点明显偏后。

结果。以公开的阀门试验数据对标：SST 偏差百分之二，RSM 偏差百分之三，epsilon 版偏差百分之八。流阻结论最终以 SST 为准，epsilon 版被标注为"趋势参考不进报告"。

解读。三个模型都没"算错"——它们只是诚实地执行了各自的假设：epsilon 的参数本来就没见过这种分离结构。教训是模型差异本身就是不确定度的信息：当两个模型打架，差别大小就是模型敏感性，写进报告比单模型独角戏更有说服力。第九章的敏感性检查会把这一步标准化。

变式。强旋流场景（旋风筒、涡流燃烧室）重复同一套横向对比时，RSM 往往反超 SST 成为主答案——地盘图上的位置要按工况挪，模型没有全场冠军。

模型参数的高频疑问

来流湍流度没有实测怎么定。按场景经验值起步：风洞或洁净来流取百分之一量级，大气边界层与工业管道取百分之五上下，燃烧室一次流取厂商推荐。取值连同出处写进报告；湍流长度尺度或黏度比配对给出，两个数是一套。

转捩修正、曲率修正这些开关默认关，要不要开。按工况对号：叶轮机与强曲率通道开曲率修正，低雷诺数翼型用转捩模型而不是随手开关。每开一个修正都要在报告里声明——修正改变的是模型性格，不是免费午餐。

两方程模型算湍流黏度比飘到上万怎么办。这是分离区失控或来流污染的信号：检查入口边界组合（强度加黏度比要自洽）、复核分离区网格分辨率、必要时切换对自由流更鲁棒的外区处理。

RANS 家族一分钟选型卡：

本节验收点

• 模型选择有书面理由，且与 4.1 的三连问答案自洽。
• 来流湍流度与湍流黏度比的取值有出处，写入边界条件说明。
• 关键结论做过至少两模型横向对比或引用同场景验证文献。
• 修正开关与近壁处理方式已记录，模型配置可复现。