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;