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



  • 各位有经验的老师好,我最近在修改rhoSimpleFoam的动量方程,目标是将其中的动量方程(UEqn.H)中turbulence->divDevRhoReff(U)在komegaSSTBase中更新,添加非线性项(红框部分)9cb1e563-2c83-48bf-9386-79bebfe65989-image.png 。目前遇到的问题是在kOmegaSSTBase.C文件的constructors中,如何能同时构造Turbulence 与 eddyViscosity?也就是说在
    47b58e64-4329-4a55-bde4-acb9fbce3a93-image.png
    冒号后面如何同时构造
    ee75edfc-029a-4fb0-988c-47b85228b336-image.png
    2909fe55-4399-4140-ae0f-b8fd91eea9e8-image.png
    可以看到内部量是一致的,因此我尝试采用如下形式
    325ef48c-17b7-4133-aa7c-8539a34ef5d1-image.png
    但是会报错,我使用的是openfoam v7。才开始尝试修改openfoam内部代码,请大家多多指教!


  • 管理员

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



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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1