Stochastic analysis of unsaturated flow: One-dimensional Monte Carlo simulations and comparisons with spectral perturbation analysis and field observations
- Univ. of California, Davis (USA)
A numerical experiment was designed to study the stochastic behavior of one-dimensional transient unsaturated flow in a Monte Carlo setting. Soil hydraulic properties, log-saturated hydraulic conductivity 1n K{sub s}, pore size distribution parameter {alpha}, and the specific water capacity C are assumed to be statistically homogeneous random fields described by exponential correlation function with identical correlation lengths. Fifty realizations of each soil hydraulic property, with statistical properties obtained from a field experiment, are generated by using a nearest-neighbor model. Numerical solutions of the one-dimensional flow equation are used repeatedly to obtain realizations of soil water pressure head {psi} and flux q corresponding to the realizations of ln K{sub s}, {alpha}, and C. The dependence of {psi} and q on the statistical properties of ln K{sub s}, {alpha}, and C and the magnitude of the hydraulic head gradient G prevailing at the lower boundary are investigated. In addition, results of Monte Carlo analysis and spectral perturbation analyses are compared with field observations. The greatest variability of {psi} and 1 are observed when soil hydraulic properties are uncorrelated and have large variances and integral scales and when G at the lower boundary is unity. The Monte Carlo and spectral perturbation analyses tend to agree reasonably well for the flow domains in which ergodicity of local soil hydraulic properties is assured. Results of the Monte Carlo and spectral perturbation analyses are also supported by field observations.
- OSTI ID:
- 5273914
- Journal Information:
- Water Resources Research; (United States), Vol. 26:9; ISSN 0043-1397
- Country of Publication:
- United States
- Language:
- English
Similar Records
Comparison of Field-Scale Effective Properties of Two-Phase Flow in Heterogeneous Porous Media Obtained by Stochastic Analysis and Numerical Experiments
Thrust at N{sup 3}LL with power corrections and a precision global fit for {alpha}{sub s}(m{sub Z})