Linking Mixing Interface Deformation to Concentration Gradients in Porous Media

We study the pore-scale transport of a conservative scalar forming an advancing mixing front, which can be re-interpreted to predict instantaneous mixing-limited bimolecular reactions. We investigate this using a set of two-dimensional, high-resolution numerical simulations within a poly-disperse granular porous medium, covering a wide range of Peclet numbers. The aim is to show and exploit the direct link between pore-scale concentration gradients and mixing interface (midpoint concentration isocontour). We believe that such a perspective provides a complementary new lens for better understanding mixing and spreading in porous media. We develop and validate a theoretical model that quantifies the temporal elongation of the mixing interface and the upscaled reaction kinetics in mixing-limited systems accounting for pore-scale concentration fluctuations. Contrary to the classical belief that, given sufficient time, pore-scale fluctuations would eventually be washed out, we show that for $Pe>1$ advection generates pore-scale concentration fluctuations more rapidly than they can be fully dissipated. For such P'eclet numbers, once incomplete mixing is established, it will persist indefinitely.