mesh neighbor list overflow

Submitted by NSalloum on Wed, 11/14/2018 - 08:59

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

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