src/engine/engineTopoChangerMesh/thoboisSliding/setThoboisSlidingConstraint.H

   1 {
   2
   3     const pointZoneMesh& pZones = pointZones();
   4
   5     label constraintSize = 0;
   6
   7     forAll(valves_, valveI)
   8     {
   9
  10         vector valveVel =
  11             valves_[valveI].curVelocity()*valves_[valveI].cs().axis();
  12
  13         if(valves_[valveI].curLift() > valves_[valveI].deformationLift())
  14         {
  15             valveVel = vector::zero;
  16         }
  17
  18
  19         {
  20
  21 //          label movingPtsIndex = pZones.findZoneID("movingPointsV"+Foam::name(valveI + 1));
  22 //          const labelList& movingPointsV = pZones[movingPtsIndex];
  23
  24             labelList movingPointsV;
  25
  26             {
  27                 label valveMovingCellsIndex = cellZones().findZoneID("movingCellsZoneV" + Foam::name(valveI + 1));
  28
  29                 if (valveMovingCellsIndex < 0)
  30                 {
  31                     FatalErrorIn("bool verticalValves::update()")
  32                         << "movingCellsZoneV"+ Foam::name(valveI+1)
  33                         << abort(FatalError);
  34                 }
  35
  36                 const labelList& movingCellsZoneV = cellZones()[valveMovingCellsIndex];
  37
  38                 const labelListList& cp = cellPoints();
  39
  40                 boolList count(newPoints.size(), false);
  41
  42                 forAll (movingCellsZoneV, cellI)
  43                 {
  44                     const labelList& curCellPoints = cp[movingCellsZoneV[cellI]];
  45
  46                     forAll (curCellPoints, i)
  47                     {
  48                         count[curCellPoints[i]] = true;
  49                     }
  50                 }
  51
  52                 // Count the points
  53                 label nCounted = 0;
  54                 forAll (count, pointI)
  55                 {
  56                     if (count[pointI] == true)
  57                     {
  58                         nCounted++;
  59                     }
  60                 }
  61
  62                 movingPointsV.setSize(nCounted);
  63
  64                 // Collect the points
  65                 nCounted = 0;
  66                 forAll (count, pointI)
  67                 {
  68                     if (count[pointI] == true)
  69                     {
  70                         movingPointsV[nCounted] = pointI;
  71                         nCounted++;
  72                     }
  73                 }
  74
  75             }
  76
  77             forAll(movingPointsV, mpI)
  78             {
  79                 constrainedPoints.append(movingPointsV[mpI]);
  80                 constrainedVelocity.append(valveVel);
  81
  82                 constraintSize++;
  83             }
  84
  85         }
  86
  87         {
  88
  89 //        label staticPtsIndex = pZones.findZoneID("staticPointsV"+Foam::name(valveI + 1));
  90 //        const labelList& staticPointsV = pZones[staticPtsIndex];
  91
  92             labelList staticPointsV;
  93
  94             {
  95                 label valveStaticCellsIndex = cellZones().findZoneID("staticCellsZoneV" + Foam::name(valveI + 1));
  96
  97                 if (valveStaticCellsIndex < 0)
  98                 {
  99                     FatalErrorIn("bool verticalValves::update()")
 100                         << "staticCellsZoneV"+ Foam::name(valveI+1)
 101                         << abort(FatalError);
 102                 }
 103
 104                 const labelList& staticCellsZoneV = cellZones()[valveStaticCellsIndex];
 105
 106                 const labelListList& cp = cellPoints();
 107
 108                 boolList count(newPoints.size(), false);
 109
 110                 forAll (staticCellsZoneV, cellI)
 111                 {
 112                     const labelList& curCellPoints = cp[staticCellsZoneV[cellI]];
 113
 114                     forAll (curCellPoints, i)
 115                     {
 116                         count[curCellPoints[i]] = true;
 117                     }
 118                 }
 119
 120                 // Count the points
 121                 label nCounted = 0;
 122                 forAll (count, pointI)
 123                 {
 124                     if (count[pointI] == true)
 125                     {
 126                         nCounted++;
 127                     }
 128                 }
 129
 130                 staticPointsV.setSize(nCounted);
 131
 132                 // Collect the points
 133                 nCounted = 0;
 134                 forAll (count, pointI)
 135                 {
 136                     if (count[pointI] == true)
 137                     {
 138                         staticPointsV[nCounted] = pointI;
 139                         nCounted++;
 140                     }
 141                 }
 142
 143             }
 144
 145
 146             forAll(staticPointsV, spI)
 147             {
 148                 constrainedPoints.append(staticPointsV[spI]);
 149                 constrainedVelocity.append(vector::zero);
 150                 constraintSize++;
 151             }
 152
 153         }
 154
 155     }
 156
 157 }