GAUSSIAN particle distribution

Hi rob1,

some time ago I tested the possibilities for particle insertion and all combinations for radius and density distributions worked fine for me. Here is my code, maybe it helps you (just replace the variables ${rho} and ${r1} with your density and radius values and comment/uncomment the regarding code lines):

# Distributions and template for insertion
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant ${rho} radius constant ${r1} # const density, const radius
#2 fix pts1 all particletemplate/sphere 1 atom_type 1 density constant ${rho} radius uniform number 0.10 0.14 # const density, uniform radius
#3 fix pts1 all particletemplate/sphere 1 atom_type 1 density constant ${rho} radius gaussian number 0.12 0.025 # const density, gaussian radius

#4 fix pts1 all particletemplate/sphere 1 atom_type 1 density uniform 200. 300. radius constant ${r1} # unitform density, const radius
#5 fix pts1 all particletemplate/sphere 1 atom_type 1 density gaussian 200. 20. radius constant ${r1} # gaussian density, const radius

#6 fix pts1 all particletemplate/sphere 1 atom_type 1 density uniform 200. 400. radius uniform number 0.08 0.1 # unitform density, uniform radius

Hi mates,

See a discussion by Christian some times ago: http://cfdem.com/forums/how-use-random-uniform-mass-particletemplatesphere
Christoph said the Gaussian distribution has not been implemented ever.

Regards,

Interesting!

The code lines I have posted above worked fine for me (using LIGGGHTS 3.3.1). Also in Paraview the results for both density and radius distributions seem to be as expected.
I have checked my simulation again, it still works fine. I have used a particledistribution with exactly one particle template, maybe this is somehow connected?

However, rob1 already gets an error message when defining the particletemplate itself ...

yes, just pick the line you want to use. These are just several runs that I have tested, in each run I comment and uncomment the corresponding lines.