Commit 6088f2a6 authored by Dan S. Bolintineanu's avatar Dan S. Bolintineanu
Browse files

One more bug in pair_granular.cpp, where JKR cutoff was incorrectly computed.... 

One more bug in pair_granular.cpp, where JKR cutoff was incorrectly computed. Thanks to Ishan Srivastava for finding this one
parent 0673f6f2

src/GRANULAR/pair_granular.cpp

+1 −1
Original line number Diff line number Diff line
@@ -1728,7 +1728,7 @@ double PairGranular::pulloff_distance(double radi, double radj, 
  if (Reff <= 0) return 0;

  coh = normal_coeffs[itype][jtype][3];

  E = normal_coeffs[itype][jtype][0]*THREEQUARTERS;

  a = cbrt(9*MY_PI*coh*Reff/(4*E));

  a = cbrt(9*MY_PI*coh*Reff*Reff/(4*E));

  return a*a/Reff - 2*sqrt(MY_PI*coh*a/E);

}

CRA Git | Maintained and supported by SUSTech CRA and CCSE