代码实现的问题

[shampoo]

大家好，我最近在修改OpenFoam里的雷诺应力模型。因为在对比OpenFOAM和Fluent之间的湍流模型的区别时，发现Fluent里面的快速压力应变项比Openfoam里多了由对流引起的一项C。

其中

所以想添加
这一项到代码里。但是改了很多种表达都通不过编译，比如说定义
volTensorField C(fvm:: div(alphaRhoPhi,R));
然后在方程里使用 dev(C)
都不行，于是想问问大家有什么代码思路。编程小白，真诚求助

[shampoo]

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

[shampoo]

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