Experimental Evidence of Helical Flow in Porous Media

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

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.

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

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.

Linking mixing and flow topology in porous media: An experimental proof

Transport processes in porous media are controlled by the characteristics of the flow field which are determined by the porous material properties and the boundary conditions of the system. This work provides experimental evidence of the relation between mixing and flow field topology in porous media at the continuum scale. The setup consists of a homogeneously packed quasi-two-dimensional flow-through chamber in which transient flow conditions, dynamically controlled by two external reservoirs, impact the transport of a dissolved tracer. The experiments were performed at two different flow velocities, corresponding to Péclet numbers of 191 and 565, respectively. The model-based interpretation of the experimental results shows that high values of the effective Okubo-Weiss parameter, driven by the changes of the boundary conditions, lead to high rates of increase of the Shannon entropy of the tracer distribution and, thus, to enhanced mixing. The comparison between a hydrodynamic dispersion model and an equivalent pore diffusion model demonstrates that despite the spatial and temporal variability in the hydrodynamic dispersion coefficients, the Shannon entropy remains almost unchanged because it is controlled by the Okubo-Weiss parameter. Overall, our work demonstrates that under highly transient boundary conditions, mixing dynamics in homogeneous porous media can also display complex patterns and is controlled by the flow topology.

Impact of solute charge and diffusion coefficient on electromigration and mixing in porous media

The application of electrokinetic techniques in porous media has great potential to enhance mass transfer rates and, thus, to mobilize contaminants and effectively deliver reactants and amendments. However, the transport mechanisms induced by the application of an external electric field are complex and entail the coupling of physical, chemical and electrostatic processes. In this study we focus on electromigration and we provide experimental evidence of the impact of compound-specific properties, such as the aqueous diffusivity and the valence of charged species, on the macroscopic electrokinetic transport. We performed a series of multidimensional experiments considering the displacement of three different tracer plumes (i.e., permanganate, allura red and new coccine) in different background electrolyte solutions. The outcomes of the experiments clearly show that both the compound-specific diffusivity and the charge of the injected and resident ions impact the transport of the selected color tracer plumes, whose evolution was monitored with image analysis. The investigated experimental scenarios led to distinct plume behavior characterized by different mass distribution, average displacement velocities, longitudinal and lateral plume spreading, shape of the invading and receding fronts, as well as dilution of the injected solutes. A numerical simulator, based on the Nernst-Planck-Poisson equations and on aqueous speciation reactions in the pore water, allowed us to quantitatively interpret the experimental results, to capture the observed patterns of plume evolution, and to illuminate the coupling between the governing physico-chemical mechanisms and the controlling role of small scale compound-specific and electrostatic properties. Finally, the model was also extended to a typical configuration of in situ electrokinetic remediation of contaminated groundwater to show the impact of such mechanisms at larger scale.

Impact of diffuse layer processes on contaminant forward and back diffusion in heterogeneous sandy-clayey domains

Low-permeability aquitards can significantly affect the transport, distribution, and persistence of contaminant plumes in subsurface systems. Although such low-permeability materials are often charged, the key role of charge-induced electrostatic processes during contaminant transport has not been extensively studied. This work presents a detailed investigation exploring the coupled effects of heterogeneous distribution of physical, chemical and electrostatic properties on reactive contaminant transport in field-scale groundwater systems including spatially distributed clay zones. We performed an extensive series of numerical experiments in three distinct heterogeneous sandy-clayey domains with different levels of complexity. The flow and reactive transport simulations were performed by explicitly resolving the complex velocity fields, the small-scale electrostatic processes, the compound-specific diffusive/dispersive fluxes and the chemical processes utilizing a multi-continua based reactive transport code (MMIT-Clay). In each particular domain, numerical experiments were performed focusing on both the forward and back diffusion through the sandy-clayey interfaces. The results illuminate the control of microscopic electrostatic mechanisms on macroscopic mass transfer. Coulombic interactions in the clay's diffuse layer can significantly accelerate or retard a particular species depending on its charge. Furthermore, the chemical heterogeneity plays a major role in mass storage and release during reactive transport. Neglecting such processes can lead to substantial over- or underestimation of the overall transport behavior, which underlines the need for integrated physical, chemical and electrostatic approaches to accurately describe mass transfer processes in systems including low-permeability inclusions.

Electroosmotic flow: From microfluidics to nanofluidics

