Hello everyone,
For my thesis, my supervisor would like me to calculate the overall coordination number (average of all coordination numbers per timestep) which should exclude the so-called "rattlers" (From C.Thorhton the average coordination number is defined as Zm = (2C - N1) / (N - N0 - N1), where N1 and N0 are the number of particles with only one or zero contacts respectively.
Does anyone have any suggestion, about how I could implement this in my code?
Greetings,
Luca
richti83 | Tue, 10/27/2015 - 10:58
you can count the number of
you can count the number of contacts per atoms with
http://www.cfdem.com/media/DEM/docu/compute_contact_atom.html
compute ca all contact/atom
than use dmp all dump custom 1000 post/dump*.liggghts id x y z c_ca
or use any compute reduce command to calculate your formula.
ZhaochenJiang | Fri, 11/06/2015 - 10:04
another approach is to use
another approach is to use "pair/gran/local" command
http://www.cfdem.com/media/DEM/docu/compute_pair_gran_local.html
for instance, compute contact all pair/gran/local id pos vel .
Then, save the contact pairs in corresponding time steps by
dump contact all local 1000 ../DEM/post dump*.contact c_contact[1] ...
This method also give you the chance to acquire velocity and force of the contact pairs.