komegaSST湍流模型中给动量方程雷诺应力增加非线性项

[fangyuanaza]

各位有经验的老师好，我最近在修改rhoSimpleFoam的动量方程，目标是将其中的动量方程(UEqn.H)中turbulence->divDevRhoReff(U)在komegaSSTBase中更新，添加非线性项（红框部分） 。目前遇到的问题是在kOmegaSSTBase.C文件的constructors中，如何能同时构造Turbulence 与 eddyViscosity？也就是说在

冒号后面如何同时构造
与

可以看到内部量是一致的，因此我尝试采用如下形式

但是会报错，我使用的是openfoam v7。才开始尝试修改openfoam内部代码，请大家多多指教！

[李东岳]

没看明白，你直接改linearViscousStress.C不就好了么？

[fangyuanaza]

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？