Table of Content

Open AccessOpen Access


Shared Memory OpenMP Parallelization of Explicit MPM and Its Application to Hypervelocity Impact

P. Huang1,2, X. Zhang1,3, S. Ma1, H.K. Wang1

School of Aerospace, Tsinghua University, Beijing 100084, China
Institute of System Engineering, China Academy of Engineering Physics, Mianyang 621900, China
Correspondence author. E-mail:

Computer Modeling in Engineering & Sciences 2008, 38(2), 119-148.


The material point method (MPM) is an extension of particle-in-cell method to solid mechanics. A parallel MPM code is developed using FORTRAN 95 and OpenMP in this study, which is designed primarily for solving impact dynamic problems. Two parallel methods, the array expansion method and the domain decomposition method, are presented to avoid data races in the nodal update stage. In the array expansion method, two-dimensional auxiliary arrays are created for nodal variables. After updating grid nodes in all threads, the auxiliary arrays are assembled to establish the global nodal array. In the domain decomposition method, the background grid is decomposed into some uniform patches, and each thread deals with a patch. The information of neighbor patches is exchanged through shared variables. After updating nodes in all patches, their nodal variables are assembled to establish the global nodal variables. The numerical tests show that the domain decomposition method has much better parallel scalability and higher parallel efficiency than the array expansion method. Therefore, a parallel computer code, MPM3DMP, is developed based on the domain decomposition method. Finally, MPM3DMP is applied to a large-scale simulation with 13,542,030 particles for obtaining the high-resolution results of debris cloud in hypervelocity impact.


Cite This Article

Huang, P., Zhang, X., Ma, S., Wang, H. (2008). Shared Memory OpenMP Parallelization of Explicit MPM and Its Application to Hypervelocity Impact. CMES-Computer Modeling in Engineering & Sciences, 38(2), 119–148.

This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • 1332


  • 725


  • 0


Share Link