Simultaneous FETI and block FETI: Robust domain decomposition with multiple search directions

被引:34
作者
Gosselet, Pierre [1 ]
Rixen, Daniel [2 ]
Roux, Francois-Xavier [3 ]
Spillane, Nicole [4 ]
机构
[1] Univ Paris Saclay, CNRS, ENS Cachan, LMT Cachan, F-94235 Cachan, France
[2] Tech Univ Munich, Fac Mech, Inst Appl Mech, D-85748 Garching, Germany
[3] Univ Paris 06, CNRS UMR 7598, Lab Jacques Louis Lions, F-75005 Paris, France
[4] Univ Chile, Fac Ciencias Fis & Matemat, Ctr Math Modeling, Santiago, Chile
关键词
domain decomposition; FETI; BDD; block Krylov methods; multiple preconditioner; heterogeneity; PARALLEL SOLUTION; DP METHODS; PART I; PRECONDITIONERS; ALGORITHM; EXTENSION; SYSTEMS; SCHWARZ; BDDC; COEFFICIENTS;
D O I
10.1002/nme.4946
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
Domain decomposition methods often exhibit very poor performance when applied to engineering problems with large heterogeneities. In particular, for heterogeneities along domain interfaces, the iterative techniques to solve the interface problem are lacking an efficient preconditioner. Recently, a robust approach, named finite element tearing and interconnection (FETI)-generalized eigenvalues in the overlaps (Geneo), was proposed where troublesome modes are precomputed and deflated from the interface problem. The cost of the FETI-Geneo is, however, high. We propose in this paper techniques that share similar ideas with FETI-Geneo but where no preprocessing is needed and that can be easily and efficiently implemented as an alternative to standard domain decomposition methods. In the block iterative approaches presented in this paper, the search space at every iteration on the interface problem contains as many directions as there are domains in the decomposition. Those search directions originate either from the domain-wise preconditioner (in the simultaneous FETI method) or from the block structure of the right-hand side of the interface problem (block FETI). We show on two-dimensional structural examples that both methods are robust and provide good convergence in the presence of high heterogeneities, even when the interface is jagged or when the domains have a bad aspect ratio. The simultaneous FETI was also efficiently implemented in an optimized parallel code and exhibited excellent performance compared with the regular FETI method. Copyright (c) 2015 John Wiley & Sons, Ltd.
引用
收藏
页码:905 / 927
页数:23
相关论文
共 57 条
[1]  
[Anonymous], 2003, ITERATIVE METHODS SP, DOI DOI 10.1137/1.9780898718003
[2]  
Bhardwaj M, 2000, INT J NUMER METH ENG, V47, P513, DOI 10.1002/(SICI)1097-0207(20000110/30)47:1/3<513::AID-NME782>3.0.CO
[3]  
2-V
[4]   Modelling and convergence in arterial wall simulations using a parallel FETI solution strategy [J].
Brands, D. ;
Klawonn, A. ;
Rheinbach, O. ;
Schroeder, J. .
COMPUTER METHODS IN BIOMECHANICS AND BIOMEDICAL ENGINEERING, 2008, 11 (05) :569-583
[5]   BOOTSTRAP AMG [J].
Brandt, A. ;
Brannick, J. ;
Kahl, K. ;
Livshits, I. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2011, 33 (02) :612-632
[6]   A multipreconditioned conjugate gradient algorithm [J].
Bridson, R ;
Greif, C .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2006, 27 (04) :1056-1068
[7]   ISOGEOMETRIC BDDC PRECONDITIONERS WITH DELUXE SCALING [J].
da Veiga, L. Beirao ;
Pavarino, L. F. ;
Scacchi, S. ;
Widlund, O. B. ;
Zampini, S. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2014, 36 (03) :A1118-A1139
[8]   Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity [J].
Dohrmann, Clark R. ;
Widlund, Olof B. .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2010, 82 (02) :157-183
[9]   AN OVERLAPPING SCHWARZ ALGORITHM FOR ALMOST INCOMPRESSIBLE ELASTICITY [J].
Dohrmann, Clark R. ;
Widlund, Olof B. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2009, 47 (04) :2897-2923
[10]   A preconditioner for substructuring based on constrained energy minimization [J].
Dohrmann, CR .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2003, 25 (01) :246-258