Is chaotic advection inherent to porous media flow?

Chaotic Transport of Solutes in Unsaturated Porous Media

Unsaturated porous media, characterized by the combined presence of several immiscible fluid phases in the pore space, are highly relevant systems in nature, because they control the fate of contaminants and the availability of nutrients in the subsoil. However, a full understanding of the mechanisms controlling solute mixing in such systems is still missing. In particular, the role of saturation in the development of chaotic solute mixing has remained unexplored. Using three-dimensional numerical simulations of flow and transport at the pore scale, built upon X-ray tomograms of a porous medium at different degrees of liquid (wetting)-phase saturation, we show the occurrence of chaotic dynamics in both the deformation of the solute plume, as characterized by computed chaos metrics (Lyapunov exponents), and the mixing of the injected solute. Our results show an enhancement of these chaotic dynamics at lower saturation and their occurrence even under diffusion-relevant conditions over the medium's length, also being strengthened by larger flow velocities. These findings highlight the dominant role of the pore-scale spatial heterogeneity of the system, enhanced by the presence of an immiscible phase (e.g., air), on the mixing efficiency. This represents a stepping stone for the assessment of mixing and reactions in unsaturated porous media.

Evolution of plume geometry, dilution and reactive mixing in porous media under highly transient flow fields at the surface water-groundwater interface

Highly transient boundary conditions affect mixing of dissolved solutes in groundwater. An example of these transient boundary conditions occurs at the surface water-groundwater interface, where the water level in rivers can change rapidly due to the operation of hydropower plants, leading to a regime known as hydropeaking. Inspired by this phenomenon, this work studies at laboratory scale the effects of fluctuating surface water bodies on solute transport in aquifers. We performed flow-through experiments at two different flow velocities and under steady and transient flow conditions where a conservative tracer was injected in the system and its concentration measured with optical imaging methods. The experimental results were quantitatively interpreted with numerical simulations implementing a non-linear velocity-dependent dispersive transport model. We estimated plume dilution by computing the dilution index and its evolution as well as two key geometrical metrics of the transient plumes: the perimeter and the area. We further investigated reactive mixing and mixing enhancement considering mixing-controlled bimolecular reactions using different critical mixing ratios. In general, highly transient boundary conditions lead to a larger area, perimeter and plume dilution and the results show greater relative enhancement for the scenarios with low groundwater flow velocity. A linear relationship was observed between the evolution of the area and the dilution index of the plumes for the transient flow scenarios investigated. Considering reactive transport and mixing-limited reactions at the surface water-groundwater interface, we identified different dilution and reaction dominated regimes, characterized, respectively, by increasing and decreasing plume entropies at different mixing ratios of the reactants. Furthermore, reactive mixing was enhanced by transient flows leading to a faster degradation of contaminant plumes compared to corresponding steady flow conditions.

Fluid-Solid Reaction in Porous Media as a Chaotic Restart Process

Chemical and biological reactions at fluid-solid interfaces are central to a broad range of porous material applications and research. Pore-scale solute transport limitations can reduce reaction rates, with marked consequences for a wide spectrum of natural and engineered processes. Recent advances show that chaotic mixing occurs spontaneously in porous media, but its impact on surface reactions is unknown. We show that pore-scale chaotic mixing significantly increases reaction efficiency compared to nonchaotic flows. We find that reaction rates are well described in terms of diffusive first-passage times of reactants to the solid interface subjected to a stochastic restart process resulting from Lagrangian chaos. Under chaotic mixing, the shear layer at no-slip interfaces sets the restart rate and leads to a characteristic scaling of reaction efficiency with Péclet number, in excellent agreement with numerical simulations. Reaction rates are insensitive to the flow topology as long as flow is chaotic, suggesting the relevance of this process to a broad range of porous materials.

Mixing Controlled Adsorption at the Liquid-Solid Interfaces in Unsaturated Porous Media

The unsaturated zone, located between the soil surface and the phreatic level, plays an important role in defining the fate of any substance entering the subsoil. In addition to the processes of flow and transport taking place in the liquid phase, surface reactions such as adsorption to the solid phase may occur and increase the residence time of the substance entering the system. In this study, we aim to understand the pore-scale mechanisms that control adsorption in unsaturated systems. We combine 2D pore-scale experimental images with numerical simulations to analyze flow, transport, and adsorption under different liquid saturation degrees. We demonstrate the role of mixing on adsorption at the liquid-solid interfaces by analyzing the deformation in time of a pulse-injected surfactant. We also analyze the impact of the isotherm functional shape and the inclusion of the liquid-gas interfaces as adsorption sites on this surface reaction. The enhancement of mixing as saturation decreases is accompanied by a reduction in the amount of adsorbed mass, located mainly along preferential flow paths, where the solute is primarily transported. For the same isotherm, a nonlinear behavior of adsorption as a function of liquid saturation has been observed. This is explained by the nonlinear variation of the volume fraction of liquid behaving as preferential path or stagnation zone as liquid saturation decreases, despite the linear decrease in the surface area of solids accessible for adsorption.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11242-022-01747-x.

