src/engine/engineTopoChangerMesh/thoboisMesh/setPistonMotionBoundaryConditionThobois.H

   1 // use tetrahedral decomposition of the engine mesh
   2
   3     vector pistonVel =
   4             piston().cs().axis()*engineTime_.pistonSpeed().value();
   5
   6 {
   7     tetPointVectorField& motionU = mSolver.motionU();
   8
   9     Info << "setting the piston velocity" << endl;
  10
  11     // Set piston velocity
  12     if (piston().patchID().active())
  13     {
  14
  15         if (debug)
  16         {
  17             Info << "Piston velocity: " << pistonVel;
  18         }
  19
  20         componentMixedTetPolyPatchVectorField& pp =
  21             refCast<componentMixedTetPolyPatchVectorField>
  22             (
  23                 motionU.boundaryField()[piston().patchID().index()]
  24             );
  25
  26         pp.refValue() = vector::zero;
  27
  28     }
  29
  30 }