Im trying to insert particles in a stream from a ternary number based distribution and the resultant distribution does not resemble the desired one.
atom_style granular
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si
region reg block -0.5 0.5 -0.5 0.5 0 2
create_box 2 reg
neighbor 0.002 bin
neigh_modify delay 0
hard_particles yes
fix m1 all property/global youngsModulus peratomtype 1.9e9 5e6
fix m2 all property/global poissonsRatio peratomtype 0.3 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.3 0.1 0.1 0.06
fix m4 all property/global coefficientFriction peratomtypepair 2 0.5 0.5 0.5 0.5
pair_style gran model hertz tangential history #Hertzian withoutcohesion
pair_coeff * *
timestep 0.00001
fix xwalls1 all wall/gran model hertz tangential history primitive type 2 xplane -0.5
fix xwalls2 all wall/gran model hertz tangential history primitive type 2 xplane 0.5
fix ywalls1 all wall/gran model hertz tangential history primitive type 2 yplane -0.5
fix ywalls2 all wall/gran model hertz tangential history primitive type 2 yplane 0.5
fix zwalls1 all wall/gran model hertz tangential history primitive type 2 zplane 0
fix zwalls2 all wall/gran model hertz tangential history primitive type 2 zplane 2
fix inface all mesh/surface file /tmp/mesh/insertion_face.stl type 1
group walls id xwalls1 xwalls2 ywalls1 ywalls2 zwalls1 zwalls2
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
fix pts1 all particletemplate/sphere 38636921 atom_type 1 density constant 7500 radius constant 0.02
fix pts2 all particletemplate/sphere 1848823 atom_type 1 density constant 7500 radius constant 0.015
fix pts3 all particletemplate/sphere 53291611 atom_type 1 density constant 7500 radius constant 0.01
fix pdd1 all particledistribution/discrete/numberbased 12980269 3 pts1 0.27 pts2 0.34 pts3 0.39
fix ins1 all insert/stream seed 92270081 distributiontemplate pdd1 nparticles INF particlerate 400000 insert_every 1000 overlapcheck yes all_in no vel constant 0.0 0.0 -10.0 insertion_face inface
fix integr all nve/sphere
run 1
dump dmp all custom/vtk 10000 /tmp/post/packing_*.vtk id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
run 100000
unfix ins1
deepakpawar.2310 | Wed, 05/19/2021 - 14:00
Hi
Hi
I guess you can use the gaussian command in LIGGGHTS, for generating differrent size particle in LIGGGHTS,
#particle distributions
fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant 2640 radius gaussian number 0.0005 0.0001 # mu (mean) and sigma (standard devation)
fix pdd1 all particledistribution/discrete 15485867 1 pts1 1.0
I hope it might useful for you
Happy learning !
Thank you
mschramm | Wed, 05/19/2021 - 21:42
Gaussian command removed in LIGGGHTS
Hello,
Gaussian is not a valid option in LIGGGHTS anymore...
I do not know exactly when it was removed but it is not an option that can be used in the public release.
Going back to the original question, there may be an issue with using the inf command.
You may want to try swapping out stream for pack and simply state that 1kg or so of material must be in your region (don't make it too large or you will run into the same issue).
Then have LIGGGHTS insert more particles into your region once the particles that were in there have been removed (you should be able to calculate an appropriate insert_every value.
This will keep inserting material into the region until there is 1kg (or so) of material, thus replicating your inf command.
Are you getting too many of pts1?