×

Physics-based balancing domain decomposition by constraints for multi-material problems. (English) Zbl 1419.65138

Summary: In this work, we present a new variant of the balancing domain decomposition by constraints preconditioner that is robust for multi-material problems. We start with a well-balanced subdomain partition, and based on an aggregation of elements according to their physical coefficients, we end up with a finer physics-based (PB) subdomain partition. Next, we define corners, edges, and faces for this PB partition, and select some of them to enforce subdomain continuity (primal faces/edges/corners). When the physical coefficient in each PB subdomain is constant and the set of selected primal faces/edges/corners satisfy a mild condition on the existence of acceptable paths, we can show both theoretically and numerically that the condition number does not depend on the contrast of the coefficient across subdomains. An extensive set of numerical experiments for 2D and 3D for the Poisson and linear elasticity problems is provided to support our findings. In particular, we show robustness and weak scalability of the new preconditioner variant up to 8232 cores when applied to 3D multi-material problems with the contrast of the physical coefficient up to \(10^8\) and more than half a billion degrees of freedom. For the scalability analysis, we have exploited a highly scalable advanced inter-level overlapped implementation of the preconditioner that deals very efficiently with the coarse problem computation. The proposed preconditioner is compared against a state-of-the-art implementation of an adaptive BDDC method in PETSc for thermal and mechanical multi-material problems.

MSC:

65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs
65F08 Preconditioners for iterative methods

Software:

