Premium
On numerical errors associated with the iterative alternating direction implicit (Iadi) finite difference solution of the two dimensional transient saturated‐unsaturated flow (Richards) equation
Author(s) -
Parissopoulos G. A.,
Wheater H. S.
Publication year - 1988
Publication title -
hydrological processes
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.222
H-Index - 161
eISSN - 1099-1085
pISSN - 0885-6087
DOI - 10.1002/hyp.3360020208
Subject(s) - richards equation , vadose zone , mathematics , flow (mathematics) , convergence (economics) , computer simulation , numerical analysis , infiltration (hvac) , finite difference , hydraulic conductivity , transient (computer programming) , sensitivity (control systems) , mechanics , water content , mathematical analysis , computer science , geometry , geotechnical engineering , soil water , geology , thermodynamics , soil science , physics , statistics , operating system , electronic engineering , engineering , economic growth , economics
Numerical models for simulation of multidimensional unsaturated flow are becoming increasingly available, but relatively little has been reported on the detailed analysis of numerical errors associated with such schemes. For unsaturated‐saturated flow, further complexity is introduced to the highly non‐linear unsaturated problem as the form of governing equation changes within the flow field. In this paper, two‐dimensional simulation of infiltration to a water table is considered. Numerical errors associated with selection of time step, grid geometry, convergence criteria, and the representation of internodal hydraulic conductivity are discussed with respect to moisture content profiles and mass balance error. Although solution sensitivity to numerical parameters is problem specific, the results presented indicate the nature and magnitude of numerical effects which should not be overlooked in model applications.