湍流模型参数（k、epsilon、omega、nut、nuTilda、I等）计算方法及相应工具。

1. 流程总览

输入参数
   │
   ├─(可选) 由 Re 估算湍流强度 I = 0.16·Re^(-1/8)
   │
   ├─▶ 湍动能    k = 3/2·(|U|·I)²
   │        │
   │        ├─▶ 混合长度  L = 0.07·Lc
   │        │       │
   │        │       ├─▶ 耗散率 ε = Cmu^0.75·k^1.5 / L
   │        │       └─▶ 耗散频率 ω = √k / (Cmu^0.25·L)   [等价: ω = ε/(Cmu·k)]
   │        │
   │        └─▶ nuTilda 初始值（SA 模型前置，可选）
   │                nuTilda ≈ nut = Cmu·k²/ε（精确）
   │                或工程估算：3ν / 5ν / 10ν / ν·√(0.01·Re)
   │
   ├─▶ 湍流粘度（按所选模型取其一）
   │        ├─ k-ε 模型: nut = Cmu·k²/ε
   │        ├─ k-ω 模型: nut = k/ω
   │        └─ SA 模型:  nut = nuTilda·fv1
   │
   └─▶ 湍流粘度比  nutRatio = nut / ν
   │
   ▼
输出参数

2. 输入参数列表

3. 计算步骤

步骤 0（可选）：估算湍流强度

已知雷诺数 Re 时（Fluent 管道充分发展流经验公式）：

若用户直接给定 I（如自由来流取 5%），跳过此步。

”

步骤 1：湍动能 k

步骤 2：混合长度 L

步骤 3：湍流耗散率 ε

步骤 4：湍流耗散频率 ω

基于混合长度理论：

等价形式（由 ε 换算，两者结果一致，可用于交叉验证）：

步骤 5：nuTilda 初始值（仅 SA 模型需要，可选）

OpenFOAM 源代码未直接提供 nuTilda 初始值公式，常用以下方法估算：

5.1 基于湍流粘度（较精确）

自由流区域 fv1 ≈ 1，因此 nuTilda ≈ nut：

需先完成步骤 1~3。

”

5.2 简化工程经验公式（最常用）

5.3 基于雷诺数

注意事项：

1. nuTilda 初始值只影响收敛初期的行为，通常会在若干时间步内调整到正确值；
2. 不确定时取 nuTilda = 3·ν 是安全的选择；
3. 壁面边界条件通常设置为 nuTilda = 0。

”

步骤 6：湍流粘度 nut（按所选湍流模型取其一）

6.1 k-ε 模型

6.2 k-ω 模型

6.3 Spalart-Allmaras 模型

其中阻尼函数：

步骤 7：湍流粘度比

4. 输出参数列表

各模型典型输出组合：

5. 数值算例

输入： |U| = 10 m/s，I = 5%，Lc = D = 0.5 m，ν = 1e-5 m²/s，Cmu = 0.09

算例验证：

• ε 与 ω 的两种换算公式结果一致，k-ε 与 k-ω 模型计算的 nut 一致；
• 步骤 5a 的 nuTilda ≈ nut，代入 SA 模型后 fv1 ≈ 1，nut 结果与 k-ε / k-ω 一致，说明整套公式自洽；
• 工程简化法（5b/5c）给出的 nuTilda 明显小于精确值（5a），属于保守初始化，收敛后会自动调整。

”

6. 反向校核（已知 k 求 I）

若已知 k 与 U，可反算湍流强度（与输入 I 相互校验）：

通常入口处 |U| 远大于 √(2k/3)，分母取 |U|，公式退化为 I = √(2k/3)/|U|。

”

注：可以根据计算流程快速编制程序。程序下载地址：https://t.zsxq.com/DqCH7

（完）