in.pour.2d
c1aebdd391f3lammm-master-deprecated

in.pour.2d
View Options

			# Pour 2d granular particles into container

			dimension 2
			atom_style sphere
			boundary f fm p
			newton off
			comm_modify vel yes

			region reg block 0 100 0 50 -0.5 0.5 units box
			create_box 1 reg

			neighbor 0.2 bin
			neigh_modify delay 0

			# IMPORTANT NOTE: these values of stiffness (4000) and timestep (0.001)
			# are used in this example file to produce a quick simulation and movie.
			# More appropriate values for realistic simulations would be
			# k = 4.0e5 and dt = 0.0001, as in bench/in.chute (if it were Hertzian).

			pair_style gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
			pair_coeff * *

			timestep 0.001

			fix 1 all nve/sphere
			fix 2 all gravity 1.0 spherical 0.0 -180.0
			fix xwalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 &
			xplane 0 100
			fix ywalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 &
			yplane 0 NULL

			region slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
			fix ins all pour 1000 1 4767548 vol 0.4 10 &
			diam range 0.5 1.0 region slab

			fix 3 all enforce2d

			compute 1 all erotate/sphere
			thermo_style custom step atoms ke c_1 vol
			thermo 1000
			thermo_modify lost ignore norm no
			compute_modify thermo_temp dynamic yes

			#dump id all atom 250 dump.pour

			#dump 2 all image 500 image.*.jpg type type &
			# zoom 1.6 adiam 1.5
			#dump_modify 2 pad 5

			#dump 3 all movie 500 movie.mpg type type &
			# zoom 1.6 adiam 1.5
			#dump_modify 3 pad 5

			run 25000