Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media: Theoretical framework

Eric Morales-Casique, Shlomo P Neuman, Alberto Guadagnini

Research output: Contribution to journalArticle

65 Citations (Scopus)

Abstract

Solute transport in randomly heterogeneous media is described by stochastic transport equations that are typically solved via Monte Carlo simulation. A promising alternative is to solve a corresponding system of statistical moment equations directly. We present exact (though not closed) implicit conditional first and second moment equations for advective-dispersive transport in finite domains. The velocity and concentration are generally non-stationary due to possible trends in heterogeneity, conditioning on data, temporal variations in velocity, fluid and/or solute sources, initial and boundary conditions. Our equations are integro-differential and include non-local parameters depending on more than one point in space-time. To allow solving these equations, we close them by perturbation and develop recursive moment equations in Laplace space for the special case of steady state flow, to second order in σY where σY2 is a measure of (natural) log hydraulic conductivity variance. We also propose a higher-order iterative closure. Our recursive equations and iterative closure are formally valid for mildly heterogeneous media, or well-conditioned strongly heterogeneous media in which the random component of heterogeneity is relatively small. The non-local moment equations suggest (and a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] demonstrates numerically) that, in general, transport cannot be validly described by means of Fick's law with a (constant or variable) macrodispersion coefficient. We show nevertheless that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that varies generally in space-time. In a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] we present a high-accuracy computational algorithm for our iterative non-local and recursive localized moment equations, assessing their accuracy and computational efficiency in comparison to unconditional and conditional Monte Carlo simulations.

Original languageEnglish (US)
Pages (from-to)1238-1255
Number of pages18
JournalAdvances in Water Resources
Volume29
Issue number8
DOIs
StatePublished - Aug 2006

Fingerprint

solute transport
porous medium
heterogeneous medium
advection
conditioning
simulation
hydraulic conductivity
solute
temporal variation
boundary condition
perturbation
water
fluid

Keywords

  • Advection
  • Dispersion
  • Moments
  • Random media
  • Stochastic
  • Transport

ASJC Scopus subject areas

  • Earth-Surface Processes

Cite this

Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media : Theoretical framework. / Morales-Casique, Eric; Neuman, Shlomo P; Guadagnini, Alberto.

In: Advances in Water Resources, Vol. 29, No. 8, 08.2006, p. 1238-1255.

Research output: Contribution to journalArticle

@article{da85df7126f8419d9906095654b3999d,
title = "Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media: Theoretical framework",
abstract = "Solute transport in randomly heterogeneous media is described by stochastic transport equations that are typically solved via Monte Carlo simulation. A promising alternative is to solve a corresponding system of statistical moment equations directly. We present exact (though not closed) implicit conditional first and second moment equations for advective-dispersive transport in finite domains. The velocity and concentration are generally non-stationary due to possible trends in heterogeneity, conditioning on data, temporal variations in velocity, fluid and/or solute sources, initial and boundary conditions. Our equations are integro-differential and include non-local parameters depending on more than one point in space-time. To allow solving these equations, we close them by perturbation and develop recursive moment equations in Laplace space for the special case of steady state flow, to second order in σY where σY2 is a measure of (natural) log hydraulic conductivity variance. We also propose a higher-order iterative closure. Our recursive equations and iterative closure are formally valid for mildly heterogeneous media, or well-conditioned strongly heterogeneous media in which the random component of heterogeneity is relatively small. The non-local moment equations suggest (and a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] demonstrates numerically) that, in general, transport cannot be validly described by means of Fick's law with a (constant or variable) macrodispersion coefficient. We show nevertheless that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that varies generally in space-time. In a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] we present a high-accuracy computational algorithm for our iterative non-local and recursive localized moment equations, assessing their accuracy and computational efficiency in comparison to unconditional and conditional Monte Carlo simulations.",
keywords = "Advection, Dispersion, Moments, Random media, Stochastic, Transport",
author = "Eric Morales-Casique and Neuman, {Shlomo P} and Alberto Guadagnini",
year = "2006",
month = "8",
doi = "10.1016/j.advwatres.2005.10.002",
language = "English (US)",
volume = "29",
pages = "1238--1255",
journal = "Advances in Water Resources",
issn = "0309-1708",
publisher = "Elsevier Limited",
number = "8",

}

