代码实现的问题
-
大家好，我最近在修改OpenFoam里的雷诺应力模型。因为在对比OpenFOAM和Fluent之间的湍流模型的区别时，发现Fluent里面的快速压力应变项比Openfoam里多了由对流引起的一项C。
其中
所以想添加
这一项到代码里。但是改了很多种表达都通不过编译，比如说定义
volTensorField C(fvm:: div(alphaRhoPhi,R));
然后在方程里使用 dev(C)
都不行，于是想问问大家有什么代码思路。编程小白，真诚求助 -
啊，我现在改成这样了：
外层定义volSymmTensorField devR(dev(R));
R方程里面用 fvm::div(alphaRhoPhi,devR)
编译可以通过，而且貌似是可以代表上述表达的，不知道这样对不对？ -
更新：使用fvm的时候可以编译但计算不行，改成fvc 可以计算接下来就是看能不能收敛了