CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    phiHbyA边界问题

    OpenFOAM
    2
    2
    654
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • V
      Vortex 最后由 编辑

                  volScalarField rAU(1.0/UEqn.A());
      
                  volVectorField HbyA("HbyA", U);
                  HbyA = rAU*UEqn.H();
                  surfaceScalarField phiHbyA
                  (
                      "phiHbyA",
                      (fvc::interpolate(HbyA) & mesh.Sf())
                     + fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)        
                  );
      

      rAU是1/Ap, Ap是矩阵的对角线元素。

      fvm::laplacian(rAU, p) == fvc::div(phiHbyA) 这一步是求解泊松方程,但需要用到phyHbyA的边界值。phiHbyA在两个网格之间的内部面是采用线性插值的interpolate(phiHbyA)这个好理解,但它是如何插值到外部边界上的?
      cfdonline有个同样的讨论,但我还是没弄清楚。https://www.cfd-online.com/Forums/openfoam-programming-development/214351-fvc-interpolate-rau-boundary-faces.html

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        phi,phiHbyA是通量,通量的边界是固定的,不需要插值

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • First post
          Last post