Block FETI–DP/BDDC preconditioners for mixed isogeometric discretizations of three-dimensional almost incompressible elasticity

Widlund, O.; Zampini, S.; Scacchi, S.; Pavarino, L.

doi:10.1090/mcom/3614

Block FETI–DP/BDDC preconditioners for mixed isogeometric discretizations of three-dimensional almost incompressible elasticity
HTML articles powered by AMS MathViewer

by O. B. Widlund, S. Zampini, S. Scacchi and L. F. Pavarino HTML | PDF

Math. Comp. 90 (2021), 1773-1797 Request permission

Abstract:

A block FETI–DP/BDDC preconditioner for mixed formulations of almost incompressible elasticity is constructed and analyzed; FETI–DP (Finite Element Tearing and Interconnecting Dual-Primal) and BDDC (Balancing Domain Decomposition by Constraints) are two very successful domain decomposition algorithms for a variety of elliptic problems. The saddle point problems of the mixed problems are discretized with mixed isogeometric analysis with continuous pressure fields. As in previous work by Tu and Li (2015) for finite element discretizations of the incompressible Stokes system, the proposed preconditioner is applied to a reduced positive definite system involving only the pressure interface variable and the Lagrange multiplier of the FETI–DP algorithm. In this work, we extend the theory to a wider class of saddle point problems and we propose a novel block-preconditioning strategy, which consists in using BDDC with deluxe scaling for the interface pressure block as well as deluxe scaling for the FETI–DP preconditioner for the Lagrange multiplier block. A convergence rate analysis is presented with a condition number bound for the preconditioned operator which depends on the inf-sup parameter of the fully assembled problem and the condition number of a closely related BDDC algorithm for compressible elasticity. This bound is scalable in the number of subdomains, poly-logarithmic in the ratio of subdomain and element sizes, and robust with respect to material incompressibility. Parallel numerical experiments validate the theory, demonstrate robustness in the presence of discontinuities of the Lamé parameters, and indicate how the rate of convergence varies with respect to the spline polynomial degree and regularity and the deformation of the domain.

References

