The bimolecular reactive transport in heterogeneous porous media: Sub-diffusion in interpretation of laboratory experiment.

This present work consists of investigating the effects of particle size heterogeneity and flow rates on transport-reaction kinetics of CuSO and NaEDTA in porous media, via the combination of a bimolecular reaction experiment and model simulations. In the early stages of transport, a peak is observed in the concentration breakthrough curve of the reactant CuSO, related to the delayed mixing and reaction of the reactants. The numerical results show that an increase in flow rate promotes the mixing processes between the reactants, resulting in a larger peak concentration and a slighter tail of breakthrough curves, while an increase in medium heterogeneity leads to a more significant heavy tail. The apparent anomalous diffusion and heavy-tailing behavior can be effectively quantified by a novel truncated fractional derivative bimolecular reaction model. The truncated fractional-order model, taking into account the incomplete mixing, offers a satisfactory reproduction of the experimental data.