LIGGGHTS® - User Forum

LIGGGHTS® related topics can be discussed here: discussion about models, installation, feature requests and general discussion

Statistics of particles

Theo_Score's picture
Submitted by Theo_Score on Mon, 01/23/2017 - 06:53

Hi All,

Find attached insertion_face.stl (to be put in folder meshes) and an input file.

I am troubled with need to count the particles of Type 1 with time in the domain and I have defined the particles of Type 1 as tracer and the output is zero, I am wondering what could be wrong with my script.

Thank you for your help.

Coordination number per atom

Theo_Score's picture
Submitted by Theo_Score on Wed, 01/11/2017 - 15:02

Hi all,

I am trying to get the coordination number of every atom during a simulation and I set out coord/atom command as

#Output settings, include total thermal energy
fix ts all check/timestep/gran 1000 0.15 0.15
compute rke all erotate/sphere
compute nconts all coord/atom 0.015 mix yes
thermo_style custom step atoms ke c_rke c_nconts f_ts[1] f_ts[2] vol
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

How to get the instant number of particles?

Submitted by zhou_yu on Tue, 01/10/2017 - 12:36

I am doing a simulation of sandpile formation. I am wodering how can I get the instant number of particles in the sandpile since I use "lost ignore" to delete the paticles that gose to far from the pile.

Actually, I want to use this number in a judgement condition of a loop to control the max number I should pour on the sandpile.

set command and region command with move keyword

Submitted by Weijing on Tue, 01/10/2017 - 01:10

Hi everyone,

I am trying to learn more about the LIGGGHTS. In this scenarios, I tried to set a moving region with constant temperature, so I wrote

variable dy equal "0.08*(step-20000)/30000"
region pos_source block -0.05 0.05 -0.05 -0.03 0.0 0.05 move NULL v_dy NULL
fix test all print 1000 "${dy}"
set region pos_source property/atom Temp 500

But the result shows the set command only affect once and give the same result as

how to make a particle spin?

Submitted by chengqiccc on Sun, 01/08/2017 - 13:21

Hi all,

I'm new to LIGGGHTS. I want make a particle spin in the simulation.
At first, the atom_style in my simulation is granular. Under these circumstances, I set the angular velocity via "velocity......setAngular......" , and the terminal give a feedback:

"set Angular requires the atom property 'angmom' as defined by atom_style ellipsoid. (../velocity.cpp:426)"

Then I use atom_style "ellipsoid", but the terminal give a feedback this time:

"Fix wall/gran (id boxwalls_x1): requires atom attributes radius, omega, torque (../fix_wall_gran.cpp:107)"

Pages

Subscribe to RSS - LIGGGHTS® - User Forum