# Conveyor
atom_style    granular
atom_modify    map array
atom_modify    sort 0 1.0
boundary    m m m
newton        off

communicate    single vel yes

units    si

region    reg block -0.1927 0.1927 0.0 1.8386 -0.1674 0.1676 units box
create_box    1 reg

neighbor    0.05 bin
neigh_modify    binsize 0.5 delay 0 every 1 check yes

# Material properties required for new pair styles
fix        m1 all property/global youngsModulus peratomtype 5.e6
fix        m2 all property/global poissonsRatio peratomtype 0.45
fix        m3 all property/global coefficientRestitution peratomtypepair 1 0.7
fix        m4 all property/global coefficientFriction peratomtypepair 1 0.5

# New pair style
pair_style gran model hertz tangential history # Hertzian without cohesion
pair_coeff    * *

timestep    0.00001

fix        gravi all gravity 9.81 vector 0.0 -1.0 0.0 # Gravity along the Y-axis

# Granular walls
fix        bx all mesh/surface file meshes/box.stl type 1 
fix        bxx all mesh/surface file meshes/auger.stl type 1 
fix        bxxx all mesh/surface file meshes/drum.stl type 1

# Move the drum mesh, rotate around the y-axis (this simulates rotation)
fix        movecad_drum all move/mesh mesh bxxx rotate origin 0.0 0.0 0.0 axis 0.0 1.0 0.0 period 250.0

# Insertion face
fix        inface all mesh/surface file meshes/insertion_face.stl type 1

# Granular wall model for the static auger and box
fix wall all wall/gran model hertz tangential history mesh n_meshes 2 meshes bx bxx

fix drum_wall all wall/gran model hertz tangential history mesh n_meshes 1 meshes bxxx


# Define a single particle template
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.007

# Define the distribution template using the single particle template
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0

# Region for insertion
group        nve_group region reg
region        bc block 0.3 0.5 -0.2 0.2 0.1 0.2 units box

# Particle insertion
fix            ins nve_group insert/stream seed 32452867 distributiontemplate pdd1 &
    maxattempt 100 mass 30. massrate 30. overlapcheck yes &
    vel constant 0.0 -1.0 0.0 insertion_face inface extrude_length 0.5

# Run particle insertion for 5000 timesteps
run            5000
unfix        ins  # Stop particle insertion after 5000 timesteps

# Apply NVE integration to all particles that are inserted as single particles
fix        integr nve_group nve/sphere

# Output settings, include total thermal energy
fix                ts all check/timestep/gran 1000 0.1 0.1
compute            rke all erotate/sphere
thermo_style    custom step atoms ke c_rke f_ts[1] f_ts[2] vol
thermo            1000
thermo_modify    lost ignore norm no

# Insert the first particles so that dump is not empty
run            1
dump        dmp all custom/vtk 400 post/conveyor_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius 

# Insert particles for 140000 timesteps
run            1400000 upto
unfix        ins