# 2-d granular flow simulation dimension 2 # 2-d granular flow simulation dimension 2 boundary p s p atom_style granular neighbor 0.3 bin neigh_modify delay 5 newton off # create geometry read_data data.granular # granular potentials pair_style gran/history 200000.0 0.5 0.5 1 # initial velocities compute thermo_rotate all rotate/gran fix 1 all nve/gran thermo_style granular # Couette flow velocity 2 set 0.0 0.0 0.0 velocity 3 set 0.1 0.0 0.0 fix 3 boundary setforce 0.0 0.0 0.0 fix 4 all enforce2d # Run timestep 0.001 thermo 100 #thermo_modify temp mobile dump 1 all atom 50 dump.flow1x dump 2 flow custom 50 dump.velocity1x tag type x y vx vy run 1000