Abstract
Numerically modeling groundwater flow on finely discretized two- and three-dimensional domains requires solution algorithms appropriate for distributed memory multiprocessor architectures. Multilevel and domain decomposition algorithms are appropriate for preconditioning or solving linear systems in parallel and have, therefore, been applied to linear models for saturated groundwater flow. These algorithms have also been incorporated into more complex nonlinear multiphase flow models in the context of a linearization procedure such as Newton's method. In this work, we study a class of parallel preconditioners based on two-level Schwarz domain decomposition applied in a nonlinear two-phase flow numerical model. The restriction and interpolation operators are based on an aggregation approach that has a straightforward implementation for a variety of applications arising in subsurface modeling: structured and unstructured discretizations, finite elements and finite differences, and multicomponent model equations. We present model formulations, results from numerical experiments, and a comparison of a standard one-level Schwarz method to three two-level aggregation-based methods.
Similar content being viewed by others
References
E. Anderson, Z. Bai, C. Bischof, J. Demmel, J. Dongarra, J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, S. Ostrouchov and D. Sorensen, LAPACK Users' Guide (SIAM, Philadelphia, PA, 1992).
O. Axelsson, Iterative Solution Methods (Cambridge Univ. Press, Cambridge, 1994).
K. Aziz and A. Settari, Petroleum Reservoir Simulation (Applied Science, London, 1979).
S. Balay, K. Buschelman, W.D. Gropp, D. Kaushik, L.C. McInnes and B.F. Smith, PETSc home page, http://www.mcs.anl.gov/petsc (2001).
S. Balay, W.D. Gropp, L.C. McInnes and B.F. Smith, Efficient management of parallelism in object oriented numerical software libraries, in: Modern Software Tools in Scientific Computing, eds. E. Arge, A.M. Bruaset and H.P. Langtangen (1997) pp. 163-202.
S. Balay, W.D. Gropp, L.C. McInnes and B.F. Smith, PETSc users manual, Technical Report ANL-95/11, Revision 2.1.3, Argonne National Laboratory (2002).
P. Bastian and R. Helmig, Efficient fully-coupled solution techniques for two-phase flow in porous media–Parallel multigrid solution and large scale computations, Adv. inWater Resourc. 23(3) (1999) 199-216.
J. Bear, Dynamics of Fluids in Porous Media (Dover, New York, 1972).
K.E. Brenan, S.L. Campbell and L.R. Petzold, The Numerical Solution of Initial Value Problems in Differential-Algebraic Equations (SIAM, Philadelphia, PA, 1996).
M. Brezina and P. Vanek, A black-box iterative solver based on a two-level Schwarz method, Computing 63(3) (1999) 233-263.
W.L. Briggs, V.E. Henson and S.F. McCormick, AMultigrid Tutorial (SIAM, Philadelphia, PA, 2000).
R.H. Brooks and A.T. Corey, Properties of porous media affecting fluid flow, J. Irrigation Drainage Division, Proc. of the American Society of Civil Engineers IR2 (1966) 61-88.
P.N. Brown, A.C. Hindmarsh and L.R. Petzold, Using Krylov methods in the solution of large-scale differential–algebraic systems, SIAM J. Scientific Comput. 15 (1994) 1467-1488.
C.N. Dawson, H. Klie, M.F. Wheeler and C. Woodward, A parallel, implicit cell-centered method for two-phase flow with a preconditioned Newton–Krylov solver, Comput. Geosci. 1 (1998) 215-250.
M. Dryja and O.B. Widlund, An additive variant of the Schwarz alternating method for the case of many subregions, Technical Report, Department of Computer Science, Courant Institute (1987).
I.G. Graham and M.J. Hagger, Unstructured additive Schwarz-conjugate gradient method for elliptic problems with highly discontinuous coefficients, SIAM J. Scientific Comput. 20(6) (1999) 2041-2066.
R. Helmig and R. Huber, Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media, Adv. in Water Resourc. 21(8) (1998) 697-711.
E.W. Jenkins, C.E. Kees, C.T. Kelley and C.T. Miller, An aggregation-based domain decomposition preconditioner for groundwater flow, SIAM J. Scientific Comput. 23(2) (2001) 430-441.
J.E. Jones and C.S. Woodward, Newton-Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems, Adv. in Water Resourc. 24(7) (2001) 763-774.
C.E. Kees and C.T. Miller, C++ implementations of numerical methods for solving differential–algebraic equations: Design and optimization considerations, ACM Trans. Math. Software 25(4) (1999) 377-403.
C.E. Kees and C.T. Miller, Higher order time integration methods for two-phase flow, Adv. in Water Resourc. 25(2) (2002) 159-177.
C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations (SIAM, Philadelphia, PA, 1995).
B.H. Kueper, W. Abbott and G. Farquhar, Experimental observations of multiphase flow in heterogeneous porous media, J. Contaminant Hydrology 5 (1989) 83-95.
B.H. Kueper and E.O. Frind, Two-phase flow in heterogeneous porous media 1. Model development, Water Resourc. Res. 27(6) (1991) 1049-1057.
C. Lawson, R. Hanson, D. Kincaid and F. Krogh, Basic linear algebra subprograms for FORTRAN usage, ACM Trans. Math. Software 5 (1979) 308-329.
J. Mandel, Hybrid domain decomposition with unstructured subdomains, in: Domain Decomposition Methods in Science and Engineering: 6th Internat. Conf. on Domain Decomposition, eds. A. Quarteroni, J. Periaux, Y.A. Kuznetsov and O.B. Widlund, Providence, RI, 1992, pp. 103-112.
Y. Saad, GMRES:A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Scientific Comput. 7 (1986) 856-869.
F. Saied and G. Mahinthakumar, Efficient parallel multigrid based solvers for large scale groundwater flow simulations, Comput. Math. Appl. 35(7) (1998) 45-54.
B. Smith, P. Bjorstad and W. Gropp, Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations (Cambridge Univ. Press, New York, 1996).
M.D. Tocci, C.T. Kelley and C.T. Miller, Accurate and economical solution of the pressure-head form of Richards' equation by the method of lines, Adv. in Water Resourc. 20(1) (1997) 1-14.
M.D. Tocci, C.T. Kelley, C.T. Miller and C.E. Kees, Inexact Newton methods and the method of lines for solving Richards' equation in two space dimensions, Comput. Geosci. 2(4) (1999) 291-309.
A.F.B. Tompson, R.D. Falgout, S.G. Smith, W.J. Bosl and S.F. Ashby, Analysis of subsurface contaminant migration and remediation using high performance computing, Adv. in Water Resourc. 22(3) (1998) 203-221.
H.A. van der Vorst, BI-CGSTAB: A fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems, SIAM J. Scientific Statist. Comput. 13(2) (1992) 631-644.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Kees, C., Miller, C., Jenkins, E. et al. Versatile Two-Level Schwarz Preconditioners for Multiphase Flow. Computational Geosciences 7, 91–114 (2003). https://doi.org/10.1023/A:1023514922877
Issue Date:
DOI: https://doi.org/10.1023/A:1023514922877