Author(s)

Carol S. Woodward

Abstract

A scalable, parallel solver for Richards' equation is described. This solver uses a globalized Newton method for linearization, a Krylov iterative method as the linear Jacobian system solver and the symmetric part of the Jacobian with a multigrid solver as a preconditioner. Scalability results are shown for both refining and enlarging the domain. 1 Introduction Accurate simulation of water resource management problems requires the solution of large problems with many spatial zones. Although par- allel processing is necessary to solve these large problems, it is not always sufficient; scalable and efficient solvers are also needed. By "scalable," we generally mean the ability to effectively use additional computational resources to solve l

Keywords