Skip to main content
SpringerLink
Log in

Inexact Newton methods and the method of lines for solving Richards' equation in two space dimensions

  • Published:
Computational Geosciences Aims and scope Submit manuscript

Abstract

Richards' equation (RE) is often used to model flow in unsaturated porous media. This model captures physical effects, such as sharp fronts in fluid pressures and saturations, which are present in more complex models of multiphase flow. The numerical solution of RE is difficult not only because of these physical effects but also because of the mathematical problems that arise in dealing with the nonlinearities. The method of lines has been shown to be very effective for solving RE in one space dimension. When solving RE in two space dimensions, direct methods for solving the linearized problem for the Newton step are impractical. In this work, we show how the method of lines and Newton-iterative methods, which solve linear equations with iterative methods, can be applied to RE in two space dimensions. We present theoretical results on convergence and use that theory to design an adaptive method for computation of the linear tolerance. Numerical results show the method to be effective and robust compared with an existing approach.

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

Access this article

Log in via an institution

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

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  1. C.H. Bischof, Incremental condition estimation, SIAM J. Matrix Anal. Appl. 11 (1990) 312–322.

    Google Scholar 

  2. K.E. Brenan, S.L. Campbell and L.R. Petzold, Numerical Solution of Initial-Value Problems in Differential–Algebraic Equations, Classics in Applied Mathematics, Vol. 14 (SIAM, Philadelphia, PA, 1996).

    Google Scholar 

  3. P.N. Brown and A.C. Hindmarsh, Matrix-free methods for stiff systems of ODE's, SIAM J. Numer. Anal. 23(3) (1986) 610–638.

    Google Scholar 

  4. P.N. Brown and A.C. Hindmarsh, Reduced storage matrix methods in stiff ODE systems, Appl. Math. Comput. 31 (1989) 40–91.

    Google Scholar 

  5. P.N. Brown, A.C. Hindmarsh and L.R. Petzold, Using Krylov methods in the solution of large-scale differential–algebraic systems, SIAM J. Sci. Comput. 15(6) (1994) 1467–1488.

    Google Scholar 

  6. P.N. Brown and H.F. Walker, GMRES on (nearly) singular systems, SIAM J. Sci. Comput. 18(1) (1997) 37–51.

    Google Scholar 

  7. R. Dembo, S. Eisenstat and T. Steihaug, Inexact Newton methods, SIAM J. Numer. Anal. 19(2) (1982) 400–408.

    Google Scholar 

  8. S.C. Eisenstat and H.F. Walker, Choosing the forcing terms in an inexact Newton method, SIAM J. Sci. Comput. 17(1) (1996) 16–32.

    Google Scholar 

  9. T. Gudmundsson, C.S. Kenney and A.J. Laub, Small-sample statistical estimates for matrix norms, SIAM J. Matrix Anal. Appl. 16(3) (1995) 776–792.

    Google Scholar 

  10. C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations, Frontiers in Applied Mathematics, Vol. 16 (SIAM, Philadelphia, PA, 1995).

    Google Scholar 

  11. C.T. Kelley, C.T. Miller and M.D. Tocci, Termination of Newton/chord iterations and the method of lines, Technical Report CRSC-TR96-19, North Carolina State University (May 1996); SIAM J. Sci. Comput., to appear.

  12. C.S. Kenney and A.J. Laub, Small-sample statistical condition estimates for general matrix functions, SIAM J. Sci. Comput. 15(1) (1994) 36–61.

    Google Scholar 

  13. C.S. Kenney, A.J. Laub and M.S. Reese, Statistical condition estimation for linear systems, SIAM J. Sci. Comput. 19(2) (1998) 566–583.

    Google Scholar 

  14. S. Kumar, M. Pernice, R. Rawat and P.J. Smith, Nonlinear GMRES acceleration of CFD codes for 3-D turbulent flow, Technical Report, University of Utah (1998); submitted to SIAM J. Sci. Comput.

  15. C.T. Miller and C.T. Kelley, A comparison of strongly convergent solution schemes for sharp front infiltration problems, in: Computational Methods in Water Resources X, Vol. 1, eds. A. Peters, G. Wittum, B. Herrling, U. Meissner, C.A. Brebbia, W.G. Gray and G.F. Pinder (Kluwer Academic, Dordrecht, 1994) pp. 325–332.

    Google Scholar 

  16. Y. Mualem, A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resources Research 12(3) (1976) 513–522.

    Google Scholar 

  17. M. Pernice and H.F. Walker, NITSOL: A Newton iterative solver for nonlinear systems, SIAM J. Sci. Comput. 19(1) (1998) 302–318.

    Google Scholar 

  18. L.A. Richards, Capillary conduction of liquids through porous media, Phys. 1 (1931) 318–333.

    Google Scholar 

  19. Y. Saad and M.H. Schultz, GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 7 (1986) 856–869.

    Google Scholar 

  20. L.F. Shampine, Implementation of implicit formulas for the solution of ODE's, SIAM J. Sci. Statist. Comput. 1(1) (1980) 103–117.

    Google Scholar 

  21. G.L.G. Sleijpen and D.R. Fokkema, BICGSTAB(l) for linear equations involving unsymmetric matrices with complex spectrum, Electron. Trans. Numer. Anal. 1 (1993) 11–32.

    Google Scholar 

  22. B. Smith, P. Bjørstad and W. Gropp, Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations (Cambridge University Press, Cambridge, 1996).

    Google Scholar 

  23. P. Sonneveld, CGS: A fast Lanczos-type solver for nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 10 (1989) 36–52.

    Google Scholar 

  24. C.A. San Soucie, Mixed finite element methods for variably saturated subsurface flow, Ph.D. thesis, Rice University (1996).

  25. M.D. Tocci, Numerical methods for variably saturated flow and transport models, Ph.D. thesis, North Carolina State University (1997).

  26. 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, Advances in Water Resources 20(1) (1997) 1–14.

    Google Scholar 

  27. 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. Sci. Statist. Comput. 13(2) (1992) 631–644.

    Google Scholar 

  28. H.A. Van der Vorst, Parallel iterative solution methods for linear systems arising from discretized PDE's, in: Special Course on Parallel Computing in CFD, AGARD-R-807 (AGARD, Neuilly-sur-Seine, 1995).

    Google Scholar 

  29. M.T. van Genuchten, Predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. America J. 44(5) (1980) 892–898.

    Google Scholar 

  30. C.S. Woodward, A Newton–Krylov-multigrid solver for variably saturated flow problems, Technical Report UCRL-JC-129371, Lawrence Livermore National Laboratory (January 1998); submitted to 12th International Conference on Computational Methods in Water Resources.

Download references

Author information

Authors and Affiliations

  1. Department of Mathematical Sciences, Worcester Polytechnic Institute, Worcester, MA, 01609, USA

    Michael D. Tocci

  2. Center for Research in Scientific Computation, Department of Mathematics, North Carolina State University, Box 8205, Raleigh, NC, 27695-8205, USA

    C.T. Kelley, Cass T. Miller & Christopher E. Kees

Authors
  1. Michael D. Tocci
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. C.T. Kelley
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Cass T. Miller
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Christopher E. Kees
    View author publications

    You can also search for this author in PubMed Google Scholar

Rights and permissions

Reprints and permissions

About this article

Cite this article

Tocci, M.D., Kelley, C., Miller, C.T. et al. Inexact Newton methods and the method of lines for solving Richards' equation in two space dimensions. Computational Geosciences 2, 291–309 (1998). https://doi.org/10.1023/A:1011562522244

Download citation

  • Issue Date:

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

Search

Navigation