Hi,
I try to count all collision events as a function of time in a granular system using contact/atom. This works fine. I get all contact events for each single particle and timestep, but in most cases the collision time is longer than one timestep. As a consequence of this multiple counts for one collision event are recorded. If I use larger timesteps I might miss some collision events.
Is there already any solution for this problem or has anyboy an idea how to fix it? I'm not a c++ expert, therefore it is difficult for me to change the source code.
Thank you in advance and best regards,
Tobias
Laure | Thu, 09/29/2016 - 11:28
pair/gran/local and postprocess
Hello,
The easiest would probably be to dump contacts with the id of both particles of each contact, then use python to load each dump in turn, see which contacts broke or where created, and increase a collision counter accordingly.
See the compute 'pair/gran/local' : http://www.cfdem.com/media/DEM/docu/compute_pair_gran_local.html
Note that pairs ordering can vary from file to file.
Best,
Laure
ckloss | Mon, 10/03/2016 - 11:15
Hi,
Hi,
we have recently developed a fix contact/atom/counter that calculates a histogram of contact events, counting each contact once. It will be in the next release!
Christoph
FelipeL | Mon, 09/04/2017 - 15:20
Hi Christoph,
Hi Christoph,
i have experienced the same problem in my simulation. I get all contacts events but i have multiples colision depending on the chosen timestep. I am trying to implement the fix that you mentioned (fix contact/atom/counter), but i cant find it in the manual. Does it has not yet been implemented or has another fix name?.
Best regards!,
Felipe