Transition from diffusion to advection controlled contaminant adsorption in saturated chemically heterogeneous porous subsurfaces

TL;DR

This study investigates how the dominant mechanism of contaminant adsorption in porous media shifts from diffusion to advection as the fraction of adsorbing particles varies, using pore-scale simulations to quantify the transition.

Contribution

It introduces a quantitative framework for understanding the transition from diffusion- to advection-controlled adsorption in heterogeneous porous media based on pore-scale dynamics.

Findings

Transition from diffusion to advection dominance depends on adsorber fraction

Critical interparticle distance scales with pore size and Peclét number

Plume deformation and growth rate are proportional to pore velocity

Abstract

We show the impact that scalar structures deformation and mixing has on the fate of plumes of waterborne contaminant transported through a chemically heterogeneous, partially adsorbing porous medium. Via pore-scale simulations, we follow the dynamic of a passive scalar injected in a packed bed consisting of a mixture of chemically inert and adsorbing spherical particles. By varying the fraction of the adsorbers $\xi$, randomly distributed in the porous volume, we find that the waterborne solute forms concentration plumes emerging between pairs of adsorbing particles. This deformation is a consequence of the different mechanisms of transport characterising the transport of molecules in the proximal and remote pores relative to the adsorbers, diffusion and advection, respectively. The resulting isoscalar surface embedding the plumes grows at a rate proportional to the average pore-scale velocity $U$ and inversely proportional to the adsorbers' interparticle dimensionless distance, i.e. $\gamma \propto U/\ell_\xi$. We provide a quantification of the characteristic diffusive time scale of the plume $t_\eta \propto \ell_\xi^2/D_m$, which dissipates the concentration differences in the vicinity of the adsorbers. Thus, by quantifying the relative importance of the advection-sustained stretching rate $\gamma$ and plume mixing rate $1/t_\eta$ for different values of fraction of adsorbers $\xi$, we establish a transition from diffusion- to advection- dominated macroscopic adsorption, whose time evolutions scale as $\propto \sqrt{t}$ and $\propto t$, respectively. Such a transition is determined by the amount of adsorbers within the medium, with diffusion and advection dominating at high and low fractions, respectively. Our numerical analysis provides $\ell_\xi/d \approx 4/\ln(2) \textit{Pe}^{-1}$ as the critical distance between adsorbers that sets the transition, being $d$ the pore size.