Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    icoFoam细节问题讨论

    OpenFOAM
    2
    2
    1472
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • M
      mushtime last edited by CFD中文网

      在icoFoam解析中,提到压力场修正时,

      “ phi = phiHbyA - pEqn.flux();// 使用求解的压力场修正通量场,在最后一次修正的时候通量守恒,Issa指出,大约需要2-3次内循环步。对应公式(25),pEqn.flux()返回公式(25)方程右边第二项,也为fvc::interpolate(rUA)*fvc::snGrad(p)*mag(mesh.Sf())。某些可压缩求解器其中的pEqn.flux()可能为+号,即为phi = phiHbyA + pEqn.flux()。这是因为pEqn中的laplacian项为−号”

      这其中,提到公式25的右边第二项也叫fvc::interpolate(rUA)*fvc::snGrad(p)*mag(mesh.Sf())。可是公式25右边第二项体积分之后,再通过高斯变换之后,就应该是p在cell面的值乘以一个系数以及mesh.Sf(),为什么这边解释是snGrad(p)。难道是我理解错了?

      1 Reply Last reply Reply Quote
      • C
        CFD中文网 last edited by CFD中文网

        公式(25)如下:
        0_1466491773445_捕获.JPG
        你说的:

        可是公式25右边第二项体积分之后,再通过高斯变换之后...

        为什么你还要对上述公式进行体积分?

        CFD中国标准用户测试帐号
        目前由徐笑笑登录

        1 Reply Last reply Reply Quote
        • First post
          Last post

        CFD中文网 | 东岳流体 | 京ICP备15017992号-2
        论坛登录问题反馈可联系 li.dy@dyfluid.com