LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Errors in lubricate/poly ?

# [lammps-users] Errors in lubricate/poly ?

 From: Ranga Radhakrishnan Date: Wed, 21 Jun 2017 11:54:02 +0100

Hi,


I think that I have a few issues with the lubricate/poly implementation in LAMMPS based on my reading of Microhydrodynamics book by Kim and Karilla [1].



2) The first term in the squeeze force seems to be missing a prefactor of 2.


3) \omega^\infty seems to have the wrong units. It is because "h_rate" has the units of length.


4) The pump term is also incorrect for particles of different sizes. Briefly, specific cases of calculation of torques in Ref. [1] cannot be used to write down a generalized version.


Please look at the attached pdf for a more detailed explanation on why I raised these concerns, and how to implement a "corrected" lubrication force if you agree with my concerns. Just to be clear, I have looked at previous messages in the mailing list before I send this message, and I don't think any of the previous messages have answered my concerns.

Regards,

Ranga



Attachment: Forces_Torques.pdf