Computer simulation of random packing of unequal particles

A Monte Carlo simulation model for the random packing of unequal spherical particles is presented in this paper. With this model, the particle radii obeying a given distribution are generated and randomly placed within a cubic packing domain with high packing density and many overlaps. Then a relaxation iteration is applied to reduce or eliminate the overlaps, while the packing space is gradually expanded. The simulation is completed once the mean overlap value falls below a preset value. To simulate the random close packing, a "vibration" process is applied after the relaxation iteration. For log-normal distributed particles, the effect of particle size standard deviation, and for bidisperse particles, the effects of particle size ratio and the volume fraction of large particles on packing density and on coordination number are investigated. Simulation results show good agreement with that obtained by experiments and by other simulations. The randomness, homogeneity, and isotropy, which have not been evaluated before for packing of distributed particles, are also examined using statistical measures.