Apologies for the same issue again.
According to error I know that I have to change my particle insertion region but I couldn't figure it out right region in my simulation. i am totally stuck here.
It would be a great help if anyone could let me to know the idea to fix my error.
Error:
ERROR: Mesh (id cad): 1 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.000000 / 0.000000 / 0.000000
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.101603 / 0.030640 / -0.093478
(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)
(/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/multi_node_mesh_parallel_I.h:605)
Code:
################################################
# Example of LIGGGHTS #
# Flow particles one by one 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
#create_box 1 reg
#####################################################################################################################
#neighbor 0.002 bin
neighbor 0.003 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 cohesio
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/inclined.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 1 overlapcheck yes all_in yes particles_in_region 1 region bc
#fix ins all insert/pack seed 32452867 distributiontemplate pdd1 vel constant 0. 0. -1. &
# insert_every once overlapcheck yes region reg ntry_mc 10000 volumefraction_region 0.4
#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 plane/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 plane/mesh_*.vtk stress wear cad
#################################################################################################
#restart 5000 restart/restart1*.bin
#unfix ins
run 1000000 upto
unfix ins
stl file:
I used from simple_chute example.
mschramm | Thu, 01/14/2021 - 19:43
size of mesh
Hello,
Your simulation box is too small and a triangular element of your mesh leaves the simulation domain.
From the error, your mesh is centered at 0,0,0. So your simulation domain needs to be at least the same size as the minmax values of your mesh in each direction.
What other changes to the mesh have you done?
I ask because I tried to recreate your error message using the chute_wear example (importing the simple_chute) and your domain, but I got 683 mesh elements leaving the domain and not 1 element leaving the domain.
laudari | Wed, 01/20/2021 - 04:28
Thanks for your reply
Really thank you for your reply.
I already fixed my problem.
Now simulation looks fine.