applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/compressibleContinuityErrors.H

   1 {
   2     dimensionedScalar totalMass = fvc::domainIntegrate(rho);
   3
   4     scalar sumLocalContErr =
   5     (
   6         fvc::domainIntegrate(mag(rho - thermo.rho()))/totalMass
   7     ).value();
   8
   9     scalar globalContErr =
  10     (
  11         fvc::domainIntegrate(rho - thermo.rho())/totalMass
  12     ).value();
  13
  14     cumulativeContErr[i] += globalContErr;
  15
  16     Info<< "time step continuity errors (" << mesh.name() << ")"
  17         << ": sum local = " << sumLocalContErr
  18         << ", global = " << globalContErr
  19         << ", cumulative = " << cumulativeContErr[i]
  20         << endl;
  21 }