Numerical modeling of chemical interaction between a fluid and rocks

A new algorithm for the numerical modeling of chemical fluid-rock interaction at the pore scale is proposed. The algorithm is based on splitting the problem into physical processes. It is assumed that the fluid rate is low and the fluid flow is stabilized almost instantly in the case of small changes in the pore space geometry. In the pore space, thus, the fluid flow is modeled using the Stokes equation for steady flows. The chemical reactant transport is described by the convection-diffusion equation with Robin boundary conditions at the fluid-rock interface. The pore space boundary changes with time and is implicitly given by a level-set function. We use finite-difference schemes with immersed boundary conditions to solve the Stokes and convection-diffusion equations