Abstract
We investigate the development of wormholes in laboratory-scale fractures using three-dimensional numerical simulations. Well-controlled initial conditions, involving a small perturbation near the inlet of an otherwise flat fracture aperture field, were used to make a systematic study of the effects of flow rate and reaction rate on the aperture evolution. We find at least two characteristic wormhole shapes, which can be grouped within a phase diagram in the space of Péclet and Damköhler numbers. We investigate how this phase diagram depends on fracture geometry, specifically the length (L) and width (W) in comparison to the initial aperture (h0). This information is used to determine an effective Damköhler number that leads to a Péclet and length-independent phase boundary between wormholes and uniform dissolution. We relate these observations to experimental studies of wormhole formation in dissolving fractures.
Original language | English |
---|---|
Pages (from-to) | 7946-7959 |
Number of pages | 14 |
Journal | Water Resources Research |
Volume | 54 |
Issue number | 10 |
DOIs | |
State | Published - Oct 2018 |
Funding
The data in this paper were generated with the OpenFOAM ® toolkit (v1706), http://www.openfoam.com/. Figures were prepared using Paraview visualization software, http://www.paraview.org/. Source codes and a sampling of the input files used to produce the data shown in this paper are provided as Supplementary Information. Additional output data is available on request.This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. A.J.C. Ladd was funded under Award Number DE-FG02-98ER14853.
Funders | Funder number |
---|---|
U.S. Department of Energy | |
Office of Science | |
Basic Energy Sciences | |
Chemical Sciences, Geosciences, and Biosciences Division | DE-FG02-98ER14853 |
Keywords
- fracture dissolution
- laboratory-scale fractures
- wormholes