A Study of Contact Detection between Noncircular Particles in Discrete Element Method

Abstract

One of the key points of modeling noncircular particles in the discrete element method is the contact detection of particles. In this study, a general contact detection algorithm of two-dimensional particles with analytic shape functions is provided. The contact detection of particles with strictly convex shape function, such as ellipses and superellipses, is solved by Newton–Raphson method, and a grid method is provided to deal with heart-shaped particles. The grid method can be generalized into a particle system, in which the shape function is not convex. The accuracy and stability of the algorithm are verified by a series of tests. For the collision of a pair of ellipses with an aspect ratio of $\alpha$  =  $a/b$  = 1000, the efficiency is not worse than the Newton–Raphson method. For random packing under gravity, no residual kinetic energy is observed, and the force that acts on the bottom is equal to the gravity, that is, the system reaches a mechanical equilibrium state. After equilibrium, the process of hopper discharge is also simulated. The present method is suitable for arbitrarily shaped particles with analytic shape functions in two-dimensional cases. In addition to superellipses, the method can also be applied to other particle systems as long as the shape functions are given.