📚 课程目录第4章 k-ε模型家族

4. k-ε模型家族

k-ε模型家族是工程湍流计算中使用最广的一类两方程模型。它们都以湍动能 k 及其耗散率 ε 作为求解变量，通过涡黏假设把雷诺应力与平均速度梯度联系起来，因此计算量适中、鲁棒性较好。但标准、RNG 与 Realizable 三个版本在方程构造、可实现性约束和壁面处理上存在实质差异，这些差异直接决定了它们对分离流、旋流和强曲率流动的预测能力。

一、三个模型的方程差异

标准 k-ε 模型由 Launder 与 Spalding 提出，k 方程与 ε 方程均为经验输运方程。其涡黏系数按 μ_t = ρC_μk²/ε 计算，C_μ 取常数 0.09。该模型在 ε 方程中采用 C_1ε=1.44、C_2ε=1.92、σ_k=1.0、σ_ε=1.3。它的主要问题在于：当平均应变率很大时，C_μ 为常数会导致雷诺正应力出现负值，即违反“可实现性”条件，从而在强旋流和强分离区产生非物理结果。

RNG k-ε 模型通过重整化群理论对 Navier-Stokes 方程进行系统粗粒化，在 ε 方程中引入额外的应变率相关项 R_ε，使耗散率在强应变区被抑制。其涡黏系数仍为 μ_t=ρC_μk²/ε，但 C_μ 在 Fluent 中默认取 0.0845，且 ε 方程系数由理论推导得到：C_1ε=1.42、C_2ε=1.68。RNG 模型还提供旋流修正选项，用于改善强旋流预测。

Realizable k-ε 模型的核心改动是让 C_μ 不再为常数，而是与平均应变率和旋转率关联：C_μ=1/(A₀+A_skU*/ε)。同时 ε 方程不再由 k 方程乘以固定系数得到，而是从涡量脉动均方方程出发重新推导，使模型在强应变、强旋转和分离流动中更符合物理约束。其常用系数为 C_1ε=1.43、C_2ε=1.9、σ_k=1.0、σ_ε=1.2。

二、壁面函数的选择

k-ε 模型属于高雷诺数模型，在近壁黏性底层内不直接求解，必须借助壁面函数把壁面剪切应力与第一层网格节点的速度、湍流量联系起来。Fluent 中常用的壁面处理包括标准壁面函数、可缩放壁面函数、增强壁面处理和非平衡壁面函数。

标准壁面函数基于对数律，要求第一层网格节点落在对数律层内，通常 y+ 在 30～300 之间。若 y+ 过小（如小于 15），节点落入黏性底层，对数律失效，壁面剪切应力会被高估；若 y+ 过大（如超过 300），节点进入尾流区，同样产生偏差。可缩放壁面函数通过引入 y+ 相关的混合函数，在 y+ 较小时自动过渡到低雷诺数形式，对 y+ 的敏感性低于标准壁面函数，适合网格 y+ 分布不均匀的复杂几何。

增强壁面处理采用两层模型，在近壁区直接求解到壁面，要求 y+≈1，能较准确地捕捉边界层分离和壁面热传导，但网格量显著增加。非平衡壁面函数考虑压力梯度和非平衡效应，适合分离流和驻点附近流动，但对 y+ 的要求仍较严格。工程上若以分离流或旋流为主，建议优先采用增强壁面处理并保证 y+≈1；若计算资源有限，可缩放壁面函数是较稳妥的折中。

三、对分离流与旋流的预测能力对比

分离流和旋流是 k-ε 模型家族差异最明显的两类流动。标准 k-ε 在逆压梯度下倾向于高估湍动能，使分离点推迟、分离区偏小；在强旋流中，由于 C_μ 为常数，涡黏被高估，旋流核心的切向速度衰减过快，回流区尺寸偏小。RNG 模型通过 ε 方程中的应变率附加项抑制强应变区的耗散，使分离区预测有所改善，开启旋流修正后对中等旋流更准确。Realizable 模型由于 C_μ 随应变率和旋转率变化，在强旋流和强分离中能保持雷诺正应力非负，对旋流核心和分离泡的预测通常最接近实验。

以一个典型旋风分离器为例：入口速度 15 m/s，筒体直径 300 mm，采用三种模型分别计算。标准 k-ε 得到的旋流数衰减最快，中心回流区长度约 1.2 倍筒径；RNG 模型约 1.5 倍筒径；Realizable 模型约 1.8 倍筒径，与激光多普勒测速实验趋势更一致。该结果属于教学示例，实际数值会随网格、离散格式和边界条件变化。

选型建议：若流动以弱分离和简单剪切为主，标准 k-ε 即可满足；若存在中等旋流或分离，优先选 RNG 并开启旋流修正；若旋流强、分离区大或对回流区尺寸敏感，优先选 Realizable k-ε，并配合增强壁面处理与 y+≈1 的网格。

四、Fluent 中的关键设置与检查

在 Fluent 中选择 k-ε 模型时，需在 Viscous 面板中指定 Standard、RNG 或 Realizable，并选择对应的壁面处理。RNG 模型可勾选 Swirl Dominated Flow 选项；Realizable 模型默认采用标准壁面函数，可改为 Enhanced Wall Treatment。近壁网格检查可通过后处理中的 Wall Yplus 云图完成，确保目标壁面 y+ 落在所选壁面函数的要求范围内。

# Fluent TUI 示例：设置 Realizable k-ε + 增强壁面处理
define/models/viscous/ke-realizable yes
define/models/viscous/near-wall-treatment/enhanced-wall-treatment yes
solve/set/expert
# 建议开启高阶离散以降低数值耗散
solve/set/discretization-scheme/momentum 1
solve/set/discretization-scheme/k 1
solve/set/discretization-scheme/epsilon 1

收敛后应重点检查三项：残差是否降至 1e-4 以下并保持平稳；壁面 y+ 云图是否落在目标区间；分离区或旋流核心的流场形态是否与物理预期一致。若 y+ 大面积偏离，应先调整边界层网格，而不是反复修改模型系数。

对比项	标准 k-ε	RNG k-ε	Realizable k-ε
Cμ	0.09 常数	0.0845 常数	随应变率/旋转率变化
ε 方程特点	经验输运	含应变率附加项	由涡量脉动方程推导
可实现性	不满足	部分改善	满足
旋流预测	偏弱	中等	较强
分离流预测	分离区偏小	优于标准	较接近实验
推荐壁面处理	标准壁面函数	标准/增强	增强/可缩放

总体而言，三个模型并非简单的“越新越好”，而是适用场景不同。标准 k-ε 计算最稳定、成本最低，适合外流和弱分离；RNG 在中等旋流和分离中提供更好的折中；Realizable 在强旋流、强分离和射流类流动中通常表现最好，但对网格和壁面处理的要求也更高。实际选型应结合流动特征、网格 y+ 分布和计算资源综合判断。