Instability of the Dissolution Front in a Smooth Fracture

In the standard model for fracture dissolution, the aperture is assumed to be constant in the direction perpendicular to the flow [1]. However, such models do not predict erosion deep into the fracture without a strong concentration dependence of the reaction rate. This kinetic "trigger" mechanism for Karst formation is well verified in calcite [1]; near saturation it is known that the dissolution rate drops by several orders of magnitude. However in other minerals, for instance Gypsum, the rate coefficient is essentially independent of concentration yet caverns develop in these formations as well.

However, we have recently found that a planar dissolution front in a smooth fracture is unstable to infinitessimally small perturbations [2]. The most striking result of the linear stability analysis is that there is a maximal growth rate at a specific wavelength. We have confirmed this prediction by numerical simulations; an example of a growing instability is illustrated in the figure below.

As a result of the instability, an individual fracture develops a strongly heterogeneous permeability during dissolution, with an inherent length scale that depends on the kinetics and flow rate, but not the initial topography. Furthermore, there is no lower limit to the reaction rate for unstable dissolution. Only the wavelength and penetration depth are affected, scaling with the inverse of the Damköhler number.

References

W. Dreybrodt, Water Resources Res., 32:2923, 1996.

P. Szymczak and A. J. C. Ladd, EPSL, 201:424-432, 2011.