On the tensorial nature of advective porosity

Research output: Contribution to journalArticle

10 Citations (Scopus)

Abstract

Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ψ of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity φ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than φ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ≤ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.

Original languageEnglish (US)
Pages (from-to)149-159
Number of pages11
JournalAdvances in Water Resources
Volume28
Issue number2
DOIs
StatePublished - Feb 2005

Fingerprint

tracer
porosity
hydraulics
material form
dispersivity
fracture network
channel flow
void
anisotropy
plume
permeability
timescale
sampling
rock
soil
effect

Keywords

  • Advective porosity
  • Directional porous media
  • Fractured rocks
  • Tensor

ASJC Scopus subject areas

  • Earth-Surface Processes

Cite this

On the tensorial nature of advective porosity. / Neuman, Shlomo P.

In: Advances in Water Resources, Vol. 28, No. 2, 02.2005, p. 149-159.

Research output: Contribution to journalArticle

@article{379114fcd67d4f7db391c933a9a26e47,
title = "On the tensorial nature of advective porosity",
abstract = "Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ψ of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity φ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than φ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ≤ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.",
keywords = "Advective porosity, Directional porous media, Fractured rocks, Tensor",
author = "Neuman, {Shlomo P}",
year = "2005",
month = "2",
doi = "10.1016/j.advwatres.2004.09.008",
language = "English (US)",
volume = "28",
pages = "149--159",
journal = "Advances in Water Resources",
issn = "0309-1708",
publisher = "Elsevier Limited",
number = "2",

}

TY - JOUR

T1 - On the tensorial nature of advective porosity

AU - Neuman, Shlomo P

PY - 2005/2

Y1 - 2005/2

N2 - Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ψ of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity φ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than φ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ≤ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.

AB - Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ψ of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity φ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than φ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ≤ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.

KW - Advective porosity

KW - Directional porous media

KW - Fractured rocks

KW - Tensor

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

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

U2 - 10.1016/j.advwatres.2004.09.008

DO - 10.1016/j.advwatres.2004.09.008

M3 - Article

AN - SCOPUS:11944262691

VL - 28

SP - 149

EP - 159

JO - Advances in Water Resources

JF - Advances in Water Resources

SN - 0309-1708

IS - 2

ER -