TRANSIENT FLOW IN BOUNDED RANDOMLY HETEROGENEOUS DOMAINS-1 - EXACT CONDITIONAL MOMENT EQUATIONS AND RECURSIVE APPROXIMATIONS

Citation
Dm. Tartakovsky et Sp. Neuman, TRANSIENT FLOW IN BOUNDED RANDOMLY HETEROGENEOUS DOMAINS-1 - EXACT CONDITIONAL MOMENT EQUATIONS AND RECURSIVE APPROXIMATIONS, Water resources research, 34(1), 1998, pp. 1-12
Citations number
12
Categorie Soggetti
Limnology,"Environmental Sciences","Water Resources
Journal title
ISSN journal
00431397
Volume
34
Issue
1
Year of publication
1998
Pages
1 - 12
Database
ISI
SICI code
0043-1397(1998)34:1<1:TFIBRH>2.0.ZU;2-T
Abstract
We consider the effect of measuring randomly varying local hydraulic c onductivities K(x) on one's ability to predict transient flow within b ounded domains, driven by random sources, initial head, and boundary c onditions. Our aim is to allow optimum unbiased prediction of local hy draulic heads h(x, t) and Darcy fluxes q(x, t) by means of their ensem ble moments, [h(x, t)](c) and [q(x, t)](c), conditioned on measurement s of K(x). We show that these predictors satisfy a compact determinist ic flow equation which contains a space-time integrodifferential ''res idual flux'' term. This term renders [q(g, t)](c) nonlocal and non-Dar cian so that the concept of effective hydraulic conductivity looses me aning in all but a few special cases, Instead, the residual flux conta ins kernels that constitute nonlocal parameters in space-time that are additionally conditional on hydraulic conductivity data and thus nonu nique. The kernels include symmetric and nonsymmetric second-rank tens ors as well as vectors, We also develop nonlocal equations for second conditional moments of head and flux which constitute measures of pred ictive uncertainty. The nonlocal expressions cannot be evaluated direc tly without either a closure approximation or high-resolution conditio nal Monte Carlo simulation. To render our theory workable, we develop recursive closure approximations for the moment equations through expa nsion in powers of a small parameter which represents the standard est imation error of natural log K(x). These approximations are valid to a rbitrary order for either mildly heterogeneous or well-conditioned str ongly heterogeneous media. They allow, in principle, evaluating the co nditional moments numerically on relatively coarse grids, without upsc aling, by standard methods such as finite elements.