applications/solvers/combustion/dieselEngineFoam/createSpray.H

   1 Info << "Constructing Spray" << endl;
   2
   3 PtrList<specieProperties> gasProperties(Y.size());
   4 forAll(gasProperties, i)
   5 {
   6     gasProperties.set
   7     (
   8         i,
   9         new specieProperties
  10         (
  11             dynamic_cast<const reactingMixture&>(thermo()).speciesData()[i]
  12         )
  13     );
  14 }
  15
  16 spray dieselSpray
  17 (
  18     U,
  19     rho,
  20     p,
  21     T,
  22     composition,
  23     gasProperties,
  24     thermo,
  25     environmentalProperties
  26 );
  27
  28 scalar gasMass0 = fvc::domainIntegrate(rho).value();
  29
  30 if (dieselSpray.twoD())
  31 {
  32     gasMass0 *= 2.0*mathematicalConstant::pi/dieselSpray.angleOfWedge();
  33 }
  34
  35 gasMass0 -=
  36     dieselSpray.injectedMass(runTime.value()) - dieselSpray.liquidMass();