代码实现的问题



  • 大家好,我最近在修改OpenFoam里的雷诺应力模型。因为在对比OpenFOAM和Fluent之间的湍流模型的区别时,发现Fluent里面的快速压力应变项比Openfoam里多了由对流引起的一项C。
    phi12.png
    其中CodeCogsEqn-6.gif
    CodeCogsEqn-7.gif
    所以想添加CodeCogsEqn-8.gif
    这一项到代码里。但是改了很多种表达都通不过编译,比如说定义
    volTensorField C(fvm:: div(alphaRhoPhi,R));
    然后在方程里使用 dev(C)
    都不行,于是想问问大家有什么代码思路。编程小白,真诚求助



  • 啊,我现在改成这样了:
    外层定义volSymmTensorField devR(dev(R));
    R方程里面用 fvm::div(alphaRhoPhi,devR)
    编译可以通过,而且貌似是可以代表上述表达的,不知道这样对不对?



  • 更新:使用fvm的时候可以编译但计算不行,改成fvc 可以计算接下来就是看能不能收敛了


Log in to reply