A theoretical and numerical framework for modeling gypsum cavity dissolution
Abstract
In this paper, a large-scale diffuse interface model is used to describe the evolution of a gypsum cavity formation induced by dissolution. The method is based upon the assumption of a pseudo-component dissolving with a thermodynamic equilibrium boundary condition. A methodology is proposed based on numerical computations with fixed boundaries in order to choose suitable parameters for the diffuse interface model, and hence predict the correct dissolution fluxes and surface recession velocity. Additional simulations were performed to check which type of momentum balance equation should be used. The numerical results did not show a strong impact of this choice for the typical initial boundary value problems under consideration. Calculations with a variable density and Boussinesq approximation were also performed to evaluate the potential for natural convection. The results showed that the impact of density driven flows was negligible in the cases under investigation. The potential of the methodology is illustrated on two large-scale configurations: one corresponding to a gypsum lens located strictly within a porous rock formation and the other to an isolated pillar in a flooded gypsum room and pillar quarry.