Hi,
I'm looking for reference(s) on the way LIGGGHTS/LAMMPS is doing the domain decomposition to track the history variables such as contact information (whether a contact existed before or if it is a new contact). This is for updating, say the incremental frictional forces. Any pointers?
I think the history dependence makes DEM codes slightly different than MD codes.
Thanks,
Joe
ckloss | Mon, 03/11/2013 - 21:31
Hi Joe,
Hi Joe,
you're right - history tracking will have subtle differences in various codes.
As far as LIGGGHTS is concerned, have a look at pair_gran_hooke_history.cpp (for p-p history) and fix_wall_gran.cpp (for p-w history) and fix_contact_history.*
History tracking and domain decomposition are two independent things. if your question is whether history tracking is correctly performed if a pair in contact is crossing domain borders, the answer is yes.
Cheers
Christoph