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



  • 在之前版本中:
    volScalarField G(GName(), nut_2magSqr(symm(fvc::grad(U_))));
    这个可以对应上fluent帮助文档上的相关公式:
    6dc47e8d-9e00-4f39-93c7-5ba4a7e391a4-image.png
    和《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®》的:
    d3a6796b-994b-422f-a963-2ebc85ad5926-image.png
    这两个是一样的,但是OF6.0表示为:

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

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



  • @羽之下 已解决。


  • OpenFOAM教授

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



  • 您好,我最近刚开始研究k-e模型的代码,想请教一下这段代码里的tgradU().v()是代表了一个什么操作呀?
    我把这段代码变成公式之后,发现openFOAM里现在的G展开比fluent里多了一项
    openfoam里的G是这样的:
    CodeCogsEqn-3.gif
    按定义CodeCogsEqn-5.gif 化成
    CodeCogsEqn-4.gif
    多了一项CodeCogsEqn.gif
    我想知道这样的展开对不对,多的这项会有什么影响吗?



  • @shampoo 没太对比你的公式,不过替代文字在不可压缩的时候是0



  • @东岳 好的,我再研究一下



  • @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,关于湍流模型这块,我也还有好多东西没搞懂,接着在看,希望共同进步!



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


Log in to reply