不可压缩湍流模型的 divDevReff 函数

OpenFOAM 中单相不可压缩求解器中，雷诺应力项调用的是湍流模型中的 divDevReff 函数。这个函数的返回值为：

• 1. 为什么是 1/3 而不是 2/3，对应到代码，即为什么用 dev 函数而不是 dev2 函数？

• 2. 根据涡粘度的 Boussinesq approximation，雷诺应力项：

中应该还包含湍动能 ，而 OpenFOAM 中的 divDefReff 函数是没有       这一项的。

所以，收敛以后， 这一项等于0。

但是在开始阶段，或者说还没有达到满足连续性的流场之前，这一项不为零。这里加上这一项是出于数值稳定性以及收敛速度的考虑，这一项不对收敛后的结果几乎没有影响。 所以1/3 或 2/3不是很重要。

但是，可压缩湍流模型里必须是 2/3，因为这个 是从 N-S 方程中严格推导而来的，而且，在可压缩的情形下，即使收敛以后，也有   。

在 OpenFOAM=3.0 以后的版本里，不可压和可压缩湍流模型纳入到一个框架下了，两种情形下，都是用 2/3      这个系数。

• 一种认为 的值相对很小，可以直接忽略不计。

• 另一种观点认为，         被放到了压力项里，即，动量方程中的压力是雷诺时均压力与雷诺应力的各向同性分量（即2/3k）之和：

不过，Tobias Holzmann 最后仍持前一种观点，即 k 项被忽略了。理由是OpenFOAM 中似乎找不到关于修改的压力场的代码，而且 OpenFOAM 那边也没见有人讨论说 OpenFOAM 中使用的是修改的压力场。