A benchmark for the numerical simulation of pore-scale mineral dissolution with fluid-solid interface evolution

International audienceThis talk presents a benchmark problem for the simulation of single-phase flow, reactive transport and solid geometry evolution at the pore scale. The problem is organized in three parts that focus on specific aspects: flow and reactive transport, dissolution-driven geometry evolution in two dimensions, and a three-dimensional dissolution-driven geometry evolution including an experimental validation. Five codes are used to obtain the solution to this benchmark problem, including Chombo-Crunch, OpenFOAM-DBS, a lattice Boltzman code, Vortex method, and dissolFoam. These codes cover a good portion of the wide range of approaches typically employed for solving pore-scale problems in the literature, including discretization methods, characterization of the fluid-solid interfaces, and methods to move these interfaces as a result of fluid-solid reactions. Among the discretization methods, one can find finite volumes, Lagrangian methods, grid based methods (stochastic and deterministic). Results from the simulations performed by the five codes show remarkable agreement both quantitatively based on upscaled parameters such as surface area, solid volume and effective reaction rate and qualitatively based on comparisons of shape evolution. This outcome is especially notable given the disparity of approaches used by the codes.References[1] S. Molins, C. Soulaine, N. I. Prasianakis, A. Abbasi, P. Poncet, A. J. C. Ladd, V.Starchenko, S. Roman, D. Trebotich, H. A. Tchelepi, and C. I. Steefel,Simulation ofmineral dissolution at the pore scale with evolving fluid-solid interfaces: Review of ap-proaches and benchmark problem set, to appear in Computational Geosciences, 2019