TY - JOUR

T1 - Non-local and localized analyses of non-reactive solute transport in bounded randomly heterogeneous porous media

T2 - Theoretical framework

AU - Morales-Casique, Eric

AU - Neuman, Shlomo P

AU - Guadagnini, Alberto

PY - 2006/8

Y1 - 2006/8

N2 - Solute transport in randomly heterogeneous media is described by stochastic transport equations that are typically solved via Monte Carlo simulation. A promising alternative is to solve a corresponding system of statistical moment equations directly. We present exact (though not closed) implicit conditional first and second moment equations for advective-dispersive transport in finite domains. The velocity and concentration are generally non-stationary due to possible trends in heterogeneity, conditioning on data, temporal variations in velocity, fluid and/or solute sources, initial and boundary conditions. Our equations are integro-differential and include non-local parameters depending on more than one point in space-time. To allow solving these equations, we close them by perturbation and develop recursive moment equations in Laplace space for the special case of steady state flow, to second order in σY where σY2 is a measure of (natural) log hydraulic conductivity variance. We also propose a higher-order iterative closure. Our recursive equations and iterative closure are formally valid for mildly heterogeneous media, or well-conditioned strongly heterogeneous media in which the random component of heterogeneity is relatively small. The non-local moment equations suggest (and a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] demonstrates numerically) that, in general, transport cannot be validly described by means of Fick's law with a (constant or variable) macrodispersion coefficient. We show nevertheless that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that varies generally in space-time. In a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] we present a high-accuracy computational algorithm for our iterative non-local and recursive localized moment equations, assessing their accuracy and computational efficiency in comparison to unconditional and conditional Monte Carlo simulations.

AB - Solute transport in randomly heterogeneous media is described by stochastic transport equations that are typically solved via Monte Carlo simulation. A promising alternative is to solve a corresponding system of statistical moment equations directly. We present exact (though not closed) implicit conditional first and second moment equations for advective-dispersive transport in finite domains. The velocity and concentration are generally non-stationary due to possible trends in heterogeneity, conditioning on data, temporal variations in velocity, fluid and/or solute sources, initial and boundary conditions. Our equations are integro-differential and include non-local parameters depending on more than one point in space-time. To allow solving these equations, we close them by perturbation and develop recursive moment equations in Laplace space for the special case of steady state flow, to second order in σY where σY2 is a measure of (natural) log hydraulic conductivity variance. We also propose a higher-order iterative closure. Our recursive equations and iterative closure are formally valid for mildly heterogeneous media, or well-conditioned strongly heterogeneous media in which the random component of heterogeneity is relatively small. The non-local moment equations suggest (and a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] demonstrates numerically) that, in general, transport cannot be validly described by means of Fick's law with a (constant or variable) macrodispersion coefficient. We show nevertheless that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that varies generally in space-time. In a companion paper [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: computational analysis. Adv Water Resour, submitted for publication] we present a high-accuracy computational algorithm for our iterative non-local and recursive localized moment equations, assessing their accuracy and computational efficiency in comparison to unconditional and conditional Monte Carlo simulations.

KW - Advection

KW - Dispersion

KW - Moments

KW - Random media

KW - Stochastic

KW - Transport

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

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

U2 - 10.1016/j.advwatres.2005.10.002

DO - 10.1016/j.advwatres.2005.10.002

M3 - Article

AN - SCOPUS:33746357914

VL - 29

SP - 1238

EP - 1255

JO - Advances in Water Resources

JF - Advances in Water Resources

SN - 0309-1708

IS - 8

ER -