applications/solvers/multiphase/cavitatingFoam/continuityErrs.H

   1 {
   2     volScalarField thermoRho = psi*p + (1.0 - gamma)*rhol0;
   3
   4     dimensionedScalar totalMass = fvc::domainIntegrate(rho);
   5
   6     scalar sumLocalContErr =
   7     (
   8         fvc::domainIntegrate(mag(rho - thermoRho))/totalMass
   9     ).value();
  10
  11     scalar globalContErr =
  12     (
  13         fvc::domainIntegrate(rho - thermoRho)/totalMass
  14     ).value();
  15
  16     cumulativeContErr += globalContErr;
  17
  18     Info<< "time step continuity errors : sum local = " << sumLocalContErr
  19         << ", global = " << globalContErr
  20         << ", cumulative = " << cumulativeContErr
  21         << endl;
  22 }