Relative rates of fluid advection, elemental diffusion and replacement govern reaction front patterns

TL;DR

This study uses numerical modeling to explore how the interplay of advection, diffusion, and reaction rates influences the formation of reaction front patterns during fluid infiltration in porous media, revealing conditions that lead to smooth or rough fronts.

Contribution

It introduces a numerical approach to quantify how relative rates of advection, diffusion, and reaction govern reaction front morphology in porous systems.

Findings

Rough fronts form when advection is dominant and reactions are slow.

No rough fronts occur without advection or with very fast reactions.

Reaction rates significantly influence pattern formation, partly independent of transport mechanisms.

Abstract

Replacement reactions during fluid infiltration into porous media, rocks and buildings are known to have important implications for reservoir development, ore formation as well as weathering. Natural observations and experiments have shown that in such systems the shape of reaction fronts can vary significantly ranging from smooth, rough to highly irregular. It remains unclear what process-related knowledge can be derived from these reaction front patterns. In this contribution we show a numerical approach to test the effect of relative rates of advection, diffusion and reaction on the development of reaction fronts patterns in granular aggregates with permeable grain boundaries. The numerical model takes (i) fluid infiltration along permeable grain boundaries, (ii) reactions and (iii) elemental diffusion into account. We monitor the change in element concentration within the fluid, while reactions occur at a pre-defined rate as a function of the local fluid concentration. In non-dimensional phase space using P{\'e}clet and Damk{\"o}hler numbers, results show that there are no rough fronts without advection (P{\'e}clet<70) nor if the reaction is too fast (Damk{\"o}hler>10-3). As advection becomes more dominant and reaction slower, roughness develops across several grains with a full microstructure mimicking replacement in the most extreme cases. The reaction front patterns show an increase in roughness with increasing P{\'e}clet number from P{\'e}clet 10 to 100 but then a decrease in roughness towards higher P{\'e}clet numbers controlled by the Damk{\"o}hler number. Our results indicate that reaction rates are crucial for pattern formation and that the shape of reaction fronts is only partly due to the underlying transport mechanism.