Electroosmotic flow (EOF), a consequence of an imposed electric field onto an electrolyte solution in the tangential direction of a charged surface, has emerged as an important phenomenon in electrokinetic transport at the micro/nanoscale. Because of their ability to efficiently pump liquids in miniaturized systems without incorporating any mechanical parts, electroosmotic methods for fluid pumping have been adopted in versatile applications—from biotechnology to environmental science. To understand the electrokinetic pumping mechanism, it is crucial to identify the role of an ionically polarized layer, the so‐called electrical double layer (EDL), which forms in the vicinity of a charged solid–liquid interface, as well as the characteristic length scale of the conducting media. Therefore, in this tutorial review, we summarize the development of electrical double layer models from a historical point of view to elucidate the interplay and configuration of water molecules and ions in the vicinity of a solid–liquid interface. Moreover, we discuss the physicochemical phenomena owing to the interaction of electrical double layer when the characteristic length of the conducting media is decreased from the microscale to the nanoscale. Finally, we highlight the pioneering studies and the most recent works on electro osmotic flow devoted to both theoretical and experimental aspects.

Plume deformation, mixing, and reaction kinetics: An analysis of interacting helical flows in three-dimensional porous media

Heterogeneity and macroscopic anisotropy of porous media play an important role for dilution and reaction enhancement of conservative and reactive plumes. In this study, we perform numerical simulations to investigate steady-state flow and transport in three-dimensional heterogeneous porous media. We consider two macroscopic anisotropic inclusions resulting in helical flows with twisting streamlines in a three-dimensional flow-through domain. The inclusions are obtained by alternating two layers of angled slices of coarse and fine porous media with different hydraulic conductivity. We investigate flow and transport scenarios considering different geometry and relative position of the two anisotropic inclusions yielding helical flow fields with different extent of interaction. We use metrics of stretching and folding to characterize the flow field and entropy-based metrics for the analysis of the conservative and reactive transport problems. The outcomes show that the two helices result in different patterns of twisting streamlines, which cause distinct deformation of the plumes. However, mixing and reaction enhancement could not be directly related to the extent of the flow field deformation: Configurations with strong deformation can result in only moderate mixing enhancement, whereas configurations with limited deformation of the flow field can lead to significant mixing of the solute plume. Finally, we explore the impact of different degradation rates on reactive transport and the role of reaction kinetics on the entropy balance for a reactant undergoing transport and mixing-controlled degradation in the twisting flow fields. The results show that strong mixing enhancement due to helical flow increases the importance of the reaction kinetics that becomes the rate-limiting process for solute reactive transport.

Model-based interpretation of methane oxidation and respiration processes in landfill biocovers: 3-D simulation of laboratory and pilot experiments

Landfill biocovers are an efficient strategy for the mitigation of greenhouse gas emissions from landfills. A complex interplay between key physical and reactive processes occurs in biocovers and affects the transport of gas components. Therefore, numerical models can greatly help the understanding of these systems, their design and optimal operation. In this study, we developed a 3-D multicomponent modeling approach to quantitatively interpret experimental datasets measured in the laboratory and in pilot-scale landfill biocovers. The proposed model is able to reproduce the observed spatial and temporal dynamics of CH₄, O₂ and CO₂ migration in biocovers under different operating conditions and demonstrates the importance of dimensionality in understanding the propagation of gas flow and migration of gas components in such porous media. The model allowed us to capture the coupled transport behavior of gas components, to evaluate the exchange of gas fluxes at the interface between the biocover surface and free air flow, and to investigate the effects of different gas injection patterns on the distribution of gas components within biocovers. The model also helps elucidating the dynamics and competition between methane oxidation and respiration processes observed in the different experimental setups. The simulation outcomes reveal that increasing availability of methane (i.e., higher injection flow rates or higher fractions of CH₄ in the injected gas composition) results in progressive dominance of methane oxidation in the biocovers and moderates the impact of respiration.

Three-Dimensional Vortex-Induced Reaction Hot Spots at Flow Intersections

We show the emergence of reaction hot spots induced by three-dimensional (3D) vortices with a simple A+B→C reaction. We conduct microfluidics experiments to visualize the spatial map of the reaction rate with a chemiluminescence reaction and cross validate the results with direct numerical simulations. 3D vortices form at spiral-saddle-type stagnation points, and the 3D vortex flow topology is essential for initiating reaction hot spots. The effect of vortices on mixing and reaction becomes more vigorous for rough-walled channels, and our findings are valid over wide ranges of channel dimensions and Damköhler numbers.

Charge interactions, reaction kinetics and dimensionality effects on electrokinetic remediation: A model-based analysis

The potential of electrokinetic remediation technologies (EKR) for the removal of different contaminants from subsurface porous media has been increasingly recognized. Despite electrokinetic applications have shown promising results, quantitative understanding of such systems is still challenging due to the complex interplay between physical transport processes, electrostatic interactions, and geochemical reactions. In this study, we perform a model-based analysis of electrokinetic transport in saturated porous media. We investigate the effects of: (i) Coulombic interactions between ions in the system mobilized by electromigration, (ii) reaction kinetics on the overall removal efficiency of a non-charged organic contaminant, and (iii) dimensionality and different electrode configurations. The results show that such effects play a major role on the performance of electrokinetic systems. The simulations illuminate the importance of microscopic processes, such as electrostatic interactions and ion-specific diffusivities, and their non-intuitive macroscopic impact on the delivery of charged amendments and on the efficiency of contaminant removal. The insights of this study are valuable to improve and optimize the design and the operational strategies of electrokinetic remediation systems.

