In this paper, a parallel overlapping domain decomposition preconditioner is proposed to solve the linear system of equations arising from the extended finite element discretization of elastic crack problems. The algorithm partitions the computational mesh into two types of subdomains: the regular subdomains and the crack tip subdomains based on the observation that the crack tips have a significant impact on the convergence of the iterative method while the impact of the crack lines is not that different from those of regular mesh points. The tip subdomains consist of mesh points at crack tips and all neighboring points where the branch enrichment functions are applied. The regular subdomains consist of all other mesh points, including those on the crack lines. To overcome the mismatch between the number of subdomains and the number of processor cores, the proposed method is divided into two steps: solve the crack tip problem and then the regular subdomain problem during each iteration. The proposed method was used to develop a parallel XFEM package which is able to test different types of iterative methods. To achieve good parallel efficiency, additional methods were introduced to reduce communication and to maintain the load balance between processors. Numerical experiments indicate that the proposed method significantly reduces the number of iterations and the total computation time compared to the classical methods. In addition, the method scales up to 8192 processor cores with over 70\% parallel efficiency to solve problems with more than $2\times10^8$ degrees of freedom.
翻译:在本文中, 提出一个平行重叠的域分解先决条件, 以解决由弹性裂缝问题的延长有限元素分解产生的等式线性系统的线性系统。 算法将计算网格分割成两种类型的子域: 常规子域和裂缝尖子域子域, 其依据的观察是, 裂缝对迭接合方法的趋同有重大影响, 而裂缝线的影响与常规网点的影响并不不同。 提示子次域由裂缝点和分支浓缩功能应用的所有相邻点的网状点组成。 常规子域由所有其他网格点组成, 包括裂缝线上的网格点。 为了克服子域数目和进程核心数之间的不匹配, 拟议的方法分为两个步骤: 解决裂隙点问题, 然后是每次断裂线的正常次域问题。 提议的方法是开发一个平行的 XFEM 包, 它可以测试不同类型的迭接方法。 为了实现良好的平行效率, 额外的子网格包括所有其他网格点, 包括裂线上的网格点 。 为了克服子域间间间断径界间数与进程之间的错错错错错,,, 提议的方法可以大大地缩小计算方法到缩缩取。 。 。 。 。 。 将计算法的计算法的计算法的计算方法比平序法 。 。