Patrick R. Amestoy, Iain S. Duff, Jean-Yves L’Excellent, and Jacko Koster, A fully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM J. Matrix Anal. Appl. 23 (2001), no. 1, 15–41. MR 1856597, DOI 10.1137/S0895479899358194
S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, A. Dener, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, D. A. May, L. Curfman McInnes, R. T. Mills, T. Munson, K. Rupp, P. Sanan, B. F. Smith, S. Zampini, and H. Zhang. PETSc Web page, http://www.mcs.anl.gov/petsc, 2018.
Y. Bazilevs, L. Beirão da Veiga, J. A. Cottrell, T. J. R. Hughes, and G. Sangalli, Isogeometric analysis: approximation, stability and error estimates for $h$-refined meshes, Math. Models Methods Appl. Sci. 16 (2006), no. 7, 1031–1090. MR 2250029, DOI 10.1142/S0218202506001455
L. Beirão da Veiga, D. Cho, L. F. Pavarino, and S. Scacchi, Overlapping Schwarz methods for isogeometric analysis, SIAM J. Numer. Anal. 50 (2012), no. 3, 1394–1416. MR 2970748, DOI 10.1137/110833476
L. Beirão da Veiga, D. Cho, L. F. Pavarino, and S. Scacchi, BDDC preconditioners for isogeometric analysis, Math. Models Methods Appl. Sci. 23 (2013), no. 6, 1099–1142. MR 3037302, DOI 10.1142/S0218202513500048
L. Beirão da Veiga, D. Cho, L. F. Pavarino, and S. Scacchi, Isogeometric Schwarz preconditioners for linear elasticity systems, Comput. Methods Appl. Mech. Engrg. 253 (2013), 439–454. MR 3002803, DOI 10.1016/j.cma.2012.10.011
L. Beirão da Veiga, L. F. Pavarino, S. Scacchi, O. B. Widlund, and S. Zampini, Isogeometric BDDC preconditioners with deluxe scaling, SIAM J. Sci. Comput. 36 (2014), no. 3, A1118–A1139. MR 3216651, DOI 10.1137/130917399
L. Beirão da Veiga, L. F. Pavarino, S. Scacchi, O. B. Widlund, and S. Zampini, Adaptive selection of primal constraints for isogeometric BDDC deluxe preconditioners, SIAM J. Sci. Comput. 39 (2017), no. 1, A281–A302. MR 3612901, DOI 10.1137/15M1054675
Hani Benhassine and Abderrahmane Bendali, A non-overlapping domain decomposition method for continuous-pressure mixed finite element approximations of the Stokes problem, ESAIM Math. Model. Numer. Anal. 45 (2011), no. 4, 675–696. MR 2804655, DOI 10.1051/m2an/2010070
Daniele Boffi, Franco Brezzi, and Michel Fortin, Mixed finite element methods and applications, Springer Series in Computational Mathematics, vol. 44, Springer, Heidelberg, 2013. MR 3097958, DOI 10.1007/978-3-642-36519-5
Andrea Bressan and Giancarlo Sangalli, Isogeometric discretizations of the Stokes problem: stability analysis by the macroelement technique, IMA J. Numer. Anal. 33 (2013), no. 2, 629–651. MR 3047946, DOI 10.1093/imanum/drr056
J. Austin Cottrell, Thomas J. R. Hughes, and Yuri Bazilevs, Isogeometric analysis, John Wiley & Sons, Ltd., Chichester, 2009. Toward integration of CAD and FEA. MR 3618875, DOI 10.1002/9780470749081
L. Dalcin, N. Collier, P. Vignal, A. M. A. Côrtes, and V. M. Calo, PetIGA: a framework for high-performance isogeometric analysis, Comput. Methods Appl. Mech. Engrg. 308 (2016), 151–181. MR 3522274, DOI 10.1016/j.cma.2016.05.011
Clark R. Dohrmann, A preconditioner for substructuring based on constrained energy minimization, SIAM J. Sci. Comput. 25 (2003), no. 1, 246–258. MR 2047204, DOI 10.1137/S1064827502412887
Clark R. Dohrmann and Olof B. Widlund, Some recent tools and a BDDC algorithm for 3D problems in $H(\textrm {curl})$, Domain decomposition methods in science and engineering XX, Lect. Notes Comput. Sci. Eng., vol. 91, Springer, Heidelberg, 2013, pp. 15–25. MR 3242973, DOI 10.1007/978-3-642-35275-1_{2}
Clark R. Dohrmann and Olof B. Widlund, A BDDC algorithm with deluxe scaling for three-dimensional $H(\textbf {curl})$ problems, Comm. Pure Appl. Math. 69 (2016), no. 4, 745–770. MR 3465088, DOI 10.1002/cpa.21574
Marco Donatelli, Carlo Garoni, Carla Manni, Stefano Serra-Capizzano, and Hendrik Speleers, Symbol-based multigrid methods for Galerkin B-spline isogeometric analysis, SIAM J. Numer. Anal. 55 (2017), no. 1, 31–62. MR 3592079, DOI 10.1137/140988590
Charbel Farhat, Michel Lesoinne, Patrick LeTallec, Kendall Pierson, and Daniel Rixen, FETI-DP: a dual-primal unified FETI method. I. A faster alternative to the two-level FETI method, Internat. J. Numer. Methods Engrg. 50 (2001), no. 7, 1523–1544. MR 1813746, DOI 10.1002/nme.76
K. P. S. Gahalaut, J. K. Kraus, and S. K. Tomar, Multigrid methods for isogeometric discretization, Comput. Methods Appl. Mech. Engrg. 253 (2013), 413–425. MR 3002801, DOI 10.1016/j.cma.2012.08.015
Paola Gervasio, Luca Dedè, Ondine Chanon, and Alfio Quarteroni, A computational comparison between isogeometric analysis and spectral element methods: accuracy and spectral properties, J. Sci. Comput. 83 (2020), no. 1, Paper No. 18, 45. MR 4081217, DOI 10.1007/s10915-020-01204-1
Paulo Goldfeld, Luca F. Pavarino, and Olof B. Widlund, Balancing Neumann-Neumann preconditioners for mixed approximations of heterogeneous problems in linear elasticity, Numer. Math. 95 (2003), no. 2, 283–324. MR 2001079, DOI 10.1007/s00211-002-0450-9
R. Haferssas, P. Jolivet, and F. Nataf, An additive Schwarz method type theory for Lions’s algorithm and a symmetrized optimized restricted additive Schwarz method, SIAM J. Sci. Comput. 39 (2017), no. 4, A1345–A1365. MR 3679920, DOI 10.1137/16M1060066
Christoph Hofer, Analysis of discontinuous Galerkin dual-primal isogeometric tearing and interconnecting methods, Math. Models Methods Appl. Sci. 28 (2018), no. 1, 131–158. MR 3737080, DOI 10.1142/S0218202518500045
Christoph Hofer, Parallelization of continuous and discontinuous Galerkin dual-primal isogeometric tearing and interconnecting methods, Comput. Math. Appl. 74 (2017), no. 7, 1607–1625. MR 3706620, DOI 10.1016/j.camwa.2017.06.051
Christoph Hofer and Ulrich Langer, Dual-primal isogeometric tearing and interconnecting solvers for multipatch dG-IgA equations, Comput. Methods Appl. Mech. Engrg. 316 (2017), 2–21. MR 3610089, DOI 10.1016/j.cma.2016.03.031
Clemens Hofreither and Stefan Takacs, Robust multigrid for isogeometric analysis based on stable splittings of spline spaces, SIAM J. Numer. Anal. 55 (2017), no. 4, 2004–2024. MR 3686804, DOI 10.1137/16M1085425
Clemens Hofreither, A black-box low-rank approximation algorithm for fast matrix assembly in isogeometric analysis, Comput. Methods Appl. Mech. Engrg. 333 (2018), 311–330. MR 3771898, DOI 10.1016/j.cma.2018.01.014
Hyea Hyun Kim and Chang-Ock Lee, A two-level nonoverlapping Schwarz algorithm for the Stokes problem: numerical study, Comput. Methods Appl. Mech. Engrg. 223/224 (2012), 153–160. MR 2917487, DOI 10.1016/j.cma.2012.02.024
Axel Klawonn and Luca F. Pavarino, Overlapping Schwarz methods for mixed linear elasticity and Stokes problems, Comput. Methods Appl. Mech. Engrg. 165 (1998), no. 1-4, 233–245. MR 1663528, DOI 10.1016/S0045-7825(98)00059-0
Axel Klawonn and Olof B. Widlund, FETI and Neumann-Neumann iterative substructuring methods: connections and new results, Comm. Pure Appl. Math. 54 (2001), no. 1, 57–90. MR 1787107, DOI 10.1002/1097-0312(200101)54:1<57::AID-CPA3>3.0.CO;2-D
Axel Klawonn and Olof B. Widlund, Dual-primal FETI methods for linear elasticity, Comm. Pure Appl. Math. 59 (2006), no. 11, 1523–1572. MR 2254444, DOI 10.1002/cpa.20156
Stefan K. Kleiss, Clemens Pechstein, Bert Jüttler, and Satyendra Tomar, IETI—isogeometric tearing and interconnecting, Comput. Methods Appl. Mech. Engrg. 247/248 (2012), 201–215. MR 2975563, DOI 10.1016/j.cma.2012.08.007
Jing Li and Xuemin Tu, A nonoverlapping domain decomposition method for incompressible Stokes equations with continuous pressures, SIAM J. Numer. Anal. 51 (2013), no. 2, 1235–1253. MR 3045654, DOI 10.1137/120861503
Jing Li and Olof B. Widlund, FETI-DP, BDDC, and block Cholesky methods, Internat. J. Numer. Methods Engrg. 66 (2006), no. 2, 250–271. MR 2224479, DOI 10.1002/nme.1553
Jan Mandel and Clark R. Dohrmann, Convergence of a balancing domain decomposition by constraints and energy minimization, Numer. Linear Algebra Appl. 10 (2003), no. 7, 639–659. Dedicated to the 70th birthday of Ivo Marek. MR 2030628, DOI 10.1002/nla.341
Duk-Soon Oh, Olof B. Widlund, Stefano Zampini, and Clark R. Dohrmann, BDDC algorithms with deluxe scaling and adaptive selection of primal constraints for Raviart-Thomas vector fields, Math. Comp. 87 (2018), no. 310, 659–692. MR 3739213, DOI 10.1090/mcom/3254
L. F. Pavarino and S. Scacchi, Isogeometric block FETI-DP preconditioners for the Stokes and mixed linear elasticity systems, Comput. Methods Appl. Mech. Engrg. 310 (2016), 694–710. MR 3548578, DOI 10.1016/j.cma.2016.07.012
L. F. Pavarino, S. Scacchi, O. B. Widlund, and S. Zampini, Isogeometric BDDC deluxe preconditioners for linear elasticity, Math. Models Methods Appl. Sci. 28 (2018), no. 7, 1337–1370. MR 3815659, DOI 10.1142/S0218202518500367
Clemens Pechstein and Clark R. Dohrmann, A unified framework for adaptive BDDC, Electron. Trans. Numer. Anal. 46 (2017), 273–336. MR 3678572
Giancarlo Sangalli and Mattia Tani, Isogeometric preconditioners based on fast solvers for the Sylvester equation, SIAM J. Sci. Comput. 38 (2016), no. 6, A3644–A3671. MR 3572372, DOI 10.1137/16M1062788
A. F. Sarmiento, A. M. Cortes, D. A. Garcia, L. Dalcin, N. Collier, V. M. Calo. PetIGA-MF: a multi-field high-performance toolbox for structure-preserving B-splines spaces. J. Comput. Sci. 18 (2017), 117–131.
Jakub Šístek, Bedřich Sousedík, Pavel Burda, Jan Mandel, and Jaroslav Novotný, Application of the parallel BDDC preconditioner to the Stokes flow, Comput. & Fluids 46 (2011), 429–435. MR 2949036, DOI 10.1016/j.compfluid.2011.01.002
Stefan Takacs, A robust multigrid method for the time-dependent Stokes problem, SIAM J. Numer. Anal. 53 (2015), no. 6, 2634–2654. MR 3432850, DOI 10.1137/140969658
Andrea Toselli and Olof Widlund, Domain decomposition methods—algorithms and theory, Springer Series in Computational Mathematics, vol. 34, Springer-Verlag, Berlin, 2005. MR 2104179, DOI 10.1007/b137868
Xuemin Tu and Jing Li, A FETI-DP type domain decomposition algorithm for three-dimensional incompressible Stokes equations, SIAM J. Numer. Anal. 53 (2015), no. 2, 720–742. MR 3317383, DOI 10.1137/13094997X
Xuemin Tu and Jing Li, A unified dual-primal finite element tearing and interconnecting approach for incompressible Stokes equations, Internat. J. Numer. Methods Engrg. 94 (2013), no. 2, 128–149. MR 3040516, DOI 10.1002/nme.4439
Olof B. Widlund, BDDC domain decomposition algorithms, 75 years of mathematics of computation, Contemp. Math., vol. 754, Amer. Math. Soc., [Providence], RI, [2020] ©2020, pp. 261–281. MR 4132126, DOI 10.1090/conm/754/15146
Olof B. Widlund and Clark R. Dohrmann, BDDC deluxe domain decomposition, Domain decomposition methods in science and engineering XXII, Lect. Notes Comput. Sci. Eng., vol. 104, Springer, Cham, 2016, pp. 93–103. MR 3643544
Ronald Haynes, Scott MacLachlan, Xiao-Chuan Cai, Laurence Halpern, Hyea Hyun Kim, Axel Klawonn, and Olof Widlund (eds.), Domain decomposition methods in science and engineering XXV, Lecture Notes in Computational Science and Engineering, vol. 138, Springer, Cham, [2020] ©2020. Selected papers from the 25th International Conference on Domain Decomposition Methods held in St. John’s, NL, July 23–27, 2018. MR 4167534, DOI 10.1007/978-3-030-56750-7
Stefano Zampini, PCBDDC: a class of robust dual-primal methods in PETSc, SIAM J. Sci. Comput. 38 (2016), no. 5, S282–S306. MR 3565564, DOI 10.1137/15M1025785
Stefano Zampini, Adaptive BDDC deluxe methods for $\rm H(curl)$, Domain decomposition methods in science and engineering XXIII, Lect. Notes Comput. Sci. Eng., vol. 116, Springer, Cham, 2017, pp. 285–292. MR 3718365, DOI 10.1007/978-3-319-52389-7_{2}9
Stefano Zampini and Xuemin Tu, Multilevel balancing domain decomposition by constraints deluxe algorithms with adaptive coarse spaces for flow in porous media, SIAM J. Sci. Comput. 39 (2017), no. 4, A1389–A1415. MR 3682184, DOI 10.1137/16M1080653
Stefano Zampini, Panayot Vassilevski, Veselin Dobrev, and Tzanio Kolev, Balancing domain decomposition by constraints algorithms for curl-conforming spaces of arbitrary order, Domain decomposition methods in science and engineering XXIV, Lect. Notes Comput. Sci. Eng., vol. 125, Springer, Cham, 2018, pp. 103–116. MR 3989859, DOI 10.1007/978-3-319-93873-8_{8}