k-ε模型 (k-epsilon Model) — CAE术语解说
k-ε模型是什么
k-ε模型是什么？我听说在RANS中尤其常用。
$k$ 是湍流的能量，$\varepsilon$ 是它消散的速度……是这样吗？
完全正确。$k = \frac{1}{2}\overline{u_i' u_i'}$，是湍流脉动速度平方平均值的一半。而 $\varepsilon$ 是该 $k$ 被分子粘性转化为热的速率。简单来说，$k$ 是"湍流的功率"，$\varepsilon$ 是"功率的衰减速率"。有了这两个量，就可以计算涡粘性系数 $\mu_t = \rho C_\mu \frac{k^2}{\varepsilon}$，从而用来近似湍流对动量输运的影响。
比如计算建筑空调风管内的流动时，就是这样把湍流的影响通过 $\mu_t$ 包含进去吧？
正是这样。风管、管路、热交换器、汽车外形空气动力学——工业中大部分内流和外流都是先尝试用k-ε模型计算。计算稳定性好，有丰富的实验验证数据。
输运方程
k-ε模型的方程是什么样的？我想看看数学表达式。
Standard k-ε模型的输运方程是这样的。首先是 $k$ 的方程：
左边是时间变化加对流，右边是扩散加生成减散逸。$P_k = \mu_t S^2$ 是平均流速度梯度产生湍流能的项。
$\varepsilon$ 的方程也是差不多的形式吗？
结构相同但系数不同：
Standard k-ε的模型常数为 $C_\mu = 0.09$、$C_{\varepsilon 1} = 1.44$、$C_{\varepsilon 2} = 1.92$、$\sigma_k = 1.0$、$\sigma_\varepsilon = 1.3$。这些值是根据简单剪切流实验数据拟合出来的。
有5个常数啊。用户可以改这些吗？
一般保持默认值。只有在燃烧或喷雾这样特殊流动中，才会微调 $C_{\varepsilon 1}$ 的研究。但通常的做法是用默认值先计算，如果与实验不符就换成其他模型（RNG、Realizable、k-ω SST等）。这是最稳妥的方法。
三种变体 — Standard、RNG、Realizable
除了Standard k-ε，还听说过RNG和Realizable，这些有什么区别？
主要有三种变体。我来整理一下各自的特点。
Standard k-ε（Launder & Spalding, 1974）：最简单，收敛性好。常数来自实验拟合。但对强旋流或大分离流动精度会下降。
RNG k-ε（Yakhot & Orszag, 1986）：使用重整化群（Renormalization Group）理论来解析推导常数。ε方程中增加了修正项 $R_\varepsilon$。这一项在速度梯度急剧变化的区域修正ε的生成，使得在中等雷诺数和轻微旋流下比Standard精度更好。
Realizable k-ε（Shih et al., 1995）：为了满足"可实现性条件"（雷诺应力张量的特征值非负），将 $C_\mu$ 从常数改为随流场变化的变量。还重新推导了ε方程。对喷流（射流）膨胀角的预测有大幅改进，对含分离的流动也更强。
工程实践中应该用哪个？因为不同求解器的默认值好像不同。
ANSYS Fluent默认推荐Realizable k-ε。OpenFOAM提供kEpsilon（Standard）、kEpsilonRNG、realizableKE可选。STAR-CCM+则以Realizable Two-Layer k-ε较受欢迎。工程上，没有特殊原因的话选Realizable就不会出大问题。Standard之所以还保留是为了兼容性，新项目很少主动用它。
比如要计算工厂通风管道的压力损失，应该选哪一个？
直管风道用Standard往往精度就足够。但如果有弯管或T形分岔，会产生分离，此时用Realizable比较保险。计算成本和Standard几乎相同，用Realizable没有坏处。
壁面函数与壁面处理
听说k-ε模型需要"壁面函数"，这是做什么用的？
壁面附近有一个叫"粘性底层"的薄层，那里是分子粘性而非湍流占主导。k-ε模型直接求解这个区域很困难，因为ε在壁面的边界条件具有奇异性。
所以，把壁面第一个网格单元的中心放在粘性底层外侧（对数律层，$y^+ \approx 30$～$300$）的相对粗网格位置，通过对数律来关联壁面的剪应力 $\tau_w$ 和壁面附近的速度。这就是壁面函数，可以大幅减少所需的网格数。
$y^+$ 是什么？做网格时总要注意它？
$y^+$ 是从壁面的无量纲距离，定义为 $y^+ = \frac{\rho u_\tau y}{\mu}$。其中 $u_\tau = \sqrt{\tau_w / \rho}$ 是摩擦速度。
使用壁面函数时应该瞄准 $y^+ \approx 30$～$300$。相反，如果要直接求解到壁面（Low-Reynolds数模型或Enhanced Wall Treatment），则需要 $y^+ \approx 1$。
如果 $y^+$ 落在5～25这个中间范围怎么办？
好问题。这叫"缓冲层"，对数律和粘性律都不完全适用。壁面第一层网格落在这里精度会大幅下降。现在新的求解器有"Enhanced Wall Treatment"（Fluent）或"All-y+ Wall Treatment"（STAR-CCM+）这样的功能，能自动根据 $y^+$ 值进行混合处理。配合k-ε模型使用这个功能是最佳实践。
适用范围和局限
k-ε模型有什么缺点或不应该用的情况吗？
主要有三大弱点。
1. 强旋流： 如旋流分离器或涡轮叶片通道这样强旋流存在的情况，k-ε模型容易过度评估涡粘性。旋流衰减往往预测得太快。
2. 大规模分离： 建筑物后的尾流或急扩管这样大规模分离支配的流动，分离点位置和再附着长度的预测精度较差。
3. 强非各向同性流动： k-ε模型基于涡粘性假设（Boussinesq假设），假定雷诺应力张量是各向同性的。应力非各向同性很重要的流动，如管道拐角处的二次流，用雷诺应力模型（RSM）更准确。
那汽车车身周围的空气动力学模拟用k-ε可以吗？背面会形成大尾流。
这恰好是"k-ε不擅长的问题"的典型。汽车车身背面形成大规模分离涡，阻力大部分来自此。所以汽车厂商的CFD团队往往采用k-ω SST或DES（Detached Eddy Simulation）。但是发动机舱内的冷却气流分配这类内流，k-ε就足够了。要根据问题性质来选择。
k-ε和k-ω SST怎么取舍，有什么规律吗？
粗略规律是这样：
- k-ε适合： 充分发展的内流（管路、风道）、自由剪切流（喷流、混合层）、初期研究需要快速出结果
- k-ω SST适合： 受逆压力梯度的边界层、有分离和再附着的外流、翼型和车体空气动力学、热传递预测
犹豫不决的话k-ω SST是更稳妥的选择，但k-ε收敛性好，在复杂几何形状的"先跑出来试试"阶段很实用。
工程实践中的选择指南
总结一下，用k-ε模型的正确步骤是什么？
工程流程是这样的：
- 掌握雷诺数 → 判断是否需要湍流模型
- 分类流动性质 → 内流还是外流？分离有多大？有旋流吗？
- 用Realizable k-ε + Enhanced Wall Treatment初始计算 → 验证收敛性和稳定性
- 验证结果 → 和实验或基准数据对比，精度不足就换k-ω SST或RSM
- 检查网格相关性 → 确认壁面 $y^+$ 值，确保在壁面函数适用范围内
按这个流程走，就既不会用高成本模型去做k-ε足够的问题，也不会为了精度不达标而苦恼。
在OpenFOAM中用k-ε模型时，入口的 $k$ 和 $\varepsilon$ 初始值怎么设？我总是随便填结果被吐槽……
入口边界条件的常用推算公式是这样：
$I$ 是湍流强度，管路取5%～10%，外部自由流取1%左右。$l$ 是湍流长度尺度，管路常用 $l \approx 0.07D$（$D$ 是管径）。OpenFOAM有 turbulentIntensityKineticEnergyInlet 和 turbulentMixingLengthDissipationRateInlet 边界条件可自动计算，用这个最方便。
相关术语
- k（湍流能）：湍流脉动速度的运动能。$k = \frac{1}{2}\overline{u_i' u_i'}$
- ε（散逸率）：湍流能被分子粘性散逸的速率
- 涡粘性（eddy viscosity）：湍流导致的表观粘性。$\mu_t = \rho C_\mu k^2 / \varepsilon$
- 壁面函数（wall function）：用对数律近似壁面附近速度和温度分布的方法
- $y^+$：从壁面的无量纲距离。判断壁面函数是否可用的指标
- Boussinesq假设：用涡粘性和平均应变速率张量的乘积来近似雷诺应力的假设
对CAE术语的准确理解是团队沟通的基础。 — Project NovaSolver也致力于工程师的学习支持。
在k-ε模型实务中遇到的问题，欢迎告诉我们
Project NovaSolver致力于解决CAE工程师日常面临的挑战——建模的繁琐、计算成本、结果的解释。你的实务经验将成为更好工具开发的驱动力。
联系方式（筹备中）相关主题
价值
详细
错误