Premium
Solute transport in three‐dimensional heterogeneous media with a Gaussian covariance of log hydraulic conductivity
Author(s) -
Zhang YouKuan,
Federico Vittorio Di
Publication year - 1998
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/98wr01142
Subject(s) - covariance function , covariance , mathematics , mathematical analysis , logarithm , hydraulic conductivity , tensor (intrinsic definition) , gaussian , law of total covariance , statistics , physics , geometry , covariance intersection , geology , quantum mechanics , soil science , soil water
The analytical solutions for the velocity covariance, u ij one particle displacement covariance X ij , and the macrodispersivity tensor α ij defined as (0.5/μ)( dX ij / dt ), were derived in three‐dimensional heterogeneous media. A Gaussian covariance function (GCF) of logarithmic hydraulic conductivity, log K , was used, assuming uniform mean flow and first‐order approximation in log‐conductivity variance, where μ is the magnitude of the mean flow velocity μ . Based on these solutions, the time‐dependent ensemble averages of the second spatial moments, Z ij ≡ 〈 A ij 〉 − A ij ∥0) = X ij − R ij and the effective dispersivity tensor γ ij defined as (0.5/μ)( d 〈 A ij 〉/ dt , were evaluated for a finite line source either normal or parallel to μ , where A ij (0) is the initial value of the second spatial moment of a plume, A ij and R ij is the plume centroid covariance. The results obtained in this study were compared with previous results for an exponential covariance function (ECF). It was found that in a stationary log K field the spreading of a solute plume depends not only on the variance and integral scale of the log K field but also on the shape of its covariance function. The more strongly correlated the hydraulic conductivities at short separation distances are, the faster Z ii and γ ii grow at early time. Also, the earlier that γ ii approaches its asymptote or peak, and the higher the peak is, the larger the negative transverse dispersivity. More importantly, the ergodic limits for GCF are reached faster than those for ECF, as the initial size of a plume increases. The ergodic limit X 11 for GCF is slightly larger than that for ECF, but X 22 and X 33 are significantly smaller than those for ECF even though the asymptotic α ii for GCF is the same with that for ECF.