HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Brunone, B. Articles by Santini, A. Search for Related Content PubMed Articles by Brunone, B. Articles by Santini, A. GeoRef GeoRef Citation Agricola Articles by Brunone, B. Articles by Santini, A. Related Collections Hydraulic Conductivity Soil Hydrology Soil Models

Numerical Simulations of One-Dimensional Infiltration into Layered Soils with the Richards Equation Using Different Estimates of the Interlayer Conductivity

^a Department of Civil and Environmental Engineering, University of Perugia, Perugia, Italy
^b Department of Agricultural Engineering, University of Naples Federico II, Portici (Naples), Italy
^* Corresponding author (nunzio.romano{at}unina.it)

Received 30 August 2002.

Transient water flow processes in unsaturated soils are usually modeled using the Richards equation. This paper compares several numerical approximations to this equation for vertical infiltration in layered soil profiles. Three formulations of the governing flow equation (i.e., the h-based, -based, and mixed forms) are compared for the critical test problem of infiltration into a layered soil profile with initially relatively low soil water contents. An efficient, yet relatively simple weighting algorithm is employed that improves the estimation of the interlayer hydraulic conductivity in an h-based finite-difference formulation. Results highlight improvements in the mass conservative properties of this model, which is termed the H-IL model. A comparison is then performed between the finite-difference H-IL model and a finite-element model for infiltration toward a water table. The H-IL model was found to be computationally very efficient for this test problem. For both illustrative flow examples, the different numerical models were evaluated in terms of their ability to reproduce an exact solution developed by Srivastava and Yeh (1991) in their work on one-dimensional, transient infiltration toward the water table in homogeneous and layered soils.

Abbreviations: FD, finite difference • FE, finite element • H-GM model, same as FD h-based model H-IL, but with geometric means (Eq. [4]) for computing the interlayer conductivities • H-IL, FD h-based model with interlayer hydraulic conductivity computed by algorithm of Romano et al. (1998) (see Eq. [1] and [5]) • MIX-GM, FD mixed model with geometric mean conductivities for all interior nodes (see Eq. [3] and [4]) • THT-GM, FD -based model with geometric mean conductivities for all interior nodes (see Eq. [2] and [4])