Assessing dual continuum method for multicomponent reactive transport

Abstract

The subsurface flow and reactive transport code PFLOTRAN can accommodate fracture-matrix interaction for multicomponent reactive transport using a Dual Continuum formulation. In this work we present two examples to assess this approach. Through these examples, we also illustrate how this Dual Continuum model can be parameterised and what kind of grid spacing is needed to accurately capture reactive transport in the secondary continuum. The results of the calculations show that when geochemical reactions in the matrix (or the secondary continuum) are included, a very fine grid spacing close to the fracture-matrix interface is required to properly capture geochemical gradients at the interface. We also show that one can reduce the computational cost by using a variable grid spacing, that limits the total number in the secondary continuum, and yet obtain accurate results. This feature makes the Dual Continuum approach suitable for large-scale reactive transport modelling in fractured crystalline rock.