Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    phiHbyA边界问题

    OpenFOAM
    2
    2
    510
    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.
    • V
      Vortex last edited by

                  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 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

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

        CFD课程 改成线上了 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • First post
          Last post

        CFD中文网 | 东岳流体 | 京ICP备15017992号-2