src/engine/engineTopoChangerMesh/accordionEngineMesh/setPistonMotionBoundaryConditionAccordionEngineMesh.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         motionU.boundaryField()[piston().patchID().index()] == pistonVel;
  21 //        motionU.boundaryField()[piston().patchID().index()] == vector::zero;
  22
  23     }
  24
  25 }