src/engine/engineTopoChangerMesh/twoStrokeEngine/addPistonLayerHrv.H

   1 if (piston().patchID().active())
   2 {
   3     // add faces for piston layering
   4
   5     faceSet pistonFaceSet(*this, piston().pistonFaceSetName());
   6
   7     boolList flipPistonFaces(pistonFaceSet.toc().size(), false);
   8     label nSet = pistonFaceSet.toc().size();
   9     label nFlip = 0;
  10
  11     forAll (flipPistonFaces, facei)
  12     {
  13         scalar scalProd =
  14             (faceAreas()[pistonFaceSet.toc()[facei]] & vector(0, 0, 1));
  15
  16         if (scalProd < 0)
  17         {
  18             flipPistonFaces[facei] = true;
  19             nFlip++;
  20         }
  21     }
  22
  23     Info << "nSet = " << nSet << endl;
  24     Info << "nFlip = " << nFlip << endl;
  25
  26     fz.append
  27     (
  28         new faceZone
  29         (
  30             "pistonLayerFaces",
  31             pistonFaceSet.toc(),
  32             flipPistonFaces,
  33             nFaceZones,
  34             faceZones()
  35         )
  36     );
  37
  38     nFaceZones++;
  39
  40 //     {
  41 //         pointSet movingPistonPoints(*this, piston().pistonPointSetName());
  42
  43 //         pz.append
  44 //         (
  45 //             new pointZone
  46 //             (
  47 //                 "pistonPoints",
  48 //                 movingPistonPoints.toc(),
  49 //                 nPointZones,
  50 //                 pointZones()
  51 //             )
  52 //         );
  53
  54 //         nPointZones++;
  55 //     }
  56
  57     cellSet movingPistonCells(*this, piston().pistonCellSetName());
  58
  59     Info<< "Adding piston cell set" << endl;
  60     cz.append
  61     (
  62         new cellZone
  63         (
  64             "pistonCells",
  65             movingPistonCells.toc(),
  66             nCellZones,
  67             cellZones()
  68         )
  69     );
  70
  71     nCellZones++;
  72 }