OpenFOAM中MRF的代码


  • OpenFOAM副教授

    MRF方法,理论上应该是在相对参考系中求解绝对速度,按照OpenFOAM wiki:http://openfoamwiki.net/index.php/See_the_MRF_development 的思路,其公式应该是:0_1510037016871_MRF.PNG , 对应的OpenFOAM代码为:

    tmp<fvVectorMatrix> tUEqn
    (
        fvm::ddt(U) + fvm::div(phi, U)
      + MRF.DDt(U)
      + turbulence->divDevReff(U)
     ==
        fvOptions(U)
    );
    

    代码实现中,MRF.DDt(U)即为omega^ UI,注意此处是显式处理,作为动量方程的源项,而通过连续性方程组建压力泊松方程时,phiHbyA只考虑了相对通量omege^r,即MRF.makeRelative(phiHbyA);;但并没有考虑omega^ UI这个源项,有人注意过这段代码实现吗?是不是这里的方程错了?


  • OpenFOAM副教授

    @dzw05 但也有可能是求HbyA时将omega^UI作为源项考虑进去了,这段代码还没有看懂,有大神指点一二吗?