Random particle packing with large particle size variations using reduced-dimension algorithms

M. D. Webb,I. L. Davis

DOI: https://doi.org/10.48550/arXiv.cond-mat/0606249

2006-06-09

Abstract:We present a reduced-dimension, ballistic deposition, Monte Carlo particle packing algorithm and discuss its application to the analysis of the microstructure of hard-sphere systems with broad particle size distributions. We extend our earlier approach (the ``central string'' algorithm) to a reduced-dimension, quasi-3D approach. Our results for monomodal hard-sphere packs exhibit a calculated packing fraction that is slightly less than the generally accepted value for a maximally random jammed state. The pair distribution functions obtained from simulations of composite structures with large particle size differences demonstrate that the algorithm provides information heretofore not attainable with existing simulation methods, and yields detailed understanding of the microstructure of these composite systems.

Materials Science

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is how to effectively simulate the packing behavior of particle systems with large particle - size differences, especially the microstructural characteristics in these systems. Traditional three - dimensional simulation methods face problems of high computational cost and difficulty in obtaining sufficient statistical accuracy when dealing with particle systems with large particle - size differences. Therefore, the authors proposed a reduced - dimension algorithm to overcome these problems and provide a detailed understanding of the microstructure of composite materials.

### Specific problems include:
1. **Computational efficiency**: For particle systems with large particle - size differences, traditional three - dimensional simulations need to handle hundreds of millions of particles, and the calculation time is extremely long. For example, in a spherical particle system with a size ratio of 100:1, 10^6 small balls may be required for each large ball, which is not feasible with existing computational resources.

2. **Microstructure analysis**: Existing simulation methods are difficult to provide detailed information about the microstructure of composite materials (such as the radial distribution function and coordination number), and this information is crucial for understanding the macroscopic properties of materials.

3. **Particle stratification phenomenon**: During the simulation process, small particles are likely to sink through the gaps between large particles, resulting in axial segregation. In addition, large particles may also be stratified in the radial direction, and these phenomena will affect the accuracy of the simulation results.

### Solutions:
The authors extended the "central string algorithm" they proposed previously and developed a quasi - 3D reduced - dimension simulation method. This method represents particles of different sizes by introducing concentric cylinders and solves the above problems through the following means:

- **Catch net**: Prevent small particles from sinking excessively and ensure that they do not pass through the gaps between large particles.
- **Synthetic surface**: Support large particles and prevent them from being stratified in the growth direction.
- **Uniform distribution**: Ensure that the number of particles in each cylinder conforms to the results of three - dimensional simulations and avoid edge effects.

### Results:
This reduced - dimension algorithm can significantly reduce the calculation time while maintaining a high degree of consistency with the three - dimensional simulation results. Especially when dealing with particle systems with large particle - size differences, the advantages of the reduced - dimension algorithm are more obvious. The experimental results show that this method can be used to study the packing characteristics of multi - modal particle systems and provides new insights into the microstructure.

In summary, the main objective of this paper is to solve the computational challenges faced by traditional three - dimensional simulations in dealing with particle systems with large particle - size differences by developing an efficient reduced - dimension simulation algorithm and to provide a new tool for the microstructure analysis of composite materials.