applications/solvers/combustion/dieselEngineFoam/createSpray.H

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