Free falling of single atoms in the box

Box : loaded from STL file

### 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
dimension                  3                                          # default value is 3

neighbor     0.002 bin
neigh_modify delay 0


#region       reg block -0.05 0.05 -0.05 0.05 0. 0.15 units box
#create_box   1 reg
#neighbor     0.002 bin
#neigh_modify delay 0

#processors # need more sudy

# Declare domain
region                     domain block -0.5 0.5 -0.5 0.5 -1.0 1.0 units box   # xmin-max ymin-max zmin-max
create_box                 1 domain


# Neighbor listing
#neighbor                   0.003 bin
#neigh_modify               every 1 check no

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.5e7
fix                         m_2 all property/global poissonsRatio peratomtype 0.25
fix                         m_3 all property/global coefficientRestitution peratomtypepair 1 0.5
fix                         m_4 all property/global coefficientFriction peratomtypepair 1 0.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 1 scale 1.0
fix                        gran_walls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad_box

#fix                        box all mesh/surface file insertion_face.stl type 1 scale 1.0

create_atoms               1 single 0.025 0.025 0.1
set                        atom 1 type 1 density 2000 diameter 0.01

fix                        integral_eq_motion all nve/sphere

timestep    0.00001
### 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 density diameter
run                        100000 upto

#fix id all property/global youngsModulus peratomtype value_1 value_2 ...
#    (value_i=value for Youngs Modulus of atom type i)
#fix id all property/global poissonsRatio peratomtype value_1 value_2 ...
#    (value_i=value for Poisson ratio of atom type i)
#fix id all property/global coefficientRestitution peratomtypepair n_atomtypes value_11 value_12 .. value_21 value_22 .. .
#    (value_ij=value for the coefficient of restitution between atom type i and j; n_atomtypes is the number of atom types you want to use in your simulation)
#fix id all property/global coefficientFriction peratomtypepair n_atomtypes value_11 value_12 .. value_21 value_22 .. .
#    (value_ij=value for the (static) coefficient of friction between atom type i and j; n_atomtypes is the number of atom types you want to use in your simulation)