Numerical Simulation of a Phase-Field Model for Reactive Transport in Porous Media

Abstract

We consider a Darcy-scale model for mineral precipitation and dissolution in a porous medium. This model is obtained by homogenization techniques starting at the scale of pores. The model is based on a phase-field approach to account for the evolution of the pore geometry and the outcome is a multi-scale strongly coupled non-linear system of equations. In this work we discuss a robust numerical scheme dealing with the scale separation in the model as well as the non-linear character of the equations. We combine mesh refinement with stable linearization techniques to illustrate the behaviour of the multi-scale iterative scheme.