### 1. Initializaion Phase
# Preliminaries
atom_style granular
atom_modify map array # necessary?
boundary f f f # default is p p p
newton off # need more information and examples
communicate single vel yes
units si # This command cannot be used after the simulation box
hard_particles yes
dimension 3 # default value is 3
# Neighbor listing
neighbor 0.003 bin
neigh_modify delay 0
#processors # need more sudy
# Declare domain
region domain_1 block -0.5 0.5 -0.5 0.5 -1.0 1.0 units box # xmin-max ymin-max zmin-max
create_box 2 domain_1
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
### 2. Setup
# Material and interaction properties required
fix m_1 all property/global youngsModulus peratomtype 2.e8 2.e8
fix m_2 all property/global poissonsRatio peratomtype 0.25 0.25
fix m_3 all property/global coefficientRestitution peratomtypepair 2 0.1 0.1 0.1 0.1
fix m_4 all property/global coefficientFriction peratomtypepair 2 0.01 0.01 0.01 0.01
#fix m_5 all property/global characteristicVelocity scalar 2.
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
### 4. Import Mesh from STL
fix cad_box all mesh/surface file aa.stl type 2 scale 1.0
fix gran_walls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad_box
### 5. region of insertion
region gen_box sphere 0.025 0.025 0.15 0.02
fix particle_template_1 all particletemplate/sphere 15485863 atom_type 1 density constant 2500 radius constant 0.001
fix particle_dist_1 all particledistribution/discrete 15485867 1 particle_template_1 1.0
fix gen_paticles all insert/pack seed 32452843 distributiontemplate particle_dist_1 &
insert_every 100 overlapcheck yes all_in yes particles_in_region 30 region gen_box
fix integral_eq_motion all nve/sphere
compute rke all erotate/sphere
#thermo_style custom step atoms ke c_rke vol
#thermo 1000
#thermo_modify lost ignore norm no
### Create Imaging informations
run 1
dump dump_stl all mesh/stl 100 post/dump*.stl
dump dump_atoms all custom/vtk 100 post/atoms*.vtk id type x y z vx vy vz omegax omegay omegaz density diameter
run 50000 upto