src/engine/engineTopoChangerMesh/layerSmooth/setPistonMotionBoundaryConditionLayerSmooth.H

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