Mixing in Porous Media: Concepts and Approaches Across Scales

This review provides an overview of concepts and approaches for the quantification of passive, non-reactive solute mixing in steady uniform porous media flows across scales. Mixing in porous media is the result of the interaction of spatial velocity fluctuations and diffusion or local-scale dispersion, which may lead to the homogenization of an initially segregated system. Velocity fluctuations are induced by spatial medium heterogeneities at the pore, Darcy or regional scales. Thus, mixing in porous media is a multiscale process, which depends on the medium structure and flow conditions. In the first part of the review, we discuss the interrelated processes of stirring, dispersion and mixing, and review approaches to quantify them that apply across scales. This implies concepts of hydrodynamic dispersion, approaches to quantify mixing state and mixing dynamics in terms of concentration statistics, and approaches to quantify the mechanisms of mixing. We review the characterization of stirring in terms of fluid deformation and folding and its relation with hydrodynamic dispersion. The integration of these dynamics to quantify the mechanisms of mixing is discussed in terms of lamellar mixing models. In the second part of this review, we discuss these concepts and approaches for the characterization of mixing in Poiseuille flow, and in porous media flows at the pore, Darcy and regional scales. Due to the fundamental nature of the mechanisms and processes of mixing, the concepts and approaches discussed in this review underpin the quantitative analysis of mixing phenomena in porous media flow systems in general.

Fluidic bacterial diodes rectify magnetotactic cell motility in porous environments

Directed motility enables swimming microbes to navigate their environment for resources via chemo-, photo-, and magneto-taxis. However, directed motility competes with fluid flow in porous microbial habitats, affecting biofilm formation and disease transmission. Despite this broad importance, a microscopic understanding of how directed motility impacts the transport of microswimmers in flows through constricted pores remains unknown. Through microfluidic experiments, we show that individual magnetotactic bacteria directed upstream through pores display three distinct regimes, whereby cells swim upstream, become trapped within a pore, or are advected downstream. These transport regimes are reminiscent of the electrical conductivity of a diode and are accurately predicted by a comprehensive Langevin model. The diode-like behavior persists at the pore scale in geometries of higher dimension, where disorder impacts conductivity at the sample scale by extending the trapping regime over a broader range of flow speeds. This work has implications for our understanding of the survival strategies of magnetotactic bacteria in sediments and for developing their use in drug delivery applications in vascular networks.

Heat flows in rock cracks naturally optimize salt compositions for ribozymes

Catalytic nucleic acids, such as ribozymes, are central to a variety of origin-of-life scenarios. Typically, they require elevated magnesium concentrations for folding and activity, but their function can be inhibited by high concentrations of monovalent salts. Here we show that geologically plausible high-sodium, low-magnesium solutions derived from leaching basalt (rock and remelted glass) inhibit ribozyme catalysis, but that this activity can be rescued by selective magnesium up-concentration by heat flow across rock fissures. In contrast to up-concentration by dehydration or freezing, this system is so far from equilibrium that it can actively alter the Mg:Na salt ratio to an extent that enables key ribozyme activities, such as self-replication and RNA extension, in otherwise challenging solution conditions. The principle demonstrated here is applicable to a broad range of salt concentrations and compositions, and, as such, highly relevant to various origin-of-life scenarios.

Scalar Signatures of Chaotic Mixing in Porous Media

Steady laminar flows through porous media spontaneously generate Lagrangian chaos at pore scale, with qualitative implications for a range of transport, reactive, and biological processes. The characterization and understanding of mixing dynamics in these opaque environments is an outstanding challenge. We address this issue by developing a novel technique based upon high-resolution imaging of the scalar signature produced by push-pull flows through porous media samples. Owing to the rapid decorrelation of particle trajectories in chaotic flows, the scalar image measured outside the porous material is representative of in situ mixing dynamics. We present a theoretical framework for estimation of the Lyapunov exponent based on extension of Lagrangian stretching theories to correlated aggregation. This method provides a full characterization of chaotic mixing dynamics in a large class of porous materials.

Stretching and folding sustain microscale chemical gradients in porous media

Fluid flow in porous media drives the transport, mixing, and reaction of molecules, particles, and microorganisms across a wide spectrum of natural and industrial processes. Current macroscopic models that average pore-scale fluctuations into an effective dispersion coefficient have shown significant limitations in the prediction of many important chemical and biological processes. Yet, it is unclear how three-dimensional flow in porous structures govern the microscale chemical gradients controlling these processes. Here, we obtain high-resolution experimental images of microscale mixing patterns in three-dimensional porous media and uncover an unexpected and general mixing mechanism that strongly enhances concentration gradients at pore-scale. Our experiments reveal that systematic stretching and folding of fluid elements are produced in the pore space by grain contacts, through a mechanism that leads to efficient microscale chaotic mixing. These insights form the basis for a general kinematic model linking chaotic-mixing rates in the fluid phase to the generic structural properties of granular matter. The model successfully predicts the resulting enhancement of pore-scale chemical gradients, which appear to be orders of magnitude larger than predicted by dispersive approaches. These findings offer perspectives for predicting and controlling the vast diversity of reactive transport processes in natural and synthetic porous materials, beyond the current dispersion paradigm.

