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, q-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 a 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 carried out 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 also 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).

Numerical simulations of one-dimensional infiltration into layered soils with the Richards equation using different estimates of the interlayer conductivity / B., Brunone; M., Ferrante; Romano, Nunzio; Santini, Alessandro. - In: VADOSE ZONE JOURNAL. - ISSN 1539-1663. - ELETTRONICO. - 2:(2003), pp. 193-200.

Numerical simulations of one-dimensional infiltration into layered soils with the Richards equation using different estimates of the interlayer conductivity.

ROMANO, NUNZIO;SANTINI, ALESSANDRO
2003

Abstract

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, q-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 a 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 carried out 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 also 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).
2003
Numerical simulations of one-dimensional infiltration into layered soils with the Richards equation using different estimates of the interlayer conductivity / B., Brunone; M., Ferrante; Romano, Nunzio; Santini, Alessandro. - In: VADOSE ZONE JOURNAL. - ISSN 1539-1663. - ELETTRONICO. - 2:(2003), pp. 193-200.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11588/372255
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 46
  • ???jsp.display-item.citation.isi??? 40
social impact