Commit 1d4372c2 authored by Dan S. Bolintineanu's avatar Dan S. Bolintineanu
Browse files

Fixes bug in calculation of relative tangential velocity for granular wall particle interactions

parent 0b24a3d8

src/GRANULAR/fix_wall_gran.cpp

+3 −3
Original line number Diff line number Diff line
@@ -1186,9 +1186,9 @@ void FixWallGran::granular(double rsq, double dx, double dy, double dz, 
  vt3 = vr3 - vn3;



  // relative rotational velocity

  wr1 = radius*omega[0] * rinv;

  wr2 = radius*omega[1] * rinv;

  wr3 = radius*omega[2] * rinv;

  wr1 = radius*omega[0];

  wr2 = radius*omega[1];

  wr3 = radius*omega[2];



  // relative tangential velocities

  vtr1 = vt1 - (nz*wr2-ny*wr3);

CRA Git | Maintained and supported by SUSTech CRA and CCSE