src/engine/engineTopoChangerMesh/engineValveSliding/addValvesFacesPointZonesEngineValveSliding.H

   1
   2     // adding face zones for the vertical valves
   3
   4     for (label valveI = 0; valveI < nValves(); valveI++)
   5     {
   6         // If both sides of the interface exist, add sliding interface
   7         // for a valve
   8         if
   9         (
  10             valves_[valveI].curtainInCylinderPatchID().active()
  11          && valves_[valveI].curtainInPortPatchID().active()
  12         )
  13         {
  14             Info<< "Adding sliding interface zones for curtain of valve "
  15                 << valveI + 1 << endl;
  16
  17             pz.append
  18             (
  19                 new pointZone
  20                 (
  21                     "cutPointsV" + Foam::name(valveI + 1),
  22                     labelList(0),
  23                     nPointZones,
  24                     pointZones()
  25                 )
  26             );
  27
  28             nPointZones++;
  29
  30             const polyPatch& cylCurtain =
  31                 boundaryMesh()
  32                     [valves_[valveI].curtainInCylinderPatchID().index()];
  33
  34             labelList cylCurtainLabels(cylCurtain.size(), cylCurtain.start());
  35
  36             forAll (cylCurtainLabels, i)
  37             {
  38                 cylCurtainLabels[i] += i;
  39             }
  40
  41             fz.append
  42             (
  43                 new faceZone
  44                 (
  45                     "curtainCylZoneV" + Foam::name(valveI + 1),
  46                     cylCurtainLabels,
  47                     boolList(cylCurtainLabels.size(), false),
  48                     nFaceZones,
  49                     faceZones()
  50                 )
  51             );
  52
  53             nFaceZones++;
  54
  55             const polyPatch& portCurtain =
  56                 boundaryMesh()
  57                     [valves_[valveI].curtainInPortPatchID().index()];
  58
  59             labelList portCurtainLabels
  60             (
  61                 portCurtain.size(),
  62                 portCurtain.start()
  63             );
  64
  65             forAll (portCurtainLabels, i)
  66             {
  67                 portCurtainLabels[i] += i;
  68             }
  69
  70             fz.append
  71             (
  72                 new faceZone
  73                 (
  74                     "curtainPortZoneV" + Foam::name(valveI + 1),
  75                     portCurtainLabels,
  76                     boolList(portCurtainLabels.size(), false),
  77                     nFaceZones,
  78                     faceZones()
  79                 )
  80             );
  81
  82             nFaceZones++;
  83
  84             // Add empty zone for cut faces
  85
  86             fz.append
  87             (
  88                 new faceZone
  89                 (
  90                     "cutFaceZoneV" + Foam::name(valveI + 1),
  91                     labelList(0),
  92                     boolList(0, false),
  93                     nFaceZones,
  94                     faceZones()
  95                 )
  96             );
  97
  98             nFaceZones++;
  99
 100         }
 101         else
 102         {
 103             Info << "No valve curtain for valve " << valveI + 1 << endl;
 104         }
 105
 106
 107
 108         // for each valve the following zones have to be created:
 109         //
 110         //  - movingFaceZone
 111         //  - movingPointZone
 112         //  - staticFaceZone
 113         //  - staticPointZone
 114
 115
 116         faceSet layeringFacesTop(*this, valves_[valveI].layeringFacesTopName());
 117
 118         fz.append
 119         (
 120             new faceZone
 121             (
 122                 "layeringFacesTopZoneV" + Foam::name(valveI + 1),
 123                 layeringFacesTop.toc(),
 124                 boolList(layeringFacesTop.size(), true),
 125                 nFaceZones,
 126                 faceZones()
 127             )
 128         );
 129
 130         nFaceZones++;
 131
 132         cellSet movingCellsTop(*this, valves_[valveI].movingCellsTopName());
 133
 134         cz.append
 135         (
 136             new cellZone
 137             (
 138                 "movingCellsTopZoneV" + Foam::name(valveI + 1),
 139                 movingCellsTop.toc(),
 140                 nCellZones,
 141                 cellZones()
 142             )
 143         );
 144
 145         nCellZones++;
 146
 147 /*
 148         pointSet movingPointsTop(*this, valves_[valveI].movingPointsTopName());
 149
 150         pz.append
 151         (
 152             new pointZone
 153             (
 154                 "movingPointsTopZoneV" + Foam::name(valveI + 1),
 155                 movingPointsTop.toc(),
 156                 nPointZones,
 157                 pointZones()
 158             )
 159         );
 160
 161         nPointZones++;
 162 */
 163
 164         cellSet movingCells(*this, valves_[valveI].movingCellsName());
 165
 166         cz.append
 167         (
 168             new cellZone
 169             (
 170                 "movingCellsZoneV" + Foam::name(valveI + 1),
 171                 movingCells.toc(),
 172                 nCellZones,
 173                 cellZones()
 174             )
 175         );
 176
 177         nCellZones++;
 178
 179         cellSet staticCells(*this, valves_[valveI].staticCellsName());
 180
 181         cz.append
 182         (
 183             new cellZone
 184             (
 185                 "staticCellsZoneV" + Foam::name(valveI + 1),
 186                 staticCells.toc(),
 187                 nCellZones,
 188                 cellZones()
 189             )
 190         );
 191
 192         nCellZones++;
 193
 194 /*
 195         pointSet movingPoints(*this, valves_[valveI].movingPointsName());
 196
 197         pz.append
 198         (
 199             new pointZone
 200             (
 201                 "movingPointsV" + Foam::name(valveI + 1),
 202                 movingPoints.toc(),
 203                 nPointZones,
 204                 pointZones()
 205             )
 206         );
 207
 208         nPointZones++;
 209 */
 210
 211 /*
 212         pointSet staticPoints(*this, valves_[valveI].staticPointsName());
 213
 214         pz.append
 215         (
 216             new pointZone
 217             (
 218                 "staticPointsV" + Foam::name(valveI + 1),
 219                 staticPoints.toc(),
 220                 nPointZones,
 221                 pointZones()
 222             )
 223         );
 224
 225         nPointZones++;
 226 */
 227
 228  }
 229
 230