Recent Advances in Experimental Studies of Steady-State Dilution and Reactive Mixing in Saturated Porous Media

Transverse dispersive mixing plays an important role in controlling natural attenuation of contaminant plumes and the performance of engineered remediation strategies. The extent of transverse mixing can be significantly affected by porous media heterogeneity and anisotropy. For instance, flow focusing in the high-permeability inclusions leads to an enhancement of dilution and reactive mixing in steady-state solute transport. Numerous modeling studies have been performed to understand the mechanism of conservative and reactive transport in homogeneous and complex heterogeneous porous media. However, experimental investigations are necessary to show an intuitive phenomenon and to validate the modeling results. This paper briefly reviews recent laboratory experimental studies on dilution and reactive mixing of steady-state transport in saturated homogeneous and heterogeneous porous media. In this context, setups and measuring techniques are described in pore-scale and Darcy-scale experiments. Parameters quantifying dilution and reactive mixing in the experiments are also introduced. Finally, we discuss the further experimental works necessary to deepen our understanding of dilution and reactive mixing in natural aquifers.

Fluid Deformation in Random Steady Three Dimensional Flow

The deformation of elementary fluid volumes by velocity gradients is a key process for scalar mixing, chemical reactions and biological processes in flows. Whilst fluid deformation in unsteady, turbulent flow has gained much attention over the past half century, deformation in steady random flows with complex structure - such as flow through heterogeneous porous media - has received significantly less attention. In contrast to turbulent flow, the steady nature of these flows constrains fluid deformation to be anisotropic with respect to the fluid velocity, with significant implications for e.g. longitudinal and transverse mixing and dispersion. In this study we derive an ab initio coupled continuous time random walk (CTRW) model of fluid deformation in random steady three-dimensional flow that is based upon a streamline coordinate transform which renders the velocity gradient and fluid deformation tensors upper-triangular. We apply this coupled CTRW model to several model flows and find these exhibit a remarkably simple deformation structure in the streamline coordinate frame, facilitating solution of the stochastic deformation tensor components. These results show that the evolution of longitudinal and transverse fluid deformation for chaotic flows is governed by both the Lyapunov exponent and power-law exponent of the velocity PDF at small velocities, whereas algebraic deformation in non-chaotic flows arises from the intermittency of shear events following similar dynamics as that for steady two-dimensional flow.

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.

Coupled continuous-time random walks for fluid stretching in two-dimensional heterogeneous media

We study the relation between flow structure and fluid deformation in steady flows through two-dimensional heterogeneous media, which are characterized by a broad spectrum of stretching behaviors, ranging from sub- to superlinear. We analyze these behaviors from first principles, which uncovers intermittent shear events to be at the origin of subexponential stretching. We derive explicit expressions for Lagrangian deformation and demonstrate that stretching obeys a coupled continuous-time random walk, which for broad distributions of flow velocities becomes a Lévy walk. The derived model provides a direct link between the flow and deformation statistics, and a natural way to quantify the impact of intermittent shear events on the stretching behavior.

Experimental investigation of transverse mixing in porous media under helical flow conditions

Plume dilution and transverse mixing can be considerably enhanced by helical flow occurring in three-dimensional heterogeneous anisotropic porous media. In this study, we perform tracer experiments in a fully three-dimensional flow-through chamber to investigate the effects of helical flow on plume spiraling and deformation, as well as on its dilution. Porous media were packed in angled stripes of materials with different grain sizes to create blocks with macroscopically anisotropic hydraulic conductivity, which caused helical flows. Steady-state transport experiments were carried out by continuously injecting dye tracers at different inlet ports. High-resolution measurements of concentration and flow rates were performed at 49 outlet ports. These measurements allowed quantifying the spreading and dilution of the solute plumes at the outlet cross section. Direct evidence of plume spiraling and visual proof of helical flow was obtained by freezing and slicing the porous media at different cross sections and observing the dye-tracer distribution. We simulated flow and transport to interpret our experimental observations and investigate the effects of helical flow on mixing-controlled reactive transport. The simulation results were evaluated using metrics of reactive mixing such as the critical dilution index and the length of continuously injected steady-state plumes. The results show considerable reaction enhancement, quantified by the remarkable decrease of reactive plume lengths (up to four times) in helical flows compared to analogous scenarios in uniform flows.