A numerical scheme for the pore scale simulation of crystal dissolution and precipitation in porous media

In this paper we analyze a numerical scheme for a dissolution and precipitation model in porous media. We focus here on the chemistry, which is modeled by a parabolic problem that is coupled through the boundary conditions to an ordinary differential inclusion defined on the boundary. We use a regularization approach for constructing a semi-implicit scheme that is stable and convergent. For dealing with the emerging time discrete nonlinear problems, we propose a simple fixed-point iterative procedure. The paper is concluded by numerical results.