Transient flow in bounded randomly heterogeneous domains 1. Exact conditional moment equations and recursive approximations

Daniel M. Tartakovsky, Shlomo P Neuman

Research output: Contribution to journalArticle

84 Citations (Scopus)

Abstract

We consider the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict transient flow within bounded domains, driven by random sources, initial head; and boundary conditions. Our aim is to allow optimum unbiased prediction of local hydraulic heads h(x, t) and Darcy fluxes q(x, t) by means of their ensemble moments, <h(x, t)>(c) and <q(x, t)>(c), conditioned on measurements of K(x). We show that these predictors satisfy a compact deterministic flow equation which contains a space-time integrodifferential 'residual flux' term. This term renders <q(x, t)>(c) nonlocal and non-Darcian so that the concept of effective hydraulic conductivity looses meaning in all but a few special cases. Instead, the residual flux contains kernels that constitute nonlocal parameters in space-time that are additionally conditional on hydraulic conductivity data and thus nonunique. The kernels include symmetric and nonsymmetric second-rank tensors as well as vectors. We also develop nonlocal equations for second conditional moments of head and flux which constitute measures of predictive uncertainty. The nonlocal expressions cannot be evaluated directly without either a closure approximation or high-resolution conditional Monte Carlo simulation. To render our theory workable, we develop recursive closure approximations for the moment equations through expansion in powers of a small parameter which represents the standard estimation error of natural log K(x). These approximations are valid to arbitrary order for either mildly heterogeneous or well-conditioned strongly heterogeneous media. They allow, in principle, evaluating the conditional moments numerically on relatively coarse grids, without upscaling, by standard methods such as finite elements.

Original languageEnglish (US)
Pages (from-to)1-12
Number of pages12
JournalWater Resources Research
Volume34
Issue number1
StatePublished - Jan 1998

Fingerprint

transient flow
hydraulic conductivity
Hydraulic conductivity
Fluxes
heterogeneous medium
upscaling
hydraulic head
seeds
fluid mechanics
boundary condition
uncertainty
Error analysis
Tensors
prediction
Hydraulics
Boundary conditions
simulation
parameter
methodology

ASJC Scopus subject areas

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

Cite this

Transient flow in bounded randomly heterogeneous domains 1. Exact conditional moment equations and recursive approximations. / Tartakovsky, Daniel M.; Neuman, Shlomo P.

In: Water Resources Research, Vol. 34, No. 1, 01.1998, p. 1-12.

Research output: Contribution to journalArticle

@article{38fd33b7ea40401c8946b4ac524f0a84,
title = "Transient flow in bounded randomly heterogeneous domains 1. Exact conditional moment equations and recursive approximations",
abstract = "We consider the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict transient flow within bounded domains, driven by random sources, initial head; and boundary conditions. Our aim is to allow optimum unbiased prediction of local hydraulic heads h(x, t) and Darcy fluxes q(x, t) by means of their ensemble moments, (c) and (c), conditioned on measurements of K(x). We show that these predictors satisfy a compact deterministic flow equation which contains a space-time integrodifferential 'residual flux' term. This term renders (c) nonlocal and non-Darcian so that the concept of effective hydraulic conductivity looses meaning in all but a few special cases. Instead, the residual flux contains kernels that constitute nonlocal parameters in space-time that are additionally conditional on hydraulic conductivity data and thus nonunique. The kernels include symmetric and nonsymmetric second-rank tensors as well as vectors. We also develop nonlocal equations for second conditional moments of head and flux which constitute measures of predictive uncertainty. The nonlocal expressions cannot be evaluated directly without either a closure approximation or high-resolution conditional Monte Carlo simulation. To render our theory workable, we develop recursive closure approximations for the moment equations through expansion in powers of a small parameter which represents the standard estimation error of natural log K(x). These approximations are valid to arbitrary order for either mildly heterogeneous or well-conditioned strongly heterogeneous media. They allow, in principle, evaluating the conditional moments numerically on relatively coarse grids, without upscaling, by standard methods such as finite elements.",
author = "Tartakovsky, {Daniel M.} and Neuman, {Shlomo P}",
year = "1998",
month = "1",
language = "English (US)",
volume = "34",
pages = "1--12",
journal = "Water Resources Research",
issn = "0043-1397",
publisher = "American Geophysical Union",
number = "1",

}

TY - JOUR

T1 - Transient flow in bounded randomly heterogeneous domains 1. Exact conditional moment equations and recursive approximations

AU - Tartakovsky, Daniel M.

AU - Neuman, Shlomo P

PY - 1998/1

Y1 - 1998/1

N2 - We consider the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict transient flow within bounded domains, driven by random sources, initial head; and boundary conditions. Our aim is to allow optimum unbiased prediction of local hydraulic heads h(x, t) and Darcy fluxes q(x, t) by means of their ensemble moments, (c) and (c), conditioned on measurements of K(x). We show that these predictors satisfy a compact deterministic flow equation which contains a space-time integrodifferential 'residual flux' term. This term renders (c) nonlocal and non-Darcian so that the concept of effective hydraulic conductivity looses meaning in all but a few special cases. Instead, the residual flux contains kernels that constitute nonlocal parameters in space-time that are additionally conditional on hydraulic conductivity data and thus nonunique. The kernels include symmetric and nonsymmetric second-rank tensors as well as vectors. We also develop nonlocal equations for second conditional moments of head and flux which constitute measures of predictive uncertainty. The nonlocal expressions cannot be evaluated directly without either a closure approximation or high-resolution conditional Monte Carlo simulation. To render our theory workable, we develop recursive closure approximations for the moment equations through expansion in powers of a small parameter which represents the standard estimation error of natural log K(x). These approximations are valid to arbitrary order for either mildly heterogeneous or well-conditioned strongly heterogeneous media. They allow, in principle, evaluating the conditional moments numerically on relatively coarse grids, without upscaling, by standard methods such as finite elements.

AB - We consider the effect of measuring randomly varying local hydraulic conductivities K(x) on one's ability to predict transient flow within bounded domains, driven by random sources, initial head; and boundary conditions. Our aim is to allow optimum unbiased prediction of local hydraulic heads h(x, t) and Darcy fluxes q(x, t) by means of their ensemble moments, (c) and (c), conditioned on measurements of K(x). We show that these predictors satisfy a compact deterministic flow equation which contains a space-time integrodifferential 'residual flux' term. This term renders (c) nonlocal and non-Darcian so that the concept of effective hydraulic conductivity looses meaning in all but a few special cases. Instead, the residual flux contains kernels that constitute nonlocal parameters in space-time that are additionally conditional on hydraulic conductivity data and thus nonunique. The kernels include symmetric and nonsymmetric second-rank tensors as well as vectors. We also develop nonlocal equations for second conditional moments of head and flux which constitute measures of predictive uncertainty. The nonlocal expressions cannot be evaluated directly without either a closure approximation or high-resolution conditional Monte Carlo simulation. To render our theory workable, we develop recursive closure approximations for the moment equations through expansion in powers of a small parameter which represents the standard estimation error of natural log K(x). These approximations are valid to arbitrary order for either mildly heterogeneous or well-conditioned strongly heterogeneous media. They allow, in principle, evaluating the conditional moments numerically on relatively coarse grids, without upscaling, by standard methods such as finite elements.

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

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

M3 - Article

VL - 34

SP - 1

EP - 12

JO - Water Resources Research

JF - Water Resources Research

SN - 0043-1397

IS - 1

ER -