U.S. flag

An official website of the United States government

Dot gov

Official websites use .gov
A .gov website belongs to an official government organization in the United States.


Secure .gov websites use HTTPS
A lock ( ) or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.


Main content area

Stochastic analysis of transient flow in heterogeneous, variably saturated porous media: the van Genuchten-Mualem constitutive model

Lu, Z., Zhang, D.
Vadose zone journal 2002 v.1 no.1 pp. -
unsaturated flow, saturated flow, porous media, vadose zone, soil water movement, mathematical models, stochastic processes, equations, Monte Carlo method, temporal variation, spatial variation
In this study, on the basis of the van Genuchten-Mualem constitutive relationship, we develop a general nonstationary stochastic model for transient, variably saturated flow in randomly heterogeneous media with the method of moment equations. We first derive partial differential equations governing the statistical moments of the flow quantities by perturbation expansions and then implement these equations under general conditions with the method of finite differences. The nonstationary stochastic model developed is applicable to the entire domain of bounded, multidimensional vadose zones or integrated unsaturated-saturated systems in the presence of random or deterministic recharge and sink¡source and in the presence of multiscale, nonstationary medium features. We demonstrate the model with some two-dimensional examples of unsaturated and integrated unsaturated-saturated flows. The validity of the developed stochastic model is confirmed through high-resolution Monte Carlo simulations. We also investigate the relative contributions of the soil variabilities (K(S), alpha, and n) as well as the variability in recharge Q to the pressure head variance. It is found that the pressure head variance is sensitive to these variabilities, in the order of n, alpha, K(S) and Q. Though the variability of and n is usually smaller than that of K(S) and Q, their effect on the pressure head variance should not be ignored.