hello,
I am trying to perform a simple simulation where I have an auger rotating inside a specified volume filled with powder, but after couple of hours of simulations, I am getting the following error:
ERROR on proc 0: mesh neighbor list overflow, boost neigh_modify one and/or page (../fix_contact_history_mesh.cpp:344)
so what is exactly the cause of this error? And how can I solve it ? Would boosting neigh_modify change the results of my simulation by in terms of particles' properties?
please find the input script I am using:
###### Auger Dosing ######
#coarse graining application
coarsegraining 6.0 model_check warn
modify_timing verbose
### Initialization
# Preliminaries
########## Standard block #############
atom_style sphere
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si
######################################
# Simulation Domain
region domain block -0.05 0.05 -0.05 0.05 -0.086 0.2 units box
create_box 2 domain
# Neighbor listing
neighbor 0.000155 bin #0.25xradius_max #particle diameter
neigh_modify delay 0
### Setup
# Define Material propertiesf
fix m1 all property/global youngsModulus peratomtype 1.0e7 1.0e7
fix m2 all property/global poissonsRatio peratomtype 0.3 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.6 0.2 0.2 0.2
fix m4 all property/global coefficientFriction peratomtypepair 2 0.61 0.115 0.115 0.1
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 1.56 0.6 0.6 0.6
# Define contact model
pair_style gran model hertz tangential history rolling_friction epsd2
pair_coeff * *
# Time step
timestep 0.0000045 #4.5e-6
# Define Gravity
fix grav all gravity 9.81 vector 0.0 0.0 -1.0
# Import mesh from cad
fix cad1 all mesh/surface file meshes/meshed_Container_1.stl type 2 rotate axis 1 0 0 angle 90 move 0.0 0.0 -80.0 scale 0.001
fix cad2 all mesh/surface file meshes/meshed_Hopper_D5_1_extension.stl type 2 rotate axis 1 0 0 angle -90 move 0.0 0.0 38.51 scale 0.001
fix cad3 all mesh/surface file meshes/meshed_Simplified_Auger_D5_ST3_KD3,1_160mm.stl element_exclusion_list read exclude3 type 2 rotate axis 1 0 0 angle 90 &
move 0.0 0.0 63.0 scale 0.001 curvature 1e-5
# Use the imported meshes as granular walls
fix geometry1 all wall/gran model hertz tangential history rolling_friction epsd2 mesh n_meshes 3 meshes cad1 cad2 cad3
#region of insertion
region factory cylinder z 0 0 0.03 0.11 0.12 units box
# Integrator
fix integrate all nve/sphere
#Defining region to measure the mass flow rate
region mass_cylinder cylinder z 0 0 0.023 -0.08 -0.02 units box
#Variables for calculating the flow and the time step that are required in the print data
variable t equal "step*dt"
variable mflow equal "mass(all,mass_cylinder)"
# Create stopper for funnel
fix stopper all wall/gran model hertz tangential history rolling_friction epsd2 primitive type 2 zplane 0.0
# Particle insertion
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 1536.1 &
radius constant 0.000062
fix pdd all particledistribution/discrete 15485867 1 pts1 1.0
fix ins all insert/rate/region seed 32452843 distributiontemplate pdd &
mass 0.04 massrate 0.04 insert_every 1000 &
overlapcheck yes vel constant 0.0 0.0 -1.0 region factory ntry_mc 10000
#40g
# Thermodynamic output settings
thermo_style custom step atoms ke cpu
thermo 10000
thermo_modify lost ignore norm no
# Check time step and initialize dump file
fix ctg all check/timestep/gran 1 0.01 0.01
run 1
#unfix ctg
# Create imaging information #make a dump of particles and the stl file
dump dmp all custom 22224 post/movingMesh_*.liggghts id type type mass x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius diameter
dump dumpstl all mesh/stl 22224 post/dump*.stl
### Execution and further settings
# Fill the hopper
run 444445 upto
# Pull the stopper
unfix stopper
#apply rotation to the auger
variable w equal -41.887902 #400 RPM
fix rotatecad3 all move/mesh mesh cad3 rotate/variable origin 0.0 0.0 0.0 axis 0.0 0.0 1.0 omega v_w
#printout mass flow
fix output all print 111111 "$t,${mflow}" file output.csv screen yes title "t,mflow"
run 13333334
#Simulation is Done
print "The Simulation is now Ready"
Fraonez | Mon, 04/26/2021 - 19:07
did you solve the problem?
hi, im Joaquin Gonzalez
im getting the same error after a while running the simulation, i want to know if you could solve the problem...i increase the number of faces of my mesh and this happend.
if you can helpme i apreciate so much.
i apologize for my english
take care