openfoam6.0版本中标准k-e模型的Gk项发生了变化？

[羽之下]

在之前版本中：
volScalarField G(GName(), nut_2magSqr(symm(fvc::grad(U_))));
这个可以对应上fluent帮助文档上的相关公式：

和《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®》的：

这两个是一样的，但是OF6.0表示为：

volScalarField::Internal G
     (
         this->GName(),
         nut.v()*(dev(twoSymm(tgradU().v())) && tgradU().v())
     );

好像跟之前的表达式并不一样啊，是我理解错了吗？

[羽之下]

@羽之下 已解决。

[wwzhao]

@羽之下 是一样的，一个是不可压缩，一个是通用形式。

[shampoo]

您好，我最近刚开始研究k-e模型的代码，想请教一下这段代码里的tgradU().v()是代表了一个什么操作呀？
我把这段代码变成公式之后，发现openFOAM里现在的G展开比fluent里多了一项
openfoam里的G是这样的：

按定义 化成

多了一项
我想知道这样的展开对不对，多的这项会有什么影响吗？

[李东岳]

@shampoo 没太对比你的公式，不过在不可压缩的时候是0

[shampoo]

@东岳 好的，我再研究一下

[羽之下]

@shampoo Fluent中给出的最后公式是针对不可压缩流。较新的Openfoam版本都将可压缩和不可压缩糅合在一起了，注意G项和P项的差异，《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®》有较为详细的推导过程，最后结论也是针对不可压缩流。可以参考下这个https://zhuanlan.zhihu.com/p/54802989，关于湍流模型这块，我也还有好多东西没搞懂，接着在看，希望共同进步！

[shampoo]

@羽之下 thank you，谢谢你的回复。我记得好像在哪看过说这里引用了一个假设，我决定把这部分系统的看一看，有什么新发现我会继续在这里交流~