Increase the radius of particle

Submitted by laudari on Fri, 04/02/2021 - 07:14

Hello LIGGGHTS forum,
Really thank you for your help.

I have one more problem in my simulation.
I am trying to increase the radius of my particle from 0.01 to 0.2 or 0.3 but I encountered following error.

"""
ERROR: Unable to calculate region volume. Possible sources of error:

(a) region volume is too small or out of domain (you may want to decrease the 'volume_limit' in the input script)
(b) particles for insertion are too large when using all_in yes
(c) region is 2d, but should be 3d
(/build/liggghts-YO7u74/liggghts-3.8.0+repack1/src/region.cpp:694)

"""

This is my LIGGGHTS script:

"""
################################################Modulus
# Flow particles one by one in a plane surface #
################################################

# Initialize

atom_style sphere
#dimension 3
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
#####################################################################################################################

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.4
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.6
fix m4 all property/global coefficientFriction peratomtypepair 1 0.4
fix m5 all property/global k_finnie peratomtypepair 1 1.0
fix m6 all property/global coefficientRollingFriction peratomtypepair 1 0.02
#######################################################################################################################

#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/new_inclined.stl type 1 wear finnie
fix cad all mesh/surface/stress file meshes/down.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 2200 radius constant 0.2

#Describe the composition of the particle distribution
fix pdd1 all particledistribution/discrete 32452843 1 pts1 1.0

##############################################################################################################

#0 nve_group region domain

group nve_group region domain
region bc cylinder z 0.0 0.0 0.015 0.05 0.12 units box
#fix 1 all nve velocity all pts1 NULL 0.6 NULL units box

#group all region reg

##############################################################################################################

#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

################################################################################################################
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

#variable t equal step*dt
#variable a equal atoms
#variable x equal

##########################################################################################################################
run 1
#########################################################################################################################
dump dmp all custom/vtk 200 same_A/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 same_A/mesh_*.vtk stress wear cad
dump dmp0 all custom 1000 same_A_dump/dump*.data id type type x y z vx vy vz fx fy fz

#dump dmp all custom/vtk 200 new/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 particleB/mesh_*.vtk stress wear cad
#fix screen all print 1 "$t,$a"
#fix file all print 1 "$t,$a" file atoms_per_time.csv title "time,atoms"
#################################################################################################

run 500000 upto

unfix ins

"""

Could you please help me to solve this problem?

Thank you

Forums:

laudari | Sat, 04/03/2021 - 12:09

Really appreciate your answer.
I tried to change all of my coordinates, it works but some of them are out side of domain.

Is there any idea to visualize the coordinates? or do I have to change all of my settings again?

Thank you so much.