Are there and commands to see how many particles are in contact with a wall?

Submitted by J_Wang on Wed, 03/14/2012 - 00:56

Hi everyone,

Are there any commands to see how many particles are in contact with a wall? Or in the code, which part can give a output for the number of particles on a wall? I see the fix mesh/gran command. It can only output the force and torque.

Thanks a lot.

J.Wang

Forums:

ckloss's picture

ckloss | Wed, 03/14/2012 - 08:08

There is a

int FixWallGran::n_contacts()

function that does this. Note that it works for mesh walls only, otherwise it will return 0.
Also note that it counts the number of contacts with the wall, not the number of particles in contact (can be different if a particle contacts twice at a corner)
If you want the number of particles in contact, you would have to implement a slightly modified version of this function

Cheers, Christoph

J_Wang | Mon, 03/19/2012 - 00:47

Hi Christoph,

I may need the contact number in the script level, are there output command available for this funtion? Or do I need to write some additional code to do that?

Thanks

J.Wang

skyopener's picture

skyopener | Mon, 03/19/2012 - 02:23

I agree with you.
currently we can dump several quantities of the mesh, such as the id vel stress wear and so on.
it will be highly appreciated if we can dump more features of mesh from the input of the script, such as the conductive heat flux or heat transfer coefficient along with the contact number.........
s.l