An adaptive hierarchical domain decomposition method for parallel contact dynamics simulations of granular materials
arXiv:1104.3516 · doi:10.1016/j.jcp.2011.09.024
Abstract
A fully parallel version of the contact dynamics (CD) method is presented in this paper. For large enough systems, 100% efficiency has been demonstrated for up to 256 processors using a hierarchical domain decomposition with dynamic load balancing. The iterative scheme to calculate the contact forces is left domain-wise sequential, with data exchange after each iteration step, which ensures its stability. The number of additional iterations required for convergence by the partially parallel updates at the domain boundaries becomes negligible with increasing number of particles, which allows for an effective parallelization. Compared to the sequential implementation, we found no influence of the parallelization on simulation results.
19 pages, 15 figures, published in Journal of Computational Physics (2011)
References in corpus (8)
- Event-Driven Molecular Dynamics in Parallel
- Force indeterminacy in the jammed state of hard disks
- Refraction of shear zones in granular materials
- Shearing behavior of polydisperse media
- Extent of force indeterminacy in packings of frictional rigid disks
- Snell's Law for Shear Zone Refraction in Granular Materials
- Generation of homogeneous granular packings: Contact dynamics method with coupling to an external pressure bath
- Macroscopic and Microscopic Investigation on the History Dependence of the Mechanical Behaviour of Powders