src/engine/engineTopoChangerMesh/verticalValves/addValvePistonCellZones.H

   1 //  Moving mask: moving cells inside the valves....
   2
   3     {
   4
   5         for(label i = 1; i <= nValves(); i++)
   6         {
   7
   8
   9             if(valves_[i-1].poppetPatchID().active())
  10             {
  11
  12                 labelList movingCells(nCells());
  13
  14                 label nMovingCells = 0;
  15
  16                 label valvePatchID = valves_[i-1].poppetPatchID().index();
  17                 scalar minPopZ = min(boundary()[valvePatchID].patch().localPoints()).z();
  18
  19                 forAll(cellCentres(),cellI)
  20                 {
  21                     const vector& v = cellCentres()[cellI];
  22
  23                     if
  24                     (
  25                         inValve(v, i-1)
  26                         &&
  27                         v.z() > minPopZ
  28                     )
  29                     {
  30                         movingCells[nMovingCells] = cellI;
  31                         nMovingCells++;
  32                     }
  33
  34                 }
  35
  36                 movingCells.setSize(nMovingCells);
  37                 Info << "Number of cells in the moving region poppet valve: " << nMovingCells << endl;
  38
  39                 cz.append
  40                 (
  41                     new cellZone
  42                     (
  43                         "movingCellsTopV"+ Foam::name(i),
  44                         movingCells,
  45                         nCellZones,
  46                         cellZones()
  47                     )
  48                 );
  49
  50                 nCellZones++;
  51
  52             }
  53
  54             if(valves_[i-1].bottomPatchID().active())
  55             {
  56
  57                 labelList movingCells(nCells());
  58
  59                 label nMovingCells = 0;
  60
  61                 label valvePatchID = valves_[i-1].bottomPatchID().index();
  62                 scalar maxPopZ = max(boundary()[valvePatchID].patch().localPoints()).z();
  63
  64
  65                 forAll(cellCentres(),cellI)
  66                 {
  67                     const vector& v = cellCentres()[cellI];
  68
  69                     if
  70                     (
  71                         inValve(v, i-1)
  72                         &&
  73                         v.z() < maxPopZ + SMALL
  74                     )
  75                     {
  76                         movingCells[nMovingCells] = cellI;
  77                         nMovingCells++;
  78                     }
  79
  80                 }
  81
  82                 movingCells.setSize(nMovingCells);
  83                 Info << "Number of cells in the moving region bottom valve: " << nMovingCells << endl;
  84
  85                 cz.append
  86                 (
  87                     new cellZone
  88                     (
  89                         "movingCellsBotV"+ Foam::name(i),
  90                         movingCells,
  91                         nCellZones,
  92                         cellZones()
  93                     )
  94                 );
  95
  96                 nCellZones++;
  97
  98             }
  99
 100
 101         }
 102
 103     }