CFD Online Discussion Forums - Confused about pEqn

Hi everyone,

I am currently doing some work on the caveatingFoam. In which the pEqn is as following:

Code:

fvScalarMatrix pEqn

        (

            fvm::ddt(psi, p)

          - (rhol0 + (psil - psiv)*pSat)*fvc::ddt(alphav) - pSat*fvc::ddt(psi)

          + fvc::div(phiv, rho)

          + fvc::div(phiGradp)

          - fvm::laplacian(rhorAUf, p)

        );

According to the literature （Phd thesis by Fabian), the pEqn should be

Code:

fvm::ddt(psi, p)

          - (rhol0 + (psil - psiv)*pSat)*fvc::ddt(alphav) - pSat*fvc::ddt(psi)

          + fvc::div(phiv, rho)

is you write every term one by one.
so what is the last two terms mean? :

Code:

+ fvc::div(phiGradp)

          - fvm::laplacian(rhorAUf, p)

And I am also confused about the last step in pimple loop:

Code:

if (pimple.finalNonOrthogonalIter())

        {

            phiv += (phiGradp + pEqn.flux())/rhof;

        }

Can anyone give me some hint here? If they are for numerical reasons, where can I find information for them? I didn't see much about this in CFD books I have read. I assume they are very important as I can find the similar term in many other solvers.

Kind regards and Merry Christmas! :D