I'm just seeing a fraction of my particles in paraview, what is the problem ?

Submitted by mahdi_roozbeh on Mon, 12/19/2011 - 10:15

The code is same as pour/dev example, but I've used nparticles 5000.

atom_style granular
atom_modify map array # for allocating the memory, IDs are stored in lookup table of length N
boundary m m m
newton off
echo both

communicate single vel yes
#processors 2 1 1

units si

region reg block -0.05 0.05 -0.05 0.05 0. 0.15 units box
create_box 1 reg

neighbor 0.002 bin
neigh_modify delay 0

#Material properties required for new pair styles

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

#New pair style
pair_style gran/hertz/history 1 0 #Hertzian without cohesion
pair_coeff * *

timestep 0.00001

fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#fix zwalls all wall/gran/hertz/history 1 0 zplane 0.0 0.15 1
#fix cylwalls all wall/gran/hertz/history 1 0 zcylinder 0.05 1

#distributions for insertion
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.0085
fix pts2 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.0025
fix pdd1 all particledistribution/discrete 1. 2 pts1 0.3 pts2 0.7

#region and insertion
group nve_group region reg
region bc cylinder z 0.0 0.0 0.015 0.05 0.12 units box
fix ins nve_group pour/dev nparticles 5000 1 distributiontemplate pdd1 vol 0.55 200 massflowrate 0.1 vel uniform 0. 0. 0. 0. -1.0 region bc

#apply nve integration to all particles that are inserted as single particles
fix integr nve_group nve/sphere

#output settings, include total thermal energy
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

#insert the first particles so that dump is not empty
run 1
dump dmp all custom 200 post1/dump1.pour id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
#insert particles
run 80000 upto
dump dmp1 all custom 200 post/dump1.pour id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
unfix ins

ckloss's picture

ckloss | Mon, 12/19/2011 - 16:36

maybe its a paraview issue, for glyphs there is a maximum # of particles to be shown (default is 5000)

Christoph

mahdi_roozbeh | Wed, 12/21/2011 - 11:36

Some of the spheres run away from the bottom of the container. I dont know what is the problem?
In some cases, I see the particle deposition in different time step, but at the last step, some particles are disappeared suddenly.
I can nearly estimate from paravierw visualization that the number of spheres are significantly lower than 5000.

Far from question of this comment, is there any limitation in the number of spheres in LIGGGHTS ?

Thank you so much,

ckloss's picture

ckloss | Wed, 12/21/2011 - 17:36

I suggest you debug your simulation step-by-step

>>Far from question of this comment, is there any limitation in the number of spheres in LIGGGHTS ?
the current limit is 2 billion

Christoph