src/engine/engineTopoChangerMesh/twoStrokeEngine/addPistonLayer.H

   1 if (piston().patchID().active())
   2 {
   3     // Add faces for piston layering
   4     // Note: because of operation of layer addition/removal
   5     // (reduce function in layer addition/removal thickness)
   6     // the layering modifier needs to be present on all processors
   7     // even if the patch size is zero
   8     // HJ, 7/Mar/2011
   9
  10     faceSet pistonFaceSet(*this, piston().pistonFaceSetName());
  11
  12     boolList flipPistonFaces(pistonFaceSet.toc().size(), false);
  13     label nSet = pistonFaceSet.toc().size();
  14     label nFlip = 0;
  15
  16     forAll (flipPistonFaces, facei)
  17     {
  18         scalar scalProd =
  19             (faceAreas()[pistonFaceSet.toc()[facei]] & vector(0, 0, 1));
  20
  21         if (scalProd < 0)
  22         {
  23             flipPistonFaces[facei] = true;
  24             nFlip++;
  25         }
  26     }
  27
  28     Info << "nSet = " << nSet << endl;
  29     Info << "nFlip = " << nFlip << endl;
  30
  31     fz[nFaceZones] =
  32     (
  33         new faceZone
  34         (
  35             "pistonLayerFaces",
  36             pistonFaceSet.toc(),
  37             flipPistonFaces,
  38             nFaceZones,
  39             faceZones()
  40         )
  41     );
  42
  43     nFaceZones++;
  44
  45     cellSet movingPistonCells(*this, piston().pistonCellSetName());
  46
  47     Info<< "Adding piston cell set" << endl;
  48     cz[nCellZones] =
  49     (
  50         new cellZone
  51         (
  52             "pistonCells",
  53             movingPistonCells.toc(),
  54             nCellZones,
  55             cellZones()
  56         )
  57     );
  58
  59     nCellZones++;
  60 }