THE INFLUENCE OF PORE-SCALE DISPERSION ON CONCENTRATION STATISTICAL MOMENTS IN TRANSPORT THROUGH HETEROGENEOUS AQUIFERS

Authors
Citation
G. Dagan et A. Fiori, THE INFLUENCE OF PORE-SCALE DISPERSION ON CONCENTRATION STATISTICAL MOMENTS IN TRANSPORT THROUGH HETEROGENEOUS AQUIFERS, Water resources research, 33(7), 1997, pp. 1595-1605
Citations number
14
Categorie Soggetti
Limnology,"Environmental Sciences","Water Resources
Journal title
ISSN journal
00431397
Volume
33
Issue
7
Year of publication
1997
Pages
1595 - 1605
Database
ISI
SICI code
0043-1397(1997)33:7<1595:TIOPDO>2.0.ZU;2-T
Abstract
Transport of an inert solute in a heterogeneous aquifer is governed by two mechanisms: advection by the random velocity field V(x) and pore- scale dispersion of coefficients D-dij. The velocity field is assumed to be stationary and of constant mean U and of correlation scale I muc h larger than the pore-scale d. It is assumed that D-dij = alpha(dij)U are constant. The relative effect of the two mechanisms is quantified by the Peclet numbers Pe(ij) = UI/D-dij = I/alpha(dij), which as a ru le are much larger than unity. The main aim of the study is to determi ne the impact of finite, though high, Pe on [C] and sigma(C)(2), the c oncentration mean and variance, respectively. The solution, derived in the past, for Pe = infinity is reconsidered first. By assuming a norm al X probability density function (p.d.f.), closed form solutions are obtained for [C] and sigma(C)(2). Recasting the problem in an Eulerian framework leads to the same results if certain closure conditions are adopted. The concentration moments for a finite Pe are derived subseq uently in a Lagrangean framework. The pore-scale dispersion is viewed as a Brownian motion type of displacement X-d of solute subparticles, of scale smaller than d, added to the advective displacements X. By ad opting again a normal p.d.f. for the latter, explicit expressions for [C] and sigma(C)(2) are obtained in terms of quadratures over the join t p.d.f. of advective two particles trajectories. While the influence of high Pe on [C] is generally small, it has a significant impact on s igma(C)(2). Simple results are obtained for a small V-0, for which tra jectories are fully correlated. In particular, the concentration coeff icient of variation at the center tends to a constant value for large time. Comparison of the present solution, obtained in terms of a quadr ature, with the Monte Carlo simulations of Graham. and McLaughlin [198 9] shows a very good agreement.