{"title":"Molecular dynamics simulations of head-on low-velocity collisions between particles","authors":"Yuki Yoshida, Eiichiro Kokubo, Hidekazu Tanaka","doi":"arxiv-2408.04164","DOIUrl":null,"url":null,"abstract":"The particle contact model is important for powder simulations. Although\nseveral contact models have been proposed, their validity has not yet been well\nestablished. Therefore, we perform molecular dynamics (MD) simulations to\nclarify the particle interaction. We simulate head-on collisions of two\nparticles with impact velocities less than a few percent of the sound velocity\nto investigate the dependence of the interparticle force and the coefficient of\nrestitution (COR) on the impact velocity and particle radius. In this study, we\ntreat particles with a radius of 10-100 nm and perform simulations. We find\nthat the interparticle force exhibits hysteresis between the loading and\nunloading phases. Larger impact velocities result in strong hysteresis and\nplastic deformation. For all impact velocities and particle radii, the\ncoefficient of restitution is smaller than that given by the\nJohnson-Kendall-Robert theory. An inelastic contact model cannot reproduce our\nMD simulations. In particular, the COR is significantly reduced when the impact\nvelocity exceeds a certain value. This significant energy dissipation cannot be\nexplained even by the contact models including plastic deformation. We also\nfind that the COR increases with increasing particle radius. We also find that\nthe previous contact models including plastic deformation cannot explain the\nstrong energy dissipation obtained in our MD simulations, although they agree\nwith the MD results for very low impact velocities. Accordingly, we have\nconstructed a new dissipative contact model in which the dissipative force\nincreases with the stress generated by collisions. The new stress dependent\nmodel successfully reproduces our MD results over a wider range of impact\nvelocities than the conventional models do. In addition, we proposed another,\nsimpler, dissipative contact model that can also reproduce the MD results.","PeriodicalId":501146,"journal":{"name":"arXiv - PHYS - Soft Condensed Matter","volume":"3 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Soft Condensed Matter","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2408.04164","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
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.