Skip to main content
Log in

Versatile Two-Level Schwarz Preconditioners for Multiphase Flow

  • Published:
Computational Geosciences Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  1. 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).

    Google Scholar 

  2. O. Axelsson, Iterative Solution Methods (Cambridge Univ. Press, Cambridge, 1994).

    Google Scholar 

  3. K. Aziz and A. Settari, Petroleum Reservoir Simulation (Applied Science, London, 1979).

    Google Scholar 

  4. 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).

  5. 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.

  6. 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).

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

    Google Scholar 

  8. J. Bear, Dynamics of Fluids in Porous Media (Dover, New York, 1972).

    Google Scholar 

  9. 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).

    Google Scholar 

  10. M. Brezina and P. Vanek, A black-box iterative solver based on a two-level Schwarz method, Computing 63(3) (1999) 233-263.

    Google Scholar 

  11. W.L. Briggs, V.E. Henson and S.F. McCormick, AMultigrid Tutorial (SIAM, Philadelphia, PA, 2000).

    Google Scholar 

  12. 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.

  13. 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.

    Google Scholar 

  14. 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.

    Google Scholar 

  15. 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).

  16. 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.

    Google Scholar 

  17. 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.

    Google Scholar 

  18. 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.

    Google Scholar 

  19. 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.

    Google Scholar 

  20. 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.

    Google Scholar 

  21. 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.

    Google Scholar 

  22. C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations (SIAM, Philadelphia, PA, 1995).

  23. B.H. Kueper, W. Abbott and G. Farquhar, Experimental observations of multiphase flow in heterogeneous porous media, J. Contaminant Hydrology 5 (1989) 83-95.

    Google Scholar 

  24. 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.

    Google Scholar 

  25. C. Lawson, R. Hanson, D. Kincaid and F. Krogh, Basic linear algebra subprograms for FORTRAN usage, ACM Trans. Math. Software 5 (1979) 308-329.

    Google Scholar 

  26. 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.

  27. Y. Saad, GMRES:A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Scientific Comput. 7 (1986) 856-869.

    Google Scholar 

  28. F. Saied and G. Mahinthakumar, Efficient parallel multigrid based solvers for large scale groundwater flow simulations, Comput. Math. Appl. 35(7) (1998) 45-54.

    Google Scholar 

  29. B. Smith, P. Bjorstad and W. Gropp, Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations (Cambridge Univ. Press, New York, 1996).

    Google Scholar 

  30. 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.

    Google Scholar 

  31. 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.

    Google Scholar 

  32. 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.

    Google Scholar 

  33. 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.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to C.T. Miller.

Rights and permissions

Reprints 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

Download citation

  • Issue Date:

  • DOI: https://doi.org/10.1023/A:1023514922877

Navigation