A terrain-following grid transform and preconditioner for parallel, large-scale, integrated hydrologic modeling

Research output: Contribution to journalArticlepeer-review

113 Scopus citations

Abstract

A terrain-following grid formulation (TFG) is presented for simulation of coupled variably-saturated subsurface and surface water flow. The TFG is introduced into the integrated hydrologic model, ParFlow, which uses an implicit, Newton Krylov solution technique. The analytical Jacobian is also formulated and presented and both the diagonal and non-symmetric terms are used to precondition the Krylov linear system. The new formulation is verified against an orthogonal stencil and is shown to provide increased accuracy at lower lateral spatial discretization for hillslope simulations. Using TFG, efficient scaling to a large number of processors (16,384) and a large domain size (8.1 Billion unknowns) is shown. This demonstrates the applicability of this formulation to high-resolution, large-spatial extent hydrology applications where topographic effects are important. Furthermore, cases where the analytical Jacobian is used for the Newton iteration and as a non-symmetric preconditioner for the linear system are shown to have faster computation times and better scaling. This demonstrates the importance of solver efficiency in parallel scaling through the use of an appropriate preconditioner.

Original languageEnglish (US)
Pages (from-to)109-117
Number of pages9
JournalAdvances in Water Resources
Volume53
DOIs
StatePublished - Mar 2013
Externally publishedYes

All Science Journal Classification (ASJC) codes

  • Water Science and Technology

Keywords

  • Integrated hydrologic modeling
  • Parallel computing
  • Solvers

Fingerprint Dive into the research topics of 'A terrain-following grid transform and preconditioner for parallel, large-scale, integrated hydrologic modeling'. Together they form a unique fingerprint.

Cite this