Unsaturated hydraulic conductivity of structured porous media: A review of liquid configuration-based models

Markus Tuller, Dani Or

Research output: Contribution to journalArticle

73 Citations (Scopus)

Abstract

Common approaches for modeling hydraulic functions of unsatu- rated structured porous media (SPM) rely on macroscopic continuum representation, where parameterization schemesand constitutive rela- tionships originally developed for homogeneous porous media are extended to represent hydraulic behavior of dual (or multi) continuum SPM. Such models often result in inconsistencies due to lack of consid- erationof structuralporespacegeometry andtheneglect ofunderlying physical processes governing liquid retention and flow under unsatu- rated conditions. We review a new framework that considers equilib- rium liquid configurations in dual continuum pore space as the basis for calculation of liquid saturation and subsequent introduction of hydrodynamic considerations. The SPM pore space is represented by a bimodal distribution of pore sizes, reflecting two disparate popula- tions of matrix and structural pores. Three steady-state and laminar flow regimes are considered to derive unsaturated hydraulic conduc- tivity functions: (i) flow in completely filled pore spaces, (ii) corner flow in partially filled pores and grooves, and (iii) film flow on solid surfaces. Two key assumptions are used in deriving the average cross- sectional flow velocities in these regimes: (i) that equilibrium liquid- vapor interfaces remain stable under slow laminar flows and (ii) that flow pathways are parallel. Liquid-vapor interfacial configurations for different matric potentials are calculated and statistically upscaled to derive sample-scale saturated and unsaturated hydraulic conductivity from velocity expressions weighted by the appropriate liquid-occupied cross-sectional areas, neglecting three-dimensional (3-D) network ef- fects. Similarly, the hydraulic functions for matrix and structural pores are derived separately and later combined by weighting the individual contributions by the porosities of the associated pore spaces. A param- eter estimation scheme was developed to calculate liquid saturation and to predict sample-scale unsaturated hydraulic conductivity. Model evaluation using measured data for homogeneous porous media, frac- tured welded tuff, and macroporous and aggregated soils shows favor- able agreement (within the limitations of model assumptions). Effects of nonequilibrium conditions between matrix and structural pore do- mains on the hydraulic conductivity and approximate consideration of 3-D network effects are discussed.

Original languageEnglish (US)
Pages (from-to)14-37
Number of pages24
JournalVadose Zone Journal
Volume1
Issue number1
StatePublished - Aug 2002
Externally publishedYes

Fingerprint

unsaturated hydraulic conductivity
porous media
porous medium
hydraulic conductivity
liquid
liquids
pore space
laminar flow
fluid mechanics
hydraulics
vapors
matrix
saturation
tuff
matric potential
saturated hydraulic conductivity
films (materials)
flow velocity
hydrodynamics
porosity

ASJC Scopus subject areas

  • Soil Science

Cite this

Unsaturated hydraulic conductivity of structured porous media : A review of liquid configuration-based models. / Tuller, Markus; Or, Dani.

In: Vadose Zone Journal, Vol. 1, No. 1, 08.2002, p. 14-37.

Research output: Contribution to journalArticle

@article{f2aa8162f05244a098b299184251fe12,
title = "Unsaturated hydraulic conductivity of structured porous media: A review of liquid configuration-based models",
abstract = "Common approaches for modeling hydraulic functions of unsatu- rated structured porous media (SPM) rely on macroscopic continuum representation, where parameterization schemesand constitutive rela- tionships originally developed for homogeneous porous media are extended to represent hydraulic behavior of dual (or multi) continuum SPM. Such models often result in inconsistencies due to lack of consid- erationof structuralporespacegeometry andtheneglect ofunderlying physical processes governing liquid retention and flow under unsatu- rated conditions. We review a new framework that considers equilib- rium liquid configurations in dual continuum pore space as the basis for calculation of liquid saturation and subsequent introduction of hydrodynamic considerations. The SPM pore space is represented by a bimodal distribution of pore sizes, reflecting two disparate popula- tions of matrix and structural pores. Three steady-state and laminar flow regimes are considered to derive unsaturated hydraulic conduc- tivity functions: (i) flow in completely filled pore spaces, (ii) corner flow in partially filled pores and grooves, and (iii) film flow on solid surfaces. Two key assumptions are used in deriving the average cross- sectional flow velocities in these regimes: (i) that equilibrium liquid- vapor interfaces remain stable under slow laminar flows and (ii) that flow pathways are parallel. Liquid-vapor interfacial configurations for different matric potentials are calculated and statistically upscaled to derive sample-scale saturated and unsaturated hydraulic conductivity from velocity expressions weighted by the appropriate liquid-occupied cross-sectional areas, neglecting three-dimensional (3-D) network ef- fects. Similarly, the hydraulic functions for matrix and structural pores are derived separately and later combined by weighting the individual contributions by the porosities of the associated pore spaces. A param- eter estimation scheme was developed to calculate liquid saturation and to predict sample-scale unsaturated hydraulic conductivity. Model evaluation using measured data for homogeneous porous media, frac- tured welded tuff, and macroporous and aggregated soils shows favor- able agreement (within the limitations of model assumptions). Effects of nonequilibrium conditions between matrix and structural pore do- mains on the hydraulic conductivity and approximate consideration of 3-D network effects are discussed.",
author = "Markus Tuller and Dani Or",
year = "2002",
month = "8",
language = "English (US)",
volume = "1",
pages = "14--37",
journal = "Vadose Zone Journal",
issn = "1539-1663",
publisher = "Soil Science Society of America",
number = "1",

}

