Numerical Simulation of the Reactive Transport at the Pore Scale

Abstract

We present a numerical algorithm for simulation of the chemical fluid-solid interaction at the pore scale, focusing on the core matrix dissolution and secondary mineral precipitation. The algorithm is based on the explicit use of chemical kinetics of the heterogeneous reactions to determine the rate of the fluid-solid interface changes. Fluid flow and chemical transport are simulated by finite differences. In contrast, level-set methods with immersed boundary conditions are applied to account for arbitrary interface position (not aligned to the grid lines).