applications/solvers/combustion/dieselEngineFoam/spraySummary.H

   1     label Nparcels = dieselSpray.size();
   2     reduce(Nparcels, sumOp<label>());
   3
   4     Info<< "\nNumber of parcels in system.... | "
   5         << Nparcels << endl
   6         << "Injected liquid mass........... | "
   7         << 1e6*dieselSpray.injectedMass(runTime.value()) << " mg" << endl
   8         << "Liquid Mass in system.......... | "
   9         << 1e6*dieselSpray.liquidMass() << " mg" << endl
  10         << "SMD, Dmax...................... | "
  11         << dieselSpray.smd()*1e6 << " mu, "
  12         << dieselSpray.maxD()*1e6 << " mu"
  13         << endl;
  14
  15     scalar evapMass =
  16         dieselSpray.injectedMass(runTime.value())
  17     - dieselSpray.liquidMass();
  18
  19     scalar gasMass = fvc::domainIntegrate(rho).value();
  20
  21     if (dieselSpray.twoD())
  22     {
  23         gasMass *= 2.0*mathematicalConstant::pi/dieselSpray.angleOfWedge();
  24     }
  25
  26     scalar addedMass = gasMass - gasMass0;
  27
  28     Info<< "Added gas mass................. | " << 1e6*addedMass << " mg"
  29         << nl << "Evaporation Continuity Error... | "
  30         << 1e6*(addedMass - evapMass) << " mg" << endl;