TY - JOUR

T1 - The impact of local diffusion upon mass arrival of a passive solute in transport through three-dimensional highly heterogeneous aquifers

AU - Fiori, A.

AU - Jankovic, I.

AU - Dagan, G.

N1 - Funding Information:
Supercomputer simulations were run at Center for Computational Research, University at Buffalo that provided CPU time for our experiments. Partial support for this work was provided by the SERDP (Strategic Environmental Research and Development Program) under project number ER-1738.

PY - 2011/12

Y1 - 2011/12

N2 - Fluid flow and solute transport take place in a porous formation. The spatial variability of the hydraulic conductivity K has a large impact on the velocity field and advective spreading of solute. The logconductivity Y=lnK is modeled as a stationary random space function, of normal distribution of mean 〈Y〉=lnKG and variance σY2, and of isotropic autocorrelation of finite integral scale I. The stationary velocity field is of uniform mean U. Transport of an ergodic plume is quantified by the mass arrival (breakthrough curve) at control planes normal to U. First-order solutions in σY2, applicable to weak heterogeneity σY2≤1, have been investigated extensively in the past. However, many natural formations are highly heterogeneous σY2>1 and the more complex flow and transport problems have been attacked only recently. The present study investigates the impact of local diffusion, quantified by the Peclet number Pe=UI/D, upon the breakthrough curve (BTC), both by numerical simulations and an analytical approach based on the self consistent approximation. The latter is an extension to finite Peclet of a method that we have developed in the past, in which we solved the 3D problem of flow and advective transport by modeling the structure as an ensemble of densely packed spherical inclusions of uniform radius R and independent random Y. Unlike the first-order solution, diffusion has a large impact on the BTC, even for the common values Pe≫1. Two mechanisms associated with finite Pe are identified and quantified: removal of the solute captured by low K blocks (characterized by "slow" flow) and exchange between the fluid flowing past inclusions (i.e. "fast" flow) and their interior. The first mechanism affects primarily the tail of the BTC, while the second one impacts both the peak and the tail. All the parameters of the model U,σY2,Pe,I are physically based and depend on the structural and flow characteristics solely.

AB - Fluid flow and solute transport take place in a porous formation. The spatial variability of the hydraulic conductivity K has a large impact on the velocity field and advective spreading of solute. The logconductivity Y=lnK is modeled as a stationary random space function, of normal distribution of mean 〈Y〉=lnKG and variance σY2, and of isotropic autocorrelation of finite integral scale I. The stationary velocity field is of uniform mean U. Transport of an ergodic plume is quantified by the mass arrival (breakthrough curve) at control planes normal to U. First-order solutions in σY2, applicable to weak heterogeneity σY2≤1, have been investigated extensively in the past. However, many natural formations are highly heterogeneous σY2>1 and the more complex flow and transport problems have been attacked only recently. The present study investigates the impact of local diffusion, quantified by the Peclet number Pe=UI/D, upon the breakthrough curve (BTC), both by numerical simulations and an analytical approach based on the self consistent approximation. The latter is an extension to finite Peclet of a method that we have developed in the past, in which we solved the 3D problem of flow and advective transport by modeling the structure as an ensemble of densely packed spherical inclusions of uniform radius R and independent random Y. Unlike the first-order solution, diffusion has a large impact on the BTC, even for the common values Pe≫1. Two mechanisms associated with finite Pe are identified and quantified: removal of the solute captured by low K blocks (characterized by "slow" flow) and exchange between the fluid flowing past inclusions (i.e. "fast" flow) and their interior. The first mechanism affects primarily the tail of the BTC, while the second one impacts both the peak and the tail. All the parameters of the model U,σY2,Pe,I are physically based and depend on the structural and flow characteristics solely.

KW - Breakthrough curve

KW - Heterogeneous media

KW - Local diffusion

KW - Porous media

KW - Self consistent approximation

KW - Solute transport

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

U2 - 10.1016/j.advwatres.2011.08.010

DO - 10.1016/j.advwatres.2011.08.010

M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???

AN - SCOPUS:80054757628

SN - 0309-1708

VL - 34

SP - 1563

EP - 1573

JO - Advances in Water Resources

JF - Advances in Water Resources

IS - 12

ER -