buoyantFoam代码问题求助

[Miraitowa]

在pEqn.H中对U进行更新时有

const surfaceScalarField phig(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf());
U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf);

我想请问其中的p_rghEqn.flux()代表什么意思，以及为什么phig中要先乘rhorAUf再在U中除以rhorAUf；
以及之后的根据密度差判断收敛性中，我发现compressibleContinuityErrs.H代码如下，只计算了误差值，但没有判断是否符合要求，想请问各位大佬是在哪里判断的呢

    dimensionedScalar totalMass = fvc::domainIntegrate(rho);

    scalar sumLocalContErr =
        (fvc::domainIntegrate(mag(rho - thermo.rho()))/totalMass).value();

    scalar globalContErr =
        (fvc::domainIntegrate(rho - thermo.rho())/totalMass).value();

    cumulativeContErr += globalContErr;

    Info<< "time step continuity errors : sum local = " << sumLocalContErr
        << ", global = " << globalContErr
        << ", cumulative = " << cumulativeContErr
        << endl;
}

谢谢大家！