## Example 1: Pour grains # Code is adapted from https://github.com/lammps/lammps/blob/master/examples/pour/in.pour # More command details in the documentation https://lammps.sandia.gov/doc/Commands.html atom_style sphere boundary p p f newton off comm_modify vel yes region reg block -3 3 -2 2 0 9 units box create_box 1 reg neighbor 0.2 bin neigh_modify delay 0 pair_style gran/hertz/history 200.0 200.0 5.0 5.0 0.5 1 pair_coeff * * timestep 0.01 fix 1 all nve/sphere fix 2 all gravity 1.0 spherical 0.0 -180.0 region reg2 block -3 3 -2 2 0 9 units box open 6 fix 3 all wall/gran/region hertz/history 200.0 200.0 5.0 5.0 0.5 1 region reg2 region reg3 block -2.0 2.0 -1.0 1.0 3.0 8.0 units box fix ins all pour 60 1 1 vol 0.13 50 region reg3 vel 0 0 0 0 -3 compute 1 all ke compute 2 all erotate/sphere thermo_style custom step atoms c_1 c_2 thermo 100 thermo_modify lost warn norm no dump 2 all image 50 pour_*.jpg type type zoom 1.8 dump_modify 2 pad 5 run 1500