湍流模型应该如何选择，才能兼顾精度、稳定性和计算成本？

湍流模型选错，计算往往照样收敛，却给出系统性偏差的结果。本文从工程指标、近壁网格和计算成本出发，梳理 RANS、RSM、DES 与 LES 的适用边界和选择流程。

同一几何、同一套网格，只换湍流模型，压降、分离位置和壁面换热系数就可能明显不同，而且每组结果看上去都“合理”。选错湍流模型通常不会让计算发散，而是给出带系统性偏差的答案。

本文围绕三个问题展开：模型要解决什么物理现象？在现有网格和算力下能否稳定收敛？多付出的计算成本是否换来了决策所需的精度？讨论以 Simcenter STAR-CCM+ 为例，方法同样适用于其他 CFD 软件。

本文要点

  • 先确定设计决策需要的指标，再选模型：时均量从 RANS 起步，脉动与噪声才需要 DES 或 LES；
  • 外流与分离问题从 SST k-ω 起步，内流与整体量从 Realizable k-ε 起步，强旋流再考虑 RSM；
  • 近壁网格必须匹配模型和指标，关心壁面摩擦与换热时按 y+≈1 布网格；
  • 关键结论应做双模型对比，并用试验数据或基准算例验证。

一、湍流模型的本质：用近似换代价

工程湍流的尺度跨度极大：大尺度涡决定流动走向，能量逐级传递到最小尺度，并在那里被黏性耗散。直接数值模拟（DNS）需要解析全部尺度，网格点数大致按雷诺数的 9/4 次方增长，对绝大多数工业问题并不可行。

所有湍流模型都在决定：哪些尺度直接计算，哪些尺度交给模型。直接计算的越多，结果越接近真实流动，网格、时间步和算力要求也越高；模型化的越多，计算越便宜，结果也越依赖模型假设是否符合你的流动。

更有用的问题不是“哪个模型最准”，而是：在给定预算和周期内，哪种近似对我关心的指标影响最小？

二、三类方法：RANS、雷诺应力输运模型与尺度解析

Siemens 官方知识库把 STAR-CCM+ 的湍流模型分为三类：涡黏模型（SA、k-ε、k-ω）、雷诺应力输运模型，以及瞬态的尺度解析方法（DES、LES）。区别在于把多少湍流交给模型。

方法 直接计算 交给模型 能给出的结果 相对成本
RANS 涡黏模型 时均流场 全部脉动，用涡黏假设封闭 压降、力系数、平均温度等时均量 低
雷诺应力输运模型（RSM） 时均流场 全部脉动，各雷诺应力分量单独求解 时均量，并能表达各向异性 中
DES 远离壁面的大尺度涡 近壁边界层（RANS） 瞬态流场、脉动载荷、尾迹 高
LES 大尺度涡 亚格子尺度的涡 瞬态流场、湍流结构、噪声源 很高

大多数工业稳态问题应从 RANS 开始。关键在于决策是否需要瞬态信息：要压降、平均温度、总阻力，RANS 通常性价比最高；要脉动压力、涡脱落频率、气动噪声源，时均方法从原理上给不出。

三、常用 RANS 模型：擅长什么，短板在哪

涡黏模型都基于 Boussinesq 假设，用一个标量涡黏系数把雷诺应力与平均应变率联系起来。这一假设隐含湍流各向同性，因此在强旋流、强曲率和二次流问题上有共同短板。

Spalart-Allmaras（SA）

单方程模型，求解修正涡黏系数的输运方程，最初为翼型、机翼等航空外流开发。计算量小、收敛好，适合附着边界层为主的外流，在高马赫数外流中验证基础较好；对大分离、自由剪切流和复杂内流适用性较弱。在 STAR-CCM+ 中，SA 也是 DES 和 SA Gamma 转捩模型的基础模型之一。

k-ε：标准与可实现（Realizable）

求解湍动能 k 和耗散率 ε，在充分发展的湍流核心区稳健，是使用最久的工业模型之一。Realizable k-ε 的 Cμ 随应变率和旋转率变化，以保证雷诺正应力为正，并采用新的 ε 方程，在旋转、分离和射流问题上通常优于标准 k-ε。

k-ε 本身是高雷诺数模型，不能直接积分到壁面。STAR-CCM+ 的两层（Two-Layer）方法在近壁层内仍求解 k 方程，ε 改由壁面距离的代数关系给定，因此具备 All y+ 能力，既可用于 y+≈1 的网格，也可用于 y+>30 的壁面函数网格。两层方法另有面向浮力驱动流动的形式，适合自然对流。

SST k-ω

