Comparison of Cohesive Models in EDEM and LIGGGHTS for Simulating Powder Compaction

The purpose of this work was to analyse the compaction of a cohesive material using different Discrete Element Method (DEM) simulators to determine the equivalent contact models and to identify how some simulation parameters affect the compaction results (maximum force and compact appearance) and computational costs. For this purpose, three cohesion contact models were tested: linear cohesion in EDEM, and simplified Johnson-Kendall-Roberts (SJKR) and modified SJKR (SJKR2) in LIGGGHTS. The influence of the particle size distribution (PSD) on the results was also investigated. Further assessments were performed on the effect of (1) selecting different timesteps, (2) using distinct conversion tolerances to export the three-dimensional models to standard triangle language (STL) files, and (3) moving the punch with different speeds. Consequently, we determined that a timestep equal to a 10% Rayleigh timestep, a conversion tolerance of 0.01 mm, and a punch speed of 0.1 m/s is adequate for simulating the compaction process using the materials and the contact models in this work. The results showed that the maximum force was influenced by the PSD due to the rearrangement of the particles. The PSD was also related to the computational cost because of the number of simulated particles and their sizes. Finally, an equivalence was found between the linear cohesion and SJKR2 contact models.