I can't completely fill the cylindrical Container

Submitted by mahdi_roozbeh on Mon, 12/26/2011 - 04:05

Hi Friends:

I don't know why I can not completely fill the container. Here is the snip of my code.

boundary p p f
region reg cylinder z 0.0 0.0 10 0.00 50 units box
create_box 1 reg
fix m1 all property/global youngsModulus peratomtype 5.e9
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 .3
fix m4 all property/global coefficientFriction peratomtypepair 1 .5

pair_style gran/hertz/history 1 0 #Hertzian without cohesion
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 50 1

fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 1
fix pdd1 all particledistribution/discrete 1 1 pts1 1

group nve_group region reg
fix ins nve_group pour/dev nparticles 5000 20 distributiontemplate pdd1 vol 0.6 10 vel uniform 0. 0. 0. 0. -5.0 region reg

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

#insert particles
run 100000 upto

Thanks, if you could help me.

Forums:

raguelmoon's picture

raguelmoon | Mon, 12/26/2011 - 06:21

Hi,

You have chosen whole simulation region as an ins region, new pouring particles dont have space to be inserted.

Best,
Ram

Ram

mahdi_roozbeh | Mon, 12/26/2011 - 08:50

I put the below command before fix pour/dev:

group nve_group region reg
region bc cylinder z 0.0 0.0 10 0 50 side in units box
fix ins nve_group pour/dev nparticles 5000 20 distributiontemplate pdd1 vol 0.6 10 vel uniform 0. 0. 0. 0. -5.0 region bc

And I changed the main region as follows:

region reg cylinder z 0.0 0.0 10 0.00 150 units box

Again, same number of particles are placed within the container.

I can change the height of cylindrical containers to 150 in both regions, but it will cost too much time.

What is your idea about that?

Thanks
Mahdi

raguelmoon's picture

raguelmoon | Thu, 12/29/2011 - 17:36

I guess you want to fill cylinder by pouring particles under gravity. Your cylinder height is 50 units, so do this:

region reg cylinder z 0.0 0.0 10 0.00 101 units box

and set pouring region to:

region bc cylinder z 0.0 0.0 10 55 100 units box

and set vol to 0.3

best,
Ram

Ram