Number of Particles versus Time

Submitted by emahvelati on Fri, 02/15/2019 - 21:52

Hi Everyone,

I would like to have the plot of # particles versus time after the end of the simulation. Do you suggest to use pizza.py file from LAMMPS or paraview with programmableFilter?

Any other comments and help would be highly appreciated.

Thanks in advance.

richti83's picture

richti83 | Mon, 02/18/2019 - 17:49

You can do this on scriptlevel

variable N equal atoms #special variable, see manual variable.html
variable t equal step*dt #special variables, see manual variable.html
fix out1 all print 10 "$t,$N" file nop.csv title "time,N" screen no #print every 10 iterations to file nop.csv

best, Christian.

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

emahvelati | Sat, 02/23/2019 - 15:59

Thanks Christian for your comment. It is exactly what I finally decided to do. Indeed the best option.

I am wondering if we can export N (number of particles) in a costume region (volume). I am currently exporting mass in a pre-defined domain to calculate N. Would it be possible to export N for a domain directly as we can have for mass?

Regards

ljmutue | Fri, 05/17/2019 - 16:49

Hello,
i have met the same situation here, which is I wanna count the number of particles leaving the system (N), and export N in the system.
So I'm wondering have you had any progress?
Thanks in advance!

Kind regards,
Lijing