Numerical analysis of one-dimensional unsaturated flow in layered soils |
| |
| Institution: | 1. Institute of Agricultural Hydraulics, University of Naples ‘Federico II’, I-80055 Portici (Napoli, Italy;2. Interfaculty Institute of Hydraulics, University of Perugia, I-06125 Perugia, Italy;1. Universidad Politécnica de Madrid, E.T.S.I. de Agrónomos, Ciudad Universitaria s/n, 28040 Madrid, Spain;2. Dept. Ingeniería Agroforestal, Universidad Politécnica de Madrid, E.T.S.I. de Agrónomos, Ciudad Universitaria s/n, 28040 Madrid, Spain;1. Department of Civil and Environmental Engineering, Mississippi State University, Mississippi State, MS 39762, USA;2. Department of Aerospace Engineering, Mississippi State University, Mississippi State, MS 39762, USA;1. Université Sidi Mohamed Ben Abdellah, Faculté des Sciences (Dhar El Mehraz), Laboratoire d''Informatique, Signaux, Automatique et Cognitivisme (LISAC), P. O. Box 1796, Fès-Atlas 30003, Fès, Morocco;2. Universite de Pau et des Pays de l’Adour, E2S UPPA, CNRS, LMAP, Pau, France |
| |
| Abstract: | The paper deals with numerical solutions to the Richards equation to simulate one-dimensional flow processes in the unsaturated zone of layered soil profiles. The equation is expressed in the pressure-based form and a finite-difference algorithm is developed for accurately estimating the values of the hydraulic conductivity between two neighboring nodes positioned in different soil layers, often referred to as the interlayer hydraulic conductivity. The algorithm is based upon flux conservation and continuity of pressure potential at the interface between two consecutive layers, and does not add significantly to simulation run time. The validity of the model is established for a number of test problems by comparing numerical results with the analytical solutions developed by Srivastava and Yeh29 which hold for vertical infiltration towards the water table in a two-layer soil profile. The results show a significant reduction in relative mass balance errors when using the proposed model. Some specific insights into its numerical performance are also gained by comparisons with a numerical model in which the more common geometric averaging operator acts on the interlayer conductivities. |
| |
| Keywords: | |
| 本文献已被 ScienceDirect 等数据库收录！ |
|
