Overview#
The CATchment HYdrology (CATHY) model couples a finite element solver for the Richards equation describing ow in variably saturated porous media and a finite difference solver for the diffusion wave equation describing surface flow propagation throughout a hill slope and stream channel network identified using terrain topography and the hydraulic geometry concept The mathematical model is described by a system of two differential equations [1]:
Subsurface flow equation:
Surface flow equation:
The 3-D Richards equation is solved numerically by Galerkin finite elements in space using tetrahedral elements and linear basis functions, and by a weighted finite difference scheme for integration in time. The nonlinear characteristic relationships \(K_{r}(\Psi)\) and \(S_{w}(\Psi)\) are specified using either van Genuchten and Nielsen [1985], Brooks and Corey [1964], or Huyakorn et al. [1984] expression. Linearization via Newton or Picard iteration is used in the solution procedure.