Menter 的 SST 模型在近壁区采用 k-ω 形式，可直接积分到壁面，在逆压梯度边界层中表现更好；远场切换为 k-ε 行为，避免对自由流湍流参数过于敏感；并基于 Bradshaw 假设限制涡黏系数，改善分离起点预测。它是外流气动、带分离的内流和旋转机械最常用的起点，也是 STAR-CCM+ 中 Gamma 系列转捩模型和 SST-DES 的基础。

使用提示分离预测依赖近壁网格，y+ 不足时优势明显打折；叶片前缘等驻点附近湍动能容易偏高，传热计算应启用产生项限制；中等旋流或曲率问题，可先试 SST 加曲率修正，再考虑 RSM。

椭圆混合 k-ε 与 Lag EB k-ε

椭圆混合（Elliptic Blending）用一个椭圆方程描述壁面对湍流的阻塞效应，改善近壁行为。Lag EB k-ε 由 Lardeau 与 Billard 于 2016 年面向工业应用提出，增加一个输运方程描述雷诺应力与应变率之间的“滞后”。这类模型面向近壁、分离和传热问题；Siemens 在第五届 AIAA 高升力预测研讨会（HLPW-5）的计算中，也把 Lag EB k-ε 与 SA、SST 并列考察。

雷诺应力输运模型（RSM）

RSM 放弃涡黏假设，对雷诺应力各分量分别求解输运方程，能表达湍流各向异性，适合旋流燃烧器、旋风分离器、强弯曲管道等问题。代价是方程多、收敛难，官方知识库也指出其稳健性不如涡黏模型。实践中通常先用 SST 或 k-ε 得到收敛解，再以它为初值切换到 RSM。如果只是方截面管道这类由各向异性驱动的二次流，可先在涡黏模型上启用二次本构关系（QCR），成本更低。

模型对照

Spalart-Allmaras（SA）

方程数（三维）
1
适合的问题
航空外流、附着边界层；高马赫数外流中验证基础较好
主要短板
大分离、射流、复杂内流
稳健性
高
说明
在 STAR-CCM+ 中也是 DES 和 SA Gamma 转捩模型的基础模型之一。

Realizable k-ε（两层）

方程数（三维）
2
适合的问题
内流、管道、混合、HVAC；浮力两层形式用于自然对流
主要短板
逆压梯度下的分离
稳健性
高
说明
两层方法在近壁层内仍求解 k 方程，ε 由壁面距离的代数关系给定，具备 All y+ 能力。

SST k-ω

方程数（三维）
2
适合的问题
外流、分离、旋转机械、低 y+ 传热
主要短板
依赖近壁网格；驻点湍动能偏高
稳健性
中高
说明
传热计算应启用产生项限制；中等旋流或曲率可先试曲率修正。

椭圆混合 / Lag EB k-ε

方程数（三维）
4（含 1 个椭圆方程）
适合的问题
近壁、分离、传热
主要短板
工程应用经验相对较少
稳健性
中
说明
Lag EB k-ε 由 Lardeau 与 Billard 于 2016 年提出；Siemens 在 HLPW-5 中将其与 SA、SST 并列考察。

雷诺应力输运模型（RSM）

方程数（三维）
7（6 个应力分量 + 1 个尺度方程）
适合的问题
旋流、强曲率、二次流
主要短板
收敛难、计算量大
稳健性
较低
说明
通常先用 SST 或 k-ε 得到收敛解，再以其为初值切换；方截面二次流可先试 QCR。

注：椭圆混合 RSM 另需 1 个椭圆方程；稳健性为工程经验判断，仅供参考。

四、近壁处理与网格：决定模型能否发挥

同一模型配上不合适的近壁网格，结果可能比简单模型更差，选模型时必须同时决定壁面处理方式。

近壁区分为黏性底层（y+ 小于约 5）、缓冲层（约 5–30）和对数律层（约 30 以上）。y+ 是壁面第一层网格中心的无量纲壁面距离。STAR-CCM+ 提供三类壁面处理：

LOW Y+

积分到壁面

要求第一层网格 y+≈1，并在边界层内布置足够的棱柱层。

HIGH Y+

壁面函数

第一层网格落在对数律层（y+>30），近壁网格量大幅减少。

ALL Y+

混合处理

按局部 y+ 混合两种处理，对 y+ 分布不均匀的网格更宽容。

关心壁面摩擦、分离点和对流换热系数，应按 Low y+ 布网格；只关心整体压降和流量分配，壁面函数网格通常足够。All y+ 不是“免检证”，大量网格落在缓冲层时精度仍会下降。做网格无关性检查时，应分别控制棱柱层和核心区网格，否则加密主体网格会连带改变 y+，难以判断结果变化的来源。

经验提示计算完成后，先查看壁面 y+ 的实际分布，再评价模型表现。

五、什么时候需要转捩模型

