I modified a seeding script of mine to allow PSD - it wasn’t too hard since Python has numpy.random.lognormal 
The attached script is able to easily reach a volume fraction of 0.15 which I believe shold be enough to get the bed height you want. If you increase the fraction to 0.2 it gets very slow as it searches for locations to add the new particles…
Here’s the script - the parameters are hard-coded near the top
psd_seeder.py (4.4 KB)
(Side note: a good chunk of the script is just writing out the fixed 35-line file header … although including the instructions in the file is perhaps a nice concession to users, it’s a bit of a pain when trying to generate particle files and this header MUST be present. I think we should reconsider this design decision - perhaps the instructions can be included as optional #comments instead of a fixed number of lines at the top of the file - this would make it easier to generate and manipulate these files.)