Numerical comparison of some contact detection algorithms

Contact interaction and contact detection (CD) remain key components of any discontinua simulations. The methods of discontinua include combined finite-discrete element method (FDEM), discrete element method, molecular dynamics, etc. In recent years, a number of CD algorithms have been developed, such as Munjiza–Rougier (MR), Munjiza–Rougier–Schiava (MR-S), Munjiza-No Binary Search (NBS), Balanced Binary Tree Schiava (BBTS), 3D Discontinuous Deformation Analysis and many others. This work aims to conduct a numerical comparison of certain algorithms often used in FDEM for bodies of the same size. These include MR, MR-S, NBS and BBTS algorithms.

Computational simulations were used in this work.

In discrete element simulations where particles are introduced randomly or in which the relative position between particles is constantly changing, the MR and MR-S algorithms present an advantage in terms of CD times.

This paper presents a detailed comparison between CD algorithms. The comparisons are performed for problem cases with different lattices and distributions of particles in discrete element simulations. The comparison includes algorithms that have not been evaluated between them. Also, two new algorithms are presented in the paper, MR-S and BBTS.