icoFoam中压力泊松方程的求解



  • 压力泊松方程的离散问题

    of中构建的压力泊松方程似乎是这样的形式:
    $\nabla\cdot(\frac{1}{a_P}\nabla p)=\nabla\cdot(HbyA)$
    (icoFoam解析中的式27、33)

    按我的理解,对上面的方程似乎应该写成下面这种离散形式?

    fvm::laplacian(rAU, p) == fvc::div(HbyA)
    

    为什么实际上icoFoam使用了phiHbyA作为中间变量呢?

    phiHybA的问题:

    phiHybA是一个surfaceScalarField 对象,即面上标量场;随后求解器执行了对这个面上标量的散度:

    fvc::div(phiHbyA)
    

    为什么可以对一个面上标量计算散度呢?是否div函数在对surfaceScalarField作用的时候采用了不同的方式?我查了div的代码,确实看到div对surfaceScalarField似乎有不同的处理,但是我没看懂。请问这是否有关系呢?

    期待各位大佬指点:135:


  • 管理员

    这个问题年年都有人讨论,是个经典问题,这样,我把这个问题更新到NS方程笔记里面,更新完告诉你:duang:



  • @李东岳 多谢多谢 嘿嘿嘿


  • 管理员


Log in to reply
 


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