An adaptive hierarchical domain decomposition method for parallel contact dynamics simulations of granular materials

JOURNAL OF COMPUTATIONAL PHYSICS

卷:231

Article

Shojaaee, Zahra¹  Shaebani, M. Reza¹  Brendel, Lothar¹  Toeroek, Janos¹  Wolf, Dietrich E.¹ 

[1] Univ Duisburg Essen, Dept Theoret Phys, Computat & Stat Phys Grp, D-47048 Duisburg, Germany

关键词: Contact dynamics method;    Granular materials;    Hierarchical domain decomposition;    Load balancing;    MPI library;   

DOI  :  10.1016/j.jcp.2011.09.024

来源: Elsevier

【 摘 要 】

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. (C) 2011 Elsevier Inc. All rights reserved.

【 授权许可】

Free