Open Access

ARTICLE

An Alternated Grid Updating Parallel Algorithm for Material Point Method Using OpenMP

by Yantao Zhang¹, Xiong Zhang¹, Yan Liu¹

1 School of Aerospace, Tsinghua University, Beijing 100084, China.
2 Corresponding author: xzhang@tsinghua.edu.cn

Computer Modeling in Engineering & Sciences 2010, 69(2), 143-166. https://doi.org/10.3970/cmes.2010.069.143

Abstract

Material point method(MPM) is a promising method in solving problems involving large deformations, especially explosion and penetration. In MPM, particles can move around the computing domain dynamically, which can result in load imbalance easily. In parallelizing MPM using OpenMP, data races will occur in the stage of grid node updating if we use loop-level parallelism for these loops. Huang et al. proposed a domain decomposition method to overcome data races [Huang, Zhang, Ma and Wang (2008)]. However, significant modifications of the original serial code are required. In this paper, we proposed a new alternated grid updating method to avoid data races in the stage of grid node updating, which only need small modifications of the original code and is much easier to achieve dynamic load balance. Test results show that our algorithm performs well.

---

Keywords

---