关于icoFOAM中通量phi的更新



  • icoFOAM中更新通量phi采用的如下语句(注:我理解这个语句)

    if (piso.finalNonOrthogonalIter())
    {
        phi = phiHbyA - pEqn.flux();
    }
    

    为什么不在下面这条语句的后面,直接用phi = fvc::flux(U)来更新通量,这样更新跟上面那条语句有什么区别吗?用phi=fvc::flux(U)会有什么问题吗?

    U = HbyA - rAU*fvc::grad(p);
    U.correctBoundaryConditions();
    

    初学icoFOAM,求各位高手解答,OpenFOAM用的版本为OpenFOAM 7.


  • 管理员

    这里面phi是严格守恒的,U不存在守恒的概念。所以你用phi = fvc::flux(U),也可以,但是没有原来那样consistent



  • @李东岳 感谢李老师的回复,谢谢


Log in to reply
 


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