Numerical assessment of 3‐D macrodispersion in heterogeneous porous media. Issue 5 (22nd May 2013)

• Record Type: Journal Article
• Authors: Beaudoin, A.; de Dreuzy, J.‐R.
• Is Part Of: Water resources research. Volume 49:Issue 5(2013:May)
• Page Start: 2489
• Page End: 2496
• Publication Date: 2013-05-22
• Subjects: 3‐D macrodispersion -- parallel computing -- solute transport; Hydrology -- Periodicals; 333.91
• DOI: 10.1002/wrcr.20206
• Languages: English
• ISSNs: 0043-1397

Abstract

Key Points: Heterogeneity‐induced dispersion is much larger in 3D than in 2D Braiding of high‐velocity channels critically enhance transverse dispersion Longitudinal dispersion also comes from high correlations of the velocity field Abstract : [1] Hydrodynamic dispersion is a key controlling factor of solute transport in heterogeneous porous media that critically depends on dimensionality. The transverse macrodispersion (asymptotic dispersion transverse to the mean velocity direction) is known to vanish only in 2‐D and not in 3‐D. Using classical Gaussian correlated permeability fields with a lognormal distribution of variance σ Y 2, we determine numerically the longitudinal and transverse dispersivities as functions of heterogeneity and dimensionality. We show that the transverse macrodispersion steeply increases with σ Y 2 underlying the essential role of flow lines braiding, a mechanism specific to 3‐D systems that we qualitatively characterize by the increasing expansion of the flow lines transversally to the flow direction. The transverse macrodispersion remains however at least two orders of magnitude smaller than the longitudinal macrodispersion, which increases even more steeply with σ Y 2 . At moderate to high levels of heterogeneity, the transverse dispersion also converges much faster to its asymptotic regime than do the longitudinal dispersion. Braiding cannot be thus taken as the sole mechanism responsible for the high longitudinal macrodispersions. It could be either supplemented or superseded by stronger velocity correlations in 3‐D than in 2‐D. This assumption is supported by the much larger longitudinal macrodispersions obtained in 3‐D than in 2‐D, up to a factor of 5 for σ Y 2 = 6.25.