Premium
Numerical Solution of Richards' Equation: A Review of Advances and Challenges
Author(s) -
Farthing Matthew W.,
Ogden Fred L.
Publication year - 2017
Publication title -
soil science society of america journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.836
H-Index - 168
eISSN - 1435-0661
pISSN - 0361-5995
DOI - 10.2136/sssaj2017.02.0058
Subject(s) - richards equation , discretization , mathematics , robustness (evolution) , numerical analysis , boundary value problem , calculus (dental) , computer science , soil water , mathematical analysis , geology , medicine , biochemistry , chemistry , dentistry , gene , soil science
Core Ideas The numerical solution of Richards' equation remains challenging. Space/time discretization affects both computational effort and accuracy. Adaption of space and time discretizations produces benefits. Dissemination of codes and improved documentation are needed. Recent reformulation of one‐dimensional Richards' equation shows promise. The flow of water in partially saturated porous media is of importance in fields such as hydrology, agriculture, environment and waste management. It is also one of the most complex flows in nature. The Richards' equation describes the flow of water in an unsaturated porous medium due to the actions of gravity and capillarity neglecting the flow of the non‐wetting phase, usually air. Analytical solutions of Richards' equation exist only for simplified cases, so most practical situations require a numerical solution in one‐ two‐ or three‐dimensions, depending on the problem and complexity of the flow situation. Despite the fact that the first reasonably complete conservative numerical solution method was published in the early 1990s, the numerical solution of the Richards' equation remains computationally expensive and in certain circumstances, unreliable. A universally robust and accurate solution methodology has not yet been identified that is applicable across the range of soils, initial and boundary conditions found in practice. Existing solution codes have been modified over years to attempt to increase robustness. Despite theoretical results on the existence of solutions given sufficiently regular data and constitutive relations, our numerical methods often fail to demonstrate reliable convergence behavior in practice, especially for higher‐order methods. Because of robustness, the lack of higher‐order accuracy and computational expense, alternative solution approaches or methods are needed. There is also a need for better documentation of improved solution methodologies and benchmark test problems to facilitate consistent advances and avoid re‐inventing of the wheel.