Upscaled flow and transport properties for heterogeneous unsaturated media

Raziuddin Khaleel, Tian-Chyi J Yeh, Zhiming Lu

Research output: Contribution to journalArticle

36 Citations (Scopus)

Abstract

To represent a heterogeneous unsaturated medium by its homogeneous equivalent, stochastic theory-based analytical formulas and numerical Monte Carlo simulations are used to obtain upscaled (effective) flow and transport properties. The Monte Carlo experiments simulate steady state flow and transient transport of nonreactive solute in two-dimensional heterogeneous media. Constitutive relations for unsaturated media at the mesh-size scale are based on the van Genuchten-Mualem relationships. A unit-mean-gradient approach is used to derive upscaled properties for flow parallel and perpendicular to bedding. Upscaled moisture retention and unsaturated conductivities (K) and longitudinal dispersivities are obtained by simulating steady gravity drainage conditions for a series of applied infiltration rates, characteristic of relatively dry conditions in coarse-textured sediments. Macrodispersivities are calculated on the basis of spatial moments of the ensemble-mean plume. Results show that the calculated effective unsaturated K based on the analytical formulas compare well with those based on Monte Carlo simulations. The macroscopic anisotropy (ratio of K parallel to bedding to K perpendicular to bedding) increases with increasing tension, although the increase is rather mild. The longitudinal macrodispersivity for the equivalent homogeneous medium also increases with increasing tension. A comparison of numerical results with stochastic solutions suggests that the computed dispersivities are of the same order of magnitude at low tensions. At higher tensions the dispersivity estimates deviate significantly. Nonetheless, both numerical and analytical results show that the longitudinal dispersivities for flow parallel to bedding are higher than those for flow perpendicular to bedding. In addition, our results suggest that the Fickian regime is reached much earlier for cases with flow perpendicular to bedding than for those with flow parallel to bedding.

Original languageEnglish (US)
Pages (from-to)111-1112
Number of pages1002
JournalWater Resources Research
Volume38
Issue number5
StatePublished - May 2002

Fingerprint

unsaturated medium
Parallel flow
Transport properties
dispersivity
Infiltration
Drainage
Gravitation
Sediments
Anisotropy
Moisture
homogeneous medium
heterogeneous medium
mesh size
gravity
simulation
Experiments
solute
drainage
anisotropy
infiltration

Keywords

  • Effective properties
  • Equivalent media
  • Numerical modeling
  • Spatial variability
  • Stochastic processes
  • Vadose zone

ASJC Scopus subject areas

  • Environmental Science(all)
  • Environmental Chemistry
  • Aquatic Science
  • Water Science and Technology

Cite this

Upscaled flow and transport properties for heterogeneous unsaturated media. / Khaleel, Raziuddin; Yeh, Tian-Chyi J; Lu, Zhiming.

In: Water Resources Research, Vol. 38, No. 5, 05.2002, p. 111-1112.

Research output: Contribution to journalArticle

Khaleel, Raziuddin ; Yeh, Tian-Chyi J ; Lu, Zhiming. / Upscaled flow and transport properties for heterogeneous unsaturated media. In: Water Resources Research. 2002 ; Vol. 38, No. 5. pp. 111-1112.
@article{934923e2f6544d798818cde1b73dc49a,
title = "Upscaled flow and transport properties for heterogeneous unsaturated media",
abstract = "To represent a heterogeneous unsaturated medium by its homogeneous equivalent, stochastic theory-based analytical formulas and numerical Monte Carlo simulations are used to obtain upscaled (effective) flow and transport properties. The Monte Carlo experiments simulate steady state flow and transient transport of nonreactive solute in two-dimensional heterogeneous media. Constitutive relations for unsaturated media at the mesh-size scale are based on the van Genuchten-Mualem relationships. A unit-mean-gradient approach is used to derive upscaled properties for flow parallel and perpendicular to bedding. Upscaled moisture retention and unsaturated conductivities (K) and longitudinal dispersivities are obtained by simulating steady gravity drainage conditions for a series of applied infiltration rates, characteristic of relatively dry conditions in coarse-textured sediments. Macrodispersivities are calculated on the basis of spatial moments of the ensemble-mean plume. Results show that the calculated effective unsaturated K based on the analytical formulas compare well with those based on Monte Carlo simulations. The macroscopic anisotropy (ratio of K parallel to bedding to K perpendicular to bedding) increases with increasing tension, although the increase is rather mild. The longitudinal macrodispersivity for the equivalent homogeneous medium also increases with increasing tension. A comparison of numerical results with stochastic solutions suggests that the computed dispersivities are of the same order of magnitude at low tensions. At higher tensions the dispersivity estimates deviate significantly. Nonetheless, both numerical and analytical results show that the longitudinal dispersivities for flow parallel to bedding are higher than those for flow perpendicular to bedding. In addition, our results suggest that the Fickian regime is reached much earlier for cases with flow perpendicular to bedding than for those with flow parallel to bedding.",
keywords = "Effective properties, Equivalent media, Numerical modeling, Spatial variability, Stochastic processes, Vadose zone",
author = "Raziuddin Khaleel and Yeh, {Tian-Chyi J} and Zhiming Lu",
year = "2002",
month = "5",
language = "English (US)",
volume = "38",
pages = "111--1112",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "5",

}

