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()) - dieselSpray.liquidMass();
  17
  18         scalar gasMass = fvc::domainIntegrate(rho).value();
  19
  20         if (dieselSpray.twoD())
  21         {
  22             gasMass *= 2.0*mathematicalConstant::pi/dieselSpray.angleOfWedge();
  23         }
  24
  25         reduce(gasMass, sumOp<scalar>());
  26
  27         scalar addedMass = gasMass - gasMass0;
  28
  29         Info<< "Added gas mass = " << 1e6*addedMass << " mg" << nl
  30             << "Evaporation Continuity Error| "
  31             << 1e6*(addedMass - evapMass) << " mg" << endl;