hello guys, currently i'm working on granular flow analysis on moon and for that this is the input file that generates output simulation, can you find major issues regarding this?
# LIGGGHTS input script for a simple granular flow (Earth gravity)
units si
atom_style granular
boundary f f f
newton off
# Define simulation box
region domain block 0 0.1 0 0.1 0 0.1 units box
create_box 1 domain
# Define material properties
fix m1 all property/global youngsModulus peratomtype 1 2.0e7
fix m2 all property/global poissonRatio peratomtype 1 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.5
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
fix m5 all property/global coefficientRollingFriction peratomtypepair 1 0.01
# Define particles
region drop_zone block 0.02 0.08 0.02 0.08 0.05 0.08 units box
create_atoms 1 random 1000 15485863 drop_zone
# Particle interaction model
pair_style gran model hertz tangential history
pair_coeff * *
# Ensure proper communication of ghost atom velocity
communicate single vel yes
# Recommended neighbor list modification for granular interactions
neigh_modify delay 0
# Gravity
fix gravity all gravity 9.81 vector 0 0 -1
# Time integration
fix integrate all nve/sphere
# Output
dump myDump all custom 100 flow_output.liggghts id type x y z vx vy vz
thermo 100
# Run simulation
timestep 1e-4
run 10000
i hope you guys reply me soon that will be so helpful