Scalable computational kernels for mortar finite element methods

被引:4
|
作者
Mayr, Matthias [1 ,2 ]
Popp, Alexander [1 ]
机构
[1] Univ Bundeswehr Munchen, Inst Math & Comp Based Simulat, Werner-Heisenberg-Weg 39, D-85577 Neubiberg, Germany
[2] Univ Bundeswehr Munchen, Data Sci & Comp Lab, Werner-Heisenberg-Weg 39, D-85577 Neubiberg, Germany
关键词
Mortar methods; Contact mechanics; Interface problems; Parallel algorithms; Finite elements; Domain decomposition; CONTACT SEARCHING ALGORITHM; FLUID-STRUCTURE INTERACTION; ACTIVE SET STRATEGY; ISOGEOMETRIC ANALYSIS; DEFORMATION CONTACT; CURVED INTERFACES; DYNAMIC CONTACT; 3D; FORMULATION; INTEGRATION;
D O I
10.1007/s00366-022-01779-3
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Targeting simulations on parallel hardware architectures, this paper presents computational kernels for efficient computations in mortar finite element methods. Mortar methods enable a variationally consistent imposition of coupling conditions at high accuracy, but come with considerable numerical effort and cost for the evaluation of the mortar integrals to compute the coupling operators. In this paper, we identify bottlenecks in parallel data layout and domain decomposition that hinder an efficient evaluation of the mortar integrals. We then propose a set of computational strategies to restore optimal parallel communication and scalability for the core kernels devoted to the evaluation of mortar terms. We exemplarily study the proposed algorithmic components in the context of three-dimensional large-deformation contact mechanics, both for cases with fixed and dynamically varying interface topology, yet these concepts can naturally and easily be transferred to other mortar applications, e.g. classical meshtying problems. To restore parallel scalability, we employ overlapping domain decompositions of the interface discretization independent from the underlying volumes and then tackle parallel communication for the mortar evaluation by a geometrically motivated reduction of ghosting data. Using three-dimensional contact examples, we demonstrate strong and weak scalability of the proposed algorithms up to 480 parallel processes as well as study and discuss improvements in parallel communication related to mortar finite element methods. For the first time, dynamic load balancing is applied to mortar contact problems with evolving contact zones, such that the computational work is well balanced among all parallel processors independent of the current state of the simulation.
引用
收藏
页码:3691 / 3720
页数:30
相关论文
共 50 条