RANS动量方程中雷诺应力表达的问题

[fangyuanaza]

RANS中张量形式动量方程为：

在rhoSimpleFoam求解器中对应为UEqn.H

由于压力梯度项一般不直接放入UEqn中，而是对动量方程进行修正，因此turbulence->divDevRhoReff(U) 对应 前一项可用拉普拉斯符号表示，后一项实际上就是雷诺应力。
问题一：我发现turbulence->divDevRhoReff在/home/dyfluid/OpenFOAM/OpenFOAM-7/src/TurbulenceModels/turbulenceModels/路径下的linearViscousStress/linearViscousStress.C及ReynoldsStress/ReynoldsStress.C中均有体现，如果我要修改动量方程中的雷诺应力项，那么我两个都需要修改？
问题二：turbulence->divDevRhoReff对应的代码是
这里的第一项也就是雷诺应力项表达式为

这为什么与Boussinesq假定（如下）不一致呢？

Boussinesq假定的第一项为什么在代码中没有体现呢？

[fangyuanaza]

@fangyuanaza https://caefn.com/openfoam/diffusion-term1 我找到了这方面的解析，说是这一项放到了压力梯度中，但是我还是没有找到

[东岳]

没有显性的操作，只不过是存在压力p里面，但是没有显性的说p=p+1/3k之类。所以，其实求解出来的p，是p+1/3k

[fangyuanaza]

@东岳 懂了~ 谢谢智慧的李老师~~