2006-09-28 03:50:23 +08:00
|
|
|
# LAMMPS benchmark of granular flow
|
|
|
|
|
# chute flow of 32000 atoms with frozen base at 26 degrees
|
|
|
|
|
|
|
|
|
|
variable x index 1
|
|
|
|
|
variable y index 1
|
|
|
|
|
|
|
|
|
|
units lj
|
2011-04-14 05:40:28 +08:00
|
|
|
atom_style sphere
|
2006-09-28 03:50:23 +08:00
|
|
|
boundary p p fs
|
|
|
|
|
newton off
|
2015-03-18 00:06:46 +08:00
|
|
|
comm_modify vel yes
|
2006-09-28 03:50:23 +08:00
|
|
|
|
|
|
|
|
read_data data.chute
|
|
|
|
|
|
|
|
|
|
replicate $x $y 1
|
|
|
|
|
|
2009-01-09 23:14:47 +08:00
|
|
|
pair_style gran/hooke/history 200000.0 NULL 50.0 NULL 0.5 0
|
2009-01-07 01:12:09 +08:00
|
|
|
pair_coeff * *
|
2006-09-28 03:50:23 +08:00
|
|
|
|
|
|
|
|
neighbor 0.1 bin
|
|
|
|
|
neigh_modify every 1 delay 0
|
|
|
|
|
|
|
|
|
|
timestep 0.0001
|
|
|
|
|
|
|
|
|
|
group bottom type 2
|
|
|
|
|
group active subtract all bottom
|
|
|
|
|
neigh_modify exclude group bottom bottom
|
|
|
|
|
|
2007-07-09 22:17:03 +08:00
|
|
|
fix 1 all gravity 1.0 chute 26.0
|
2006-09-28 03:50:23 +08:00
|
|
|
fix 2 bottom freeze
|
2008-05-17 04:30:53 +08:00
|
|
|
fix 3 active nve/sphere
|
2006-09-28 03:50:23 +08:00
|
|
|
|
2008-05-17 04:30:53 +08:00
|
|
|
compute 1 all erotate/sphere
|
|
|
|
|
thermo_style custom step atoms ke c_1 vol
|
|
|
|
|
thermo_modify norm no
|
2006-09-28 03:50:23 +08:00
|
|
|
thermo 100
|
|
|
|
|
|
|
|
|
|
run 100
|
