print density in a region

Submitted by Rahul Kumar on Wed, 03/13/2019 - 08:19

Hello everyone
I am doing a simulation to check the bulk density of the particles, which is filled in a box. I want to know the bulk density or mass of the particles within the box.
How can I do this?

PPhurit | Thu, 03/14/2019 - 02:08

1) Simply use [fix property/atom/tracer] to trace particles, then [compute nparticles/tracer/region] to count them.
2) In your [thermo_style] print the [n_'name of the nparticles'[3]] to print the mass of the particles in the trace region.
I'm sure this approach won't count the overlapping like @richti83 stated above.

(I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®)

Rahul Kumar | Thu, 03/14/2019 - 12:27

Sir
I have a doubt that whether the density entered in liggghts is solid density or bulk density?
------------------------------------------------------
fix pts1 all particletemplate/multisphere 15485863 atom_type 1 density constant 2280 nspheres 4 ntry 1000000 spheres file data/particle1.multisphere scale 0.015 type 1
-------------------------------------------------------------