New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards' equation. (August 2016) - Lipnikov, Konstantin
Highlights: Preconditioning strategy for Jacobian-free solvers for Richards' equation is proposed. The strategy is based on a stable discretization of the continuum Jacobian. It simplifies the use of sophisticated advanced discretizations on polytopal meshes. The performance was improved in 3–20 times for challenging infiltration problems. 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. We also show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers.