Numerical Simulation of the Reactive Transport at Pore Scale in 3D.

The paper presents a finite-difference algorithm for reactive transport simulation at the pore scale in 3D. We simulate the matrix dissolution and crystal precipitation due to heterogeneous reactions, acquired at the fluid-solid interface. The fluid flow and the reactive transport are computed using finite difference method on a regular rectangular mesh, whereas the immersed boundary conditions are applied to account for irregular interface geometry. The time-evolving pore-to-matrix surface is defined implicitly by the level-set method. The algorithm is implemented using Graphic Processor Units.