基于静态双重区域分解的两种接触并行算法

CHAP3D是北京应用物理与计算数学研究所自主研发的Lagrange通用弹塑性流体力学分析程序.文章介绍了在CHAP3D程序中使用的、针对多处理器集群的、基于静态双重区域分解的两种接触并行算法.第一种是分配单个完整接触面的接触并行算法,此算法将一对完整的接触面分配到一个处理器上,并建立计算域与接触域的通信关系.此接触并行算法的优点是简单,在具有接触面的处理器上可以直接使用串行的接触搜索算法和接触力耦合计算算法.另一种是主面剖分区域分解的接触并行算法,此算法将所有接触面的主面区域分解到所有处理器上.须建立计算域与接触域以及接触域内各处理器间的两种通信关系.该接触并行算法是一个负载平衡的并行算法,具有很好的并行效率和可扩展性.数值算例显示,这两种接触并行算法都能够很好地模拟多种不同类型的接触问题.

Two Types of Contact Parallel Algorithms Based on the Dual Static Domain Decomposition

CHAP3D code is a three-dimensional elastic-plastic Lagrangian program developed by IAPCM. The contact parallel algorithms used in CHAP3D code for the clusters of multi-processors were introduced in this paper. The dual domain decomposition method and two strategies of contact domain decomposition were adopted in these parallel contact algorithms. One is to distribute a pair of whole contact surfaces to a processor, and a communication between the computational domain and the contact domain was established. The serial contact search algorithm and the contact forces algorithm can be applied directly. The other is to divide the master surface mesh into every processor. Two types of communications, the communication between the computational domain and the contact domain and the communication among all the processors in the contact domain, should be constructed in this contact parallel algorithm. The contact parallel algorithm is a scalable parallel algorithm with load-balance. The numerical examples illustrate that these two types of contact parallel algorithms have the abilities to simulate different types of contact problems.