phiHbyA边界问题

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

李东岳

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