标准 RANS 默认边界层全湍流。中低雷诺数的翼型和叶片、小型风扇与无人机螺旋桨、涡轮叶片换热，以及来流湍流度很低的光滑外形，层流区和层流分离泡会显著影响阻力、换热和分离，需要转捩模型。

STAR-CCM+ 提供三种方法，预测能力和成本依次上升：湍流抑制模型是零方程方法，需预先指定转捩位置，适合与试验对比；Gamma 模型与 SST k-ω 耦合，额外求解 1 个间歇因子方程，设置简单；Gamma-ReTheta 模型额外求解 2 个方程，基于经验关联式，预测更完整。SA 模型另有 SA Gamma 转捩模型可用。

转捩模型对网格敏感，除 y+≈1 外，流向也需要足够细的网格分辨转捩区；来流湍流强度及其衰减对转捩位置影响很大，不能随意给定。

六、什么时候需要 DES 或 LES

判断标准是：决策是否依赖瞬态信息，或 RANS 在大分离区的预测已不可信。典型情况包括脉动压力与涡脱落频率、气动噪声源、钝体尾迹与大范围分离，以及瞬时混合与热疲劳。

DES

分离涡模拟

近壁边界层用 RANS，远离壁面的分离区用 LES，基础模型可以是 SA 或 SST k-ω，避免了 LES 在壁面附近所需的极细网格。DDES 用于避免网格诱导分离（GIS），IDDES 在此基础上兼具壁面模化 LES 能力。

LES

大涡模拟

直接计算大尺度涡，只对亚格子尺度使用模型。STAR-CCM+ 提供 WALE 等亚格子模型，WALE 无需阻尼函数即可给出正确的近壁涡黏标度，对模型系数也不如 Smagorinsky 敏感。

采用尺度解析方法前，先估算三项成本：

  1. 空间分辨率：官方知识库指出，LES 和 DES 需要比 RANS 更细的空间和时间分辨率，分离区和剪切层网格应尽量各向同性。
  2. 时间步长：需分辨关心的最高频率，LES 区域对流库朗数通常控制在 1 左右。
  3. 统计时长：越过初始过渡后，还需采样足够长的时间，才能得到稳定的平均值和频谱。

只关心规则的周期性涡脱落时，非定常 RANS（URANS）有时能以低得多的成本给出主频，但无法表达宽频湍流和噪声。

七、一套可执行的选择流程

前四步决定用什么模型，后两步决定结果能不能用。下面的选择器按同样的顺序给出起步建议。

交互工具

湍流模型选择器

按顺序回答，结果是起步建议，最终仍需对比与验证。

1　需要瞬态信息吗？
查看完整 6 步流程（文字版）
  1. 明确决策指标压降、力、分离、换热，还是脉动与噪声。
  2. 判断是否需要瞬态信息需要则进入 DES 或 LES；规则涡脱落可先试 URANS。
  3. 判断是否存在明显层流区存在则采用 SST k-ω 加 Gamma 或 Gamma-ReTheta 转捩模型。
  4. 按主导物理选择 RANS 模型旋流、强曲率、二次流选 RSM；自然对流选浮力两层 Realizable k-ε；外流分离与近壁传热选 SST k-ω；内流整体量与混合选 Realizable k-ε；航空附着外流选 SA。
  5. 按指标准备近壁网格关心壁面摩擦和换热按 y+≈1，只关心整体量可用壁面函数。
  6. 双模型对比与试验验证确认结论对模型选择不敏感。
经验提示流程中最容易被省略的是第 6 步：用两个合理的模型计算同一工况，关键指标差别很小，说明结论对模型不敏感，可以推进；差别明显，说明问题本身对湍流模型敏感，需要试验数据或更高保真度的计算确认。

八、三个常见误区

  • 把收敛当成正确：残差下降只说明数值稳定，不说明模型适合这个流动。
  • 只换模型，不改网格：从壁面函数网格直接换成要求 y+≈1 的模型，或用 RANS 网格直接跑 LES，对比没有意义。
  • 随意给定湍流入口条件：湍流强度与长度尺度（或湍流黏度比）会影响 SST 和转捩模型的结果，应依据上游条件或试验数据设定。

湍流模型选择没有万能答案，但有可重复的方法：先定指标，再判断瞬态、转捩和主导物理，然后让近壁网格匹配模型，最后用对比和试验确认。对多数工业问题，SST k-ω 或 Realizable k-ε 是合理的起点，真正决定精度的，往往是之后的网格、边界条件和验证工作。

如果您正在使用 Simcenter STAR-CCM+ 处理流动、传热或多物理场问题，希望评估现有模型设置与网格策略，可以了解我们的 工程仿真与物理测试 服务。

相关文章

继续阅读