#2D with fixed y wall Vibration, Amp 12.5 Tp 10 units nano atom_style sphere dimension 2 boundary p f p #fixed wall delete the particles crossing## #### INITIALISATION #### lattice sq 90 origin 0.5 0.5 0.0 region box block 0 7 0 7 -0.5 0.5 create_box 1 box create_atoms 1 box set type 1 diameter 50 set type 1 mass 276.84 velocity all create 80 87287 dist gaussian ### multi neighbor and comm for efficiency### comm_modify vel yes ### Hybrid Pair Potential ### # kn kt gn gt xmu damp pair_style hybrid/overlay gran/hertz/history 0.0108 NULL 0.0488 NULL 0.1 0 colloid 150.0 pair_coeff * * gran/hertz/history pair_coeff * * colloid 60 25.0 25.0 25.0 150.0 ### Fixed Wall Vibration ### fix 1 all wall/gran hertz/history 0.0108 NULL 0.0488 NULL 0.1 0 yplane 0.0 630.0 wiggle y 48.0 10.0 fix 2 all nve/sphere ### Dump ### dump 1 all custom 500 dump.txt id x y vx vy diameter mass ### Post Processing ### compute 1 all temp/sphere thermo_style custom step ke evdwl thermo_modify temp 1 thermo 100 thermo_modify lost ignore timestep 0.01 run 100000