Numerical study of geometrical dispersion in self-affine rough fractures

We report a numerical study of passive tracer dispersion in fractures with rough walls modeled as the space between two complementary self-affine surfaces rigidly translated with respect to each other. Geometrical dispersion due to the disorder of the velocity distribution is computed using the lubrication approximation. Using a spectral perturbative scheme to solve the flow problem and a mapping coordinate method to compute dispersion, we perform extensive ensemble averaged simulations to test theoretical predictions on the dispersion dependence on simple geometrical parameters. We observe the expected quadratic dispersion coefficient dependence on both the mean aperture and the relative shift of the crack as of well as the anomalous dispersion dependence on tracer traveling distance. We also characterize the anisotropy of the dispersion front, which progressively wrinkles into a self-affine curve whose exponent is equal to that of the fracture surface.