Subterranean treatment formation using a model which takes into account the pore level physics by coupling the local pore scale phenomena to the macroscopic variables (Darcy velocity, pressure and reactant cup-mixing concentration) through the structure-property relationships (permeability-porosity, average pore size-porosity and interfacial area-porosity) and the dependence of the fluid-solid mass transfer coefficient and fluid phase dispersion coefficient on the evolving pore scale variables (average pore size, local Reynolds and Schmidt numbers).