A mathematical and computational study of the dispersivity tensor in anisotropic porous media

Dispersive transport in porous media is usually described through a Fickian model, in which the flux is the product of a dispersion tensor times the concentration gradient. This model is based on certain implicit assumptions, including slowly varying conditions. About fifty years ago, it was first suggested that the parameterization of the second-order dispersion tensor for anisotropic porous media involves a fourth-order dispersivity tensor. However, the properties of the dispersivity tensor have not been adequately studied. This work contributes to achieving a better grasp of dispersion in anisotropic porous media through a number of ways. First, with clearly stated assumptions and from first principles, we use the method of moments to derive a mathematical formula for the fourth-order dispersivity tensor, and show that it is a function of pore geometry, fluid velocity, and pore diffusion. Second, by using pore-scale flow and transport simulations through orderly and randomly packed 2-D and 3-D porous media, we evaluate the effects of the three factors on dispersivity. Different relationships with the Peclet number are observed for the longitudinal and transverse dispersivities and for orderly and randomly packed media. Third, we discuss the limitations of 2-D periodic media with simple structures in computing transverse dispersivity, which is more accurately predicted in the 3-D periodic media and 2-D randomly packed media. Fourth, we exhibit through numerical simulations that the method of moments can, computational limitations notwithstanding, be extended to stationary porous media. (C) 2013 Elsevier Ltd. All rights reserved.