37 dimensionedScalar totalMass = fvc::domainIntegrate(rho*voidfraction);
39 scalar sumLocalContErr =
40 (fvc::domainIntegrate(mag(rho - thermo.rho())*voidfraction)/totalMass).value();
42 scalar globalContErr =
43 (fvc::domainIntegrate((rho - thermo.rho())*voidfraction)/totalMass).value();
45 cumulativeContErr += globalContErr;
47 Info<<
"time step continuity errors : sum local = " << sumLocalContErr
48 <<
", global = " << globalContErr
49 <<
", cumulative = " << cumulativeContErr