FEMPAR; PCBDDC; PETSc
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Badia, S., Martín, A.F., Principe, J.: Enhanced balancing Neumann-Neumann preconditioning in computational fluid and solid mechanics. Int. J. Numer. Methods Eng. 96, 203-230 (2013) · Zbl 1352.74322
[2] Badia, S., Martín, A.F., Principe, J.: Implementation and scalability analysis of balancing domain decomposition methods. Arch. Comput. Methods Eng. 20, 239-262 (2013) · Zbl 1354.65261
[3] Badia, S., Martín, A.F., Principe, J.: A highly scalable parallel implementation of balancing domain decomposition by constraints. SIAM J. Sci. Comput. 36, C190-C218 (2014) · Zbl 1296.65177
[4] Badia, S., Martín, A.F., Principe, J.: On the scalability of inexact balancing domain decomposition by constraints with overlapped coarse/fine corrections. Parallel Comput. 50, 1-24 (2015)
[5] Badia, S., Martín, A.F., Principe, J.: Multilevel balancing domain decomposition at extreme scales. SIAM J. Sci. Comput. 38, C22-C52 (2016) · Zbl 1334.65217
[6] Badia, S., Martín, A. F., Principe, J.: FEMPAR Web page. http://www.fempar.org, 2017
[7] Badia, S., Martín, A.F., Nguyen, H.: Balancing domain decomposition by constraints associated with subobjects. Appl. Math. Lett. 87, 93-100 (2019) · Zbl 1410.65472
[8] Badia, S., Martín, A.F., Principe, J.: FEMPAR: an object-oriented parallel finite element framework. Arch. Comput. Methods Eng. 25, 195-271 (2018) · Zbl 1392.65005
[9] Badia, S., Nguyen, H.: Balancing domain decomposition by constraints and perturbation. SIAM J. Numer. Anal. 54, 3436-3464 (2016) · Zbl 1354.65262
[10] Badia, S., Nguyen, H.: Relaxing the roles of corners in BDDC by perturbed formulation. In: Domain Decomposition Methods in Science and Engineering XXIII. Springer, pp. 397-405 (2017) · Zbl 1367.65047
[11] Balay, S., Abhyankar, S., Adams, M. F., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M.G., May, D.A., McInnes, L.C., Rupp, K., Sanan, P., Smith, B.F., Zampini, S., Zhang, H., Zhang, H.: PETSc users manual, Tech. Report ANL-95/11 - Revision 3.8, Argonne National Laboratory (2017)
[12] Brenner, S.C., Scott, L.R.: The Mathematical Theory of Finite Element Methods. Texts in Applied Mathematics, vol. 15, 3rd edn. Springer, New York (2008) · Zbl 1135.65042
[13] Brenner, S.C., Sung, L.-Y.: BDDC and FETI-DP without matrices or vectors. Comput. Methods Appl. Mech. Eng. 196, 1429-1435 (2007) · Zbl 1173.65363
[14] Brömmel, D., Wylie, B.J.N., Frings, W.: JUQUEEN Extreme Scaling Workshop 2015, Tech. Report FZJ-2015-01645, Jülich Supercomputing Center (2015)
[15] Chan, T.F., Mathew, T.P.: Domain decomposition algorithms. Acta Numer. 3, 61-143 (1994) · Zbl 0809.65112
[16] Dohrmann, C.R.: A preconditioner for substructuring based on constrained energy minimization. SIAM J. Sci. Comput. 25, 246-258 (2003) · Zbl 1038.65039
[17] Dohrmann, C.R.: An approximate BDDC preconditioner. Numer. Linear Algebra Appl. 14, 149-168 (2007) · Zbl 1199.65088
[18] Dolean, V., Nataf, F., Scheichl, R., Spillane, N.: Analysis of a two-level Schwarz method with coarse spaces based on local Dirichlet-to-Neumann maps. Comput. Methods Appl. Math. 12, 391-414 (2012) · Zbl 1284.65050
[19] Dryja, M., Galvis, J., Sarkis, M.: BDDC methods for discontinuous Galerkin discretization of elliptic problems. J. Complex. 23, 715-739 (2007) · Zbl 1133.65097
[20] Dryja, M.; Sarkis, M.; Huang, Y. (ed.); Kornhuber, R. (ed.); Widlund, O. (ed.); Xu, J. (ed.), Technical tools for boundary layers and applications to heterogeneous coefficients, 205-212 (2011), Berlin · Zbl 1217.65220
[21] Dryja, M., Sarkis, M.V., Widlund, O.B.: Multilevel Schwarz methods for elliptic problems with discontinuous coefficients in three dimensions. Numer. Math. 72, 313-348 (1996) · Zbl 0857.65131
[22] Dryja, M., Smith, B.F., Widlund, O.B.: Schwarz analysis of iterative substructuring algorithms for elliptic problems in three dimensions. SIAM J. Numer. Anal. 31, 1662-1694 (1994) · Zbl 0818.65114
[23] Farhat, C., Lesoinne, M., Le Tallec, P., Pierson, K., Rixen, D.: FETI-DP: a dual – primal unified FETI method-part I: a faster alternative to the two-level FETI method. Int. J. Numer. Methods Eng. 50, 1523-1544 (2001) · Zbl 1008.74076
[24] Farhat, C., Lesoinne, M., Pierson, K.: A scalable dual – primal domain decomposition method. Numer. Linear Algebra Appl. 7, 687-714 (2000) · Zbl 1051.65119
[25] Galvis, J., Efendiev, Y.: Domain decomposition preconditioners for multiscale flows in high-contrast media. Multiscale Model. Simul. 8, 1461-1483 (2010) · Zbl 1206.76042
[26] Galvis, J., Efendiev, Y.: Domain decomposition preconditioners for multiscale flows in high contrast media: reduced dimension coarse spaces. Multiscale Model. Simul. 8, 1621-1644 (2010) · Zbl 1381.65029
[27] Gander, M.J., Loneland, A., Rahman, T.: Analysis of a new harmonically enriched multiscale coarse space for domain decomposition methods (2015). arXiv:1512.05285
[28] Gippert, S., Klawonn, A., Rheinbach, O.: Analysis of FETI-DP and BDDC for linear elasticity in 3D with almost incompressible components and varying coefficients inside subdomains. SIAM J. Numer. Anal. 50, 2208-2236 (2012) · Zbl 1255.74063
[29] Graham, I.G., Lechner, P.O., Scheichl, R.: Domain decomposition for multiscale PDEs. Numer. Math. 106, 589-626 (2007) · Zbl 1141.65084
[30] Karypis, G., Kumar, V.: A fast and high quality multilevel scheme for partitioning irregular graphs. SIAM J. Sci. Comput 20, 359-392 (1998). (electronic) · Zbl 0915.68129
[31] Kim, H.H., Chung, E.T.: A BDDC algorithm with enriched coarse spaces for two-dimensional elliptic problems with oscillatory and high contrast coefficients. Multiscale Model. Simul. 13, 571-593 (2015) · Zbl 1317.65090
[32] Klawonn, A., Radtke, P., Rheinbach, O.: FETI-DP methods with an adaptive coarse space. SIAM J. Numer. Anal. 53, 297-320 (2015) · Zbl 1327.65063
[33] Klawonn, A.; Radtke, P.; Rheinbach, O.; Dickopf, T. (ed.); Gander, JM (ed.); Halpern, L. (ed.); Krause, R. (ed.); Pavarino, FL (ed.), Adaptive coarse spaces for BDDC with a transformation of basis, 301-309 (2016), Berlin · Zbl 1339.65238
[34] Klawonn, A., Radtke, P., Rheinbach, O.: A comparison of adaptive coarse spaces for iterative substructuring in two dimensions. Electron. Trans. Numer. Anal. 45, 75-106 (2016) · Zbl 1338.65084
[35] Klawonn, A., Rheinbach, O.: Robust FETI-DP methods for heterogeneous three dimensional elasticity problems. Comput. Methods Appl. Mech. Eng. 196, 1400-1414 (2007) · Zbl 1173.74428
[36] Klawonn, A., Widlund, O.B.: Dual – primal FETI methods for linear elasticity. Commun. Pure Appl. Math. 59, 1523-1572 (2006) · Zbl 1110.74053
[37] Klawonn, A., Widlund, O.B., Dryja, M.: Dual – primal FETI methods for three-dimensional elliptic problems with heterogeneous coefficients. SIAM J. Numer. Anal. 40, 159-179 (2002). (electronic) · Zbl 1032.65031
[38] Li, J., Widlund, O.B.: FETI-DP, BDDC, and block Cholesky methods. Int. J. Numer. Methods Eng. 66, 250-271 (2006) · Zbl 1114.65142
[39] Li, J., Widlund, O.B.: On the use of inexact subdomain solvers for BDDC algorithms. Comput. Methods Appl. Mech. Eng. 196, 1415-1428 (2007) · Zbl 1173.65365
[40] Loisel, S., Nguyen, H., Scheichl, R.: Optimized Schwarz and 2-Lagrange multiplier methods for multiscale elliptic PDEs. SIAM J. Sci. Comput. 37, A2896-A2923 (2015) · Zbl 1329.65296
[41] Mandel, J.: Balancing domain decomposition. Commun. Numer. Methods Eng. 9, 233-241 (1993) · Zbl 0796.65126
[42] Mandel, J., Dohrmann, C.R.: Convergence of a balancing domain decomposition by constraints and energy minimization. Numer. Linear Algebra Appl. 10, 639-659 (2003) · Zbl 1071.65558
[43] Mandel, J., Dohrmann, C.R., Tezaur, R.: An algebraic theory for primal and dual substructuring methods by constraints. Appl. Numer. Math. 54, 167-193 (2005) · Zbl 1076.65100
[44] Mandel, J., Sousedík, B., Dohrmann, C.: Multispace and multilevel BDDC. Computing 83, 55-85 (2008) · Zbl 1163.65091
[45] Mathew, T.P.A.: Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations, vol. 61. Springer, Berlin (2008) · Zbl 1147.65101
[46] Nataf, F., Xiang, H., Dolean, V., Spillane, N.: A coarse space construction based on local Dirichlet-to-Neumann maps. SIAM J. Sci. Comput. 33, 1623-1642 (2011) · Zbl 1230.65134
[47] Pechstein, C.; Erhel, J. (ed.); Gander, MJ (ed.); Halpern, L. (ed.); Pichot, G. (ed.); Sassi, T. (ed.); Widlund, O. (ed.), On iterative substructuring methods for multiscale problems, 85-98 (2014), Berlin
[48] Pechstein, C., Dohrmann, C.R.: A unified framework for adaptive BDDC. Electron. Trans. Numer. Anal. 46, 273-336 (2017) · Zbl 1368.65043
[49] Pechstein, C., Scheichl, R.: Analysis of FETI methods for multiscale PDEs. Numer. Math. 111, 293-333 (2008) · Zbl 1170.65097
[50] Pechstein, C., Scheichl, R.: Analysis of FETI methods for multiscale PDEs. Part II: Interface variation. Numer. Math. 118, 485-529 (2011) · Zbl 1380.65388
[51] Scheichl, R., Vainikko, E.: Additive Schwarz with aggregation-based coarsening for elliptic problems with highly variable coefficients. Computing 80, 319-343 (2007) · Zbl 1171.65372
[52] Šístek, J., Čertíková, M., Burda, P., Novotný, J.: Face-based selection of corners in 3D substructuring. Math. Comput. Simul. 82, 1799-1811 (2012)
[53] Spillane, N., Dolean, V., Hauret, P., Nataf, F., Pechstein, C., Scheichl, R.: A robust two-level domain decomposition preconditioner for systems of PDEs. C. R. Math. Acad. Sci. Paris 349, 1255-1259 (2011) · Zbl 1252.65201
[54] Spillane, N., Dolean, V., Hauret, P., Nataf, F., Pechstein, C., Scheichl, R.: Abstract robust coarse spaces for systems of PDEs via generalized eigenproblems in the overlaps. Numer. Math. 126, 741-770 (2014) · Zbl 1291.65109
[55] Spillane, N., Dolean, V., Hauret, P., Nataf, F., Rixen, D.J.: Solving generalized eigenvalue problems on the interfaces to build a robust two-level FETI method. C. R. Math. Acad. Sci. Paris 351, 197-201 (2013) · Zbl 1269.65037
[56] Toselli, A., Widlund, O.: Domain Decomposition Methods—Algorithms and Theory. Springer Series in Computational Mathematics, vol. 34. Springer, Berlin (2005) · Zbl 1069.65138
[57] Tu, X.: Three-level BDDC in three dimensions. SIAM J. Sci. Comput. 29, 1759-1780 (2007) · Zbl 1163.65094
[58] Widlund, OB; Dohrmann, CR; Dickopf, T. (ed.); Gander, JM (ed.); Halpern, L. (ed.); Krause, R. (ed.); Pavarino, FL (ed.), BDDC deluxe domain decomposition, 93-103 (2016), Cham · Zbl 1339.65242
[59] Zampini, S.: PCBDDC: a class of robust dual – primal methods in PETSc. SIAM J. Sci. Comput. 38, S282-S306 (2016) · Zbl 1352.65632
[60] Zhu, Y.: Domain decomposition preconditioners for elliptic equations with jump coefficients. Numer. Linear Algebra Appl. 15, 271-289 (2008) · Zbl 1212.65502
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.