TY - JOUR

T1 - Unsaturated hydraulic conductivity of structured porous media

T2 - A review of liquid configuration-based models

AU - Tuller, Markus

AU - Or, Dani

PY - 2002/8

Y1 - 2002/8

N2 - Common approaches for modeling hydraulic functions of unsatu- rated structured porous media (SPM) rely on macroscopic continuum representation, where parameterization schemesand constitutive rela- tionships originally developed for homogeneous porous media are extended to represent hydraulic behavior of dual (or multi) continuum SPM. Such models often result in inconsistencies due to lack of consid- erationof structuralporespacegeometry andtheneglect ofunderlying physical processes governing liquid retention and flow under unsatu- rated conditions. We review a new framework that considers equilib- rium liquid configurations in dual continuum pore space as the basis for calculation of liquid saturation and subsequent introduction of hydrodynamic considerations. The SPM pore space is represented by a bimodal distribution of pore sizes, reflecting two disparate popula- tions of matrix and structural pores. Three steady-state and laminar flow regimes are considered to derive unsaturated hydraulic conduc- tivity functions: (i) flow in completely filled pore spaces, (ii) corner flow in partially filled pores and grooves, and (iii) film flow on solid surfaces. Two key assumptions are used in deriving the average cross- sectional flow velocities in these regimes: (i) that equilibrium liquid- vapor interfaces remain stable under slow laminar flows and (ii) that flow pathways are parallel. Liquid-vapor interfacial configurations for different matric potentials are calculated and statistically upscaled to derive sample-scale saturated and unsaturated hydraulic conductivity from velocity expressions weighted by the appropriate liquid-occupied cross-sectional areas, neglecting three-dimensional (3-D) network ef- fects. Similarly, the hydraulic functions for matrix and structural pores are derived separately and later combined by weighting the individual contributions by the porosities of the associated pore spaces. A param- eter estimation scheme was developed to calculate liquid saturation and to predict sample-scale unsaturated hydraulic conductivity. Model evaluation using measured data for homogeneous porous media, frac- tured welded tuff, and macroporous and aggregated soils shows favor- able agreement (within the limitations of model assumptions). Effects of nonequilibrium conditions between matrix and structural pore do- mains on the hydraulic conductivity and approximate consideration of 3-D network effects are discussed.

AB - Common approaches for modeling hydraulic functions of unsatu- rated structured porous media (SPM) rely on macroscopic continuum representation, where parameterization schemesand constitutive rela- tionships originally developed for homogeneous porous media are extended to represent hydraulic behavior of dual (or multi) continuum SPM. Such models often result in inconsistencies due to lack of consid- erationof structuralporespacegeometry andtheneglect ofunderlying physical processes governing liquid retention and flow under unsatu- rated conditions. We review a new framework that considers equilib- rium liquid configurations in dual continuum pore space as the basis for calculation of liquid saturation and subsequent introduction of hydrodynamic considerations. The SPM pore space is represented by a bimodal distribution of pore sizes, reflecting two disparate popula- tions of matrix and structural pores. Three steady-state and laminar flow regimes are considered to derive unsaturated hydraulic conduc- tivity functions: (i) flow in completely filled pore spaces, (ii) corner flow in partially filled pores and grooves, and (iii) film flow on solid surfaces. Two key assumptions are used in deriving the average cross- sectional flow velocities in these regimes: (i) that equilibrium liquid- vapor interfaces remain stable under slow laminar flows and (ii) that flow pathways are parallel. Liquid-vapor interfacial configurations for different matric potentials are calculated and statistically upscaled to derive sample-scale saturated and unsaturated hydraulic conductivity from velocity expressions weighted by the appropriate liquid-occupied cross-sectional areas, neglecting three-dimensional (3-D) network ef- fects. Similarly, the hydraulic functions for matrix and structural pores are derived separately and later combined by weighting the individual contributions by the porosities of the associated pore spaces. A param- eter estimation scheme was developed to calculate liquid saturation and to predict sample-scale unsaturated hydraulic conductivity. Model evaluation using measured data for homogeneous porous media, frac- tured welded tuff, and macroporous and aggregated soils shows favor- able agreement (within the limitations of model assumptions). Effects of nonequilibrium conditions between matrix and structural pore do- mains on the hydraulic conductivity and approximate consideration of 3-D network effects are discussed.

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

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

M3 - Article

AN - SCOPUS:0013213650

VL - 1

SP - 14

EP - 37

JO - Vadose Zone Journal

JF - Vadose Zone Journal

SN - 1539-1663

IS - 1

ER -