00001 label Nparcels = dieselSpray.size(); 00002 reduce(Nparcels, sumOp<label>()); 00003 00004 Info<< "\nNumber of parcels in system.... | " 00005 << Nparcels << endl 00006 << "Injected liquid mass........... | " 00007 << 1e6*dieselSpray.injectedMass(runTime.value()) << " mg" << endl 00008 << "Liquid Mass in system.......... | " 00009 << 1e6*dieselSpray.liquidMass() << " mg" << endl 00010 << "SMD, Dmax...................... | " 00011 << dieselSpray.smd()*1e6 << " mu, " 00012 << dieselSpray.maxD()*1e6 << " mu" 00013 << endl; 00014 00015 scalar evapMass = 00016 dieselSpray.injectedMass(runTime.value()) 00017 - dieselSpray.liquidMass(); 00018 00019 scalar gasMass = fvc::domainIntegrate(rho).value(); 00020 00021 if (dieselSpray.twoD()) 00022 { 00023 gasMass *= 2.0*mathematicalConstant::pi/dieselSpray.angleOfWedge(); 00024 } 00025 00026 scalar addedMass = gasMass - gasMass0; 00027 00028 Info<< "Added gas mass................. | " << 1e6*addedMass << " mg" 00029 << nl << "Evaporation Continuity Error... | " 00030 << 1e6*(addedMass - evapMass) << " mg" << endl;
