00001 Info << "Constructing Spray" << endl; 00002 00003 PtrList<gasThermoPhysics> gasProperties(Y.size()); 00004 forAll(gasProperties, i) 00005 { 00006 gasProperties.set 00007 ( 00008 i, 00009 new gasThermoPhysics 00010 ( 00011 dynamic_cast<const reactingMixture<gasThermoPhysics>&> 00012 (thermo).speciesData()[i] 00013 ) 00014 ); 00015 } 00016 00017 spray dieselSpray 00018 ( 00019 U, 00020 rho, 00021 p, 00022 T, 00023 composition, 00024 gasProperties, 00025 thermo, 00026 g 00027 ); 00028 00029 scalar gasMass0 = fvc::domainIntegrate(rho).value(); 00030 00031 if (dieselSpray.twoD()) 00032 { 00033 gasMass0 *= 2.0*mathematicalConstant::pi/dieselSpray.angleOfWedge(); 00034 } 00035 00036 gasMass0 -= 00037 dieselSpray.injectedMass(runTime.value()) - dieselSpray.liquidMass(); 00038 00039 // ************************ vim: set sw=4 sts=4 et: ************************ //
