关于sonicLiquidFoam中pEqn.H的一个疑问
-
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)////??????