关于sonicLiquidFoam中pEqn.H的一个疑问

[yuan_neu]

            U = rAU*UEqn.H();

            surfaceScalarField phid
            (
                "phid",
                psi
               *(
                   (fvc::interpolate(U) & mesh.Sf())
                 + rhorAUf*fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho)
                )
            );

            phi = (rhoO/psi)*phid;

            fvScalarMatrix pEqn//for solution of pEqn, construction with respect to p0 and rho0
            (
                fvm::ddt(psi, p)
              + fvc::div(phi)////？？？？？？
              + fvm::div(phid, p)//=0, if psi=0  
              - fvm::laplacian(rhorAUf, p)
            );//

与rhoPimpleFoam相比，为什么多出了第二项 fvc::div(phi)////？？？？？？

[yuan_neu]

是否与readThermodynamicProperties.H里面的
dimensionedScalar rhoO("rhoO", rho0 - psi*p0);有关？

[李东岳]

咦？看看这个：sonicLiquidFoam解析

[yuan_neu]

@李东岳 :kiss: