Particle insertion limit in 1 time step

Hi Christoph,

thanks for your reply.

I have been giving it more work and found out that my simulation is not even working for 1 single imported particle.

I am copying here my input and csv file. If any of you guys can take a look and spot what my error is, I would be very grateful as this is the first time I work with the read_data command.

Thanks in advance!!

#Mill performance

atom_style granular

atom_modify map array

boundary m m m

newton off

communicate single vel yes

units si

region reg block -2.66775 2.66775 -0.245 0.245 -2.66775 2.66775 units box

#create_box 2 reg

# Particle position import

read_data positions_granular2.csv

neighbor 0.002 bin

neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 5.2e6 182e9

fix m2 all property/global poissonsRatio peratomtype 0.3 0.3

fix m3 all property/global coefficientRestitution peratomtypepair 2 0.1 0.2 0.2 1

fix m4 all property/global coefficientFriction peratomtypepair 2 0.545 0.5 0.5 1

#fix m5 all property/global k_finnie peratomtypepair 1 1.0

fix m6 all property/global coefficientRollingFriction peratomtypepair 2 0.01 0.01 0.01 1

#New pair style
pair_style gran model hertz tangential history rolling_friction cdt #Hertzian without cohesion

pair_coeff * *

timestep 0.0003

fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#the CAD mesh

fix inface all mesh/surface file mesh/insertion_face.stl type 1

fix ywalls1 all wall/gran model hertz tangential history primitive type 2 yplane 0.245

fix ywalls2 all wall/gran model hertz tangential history primitive type 2 yplane -0.245

fix drum all mesh/surface file mesh/UntitledM.stl type 2

fix wall all wall/gran model hertz tangential history mesh n_meshes 1 meshes drum

variable omega1 equal 1.3

fix rotate all move/mesh mesh drum rotate/variable origin 0. 0. 0. axis 0. 1. 0. omega v_omega1

#distributions for insertion

#fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2678 radius constant 0.05

#fix pdd1 all particledistribution/discrete 1. 1 pts1 1

#region and insertion

group nve_group region reg

#region bc cylinder z 0.0 0.0 0.015 0.05 0.12 units box

#particle insertion

#fix ins nve_group insert/stream seed 10001 distributiontemplate pdd1 maxattempt 500 nparticles 10000 particlerate 700 start 1 overlapcheck yes vel constant 0 0 -5 insertion_face inface extrude_length 0.5

#apply nve integration to all particles that are inserted as single particles

fix integr nve_group nve/sphere

#output settings
compute 1 all erotate/sphere

#insert the first particles so that dump is not empty

run 1

dump dmp all custom 200 post2/dump*.drum id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius

dump dumpmesh all mesh/gran/VTK 200 post2/dump*.vtk id drum
#dump dumpinface all mesh/gran/VTK 20 post/dumpinface*.vtk id inface

#insert particles
run 1000
upto
unfix ins

"#Particle positions "
"10000 atoms "
"2 atom types "
"-2.66775 2.66775 xlo xhi "
"-0.245 0.245 ylo yhi "
"-2.66775 2.66775 zlo zhi "
" "
"Atoms "
" "
"1 1 0.005 2678 -2.23638 0.033338 0.0914994"
"2 1 0.005 2678 -1.18198 -0.013711 -1.94566"
"3 1 0.005 2678 -2.05143 -0.218333 -0.974516"
"4 1 0.005 2678 -2.24308 -0.215838 -0.988424"
"5 1 0.005 2678 -1.45497 0.217887 -0.426556"
"6 1 0.005 2678 -1.09339 0.0863956 -1.71815"
"7 1 0.005 2678 -1.42043 -0.0644707 -1.07328"
"8 1 0.005 2678 -2.3322 0.0048294 0.0899998"
"9 1 0.005 2678 -2.35614 -0.142244 -0.212717"
"10 1 0.005 2678 -1.24653 0.208074 -0.849847"
"11 1 0.005 2678 -2.18464 -0.170766 0.174428"
"12 1 0.005 2678 -1.45508 0.236039 -1.49035"
"13 1 0.005 2678 -1.78912 0.144933 1.35296"
"14 1 0.005 2678 -1.63818 -0.151748 -0.973731"
"15 1 0.005 2678 -1.24386 -0.127164 -0.290288"

Hi Christoph, thank you so much for the help and attention. I am attaching both files in a tarball.

Thanks again,

Carles

Thanks Christoph!!