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



  • RANS中张量形式动量方程为:
    798afe64-8bcc-426b-80a0-a9086108325d-image.png
    在rhoSimpleFoam求解器中对应为UEqn.H
    368d0d2c-93e1-428a-9448-a00f28ef6d11-image.png
    由于压力梯度项一般不直接放入UEqn中,而是对动量方程进行修正,因此turbulence->divDevRhoReff(U) 对应07bd23db-8f4c-4d77-b1e5-afb046e5e712-image.png 前一项可用拉普拉斯符号表示,后一项实际上就是雷诺应力。
    问题一:我发现turbulence->divDevRhoReff在/home/dyfluid/OpenFOAM/OpenFOAM-7/src/TurbulenceModels/turbulenceModels/路径下的linearViscousStress/linearViscousStress.C及ReynoldsStress/ReynoldsStress.C中均有体现,如果我要修改动量方程中的雷诺应力项,那么我两个都需要修改?
    问题二:turbulence->divDevRhoReff对应的代码是
    335942ce-0263-449f-a98d-54fbf388ee5b-image.png 这里的第一项也就是雷诺应力项表达式为
    04392fb8-9aaf-4137-9af7-8349a57fee31-image.png
    这为什么与Boussinesq假定(如下)不一致呢?
    a574197b-04cf-4ac4-bbc2-2f7402831d04-image.png
    Boussinesq假定的第一项b2d0bd9f-1d34-4c49-8230-66db725d6e7e-image.png为什么在代码中没有体现呢?



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



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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2