mixing is not rightly taking place between two different size particle

Submitted by umair0737 on Mon, 09/25/2017 - 21:47

aoa

particle are not mixing like discrete . they move like chunk and do not mix properly ... i m uploading my input file plz help me because i have to submit my work this week. i m uploading a pics and code

#Poly dispersed flow

atom_style sphere
atom_modify map array
boundary f f f
newton off

communicate single vel yes

units si

region cyl cylinder z 0. 0. 0.075 0. 0.320 units box
create_box 2 cyl

neighbor 0.003 bin
neigh_modify delay 0

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 5.e6 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.3 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.4 0.4 0.4 0.4
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.01 0.01 0.01 0.01

#New pair style
pair_style gran model hertz tangential history #Hertzian without cohesion
pair_coeff * *
#physical characteristics of particles:
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2456 radius constant 0.0015
fix pdd1 all particledistribution/discrete 49979693 1 pts1 1.0

fix pts2 all particletemplate/sphere 15485867 atom_type 2 density constant 2456 radius constant 0.002
fix pdd2 all particledistribution/discrete 67867967 1 pts2 1.0

#import mesh from cad:
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.020
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.33
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.072 0. 0.
fix cad2 all mesh/surface file meshes/impeller.stl type 1 curvature 1e-5
fix mixer all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad2

#integrator
fix integrate all nve/sphere

fix grav all gravity 9.81 vector 0.0 0.0 -1.0

timestep 0.00005

#Area where the particle are introduced from is set:
region factory1 cylinder z 0. 0. 0.07 0. 0.310 units box

thermo_style custom step atoms ke cpu
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
#fix ctg all check/timestep/gran 1 0.01 0.01
#run 1
#unfix ctg

#defines how the insertion is going to be done:
#particle insertion-large particles
fix ins1 all insert/rate/region seed 86028157 distributiontemplate pdd1 &
nparticles 4000 particlerate 8000 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -2.0 region factory1 ntry_mc 10000
run 40000
unfix ins1

#particle insertion small particles
fix ins2 all insert/rate/region seed 15485863 distributiontemplate pdd2 &
nparticles 3000 particlerate 60000 insert_every 1000 &
overlapcheck yes vel constant 0. 0. -2.0 region factory1 ntry_mc 10000
run 40000
unfix ins2

#make a dump of particles and the stl file
dump dumpstl all mesh/stl 1000 post1/dump*.stl

dump dmp all custom 1000 post1/dump.1 id type type x y z ix iy iz vx vy vz &
fx fy fz omegax omegay omegaz radius
dump dmpvtk all custom/vtk 1000 post1/dump*.poly.vtk id type x y z ix iy iz vx vy vz &
fx fy fz omegax omegay omegaz radius

#rotate central shaft
fix movecad2 all move/mesh mesh cad2 rotate origin 0. 0. 0. axis 0. 0. 1. period -1

#run 5 seconds
run 1000000
unfix movecad2
run 100000

AttachmentSize
Image icon side_view.png55.22 KB
Image icon top_view.png160.26 KB
Image icon bottom_view.png138.05 KB