Molecular dynamics simulations of head-on low-velocity collisions between particles

Yuki Yoshida,Eiichiro Kokubo,Hidekazu Tanaka

DOI: https://doi.org/10.1103/PhysRevE.110.015001

2024-08-08

Abstract:The particle contact model is important for powder simulations. Although several contact models have been proposed, their validity has not yet been well established. Therefore, we perform molecular dynamics (MD) simulations to clarify the particle interaction. We simulate head-on collisions of two particles with impact velocities less than a few percent of the sound velocity to investigate the dependence of the interparticle force and the coefficient of restitution (COR) on the impact velocity and particle radius. In this study, we treat particles with a radius of 10-100 nm and perform simulations. We find that the interparticle force exhibits hysteresis between the loading and unloading phases. Larger impact velocities result in strong hysteresis and plastic deformation. For all impact velocities and particle radii, the coefficient of restitution is smaller than that given by the Johnson-Kendall-Robert theory. An inelastic contact model cannot reproduce our MD simulations. In particular, the COR is significantly reduced when the impact velocity exceeds a certain value. This significant energy dissipation cannot be explained even by the contact models including plastic deformation. We also find that the COR increases with increasing particle radius. We also find that the previous contact models including plastic deformation cannot explain the strong energy dissipation obtained in our MD simulations, although they agree with the MD results for very low impact velocities. Accordingly, we have constructed a new dissipative contact model in which the dissipative force increases with the stress generated by collisions. The new stress dependent model successfully reproduces our MD results over a wider range of impact velocities than the conventional models do. In addition, we proposed another, simpler, dissipative contact model that can also reproduce the MD results.

Soft Condensed Matter,Earth and Planetary Astrophysics,Astrophysics of Galaxies

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that the validity of the existing particle contact models has not been fully verified. Specifically, through molecular dynamics (MD) simulations, the paper studies the dependence of the inter - particle force and the coefficient of restitution (COR) on the impact velocity and particle radius during the low - speed frontal collision of two particles. The particle radii used in the study are between 10 - 100 nanometers, and simulations with up to 20 million atoms were carried out. The study found that the inter - particle force exhibits hysteresis during the loading and unloading stages, and a higher impact velocity will lead to stronger hysteresis and plastic deformation. For all impact velocities and particle radii, the experimentally measured coefficient of restitution is less than the value predicted by the Johnson - Kendall - Robert (JKR) theory, and the existing elastic particle contact models cannot reproduce the MD simulation results. In particular, when the impact velocity exceeds a certain threshold, the coefficient of restitution decreases significantly, and this significant energy dissipation cannot be explained by the contact models including plastic deformation either. In addition, the coefficient of restitution increases as the particle radius increases. Based on these findings, the author constructs a new dissipative contact model in which the dissipative force increases with the stress generated by the collision, and this model can successfully reproduce the MD simulation results in a wider range of impact velocities. In addition, another simpler dissipative contact model is also proposed, which can also reproduce the MD simulation results.

The main objective of the paper is to explore and verify the dynamic behaviors of particle collisions under different conditions, especially the energy dissipation mechanism, through detailed MD simulations, and to propose improved contact models to better describe these processes. This helps to improve the accuracy of powder simulations, especially in the research of particle aggregation, growth and fragmentation processes in fields such as mineral processing engineering, fertilizer industry and astrophysics.