Hello liggghts forum, I am trying to run a simple simulation after following an example which is available in liggghts. I am just dropping single particles on by one on the region. Some how I did this simulations after changing 'boundary m m m'. But the problem is that particles are dropping outside of my mesh domain.
I have included my script here.
mesh that I am using is from example of liggghts calles 'simple_chute.stl'.
********script*******
###############################################
# Example of LIGGGHTS #
# Flow only two particles in a plane surface #
###############################################
# Initialize
atom_style granular
atom_modify map array
###########################################################################################
boundary f f f
##################################################################################################################
newton off
communicate single vel yes
units si
####################################################################################################################
#region domain block -0.5 0.1 -0.2 0.2 -0.4 0.15 units box
region domain block -0.2 0.0 -0.2 0.2 -0.4 0.3 units box
#region domain block 0.1 0.1 0.2 0.2 0.15 0.15 units box
#####################################################################################################################
create_box 1 domain
#####################################################################################################################
#neighbor 0.002 bin
neighbor 0.001 bin
neigh_modify delay 0
#####################################################################################################################
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.3
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5
fix m5 all property/global k_finnie peratomtypepair 1 1.0
#######################################################################################################################
#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
#############################################################################################
timestep 0.00001
##############################################################################################
fix integrate all nve/sphere
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#############################################################################################
fix cad all mesh/surface/stress file meshes/simple_chute.stl type 1 wear finnie
fix inface all mesh/surface file meshes/insertion_face.stl type 1
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad
fix mass all massflow/mesh mesh inface count once point_at_outlet 0. 0. 0.
#fix granwalls all wall/gran model hertz tangential history rolling_friction epsd2 mesh n_meshes 1 meshes cad
#######################################################################################################
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 200 radius constant 0.0025
fix pts2 all particletemplate/sphere 15485867 atom_type 1 density constant 200 radius constant 0.0025
#Describe the composition of the particle distribution
fix pdd1 all particledistribution/discrete 32452843 2 pts1 0.3 pts2 0.7
##############################################################################################################
group nve_group region domain
#region bc block -0.5 0.1 -0.2 0.2 -0.4 0.15 units box
region bc block -0.2 0.0 -0.2 0.2 -0.4 0.3 units box
#region bc block 0.0 0.0 0.2 0.2 0.3 0.3 units box
##############################################################################################################
#particle insertion
fix ins nve_group insert/pack seed 32452867 distributiontemplate pdd1 insert_every 2 overlapcheck yes all_in yes vel constant -0.1 0.0 -0.1 particles_in_region 1 region bc
#run 5000
#unfix ins1
################################################################################################################
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no
##########################################################################################################################
run 1
#########################################################################################################################
dump dmp all custom/vtk 200 result/chute_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dumpstress all mesh/gran/VTK 200 result/mesh_*.vtk stress wear cad
#################################################################################################
#restart 5000 restart/restart1*.bin
run 100000 upto
unfix ins
when I use 'boundary f f f', I get following error,
ERROR: Mesh (id cad): 683 mesh elements have been lost / left the domain.
Please use 'boundary m m m' or scale/translate/rotate the mesh or change its dynamics
FYI: center of mass of mesh including scale/tranlate/rotate before cutting out elements is -0.148131 / 0.003869 / -0.127528
simulation box x from -0.200000 to 0.000000 y from -0.200000 to 0.200000 z from -0.400000 to 0.300000
center of mass of mesh after cutting out elements outside simulation box is is -0.098473 / 0.011222 / -0.102990
(gives indication about changes in scale/tranlate/rotate necessary to make simulation run)
(/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/multi_node_mesh_parallel_I.h:605)
Could someone help me to fix my region, I am totally stuck here.
I am a newbie in LIGGGHTS.
Thank you so much .