applications/solvers/combustion/dieselEngineFoam/createSpray.H

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