37 volScalarField contErr( fvc::div(phiGes) + fvc::ddt(voidfraction) );
39 scalar sumLocalContErr = runTime.deltaTValue()*
40 mag(contErr)().weightedAverage(mesh.V()).value();
42 scalar globalContErr = runTime.deltaTValue()*
43 contErr.weightedAverage(mesh.V()).value();
44 cumulativeContErr += globalContErr;
46 Info<<
"time step continuity errors : sum local = " << sumLocalContErr
47 <<
", global = " << globalContErr
48 <<
", cumulative = " << cumulativeContErr