TY - JOUR

T1 - Upscaled flow and transport properties for heterogeneous unsaturated media

AU - Khaleel, Raziuddin

AU - Yeh, Tian-Chyi J

AU - Lu, Zhiming

PY - 2002/5

Y1 - 2002/5

N2 - To represent a heterogeneous unsaturated medium by its homogeneous equivalent, stochastic theory-based analytical formulas and numerical Monte Carlo simulations are used to obtain upscaled (effective) flow and transport properties. The Monte Carlo experiments simulate steady state flow and transient transport of nonreactive solute in two-dimensional heterogeneous media. Constitutive relations for unsaturated media at the mesh-size scale are based on the van Genuchten-Mualem relationships. A unit-mean-gradient approach is used to derive upscaled properties for flow parallel and perpendicular to bedding. Upscaled moisture retention and unsaturated conductivities (K) and longitudinal dispersivities are obtained by simulating steady gravity drainage conditions for a series of applied infiltration rates, characteristic of relatively dry conditions in coarse-textured sediments. Macrodispersivities are calculated on the basis of spatial moments of the ensemble-mean plume. Results show that the calculated effective unsaturated K based on the analytical formulas compare well with those based on Monte Carlo simulations. The macroscopic anisotropy (ratio of K parallel to bedding to K perpendicular to bedding) increases with increasing tension, although the increase is rather mild. The longitudinal macrodispersivity for the equivalent homogeneous medium also increases with increasing tension. A comparison of numerical results with stochastic solutions suggests that the computed dispersivities are of the same order of magnitude at low tensions. At higher tensions the dispersivity estimates deviate significantly. Nonetheless, both numerical and analytical results show that the longitudinal dispersivities for flow parallel to bedding are higher than those for flow perpendicular to bedding. In addition, our results suggest that the Fickian regime is reached much earlier for cases with flow perpendicular to bedding than for those with flow parallel to bedding.

AB - To represent a heterogeneous unsaturated medium by its homogeneous equivalent, stochastic theory-based analytical formulas and numerical Monte Carlo simulations are used to obtain upscaled (effective) flow and transport properties. The Monte Carlo experiments simulate steady state flow and transient transport of nonreactive solute in two-dimensional heterogeneous media. Constitutive relations for unsaturated media at the mesh-size scale are based on the van Genuchten-Mualem relationships. A unit-mean-gradient approach is used to derive upscaled properties for flow parallel and perpendicular to bedding. Upscaled moisture retention and unsaturated conductivities (K) and longitudinal dispersivities are obtained by simulating steady gravity drainage conditions for a series of applied infiltration rates, characteristic of relatively dry conditions in coarse-textured sediments. Macrodispersivities are calculated on the basis of spatial moments of the ensemble-mean plume. Results show that the calculated effective unsaturated K based on the analytical formulas compare well with those based on Monte Carlo simulations. The macroscopic anisotropy (ratio of K parallel to bedding to K perpendicular to bedding) increases with increasing tension, although the increase is rather mild. The longitudinal macrodispersivity for the equivalent homogeneous medium also increases with increasing tension. A comparison of numerical results with stochastic solutions suggests that the computed dispersivities are of the same order of magnitude at low tensions. At higher tensions the dispersivity estimates deviate significantly. Nonetheless, both numerical and analytical results show that the longitudinal dispersivities for flow parallel to bedding are higher than those for flow perpendicular to bedding. In addition, our results suggest that the Fickian regime is reached much earlier for cases with flow perpendicular to bedding than for those with flow parallel to bedding.

KW - Effective properties

KW - Equivalent media

KW - Numerical modeling

KW - Spatial variability

KW - Stochastic processes

KW - Vadose zone

UR - http://www.scopus.com/inward/record.url?scp=0036563020&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0036563020&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0036563020

VL - 38

SP - 111

EP - 1112

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 5

ER -