Space-Group Symmetries Generate Chaotic Fluid Advection in Crystalline Granular Media

The classical connection between symmetry breaking and the onset of chaos in dynamical systems harks back to the seminal theory of Noether [Transp. Theory Statist. Phys. 1, 186 (1918)10.1080/00411457108231446]. We study the Lagrangian kinematics of steady 3D Stokes flow through simple cubic and body-centered cubic (bcc) crystalline lattices of close-packed spheres, and uncover an important exception. While breaking of point-group symmetries is a necessary condition for chaotic mixing in both lattices, a further space-group (glide) symmetry of the bcc lattice generates a transition from globally regular to globally chaotic dynamics. This finding provides new insights into chaotic mixing in porous media and has significant implications for understanding the impact of symmetries upon generic dynamical systems.

Solvable continuous-time random walk model of the motion of tracer particles through porous media

We consider the continuous-time random walk (CTRW) model of tracer motion in porous medium flows based on the experimentally determined distributions of pore velocity and pore size reported by Holzner et al. [M. Holzner et al., Phys. Rev. E 92, 013015 (2015)PLEEE81539-375510.1103/PhysRevE.92.013015]. The particle's passing through one channel is modeled as one step of the walk. The step (channel) length is random and the walker's velocity at consecutive steps of the walk is conserved with finite probability, mimicking that at the turning point there could be no abrupt change of velocity. We provide the Laplace transform of the characteristic function of the walker's position and reductions for different cases of independence of the CTRW's step duration τ, length l, and velocity v. We solve our model with independent l and v. The model incorporates different forms of the tail of the probability density of small velocities that vary with the model parameter α. Depending on that parameter, all types of anomalous diffusion can hold, from super- to subdiffusion. In a finite interval of α, ballistic behavior with logarithmic corrections holds, which was observed in a previously introduced CTRW model with independent l and τ. Universality of tracer diffusion in the porous medium is considered.

Statistics of highly heterogeneous flow fields confined to three-dimensional random porous media

We present a strong relationship between the microstructural characteristics of, and the fluid velocity fields confined to, three-dimensional random porous materials. The relationship is revealed through simultaneously extracting correlation functions R_{uu}(r) of the spatial (Eulerian) velocity fields and microstructural two-point correlation functions S_{2}(r) of the random porous heterogeneous materials. This demonstrates that the effective physical transport properties depend on the characteristics of complex pore structure owing to the relationship between R_{uu}(r) and S_{2}(r) revealed in this study. Further, the mean excess plot was used to investigate the right tail of the streamwise velocity component that was found to obey light-tail distributions. Based on the mean excess plot, a generalized Pareto distribution can be used to approximate the positive streamwise velocity distribution.

Intermittent Lagrangian velocities and accelerations in three-dimensional porous medium flow

Intermittency of Lagrangian velocity and acceleration is a key to understanding transport in complex systems ranging from fluid turbulence to flow in porous media. High-resolution optical particle tracking in a three-dimensional (3D) porous medium provides detailed 3D information on Lagrangian velocities and accelerations. We find sharp transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity, superlinear evolution of particle dispersion, and double-peak behavior in the propagators. The velocity distribution is quantified in terms of pore geometry and flow connectivity, which forms the basis for a continuous-time random-walk model that sheds light on the observed Lagrangian flow and transport behaviors.

Anomalous transport and chaotic advection in homogeneous porous media

The topological complexity inherent to all porous media imparts persistent chaotic advection under steady flow conditions, which, in concert with the no-slip boundary condition, generates anomalous transport. We explore the impact of this mechanism upon longitudinal dispersion via a model random porous network and develop a continuous-time random walk that predicts both preasymptotic and asymptotic transport. In the absence of diffusion, the ergodicity of chaotic fluid orbits acts to suppress longitudinal dispersion from ballistic to superdiffusive transport, with asymptotic variance scaling as σ(L)(2)(t)∼t(2)/(ln t)(3). These results demonstrate that anomalous transport is inherent to homogeneous porous media and has significant implications for macrodispersion.

Predicting the enhancement of mixing-driven reactions in nonuniform flows using measures of flow topology

The ability for reactive constituents to mix is often the key limiting factor for the completion of reactions across a huge range of scales in a variety of media. In flowing systems, deformation and shear enhance mixing by bringing constituents into closer proximity, thus increasing reaction potential. Accurately quantifying this enhanced mixing is key to predicting reactions and typically is done by observing or simulating scalar transport. To eliminate this computationally expensive step, we use a Lagrangian stochastic framework to derive the enhancement to reaction potential by calculating the collocation probability of particle pairs in a heterogeneous flow field accounting for deformations. We relate the enhanced reaction potential to three well known flow topology metrics and demonstrate that it is best correlated to (and asymptotically linear with) one: the largest eigenvalue of the (right) Cauchy-Green tensor.