×

Simultaneous FETI and block FETI: robust domain decomposition with multiple search directions. (English) Zbl 1352.65625

Summary: 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.

MSC:

65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs
65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs

Software:

METIS; BDDCML; BDDC; PARDISO

References:

[1] FarhatC, RouxFX. A method of finite element tearing and interconnecting and its parallel solution algorithm. International Journal for Numerical Methods in Engineering1991; 32(6):1205. · Zbl 0758.65075
[2] MandelJ. Balancing domain decomposition. Communications in Numerical Methods in Engineering1993; 9(3):233-241. · Zbl 0796.65126
[3] BhardwajM, DayD, FarhatC, LesoinneM, PiersonK, RixenD. Application of the FETI method to ASCI problems: scalability results on a thousand‐processor and discussion of highly heterogeneous problems. International Journal for Numerical Methods in Engineering2000; 47(1‐3):513-536. · Zbl 0970.74069
[4] BrandsD, KlawonnA, RheinbachO, SchröderJ. Modelling and convergence in arterial wall simulations using a parallel FETI solution strategy. Computer Methods in Applied Mechanics and Engineering2008; 11(5):569-583.
[5] GosseletP, ReyC. Non‐overlapping domain decomposition methods in structural mechanics. Archives of Computational Methods in Engineering2006; 13(4):515-572. · Zbl 1171.74041
[6] FarhatC, RixenD. A new coarsening operator for the optimal preconditioning of the dual and primal domain decomposition methods: application to problems with severe coefficient jumps. In Proceedings of the Seventh Copper Mountain Conference on Multigrid Methods, Duane MelsonN (ed.), McCormickSF (ed.), ManteuffelTA (ed.), DouglasMCC (ed.) (eds).: NASA, Langley Virginia US, 1995; 301-316.
[7] RixenD, FarhatC. A simple and efficient extension of a class of substructure based preconditioners to heterogeneous structural mechanics problems. International Journal for Numerical Methods in Engineering1999; 44(4):489-516. · Zbl 0940.74067
[8] KlawonnA, RheinbachO. Robust FETI‐DP methods for heterogeneous three dimensional elasticity problems. Computer Methods in Applied Mechanics and Engineering2007; 196(8):1400-1414. · Zbl 1173.74428
[9] KlawonnA, RheinbachO, WidlundOB. An analysis of a FETI-DP algorithm on irregular subdomains in the plane. SIAM Journal on Numerical Analysis2008; 46(5):2484-2504. · Zbl 1176.65135
[10] DryjaM, SarkisMV, WidlundOB. Multilevel Schwarz methods for elliptic problems with discontinuous coefficients in three dimensions. Numerische Mathematik1996; 72(3):313-348. · Zbl 0857.65131
[11] MandelJ, BrezinaM. Balancing domain decomposition for problems with large jumps in coefficients. Mathematics of Computation1996; 65(216):1387-1401. · Zbl 0853.65129
[12] DohrmannCR, WidlundOB. An overlapping Schwarz algorithm for almost incompressible elasticity. SIAM Journal on Numerical Analysis2009; 47(4):2897-2923. · Zbl 1410.74064
[13] DohrmannCR, WidlundOB. Hybrid domain decomposition algorithms for compressible and almost incompressible elasticity. International Journal for Numerical Methods in Engineering2010; 82(2):157-183. · Zbl 1188.74053
[14] PechsteinC, ScheichlR. Scaling up through domain decomposition. Applicable Analysis2009; 88(10‐11):1589-1608. · Zbl 1281.76040
[15] PechsteinC, ScheichlR. Analysis of FETI methods for multiscale PDEs. Numerische Mathematik2008; 111(2):293-333. · Zbl 1170.65097
[16] GippertS, KlawonnA, RheinbachO. Analysis of FETI‐DP and BDDC for linear elasticity in 3D with almost incompressible components and varying coefficients inside subdomains. SIAM Journal on Numerical Analysis2012; 50(5):2208-2236. · Zbl 1255.74063
[17] PechsteinC, ScheichlR. Analysis of FETI methods for multiscale PDEs. Part II: interface variation. Numerische Mathematik2011; 118(3):485-529. · Zbl 1380.65388
[18] PechsteinC, SarkisM, ScheichlR. New theoretical coefficient robustness results for FETI‐DP. In Domain Decomposition Methods in Science and Engineering XX, Lecture Notes in Computational Science and Engineering. Springer‐Verlag, Heidelberg, 2011; 327-335.
[19] PechsteinC. Finite and Boundary Element Tearing and Interconnecting Solvers for Multiscale Problems, Lecture Notes in Computational Science and Engineering, vol. 90. Springer‐Verlag: Berlin Heidelberg, 2013. · Zbl 1272.65100
[20] PechsteinC, ScheichlR. Weighted Poincaré inequalities. IMA Journal on Numerical Analysis2013; 33(2):652-686. · Zbl 1275.26030
[21] Beirão da VeigaL, PavarinoLF, ScacchiS, WidlundOB, ZampiniS. Isogeometric BDDC preconditioners with deluxe scaling. SIAM Journal of Scientific Computing2014; 36(3):A1118-A1139. · Zbl 1320.65047
[22] NatafF, XiangH, DoleanV, SpillaneN. A coarse space construction based on local Dirichlet‐to‐Neumann maps. SIAM Journal of Scientific Computing2011; 33(4):1623-1642. · Zbl 1230.65134
[23] DoleanV, NatafF, ScheichlR, SpillaneN. Analysis of a two‐level Schwarz method with coarse spaces based on local Dirichlet‐to‐Neumann maps. Computational Methods in Applied Mathematics2012; 12(4):391-414. · Zbl 1284.65050
[24] SpillaneN, DoleanV, HauretP, NatafF, PechsteinC, ScheichlR. Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps. Numerische Mathematik2014; 126(4):741-770. · Zbl 1291.65109
[25] GalvisJ, EfendievY. Domain decomposition preconditioners for multiscale flows in high contrast media: reduced dimension coarse spaces. Multiscale Modeling and Simulation2010; 8(5):1621-1644. · Zbl 1381.65029
[26] EfendievY, GalvisJ, LazarovR, WillemsJ. Robust domain decomposition preconditioners for abstract symmetric positive definite bilinear forms. ESAIM Mathematical Modelling and Numerical Analysis2012; 46(5):1175-1199. · Zbl 1272.65098
[27] MandelJ, SousedíkB. Adaptive selection of face coarse degrees of freedom in the BDDC and the FETI‐DP iterative substructuring methods. Computer Methods in Applied Mechanics and Engineering2007; 196(8):1389-1399. · Zbl 1173.74435
[28] SousedíkB, ŠístekJ, MandelJ. Adaptive‐multilevel BDDC and its parallel implementation. Computing2013; 95(12):1087-1119. · Zbl 1307.65175
[29] SpillaneN, RixenDJ. Automatic spectral coarse spaces for robust FETI and BDD algorithms. International Journal for Numerical Methods in Engineering2013; 95(11):953-990. · Zbl 1352.65553
[30] KlawonnA, RadtkeP, RheinbachO. Adaptive coarse spaces for BDDC with a transformation of basis. Twenty Second International Conference on Domain Decomposition Methods, Università della Svizzera Italiana, Lugano, Switzerland, 2014.
[31] KlawonnA, RadtkeP, RheinbachO. FETI‐DP methods with an adaptive coarse space. SIAM Journal on Numerical Analysis2015; 53(1):297-320. · Zbl 1327.65063
[32] SpillaneN. Robust domain decomposition methods for symmetric positive definite problems. Ph.D. Thesis, Thèse de l’Ecole doctorale de Mathmatiques de Paris centre, Laboratoire Jacques Louis Lions, Université Pierre et Marie Curie, Paris, 2014.
[33] RixenD. Substructuring and dual methods in structural analysis. Ph.D. Thesis, Université de Liège, Belgium, Collection des Publications de la Faculté des Sciences appliquées, n.175, 1997.
[34] RixenDJ. A domain decomposition interface solver with multiple direction of descent for heterogeneous problems. 84th Annual Meeting of the International Association of Applied Mathematics and Mechanics, GAMM, Novi‐Sad, Serbia (oral presentation only), University of Novi Sad, Serbia, 2013.
[35] GosseletP, RixenDJ, ReyC. A domain decomposition strategy to efficiently solve structures containing repeated patterns. International Journal for Numerical Methods in Engineering2009; 78(7):828-842. · Zbl 1183.74277
[36] GreifC, ReesT, SzyldDB. MPGMRES: a generalized minimum residual method with multiple preconditioners. Technical Report 11‐12‐23, Department of Mathematics, Temple University, 2011. Revised September 2012 and January 2014. Also available as Technical Report TR‐2011‐12, Department of Computer Science, University of British Columbia.
[37] O’LearyDP. The block conjugate gradient algorithm and related methods. Linear Algebra and its Applications1980; 29:293-322. · Zbl 0426.65011
[38] FarhatC, PiersonK, LesoineM. The second generation FETI methods and their application to the parallel solution of large‐scale linear and geometrically non‐linear structural analysis problems. Computer Methods in Applied Mechanics and Engineering2000; 184(2‐4):333-374. · Zbl 0981.74064
[39] FarhatC, LesoinneM, LeTallecP, PiersonK, RixenD. FETI‐DP: a dual‐primal unified FETI method ‐ part I: a faster alternative to the two‐level FETI method. International Journal for Numerical Methods in Engineering2001; 50(7):1523-1544. · Zbl 1008.74076
[40] DohrmannCR. A preconditioner for substructuring based on constrained energy minimization. SIAM Journal of Scientific Computing2003; 25(1):246-258 (electronic). · Zbl 1038.65039
[41] RixenDJ, FarhatC, TezaurR, MandelJ. Theoretical comparison of the FETI and algebraically partitioned FETI methods, and performance comparisons with a direct sparse solver. International Journal for Numerical Methods in Engineering1999; 46(4):501-533. · Zbl 0977.74065
[42] FarhatC, MandelJ. The two‐level FETI method for static and dynamic plate problems. I. An optimal iterative solver for biharmonic systems. Computer Methods in Applied Mechanics and Engineering1998; 155(1‐2):129-151. · Zbl 0964.74062
[43] FarhatC, ChenPS, MandelJ, RouxFX. The two‐level FETI method. II. Extension to shell problems, parallel implementation and performance results. Computer Methods in Applied Mechanics and Engineering1998; 155(1‐2):153-179. · Zbl 1040.74513
[44] GosseletP, ReyC, PebrelJ. Total and selective reuse of Krylov subspaces for the solution to a sequence of nonlinear structural problems. International Journal for Numerical Methods in Engineering2013; 94(1):60-83. · Zbl 1352.65105
[45] ReyV, ReyC, GosseletP. A strict error bound with separated contributions of the discretization and of the iterative solver in non‐overlapping domain decomposition methods. Computer Methods in Applied Mechanics and Engineering2013; 270(1):293-303. · Zbl 1296.65167
[46] BridsonR, GreifC. A multipreconditioned conjugate gradient algorithm. SIAM Journal on Matrix Analysis and Applications2006; 27(4):1056-1068 (electronic). · Zbl 1104.65027
[47] SaadY. Iterative Methods for Sparse Linear Systems (2nd edn). Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, 2003. · Zbl 1031.65046
[48] GreifC, ReesT, SzyldDB. Additive Schwarz with variable weights. In Domain Decomposition Methods in Science and Engineering XXI, vol. 98, ErhelJ (ed.), GanderMJ (ed.), HalpernL (ed.), PichotG (ed.), SassiT (ed.), WidlundO (ed.) (eds)., Lecture Notes in Computational Science and Engineering. Springer International Publishing: Cham, Switzerland, 2014; 779-787. · Zbl 1382.65085
[49] NikishinAA, YereminAY. Variable block CG algorithms for solving large sparse symmetric positive definite linear systems on parallel computers. I. General iterative scheme. SIAM Journal on Matrix Analysis and Applications1995; 16(4):1135-1153. · Zbl 0837.65029
[50] BrandtA, BrannickJ, KahlK, LivshitsI. Bootstrap AMG. SIAM Journal of Scientific Computing2011; 33(2):612-632. · Zbl 1227.65120
[51] GosseletP, ReyC, RixenD. On the initial estimate of interface forces in FETI methods. Computer Methods in Applied Mechanics and Engineering2003; 192(25):2749-2764. · Zbl 1054.74728
[52] KarypisG, KumarV. METIS: a software package for partitioning unstructured graphs, partitioning meshes, and computing fill‐reducing orderings of sparse matrices, Department of Computer Science, University of Minnesota, 1998. http://glaros.dtc.umn.edu/gkhome/views/metis.
[53] VereeckeB, BavestrelloH, DureisseixD. An extension of the FETI domain decomposition method for incompressible and nearly incompressible problems. Computer Methods in Applied Mechanics and Engineering2003; 192(31):3409-3429. · Zbl 1054.74739
[54] ZsakiA, RixenDJ, ParaschivoiuM. A substructured based iterative inner solver coupled with Uzawa’s algorithm for the Stokes problem. International Journal for Numerical Methods in Fluids2003; 43:215-230. · Zbl 1032.76608
[55] SchenkO, GartnerK. Solving unsymmetric sparse systems of linear equations with PARDISO. Journal of Future Generation Computer Systems2004; 20(3):475-487.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.