Is there a possible way to know the number of contacts of each and every particle(may be with the wall or with the other atoms) and write the data(total number of contacts) into a new file?
I tried to follow this FORUM TOPIC but I couldnt figure it out : http://www.cfdem.com/forums/contactatom
It was printing only zeros for every ntimesteps.
Thank you.
ckloss | Mon, 03/02/2015 - 09:21
Hi,
Hi,
using compute contact atom and a dump local command to output it should do it!
Best wishes
Christoph
ravirasoon1991 | Fri, 03/13/2015 - 08:47
I actually used compute
I actually used compute contact/atom but I was not successfull in achieving the exact number of contacts. (I have 180 particles of two types in my domain)
Might be little stupid to ask you this way.
Can you tell me the way I am supposed to define this "compute" and "dump local" command for the number of contacts with the wall and particle-particle?
Thank you.
richti83 | Fri, 03/13/2015 - 09:36
howto
compute ppc all contact/atom #particle particle contacts,
compute allcontactsinsystem all reduce sum c_ppc #sum all contacts over all procs
variable number_of_contacts equal c_allcontactsinsystem/2 #half bc every atom see the contact AND route to a variable for which can be outputted to file
variable N_atoms equal atoms #number of atoms in system
#screen output
fix screen all print 1 "There are ${number_of_contacts} contacts and ${N_atoms} atoms in system" screen yes
variable t equal step*dt
#file output
fix file all print 1 "$t,${N_atoms},${number_of_contacts} " screen no file collisions.csv title "time,Atoms,Contacts"
[...]
#alternativly dump per atom to postprozess in PV
dump dmp all custom 1000 post/dump*.liggghts id type x y z vx vy vz fx fy fz radius mass c_ppc
ravirasoon1991 | Sat, 03/14/2015 - 05:20
Thanks alot
Thanks alot