Realizable k-ε模型中的'Realizable'是指'可实现'的意思吗？具体是指什么可实现？

该模型由Shih等人（1995）提出。'Realizable'意味着模型满足雷诺应力张量的物理可实现性条件。具体来说，是指法向雷诺应力非负（$\overline{u_\alpha'^2} \geq 0$），并且满足施瓦茨不等式（$\overline{u_\alpha' u_\beta'}^2 \leq \overline{u_\alpha'^2}\cdot\overline{u_\beta'^2}$）。

标准k-ε模型会违反这个条件吗？

会的。在标准k-ε模型中，由于$C_\mu = 0.09$是常数，在高应变率区域（例如旋流中心），法向应力可能出现负值。这在物理上是不可能的。

其最大的特点就是$C_\mu$变成了变量，对吗？

没错。在Realizable模型中，涡粘系数由以下公式定义： $$ \mu_t = \rho C_\mu \frac{k^2}{\varepsilon}, \quad C_\mu = \frac{1}{A_0 + A_s \frac{kU^*}{\varepsilon}} $$

可实现k-ε模型

分类: 流体解析（CFD） | 综合版 2026-04-06

CAE visualization for k epsilon realizable theory - technical simulation diagram

Realizable k-εモデル

理论与物理

概述

🧑‍🎓

Realizable k-ε模型中的"Realizable"是"可实现"的意思吗？具体是什么"可实现"呢？

🎓

这是由Shih等人（1995）提出的模型。"Realizable"意味着满足雷诺应力张量的物理一致性条件（realizability constraints）。具体来说，是指法向雷诺应力非负（$\overline{u_\alpha'^2} \geq 0$），以及满足施瓦茨不等式（$\overline{u_\alpha' u_\beta'}^2 \leq \overline{u_\alpha'^2}\cdot\overline{u_\beta'^2}$）。

🧑‍🎓

标准k-ε模型有时会违反这个条件吗？

🎓

是的。在标准k-ε模型中，由于 $C_\mu = 0.09$ 是常数，在应变率非常大的区域（例如旋流的中心），法向应力有时会变为负值。这在物理上是不可能的。

输运方程

🎓

$k$ 方程与标准k-ε相同，但 $\varepsilon$ 方程有很大不同。

$$ \frac{\partial(\rho k)}{\partial t}+\frac{\partial(\rho U_j k)}{\partial x_j}=\frac{\partial}{\partial x_j}\left[\left(\mu+\frac{\mu_t}{\sigma_k}\right)\frac{\partial k}{\partial x_j}\right]+P_k-\rho\varepsilon $$

$$ \frac{\partial(\rho\varepsilon)}{\partial t}+\frac{\partial(\rho U_j\varepsilon)}{\partial x_j}=\frac{\partial}{\partial x_j}\left[\left(\mu+\frac{\mu_t}{\sigma_\varepsilon}\right)\frac{\partial\varepsilon}{\partial x_j}\right]+\rho C_1 S\varepsilon - \rho C_2\frac{\varepsilon^2}{k+\sqrt{\nu\varepsilon}} $$

🎓

这里重要的是 $C_1$ 的定义。

$$ C_1 = \max\left(0.43, \frac{\eta}{\eta+5}\right), \quad \eta = S\frac{k}{\varepsilon} $$

可变 $C_\mu$

🧑‍🎓

$C_\mu$ 变成变量是最大的特点吧？

🎓

没错。在Realizable模型中，涡粘性定义如下。

$$ \mu_t = \rho C_\mu \frac{k^2}{\varepsilon}, \quad C_\mu = \frac{1}{A_0 + A_s \frac{kU^*}{\varepsilon}} $$

🎓

这里 $U^* = \sqrt{S_{ij}S_{ij} + \tilde{\Omega}_{ij}\tilde{\Omega}_{ij}}$，$A_0 = 4.04$，$A_s = \sqrt{6}\cos\phi$（$\phi$ 由应变率张量的第三不变量计算得出）。

🧑‍🎓

也就是说，当应变率变大时，$C_\mu$ 会自动变小，对吗？

🎓

正是如此。在旋流中心，$S$ 很大，因此 $C_\mu$ 减小，抑制了湍流动能的过度预测。在静止流（$S k/\varepsilon \to 0$）中，$C_\mu \to 1/A_0 \approx 0.25$，但在通常的剪切流中，其值接近约0.09，与标准k-ε模型保持一致。

模型常数

常数	值
$C_2$	1.9
$\sigma_k$	1.0
$\sigma_\varepsilon$	1.2
$A_0$	4.04

🧑‍🎓

Realizable k-ε模型有什么弱点吗？

🎓

由于 $\varepsilon$ 方程的源项包含 $\sqrt{\nu\varepsilon}$，即使 $\varepsilon \to 0$ 也能避免奇异性。但是，在使用多个旋转参考系的多区域计算中，有时会产生非物理的湍流粘度。Fluent中已通过修正版（用于滑移网格的修正）对此进行了处理。

Coffee Break 闲谈

"Realizable"一词的含义——守护"理所当然"的艰辛

Realizable k-ε中的"Realizable（可实现）"这个名字，意味着满足湍流的数学约束——法向应力的非负性和施瓦茨不等式。实际上，标准k-ε并不总是满足这些约束，在强应变流中可能会预测出非物理的负法向应力。Shih等人（1994）展示了"仅仅通过遵守这些理所当然的物理约束，模型就能得到如此大的改进"。Realizable k-ε之所以成为k-ε家族在压缩、膨胀、分离流中的标准选择，原因就在于这个"理所当然的保证"。

各项的物理意义

时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才会变成稳定的水流，对吧？描述这个"变化过程中"的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭导致流动变化，这些都是非定常现象。那么定常分析是什么？就是只看"经过足够长时间、流动稳定之后"的状态——也就是将此项设为零。由于计算成本大幅降低，先用定常分析求解是CFD的基本策略。
对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？它会随水流被带到下游，对吧？这就是"对流"——流体运动携带物质的效果。暖风的暖气能到达房间的另一端，也是因为空气这个"搬运工"通过对流输送热量。这里有趣的是——这项包含"速度×速度"，因此是非线性的。也就是说，流速变快时，这项会急剧增强，变得难以控制。这就是湍流的根本原因。一个常见的误解："对流和传导差不多"→ 完全不一样！对流是流动携带，传导是分子传递。两者效率有天壤之别。
扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘性（$\mu$）高，所以不易流动。粘性越大，扩散项越强，流体的运动就变得"粘稠"。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，在Re数大的流动中，对流占压倒性优势，扩散则成为配角。
压力项 $-\nabla p$：按下注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差产生了推动流体的力。大坝放水也是同样的原理。天气图中等压线密集的地方会怎样？没错，会刮强风。"有压力差的地方就会产生流动"——这就是纳维-斯托克斯方程压力项的物理意义。这里容易误解的点是：CFD中的"压力"通常指表压而非绝对压力。切换到可压缩分析时，如果结果变得奇怪，原因可能就是混淆了绝对压力/表压。
源项 $S_\phi$：受热的空气会上升——为什么呢？因为变得比周围轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶的火焰产生化学反应热，工厂的电磁泵对金属熔液施加洛伦兹力……这些都是"从外部向流体注入能量或力"的作用，都用源项来表示。如果忘记源项会怎样？在自然对流分析中忘记加入浮力，流体就完全不动——就像冬天在房间里开了暖气，暖空气却不上浮一样，会得到物理上不可能的结果。

假设条件与适用范围

连续介质假设：在克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
牛顿流体假设：剪切应力与应变率呈线性关系（非牛顿流体需要粘度模型）
不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项中使用恒定密度
不适用的情况：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）

量纲分析与单位制

变量	SI单位	注意事项·换算备忘
速度 $u$	m/s	从入口条件的体积流量换算时，注意截面积的单位
压力 $p$	Pa	区分表压和绝对压力。可压缩分析中使用绝对压力
密度 $\rho$	kg/m³	空气: 约1.225 kg/m³@20°C，水: 约998 kg/m³@20°C
粘性系数 $\mu$	Pa·s	注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆
雷诺数 $Re$	无量纲	$Re = \rho u L / \mu$。层流/湍流转换的判断指标
CFL数	无量纲	$CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性

数值解法与实现

离散化的注意事项

🧑‍🎓

Realizable k-ε在数值处理上与标准k-ε有什么不同吗？

🎓

$\varepsilon$ 方程的形式不同，因此需要注意源项的线性化。特别是分母 $k + \sqrt{\nu\varepsilon}$ 很重要，其设计确保了即使在 $k \to 0$ 的区域（壁面近壁区或非湍流区域）也不会发生除零错误。

🎓

离散化方案的推荐与标准k-ε相同。

项	推荐方案
对流项	Second Order Upwind
扩散项	Central Differencing
时间项	Second Order Implicit

$C_\mu$ 的数值计算

🧑‍🎓

$C_\mu$ 依赖于流场，这意味着每次迭代都要重新计算吗？

🎓

是的。在每个单元计算 $S_{ij}$ 和 $\Omega_{ij}$，算出 $U^*$ 和 $\phi$ 后更新 $C_\mu$。计算成本比标准k-ε略有增加，但由于不需要求解额外的输运方程，所以差异不大。

Ansys Fluent

🎓

```

Models → Viscous → k-epsilon → Realizable

Near-Wall Treatment → Enhanced Wall Treatment（推荐）

```

在Fluent中，Realizable已成为k-ε的默认推荐模型。与Enhanced Wall Treatment结合使用时，对 $y^+$ 的限制会放宽。

OpenFOAM

🎓

在 constant/turbulenceProperties 中指定如下。

```

RAS

{

RASModel realizableKE;

turbulence on;

printCoeffs on;

}

```

STAR-CCM+

🎓

选择 K-Epsilon Turbulence → Realizable K-Epsilon Two-Layer。Two-Layer模型通常与All y+ Wall Treatment结合使用。

收敛性比较

🧑‍🎓

与标准k-ε相比，更难收敛吗？

🎓

一般来说，Realizable k-ε的收敛性更好。这是因为 $C_\mu$ 可变，抑制了非物理的 $\mu_t$ 爆发。不过，在最初的几十次迭代中，如果 $C_\mu$ 变化剧烈，可能会变得不稳定。这种情况下，将湍流粘度比的URF降低到0.8左右会比较好。

Coffee Break 闲谈

Realizable k-ε成为汽车外气动分析主角的经过

2000年代以后，Realizable k-ε在汽车制造商的CFD部门中作为外气动分析的标准模型普及开来。原因之一是其对A柱附近的分离/再附着以及后视镜尾流的再现性比标准k-ε更高。有多个案例报告称，某OEM设计团队使用相同网格和设置比较了4个模型，结果Realizable k-ε的阻力系数与风洞实验值的差异最小，因此被后续的设计周期采用。即使在今天，ANSYS Fluent入门研讨会上也经常教导"外气动分析就从Realizable k-ε开始"。

迎风格式（Upwind）

一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。

中心差分（Central Differencing）

二阶精度，但当Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。

TVD格式（MUSCL、QUICK等）

通过限制器函数抑制数值振荡，同时保持高精度。对捕捉激波或陡峭梯度有效。

有限体积法 vs 有限元法

FVM：自然满足守恒定律。CFD的主流。FEM：对复杂形状/多物理场有利。SPH等无网格法也在发展中。

CFL条件（库朗数）

显式方法：CFL ≤ 1为稳定条件。隐式方法：即使CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：一个时间步内信息传播不超过一个网格。