Applications of Mathematics, Vol. 64, No. 3, pp. 309-334, 2019


On adaptive BDDC for the flow in heterogeneous porous media

Bedřich Sousedík

Received August 15, 2018.   Published online April 23, 2019.

Abstract:  We study a method based on Balancing Domain Decomposition by Constraints (BDDC) for numerical solution of a single-phase flow in heterogeneous porous media. The method solves for both flux and pressure variables. The fluxes are resolved in three steps: the coarse solve is followed by subdomain solves and last we look for a divergence-free flux correction and pressures using conjugate gradients with the BDDC preconditioner. Our main contribution is an application of the adaptive algorithm for selection of flux constraints. Performance of the method is illustrated on the benchmark problem from the 10th SPE Comparative Solution Project (SPE 10). Numerical experiments in both 2D and 3D demonstrate that the first two steps of the method exhibit some numerical upscaling properties, and the adaptive preconditioner in the last step allows a significant decrease in the number of iterations of conjugate gradients at a small additional cost.
Keywords:  iterative substructuring; balancing domain decomposition; BDDC; multiscale methods; adaptive methods, flow in porous media; reservoir simulation; SPE 10 benchmark
Classification MSC:  65F08, 65F10, 65M55, 65N55


References:
[1] J. E. Aarnes, T. Gimse, K.-A. Lie: An introduction to the numerics of flow in porous media using Matlab. Geometric Modelling, Numerical Simulation, and Optimization: Applied Mathematics at SINTEF. Springer, Berlin 2007, pp. 265-306. DOI 10.1007/978-3-540-68783-2_9 | MR 2348925 | Zbl 1330.76004
[2] J. E. Aarnes, S. Krogstad, K.-A. Lie: A hierarchical multiscale method for two-phase flow based upon mixed finite elements and nonuniform coarse grids. Multiscale Model. Simul. 5 (2006), 337-363. DOI 10.1137/050634566 | MR 2247754 | Zbl 1124.76022
[3] J. E. Aarnes, S. Krogstad, K.-A. Lie: Multiscale mixed/mimetic methods on corner-point grids. Comput. Geosci. 12 (2008), 297-315. DOI 10.1007/s10596-007-9072-8 | MR 2434946 | Zbl 1259.76065
[4] F. Brezzi, M. Fortin: Mixed and Hybrid Finite Element Methods. Springer Series in Computational Mathematics 15, Springer, New York (1991). DOI 10.1007/978-1-4612-3172-1 | MR 1115205 | Zbl 0788.73002
[5] M. A. Christie, M. J. Blunt: Tenth SPE comparative solution project: A comparison of upscaling techniques. SPE Reservoir Eval. Eng. 4 (2001), 308-317. DOI 10.2118/72469-pa
[6] L. C. Cowsar, J. Mandel, M. F. Wheeler: Balancing domain decomposition for mixed finite elements. Math. Comput. 64 (1995), 989-1015. DOI 10.2307/2153480 | MR 1297465 | Zbl 0828.65135
[7] J.-M. Cros: A preconditioner for the Schur complement domain decomposition method. 14th Int. Conf. on Domain Decomposition Methods in Science and Engineering (I. Herrera et al., eds.). National Autonomous University of Mexico (UNAM), México, 2003, pp. 373-380. MR 2093729 | Zbl 1103.65004
[8] J. W. Demmel: Applied Numerical Linear Algebra. Society for Industrial and Applied Mathematics, Philadelphia (1997). DOI 10.1137/1.9781611971446 | MR 1463942 | Zbl 0879.65017
[9] C. R. Dohrmann: A preconditioner for substructuring based on constrained energy minimization. SIAM J. Sci. Comput. 25 (2003), 246-258. DOI 10.1137/S1064827502412887 | MR 2047204 | Zbl 1038.65039
[10] C. R. Dohrmann: A substructuring preconditioner for nearly incompressible elasticity problems. Technical report SAND 2004-5393, Sandia National Laboratories (2004).
[11] C. R. Dohrmann, O. B. Widlund: Some recent tools and a BDDC algorithm for 3D problems in $H(curl)$. Domain Decomposition Methods in Science and Engineering XX. Lecture Notes Computational Science and Engineering 91, Springer, Heidelberg, 2013, pp. 15-25. DOI 10.1007/978-3-642-35275-1_2 | MR 3242973 | Zbl 06125818
[12] Y. Efendiev, T. Y. Hou: Multiscale Finite Element Methods: Theory and Applications. Surveys and Tutorials in the Applied Mathematical Sciences 4, Springer, New York (2009). DOI 10.1007/978-0-387-09496-0 | MR 2477579 | Zbl 1163.65080
[13] R. E. Ewing, J. Wang: Analysis of the Schwarz algorithm for mixed finite elements methods. RAIRO, Modélisation Math. Anal. Numér. 26 (1992), 739-756. DOI 10.1051/m2an/1992260607391 | MR 1183415 | Zbl 0765.65104
[14] C. Farhat, M. Lesoinne, P. LeTallec, K. Pierson, D. Rixen: FETI-DP: a dual-primal unified FETI method. I. A faster alternative to the two-level FETI method. Int. J. Numer. Methods Eng. 50 (2001), 1523-1544. DOI 10.1002/nme.76 | MR 1813746 | Zbl 1008.74076
[15] C. Farhat, M. Lesoinne, K. Pierson: A scalable dual-primal domain decomposition method. Numer. Linear Algebra Appl. 7 (2000), 687-714. DOI 10.1002/1099-1506(200010/12)7:7/8<687::AID-NLA219>3.0.CO;2-S | MR 1802366 | Zbl 1051.65119
[16] Y. Fragakis, M. Papadrakakis: The mosaic of high performance domain decomposition methods for structural mechanics: formulation, interrelation and numerical efficiency of primal and dual methods. Comput. Methods Appl. Mech. Eng. 192 (2003), 3799-3830. DOI 10.1016/S0045-7825(03)00374-8 | Zbl 1054.74069
[17] R. Glowinski, M. F. Wheeler: Domain decomposition and mixed finite element methods for elliptic problems. First International Symposium on Domain Decomposition Methods for Partial Differential Equations. SIAM, Philadelphia, 1988, pp. 144-172. MR 0972516 | Zbl 0661.65105
[18] G. H. Golub, C. F. Van Loan: Matrix Computations. Johns Hopkins Studies in the Mathematical Sciences, Johns Hopkins University Press, Baltimore (1996). MR 1417720 | Zbl 0865.65009
[19] M. Hanek, J. Šístek, P. Burda: The effect of irregular interfaces on the BDDC method for the Navier-Stokes equations. Proc. Int. Conf. Domain Decomposition Methods in Science and Engineering XXIII. Lecture Notes Computational Science and Engineering 116, Springer, Cham, 2017, pp. 171-178. DOI 10.1007/978-3-319-52389-7_16 | MR 3718352 | Zbl 06747817
[20] G. Karypis, V. Kumar: METIS: A software package for partitioning unstructured graphs, partitioning meshes, and computing fill-reducing orderings of sparse matrices, version 4.0. Technical report, Department of Computer Science, University of Minnesota (1998).
[21] A. Klawonn, M. Kühn, O. Rheinbach: Adaptive coarse spaces for FETI-DP in three dimensions. SIAM J. Sci. Comput. 38 (2016), A2880-A2911 (2016). DOI 10.1137/15m1049610 | MR 3546980 | Zbl 1346.74168
[22] A. Klawonn, M. Kühn, O. Rheinbach: A closer look at local eigenvalue solvers for adaptive FETI-DP and BDDC. Technical report, Universität zu Köln (2018). Available at https://kups.ub.uni-koeln.de/9020/.
[23] A. Klawonn, O. Rheinbach, O. B. Widlund: An analysis of a FETI-DP algorithm on irregular subdomains in the plane. SIAM J. Numer. Anal. 46 (2008), 2484-2504. DOI 10.1137/070688675 | MR 2421044 | Zbl 1176.65135
[24] A. V. Knyazev: Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method. SIAM J. Sci. Comput. 23 (2001), 517-541. DOI 10.1137/S1064827500366124 | MR 1861263 | Zbl 0992.65028
[25] M. la Cour Christensen, U. Villa, A. P. Engsig-Karup, P. S. Vassilevski: Numerical multilevel upscaling for incompressible flow in reservoir simulation: an element-based algebraic multigrid (AMGe) approach. SIAM J. Sci. Comput. 39 (2017), B102-B137. DOI 10.1137/140988991 | MR 3612903 | Zbl 1360.65247
[26] J. Li, X. Tu: Convergence analysis of a balancing domain decomposition method for solving a class of indefinite linear systems. Numer. Linear Algebra Appl. 16 (2009), 745-773. DOI 10.1002/nla.639 | MR 2554500 | Zbl 1224.65248
[27] J. Li, O. B. Widlund: BDDC algorithms for incompressible Stokes equations. SIAM J. Numer. Anal. 44 (2006), 2432-2455. DOI 10.1137/050628556 | MR 2272601 | Zbl 1233.76077
[28] J. Li, O. B. Widlund: FETI-DP, BDDC, and block Cholesky methods. Int. J. Numer. Methods Eng. 66 (2006), 250-271. DOI 10.1002/nme.1553 | MR 2224479 | Zbl 1114.65142
[29] J. Mandel, B. Sousedík: Adaptive selection of face coarse degrees of freedom in the BDDC and the FETI-DP iterative substructuring methods. Comput. Methods Appl. Mech. Eng. 196 (2007), 1389-1399. DOI 10.1016/j.cma.2006.03.010 | MR 2277024 | Zbl 1173.74435
[30] J. Mandel, B. Sousedík, C. R. Dohrmann: Multispace and multilevel BDDC. Computing 83 (2008), 55-85. DOI 10.1007/s00607-008-0014-7 | MR 2457352 | Zbl 1163.65091
[31] J. Mandel, B. Sousedík, J. Šístek: Adaptive BDDC in three dimensions. Math. Comput. Simul. 82 (2012), 1812-1831. DOI 10.1016/j.matcom.2011.03.014 | MR 2967935 | Zbl 1255.65225
[32] T. P. Mathew: Schwarz alternating and iterative refinement methods for mixed formulations of elliptic problems. I. Algorithms and numerical results. Numer. Math. 65 (1993), 445-468. DOI 10.1007/BF01385762 | MR 1231895 | Zbl 0801.65106
[33] D.-S. Oh, O. B. Widlund, S. Zampini, C. R. Dohrmann: BDDC algorithms with deluxe scaling and adaptive selection of primal constraints for Raviart-Thomas vector fields. Math. Comput. 87 (2018), 659-692. DOI 10.1090/mcom/3254 | MR 3739213 | Zbl 1380.65065
[34] C. Pechstein: Finite and Boundary Element Tearing and Interconnecting Solvers for Multiscale Problems. Lecture Notes in Computational Science and Engineering 90, Springer, Berlin (2013). DOI 10.1007/978-3-642-23588-7 | MR 3013465 | Zbl 1272.65100
[35] C. Pechstein, C. R. Dohrmann: A unified framework for adaptive BDDC. ETNA, Electron. Trans. Numer. Anal. 46 (2017), 273-336. MR 3678572 | Zbl 1368.65043
[36] C. Pechstein, R. Scheichl: Analysis of FETI methods for multiscale PDEs. Part II: interface variation. Numer. Math. 118 (2011), 485-529. DOI 10.1007/s00211-011-0359-2 | MR 2810804 | Zbl 1380.65388
[37] J. Šístek, J. Březina, B. Sousedík: BDDC for mixed-hybrid formulation of flow in porous media with combined mesh dimensions. Numer. Linear Algebra Appl. 22 (2015), 903-929. DOI 10.1002/nla.1991 | MR 3426321 | Zbl 1389.76057
[38] B. Sousedík: Nested BDDC for a saddle-point problem. Numer. Math. 125 (2013), 761-783. DOI 10.1007/s00211-013-0548-2 | MR 3127330 | Zbl 1282.65167
[39] B. Sousedík, J. Šístek, J. Mandel: Adaptive-multilevel BDDC and its parallel implementation. Computing 95 (2013), 1087-1119. DOI 10.1007/s00607-013-0293-5 | MR 3125603 | Zbl 1307.65175
[40] N. Spillane, D. J. Rixen: Automatic spectral coarse spaces for robust finite element tearing and interconnecting and balanced domain decomposition algorithms. Int. J. Numer. Methods Eng. 95 (2013), 953-990. DOI 10.1002/nme.4534 | MR 3093793 | Zbl 1352.65553
[41] A. Toselli, O. Widlund: Domain Decomposition Methods - Algorithms and Theory. Springer Series in Computational Mathematics 34, Springer, Berlin (2005). DOI 10.1007/b137868 | MR 2104179 | Zbl 1069.65138
[42] X. Tu: A BDDC algorithm for a mixed formulation of flow in porous media. ETNA, Electron. Trans. Numer. Anal. 20 (2005), 164-179. MR 2175341 | Zbl 1160.76368
[43] X. Tu: A BDDC algorithm for flow in porous media with a hybrid finite element discretization. ETNA, Electron. Trans. Numer. Anal. 26 (2007), 146-160. MR 2366094 | Zbl 1170.76034
[44] X. Tu: Three-level BDDC in three dimensions. SIAM J. Sci. Comput. 29 (2007), 1759-1780. DOI 10.1137/050629902 | MR 2341811 | Zbl 1163.65094
[45] X. Tu: Three-level BDDC in two dimensions. Int. J. Numer. Methods Eng. 69 (2007), 33-59. DOI 10.1002/nme.1753 | MR 2282536 | Zbl 1134.65087
[46] X. Tu: A three-level BDDC algorithm for a saddle point problem. Numer. Math. 119 (2011), 189-217. DOI 10.1007/s00211-011-0375-2 | MR 2824859 | Zbl 1230.65136
[47] X. Tu, J. Li: A balancing domain decomposition method by constraints for advection-diffusion problems. Commun. Appl. Math. Comput. Sci. 3 (2008), 25-60. DOI 10.2140/camcos.2008.3.25 | MR 2425545 | Zbl 1165.65402
[48] E. Vecharynski, Y. Saad, M. Sosonkina: Graph partitioning using matrix values for preconditioning symmetric positive definite systems. SIAM J. Sci. Comput. 36 (2014), A63-A87. DOI 10.1137/120898760 | MR 3151390 | Zbl 1290.65025
[49] Y. Yang, S. Fu, E. T. Chung: A two-grid preconditioner with an adaptive coarse space for flow simulations in highly heterogeneous media. Available at https://arxiv.org/abs/1807.07220 (2018), 17 pages.
[50] S. Zampini, X. Tu: Multilevel balancing domain decomposition by constraints deluxe algorithms with adaptive coarse spaces for flow in porous media. SIAM J. Sci. Comput. 39 (2017), A1389-A1415. DOI 10.1137/16M1080653 | MR 3682184 | Zbl 06760251

Affiliations:   Bedřich Sousedík, Department of Mathematics and Statistics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA, e-mail: sousedik@umbc.edu


 
PDF available at: