Models of sub-grid variability in numerical simulations of solute transport in heterogeneous porous formations: three-dimensional flow and effect of pore-scale dispersion

Modeling transport of contaminants in the earths subsurface relies on numerical solutions over grids with blocks larger than Darcys scale. The hydraulic conductivity is homogenized over the grid blocks and the plumes spreading is reduced as a consequence of the wiped-out variability. To compensate for this loss Rubin et al. (1999) proposed to augment mixing by block-effective dispersion coefficients, and Rubin et al. (2003) showed, by means of two dimensional simulations, how this concept can be applied in practice. In this paper, we present new solutions of the block-effective dispersion tensor for an axisymmetric exponential covariance model. In addition, we discuss the influence of pore-scale dispersion in both two- and three-dimensional applications.