Dissolution-seepage coupled analysis through formations containing soluble materials

Seepage flow can dissolve particulate soluble materials contained in soil layers and rock formations. The above-mentioned dissolution increases the porosity of the formation and hence seepage flow, which in turn progressively increases the dissolution rate. Due to progressive dissolution, several dams around the world have lost functionality or even failed. Dissolution propagation can be modeled as progress of a solution front, with its progression and resulting excess seepage coupled in the analysis. This is made possible in this paper by simultaneously solving the governing differential equation of seepage and the equation expressing progress of the solution front. The outcome (coupled differential equation) is nonlinear and transient, since both porosity and coefficient of permeability vary with the advancement of the solution front through the medium. The finite-element method is used to solve the resulting nonlinear partial differential equation. Using several examples, influence of material properties and geometry characteristics on the solution front progress and the resulting excess seepage loss is evaluated. Furthermore, effectiveness of different countermeasures (e.g., positive cutoffs and their positions) in dam foundations are studied. Contaminant transport can also be easily modeled and analyzed after applying some modifications into the approach.