src/engine/engineTopoChangerMesh/engineValveSliding/moveValvePointsEngineValveSliding.H

   1     const labelListList& cp = cellPoints();
   2     // valves
   3     forAll(valves_,valveI)
   4     {
   5         vector valveVel =
   6             valves_[valveI].curVelocity()*valves_[valveI].cs().axis();
   7
   8 //         scalar valveDeltaX = valves_[valveI].curVelocity()*
   9 //             valves_[valveI].cs().axis().x()*engTime().deltaT().value() ;
  10 //         scalar valveDeltaY = valves_[valveI].curVelocity()*
  11 //             valves_[valveI].cs().axis().y()*engTime().deltaT().value() ;
  12         scalar valveDeltaZ = valves_[valveI].curVelocity()*
  13             valves_[valveI].cs().axis().z()*engTime().deltaT().value() ;
  14
  15         if(valves_[valveI].curLift() >= valves_[valveI].deformationLift())
  16         {
  17             if(mag(valves_[valveI].curVelocity()) > 0)
  18             {
  19                 Info<< "Valve n. " << valveI << " is moving with velocity = "
  20                     << valves_[valveI].curVelocity() << endl;
  21             }
  22
  23             Info<< "Valve displacement for valve " << valveI << " = "
  24                 << valveDeltaZ << endl;
  25
  26             if
  27             (
  28                 valves_[valveI].poppetPatchID().active()
  29              && valves_[valveI].isOpen()
  30             )
  31             {
  32                 List<bool> valveTopPoint(newPoints.size(), false);
  33
  34 /*
  35                 label layeringVPtsIndex =
  36                     pointZones().findZoneID
  37                     (
  38                         "movingPointsTopZoneV" + Foam::name(valveI+1)
  39                     );
  40
  41                 const labelList& layeringVPoints =
  42                     pointZones()[layeringVPtsIndex];
  43
  44                 forAll(layeringVPoints, i)
  45                 {
  46                     point& p = newPoints[layeringVPoints[i]];
  47
  48                     p += valveVel * engTime().deltaT().value();
  49
  50                     valveTopPoint[layeringVPoints[i]] = true;
  51                 }
  52 */
  53
  54                 label layeringVCellsIndex =
  55                     cellZones().findZoneID
  56                     (
  57                         "movingCellsTopZoneV" + Foam::name(valveI+1)
  58                     );
  59
  60                 const labelList& layeringVCells =
  61                     cellZones()[layeringVCellsIndex];
  62
  63                 labelList layeringVPoints;
  64
  65                 boolList count(newPoints.size(), false);
  66
  67                 forAll(layeringVCells, cellI)
  68                 {
  69                     const labelList& curCellPoints = cp[layeringVCells[cellI]];
  70
  71                     forAll (curCellPoints, i)
  72                     {
  73                         count[curCellPoints[i]] = true;
  74                     }
  75                 }
  76
  77                 // Count the points
  78                 label nCounted = 0;
  79                 forAll (count, pointI)
  80                 {
  81                     if (count[pointI] == true)
  82                     {
  83                         nCounted++;
  84                     }
  85                 }
  86
  87                 layeringVPoints.setSize(nCounted);
  88
  89                 // Collect the points
  90                 nCounted = 0;
  91                 forAll (count, pointI)
  92                 {
  93                     if (count[pointI] == true)
  94                     {
  95                         layeringVPoints[nCounted] = pointI;
  96                         nCounted++;
  97                     }
  98                 }
  99
 100
 101                 forAll(layeringVPoints, i)
 102                 {
 103                     point& p = newPoints[layeringVPoints[i]];
 104
 105                     p += valveVel * engTime().deltaT().value();
 106
 107                     valveTopPoint[layeringVPoints[i]] = true;
 108                 }
 109
 110 /*
 111                 label movingTopPtsVIndex =
 112                     pointZones().findZoneID
 113                     (
 114                         "movingPointsV" + Foam::name(valveI+1)
 115                     );
 116
 117                 const labelList& movingTopPointsV =
 118                     pointZones()[movingTopPtsVIndex];
 119
 120                 forAll(movingTopPointsV, i)
 121                 {
 122                     if(!valveTopPoint[movingTopPointsV[i]])
 123                     {
 124
 125                         point& p = newPoints[movingTopPointsV[i]];
 126                         p += valveVel * engTime().deltaT().value();
 127
 128                     }
 129                 }
 130 */
 131
 132                 label movingTopCellsVIndex =
 133                     cellZones().findZoneID
 134                     (
 135                         "movingCellsZoneV" + Foam::name(valveI+1)
 136                     );
 137
 138                 const labelList& movingTopCellsV =
 139                     cellZones()[movingTopCellsVIndex];
 140
 141                 labelList movingTopPointsV;
 142
 143                 boolList count2(newPoints.size(), false);
 144
 145                 forAll(movingTopCellsV, cellI)
 146                 {
 147                     const labelList& curCellPoints = cp[movingTopCellsV[cellI]];
 148
 149                     forAll (curCellPoints, i)
 150                     {
 151                         count2[curCellPoints[i]] = true;
 152                     }
 153                 }
 154
 155                 // Count the points
 156                 nCounted = 0;
 157                 forAll (count2, pointI)
 158                 {
 159                     if (count2[pointI] == true)
 160                     {
 161                         nCounted++;
 162                     }
 163                 }
 164
 165                 movingTopPointsV.setSize(nCounted);
 166
 167                 // Collect the points
 168                 nCounted = 0;
 169                 forAll (count2, pointI)
 170                 {
 171                     if (count2[pointI] == true)
 172                     {
 173                         movingTopPointsV[nCounted] = pointI;
 174                         nCounted++;
 175                     }
 176                 }
 177
 178                 forAll(movingTopPointsV, i)
 179                 {
 180                     if(!valveTopPoint[movingTopPointsV[i]])
 181                     {
 182
 183                         point& p = newPoints[movingTopPointsV[i]];
 184                         p += valveVel * engTime().deltaT().value();
 185
 186                     }
 187                 }
 188
 189
 190             }
 191             else
 192             {
 193                 Info << "Valve " << valveI << " is CLOSED!!!" << endl;
 194             }
 195
 196         }
 197     }