Random insertion or particle's properties wise in Silo simulation

Submitted by gibel on Mon, 01/17/2022 - 04:54

Hello Mschramm,

Apologies for many questions.

I have been doing simulation in Silo with liggghts. I have one more question about particles filling.

This is how I have been generating particles and inserting :
```
###################################Generation and Insertion of the particles#########################################

#distributions for insertion

fix pts1 all particletemplate/sphere 15485863 atom_type 1 density constant ${density} radius constant ${radius1}
fix pts2 all particletemplate/sphere 15485867 atom_type 2 density constant ${density} radius constant ${radius2}
fix pdd1 all particledistribution/discrete 32452843 ${nradii} pts1 ${frac1} pts2 ${frac2}

fix ins_mesh all mesh/surface/planar file Insertionsface.stl type 1 scale 0.001

fix ins all insert/stream seed 32452867 distributiontemplate pdd1 &
mass ${fillmass} massrate ${fillmassrate} insert_every 3000 overlapcheck yes all_in no vel constant 0 0 -0.75 &
insertion_face ins_mesh
```

1) Does this set up fills silo randomly with particles 1 and 2 or with some pattern(like bigger particles first and smaller second)?

If it doesn't fill randomly, is there any idea to fill silo randomly with two types of particles?

THank you so much.

AttachmentSize
Plain text icon Main simulation file7.61 KB

mschramm | Mon, 01/17/2022 - 18:14

Hello,
During the insertion process, the particledistribution/discrete takes your user supplied mass fractions and converts them to number fractions.
Then given the number of spheres/mass you would like in your simulation, calculates the number of spheres of each type of particle to insert.
LIGGGHTS' then takes the largest particle and inserts them all randomly into the simulation first.
Then the smaller particle is selected and is inserted randomly around the larger particles.

This is why you may try inserting too many particles into too small of a volume and only one particle type is inserted.

I do not know of any user created insertion type that randomly selects a type of particle and then tries to randomly insert it into the simulation.

You could do something like this using python and a "create_atom" statement.