New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation
Abstract
We develop a new approach for solving the nonlinear Richards’ equation arising in variably saturated flow modeling. The growing complexity of geometric models for simulation of subsurface flows leads to the necessity of using unstructured meshes and advanced discretization methods. Typically, a numerical solution is obtained by first discretizing PDEs and then solving the resulting system of nonlinear discrete equations with a Newton-Raphson-type method. Efficiency and robustness of the existing solvers rely on many factors, including an empiric quality control of intermediate iterates, complexity of the employed discretization method and a customized preconditioner. We propose and analyze a new preconditioning strategy that is based on a stable discretization of the continuum Jacobian. We will show with numerical experiments for challenging problems in subsurface hydrology that this new preconditioner improves convergence of the existing Jacobian-free solvers 3-20 times. Furthermore, we show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers.
- Authors:
-
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Publication Date:
- Research Org.:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Office of Environmental Management (EM)
- OSTI Identifier:
- 1256101
- Alternate Identifier(s):
- OSTI ID: 1325352
- Report Number(s):
- LA-UR-15-27929
Journal ID: ISSN 0309-1708
- Grant/Contract Number:
- AC52-06NA25396
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Advances in Water Resources
- Additional Journal Information:
- Journal Volume: 94; Journal Issue: C; Journal ID: ISSN 0309-1708
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 58 GEOSCIENCES; 97 MATHEMATICS AND COMPUTING; Mathematics; Planetary Sciences; Richards' equation, variable saturated flow, nonlinear solver, Picard method
Citation Formats
Lipnikov, Konstantin, Moulton, David, and Svyatskiy, Daniil. New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation. United States: N. p., 2016.
Web. doi:10.1016/j.advwatres.2016.04.016.
Lipnikov, Konstantin, Moulton, David, & Svyatskiy, Daniil. New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation. United States. https://doi.org/10.1016/j.advwatres.2016.04.016
Lipnikov, Konstantin, Moulton, David, and Svyatskiy, Daniil. Fri .
"New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation". United States. https://doi.org/10.1016/j.advwatres.2016.04.016. https://www.osti.gov/servlets/purl/1256101.
@article{osti_1256101,
title = {New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation},
author = {Lipnikov, Konstantin and Moulton, David and Svyatskiy, Daniil},
abstractNote = {We develop a new approach for solving the nonlinear Richards’ equation arising in variably saturated flow modeling. The growing complexity of geometric models for simulation of subsurface flows leads to the necessity of using unstructured meshes and advanced discretization methods. Typically, a numerical solution is obtained by first discretizing PDEs and then solving the resulting system of nonlinear discrete equations with a Newton-Raphson-type method. Efficiency and robustness of the existing solvers rely on many factors, including an empiric quality control of intermediate iterates, complexity of the employed discretization method and a customized preconditioner. We propose and analyze a new preconditioning strategy that is based on a stable discretization of the continuum Jacobian. We will show with numerical experiments for challenging problems in subsurface hydrology that this new preconditioner improves convergence of the existing Jacobian-free solvers 3-20 times. Furthermore, we show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers.},
doi = {10.1016/j.advwatres.2016.04.016},
journal = {Advances in Water Resources},
number = C,
volume = 94,
place = {United States},
year = {Fri Apr 29 00:00:00 EDT 2016},
month = {Fri Apr 29 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Iterative Procedures for Nonlinear Integral Equations
journal, October 1965
- Anderson, Donald G.
- Journal of the ACM, Vol. 12, Issue 4
On equivalent hydraulic conductivity for oscillation–free solutions of Richard’s equation
journal, November 2013
- Belfort, Benjamin; Younes, Anis; Fahs, Marwan
- Journal of Hydrology, Vol. 505
Mixed finite elements and Newton-type linearizations for the solution of Richards' equation
journal, July 1999
- Bergamaschi, Luca; Putti, Mario
- International Journal for Numerical Methods in Engineering, Vol. 45, Issue 8
A class of methods for solving nonlinear simultaneous equations
journal, January 1965
- Broyden, C. G.
- Mathematics of Computation, Vol. 19, Issue 92
Nonlinear Krylov acceleration applied to a discrete ordinates formulation of the k-eigenvalue problem
journal, April 2013
- Calef, Matthew T.; Fichtl, Erin D.; Warsa, James S.
- Journal of Computational Physics, Vol. 238
Design and Application of a Gradient-Weighted Moving Finite Element Code I: in One Dimension
journal, May 1998
- Carlson, Neil N.; Miller, Keith
- SIAM Journal on Scientific Computing, Vol. 19, Issue 3
Orthogonal collocation and alternating-direction procedures for unsaturated flow problems
journal, December 1987
- Celia, Michael A.; Ahuja, Lajpat R.; Pinder, George F.
- Advances in Water Resources, Vol. 10, Issue 4
A general mass-conservative numerical solution for the unsaturated flow equation
journal, July 1990
- Celia, Michael A.; Bouloutas, Efthimios T.; Zarba, Rebecca L.
- Water Resources Research, Vol. 26, Issue 7
A 2D high-β Hall MHD implicit nonlinear solver
journal, July 2003
- Chacón, L.; Knoll, D. A.
- Journal of Computational Physics, Vol. 188, Issue 2
Inexact Newton Methods
journal, April 1982
- Dembo, Ron S.; Eisenstat, Stanley C.; Steihaug, Trond
- SIAM Journal on Numerical Analysis, Vol. 19, Issue 2
A Unified Approach to Mimetic Finite Difference, Hybrid Finite Volume and Mixed Finite Volume Methods
journal, February 2010
- Droniou, JÉRÔMe; Eymard, Robert; GallouËT, Thierry
- Mathematical Models and Methods in Applied Sciences, Vol. 20, Issue 02
A Comparative Study on Methods for Convergence Acceleration of Iterative Vector Sequences
journal, March 1996
- Eyert, V.
- Journal of Computational Physics, Vol. 124, Issue 2
A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils1
journal, January 1980
- van Genuchten, M. Th.
- Soil Science Society of America Journal, Vol. 44, Issue 5
M-Adaptation in the mimetic finite difference method
journal, May 2014
- Gyrya, Vitaliy; Lipnikov, Konstantin; Manzini, Gianmarco
- Mathematical Models and Methods in Applied Sciences, Vol. 24, Issue 08
Methods of conjugate gradients for solving linear systems
journal, December 1952
- Hestenes, M. R.; Stiefel, E.
- Journal of Research of the National Bureau of Standards, Vol. 49, Issue 6
Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)
journal, January 2014
- Karra, S.; Painter, S. L.; Lichtner, P. C.
- The Cryosphere, Vol. 8, Issue 5
Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004
- Knoll, D. A.; Keyes, D. E.
- Journal of Computational Physics, Vol. 193, Issue 2
A Finite-Volume Method with Hexahedral Multiblock Grids for Modeling Flow in Porous Media
journal, January 2002
- Lee, S. H.; Jenny, P.; Tchelepi, H. A.
- Computational Geosciences, Vol. 6, Issue 3/4, p. 353-379
The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient
journal, January 2016
- Lipnikov, Konstantin; Manzini, Gianmarco; Moulton, J. David
- Journal of Computational Physics, Vol. 305
A monotone finite volume method for advection–diffusion equations on unstructured polygonal meshes
journal, June 2010
- Lipnikov, K.; Svyatskiy, D.; Vassilevski, Y.
- Journal of Computational Physics, Vol. 229, Issue 11
Minimal stencil finite volume scheme with the discrete maximum principle
journal, January 2012
- Lipnikov, K.; Svyatskiy, D.; Vassilevski, Yu.
- Russian Journal of Numerical Analysis and Mathematical Modelling, Vol. 27, Issue 4
An accelerated Picard method for nonlinear systems related to variably saturated flow
journal, March 2012
- Lott, P. A.; Walker, H. F.; Woodward, C. S.
- Advances in Water Resources, Vol. 38
Second-order accurate monotone finite volume scheme for Richards’ equation
journal, April 2013
- Misiats, Oleksandr; Lipnikov, Konstantin
- Journal of Computational Physics, Vol. 239
NITSOL: A Newton Iterative Solver for Nonlinear Systems
journal, January 1998
- Pernice, Michael; Walker, Homer F.
- SIAM Journal on Scientific Computing, Vol. 19, Issue 1
Capillary Conduction of Liquids Through Porous Mediums
journal, November 1931
- Richards, L. A.
- Physics, Vol. 1, Issue 5
Database-Related Accuracy and Uncertainty of Pedotransfer Functions
journal, January 1998
- Schaap, Marcel G.; Leij, Feike J.
- Soil Science, Vol. 163, Issue 10
Accurate and economical solution of the pressure-head form of Richards' equation by the method of lines
journal, February 1997
- Tocci, Michael D.; Kelley, C. T.; Miller, Cass T.
- Advances in Water Resources, Vol. 20, Issue 1
Anderson Acceleration for Fixed-Point Iterations
journal, January 2011
- Walker, Homer F.; Ni, Peng
- SIAM Journal on Numerical Analysis, Vol. 49, Issue 4
Monotonicity of the cell-centred triangular MPFA method for saturated and unsaturated flow in heterogeneous porous media
journal, November 2013
- Younes, Anis; Fahs, Marwan; Belfort, Benjamin
- Journal of Hydrology, Vol. 504
Picard-Newton Iterative Method with Time Step Control for Multimaterial Non-Equilibrium Radiation Diffusion Problem
journal, October 2011
- Yue, Jingyan; Yuan, Guangwei
- Communications in Computational Physics, Vol. 10, Issue 4
Works referencing / citing this record:
Numerical Solution of Richards' Equation: A Review of Advances and Challenges
journal, October 2017
- Farthing, Matthew W.; Ogden, Fred L.
- Soil Science Society of America Journal, Vol. 81, Issue 6