src/engine/engineTopoChangerMesh/thoboisMesh/setThoboisMeshConstraints.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         {
  14
  15             label movingPtsIndex = pZones.findZoneID("movingPointsV"+Foam::name(valveI + 1));
  16             const labelList& movingPointsV = pZones[movingPtsIndex];
  17
  18             forAll(movingPointsV, mpI)
  19             {
  20                 constrainedPoints.append(movingPointsV[mpI]);
  21                 constrainedVelocity.append(valveVel);
  22
  23                 constraintSize++;
  24             }
  25
  26         }
  27
  28         {
  29
  30             label movingInternalPtsIndex = pZones.findZoneID("movingInternalPointsV"+Foam::name(valveI + 1));
  31             const labelList& movingInternalPointsV = pZones[movingInternalPtsIndex];
  32
  33             forAll(movingInternalPointsV, mipI)
  34             {
  35                 constrainedPoints.append(movingInternalPointsV[mipI]);
  36                 constrainedVelocity.append(valveVel);
  37
  38                 constraintSize++;
  39             }
  40
  41         }
  42
  43         {
  44
  45             label staticPtsIndex = pZones.findZoneID("staticPointsV"+Foam::name(valveI + 1));
  46             const labelList& staticPointsV = pZones[staticPtsIndex];
  47
  48             forAll(staticPointsV, spI)
  49             {
  50                 constrainedPoints.append(staticPointsV[spI]);
  51                 constrainedVelocity.append(vector::zero);
  52                 constraintSize++;
  53
  54             }
  55
  56         }
  57
  58     }
  59
  60 }