Phase Saturation Control on Mixing-Driven Reactions in 3D Porous Media

Enhanced Reaction Kinetics in Stationary Two-Phase Flow through Porous Media

Understanding the interaction between multiphase flow and reactive transport in porous media is critical for many environmental and industrial applications. When a nonwetting immiscible phase is present within the pore space, it can remain immobile, which we call unsaturated flow, or move, resulting in multiphase flow. Previous studies under unsaturated flow conditions have shown that, for a given flow rate, the product of a mixing-driven reaction increases as wetting phase saturation decreases. Conversely, the opposite effect is observed for a given Péclet number (i.e., the flow rate is adapted depending on the wetting phase saturation). However, the impact of multiphase flow dynamics on mixing-driven reactions is poorly understood due to experimental and numerical challenges. To assess the impact of multiphase flow conditions on product formation, we use an optimized chemiluminescence reaction and an experimental setup that allows the separate injection of reactants along with a stationary two-phase flow. In our experiments, the mass of the reaction product under stationary two-phase flow conditions increases faster than Fickian beyond the diffusive time. The global kinetics initially increase before experiencing a monotonic decrease with significant fluctuations caused by the displacement of the nonwetting phase. For a given flow rate of the wetting phase, product formation depends on the flow rate of the nonwetting immiscible phase.

Experimental and modeling insights into mixing-limited reactive transport in heterogeneous porous media: Role of stagnant zones

The understanding of mixing-controlled reactive dynamics in heterogeneous porous media remains limited, presenting significant challenges for modeling subsurface contaminant transport processes and for designing cost-effective environmental remedial efforts. The complexity of accurately observing, measuring, and modeling mixing-limited reactive transport has led to inadequate exploration of these critical processes. This study investigates the mixing and reaction kinetics affected by stagnant zones, which are commonly found in alluvial aquifers-aquitards and fracture-matrix systems. By conducting experiments involving conservative and bimolecular reactive transport through porous media within translucent chambers filled with two sizes of glass beads and under varying flow rates, we explored the effects of grain size and hydrodynamic conditions. Using a high-resolution camera, we monitored the concentration changes of conservative and reactive tracers, with subsequent interpretation through three-dimensional numerical simulations. The outcomes revealed the emergence of distinct mixing interfaces within both mobile and stagnant zones, culminating in a bi-peaked plume formation. Notably, the mixing and reaction times in media containing stagnant zones were found to be approximately 10 times longer than in homogeneous media. These findings, through experimental and modeling efforts, advance our understanding of mixing-limited reactive transport phenomena within heterogeneous media, underscoring the significant role of stagnant zones-a topic previously underexplored.

Internal Biofilm Heterogeneities Enhance Solute Mixing and Chemical Reactions in Porous Media

Bacterial biofilms can form in porous media that are of interest in industrial applications ranging from medical implants to biofilters as well as in environmental applications such as in situ groundwater remediation, where they can be critical locations for biogeochemical reactions. The presence of biofilms modifies porous media topology and hydrodynamics by clogging pores and consequently solutes transport and reactions kinetics. The interplay between highly heterogeneous flow fields found in porous media and microbial behavior, including biofilm growth, results in a spatially heterogeneous biofilm distribution in the porous media as well as internal heterogeneity across the thickness of the biofilm. Our study leverages highly resolved three-dimensional X-ray computed microtomography images of bacterial biofilms in a tubular reactor to numerically compute pore-scale fluid flow and solute transport by considering multiple equivalent stochastically generated internal permeability fields for the biofilm. We show that the internal heterogeneous permeability mainly impacts intermediate velocities when compared with homogeneous biofilm permeability. While the equivalent internal permeability fields of the biofilm do not impact fluid-fluid mixing, they significantly control a fast reaction. For biologically driven reactions such as nutrient or contaminant uptake by the biofilm, its internal permeability field controls the efficiency of the process. This study highlights the importance of considering the internal heterogeneity of biofilms to better predict reactivity in industrial and environmental bioclogged porous systems.

Morphogenesis of Biofilms in Porous Media and Control on Hydrodynamics

The functioning of natural and engineered porous media, like soils and filters, depends in many cases on the interplay between biochemical processes and hydrodynamics. In such complex environments, microorganisms often form surface-attached communities known as biofilms. Biofilms can take the shape of clusters, which alter the distribution of fluid flow velocities within the porous medium, subsequently influencing biofilm growth. Despite numerous experimental and numerical efforts, the control of the biofilm clustering process and the resulting heterogeneity in biofilm permeability is not well understood, limiting our predictive abilities for biofilm-porous medium systems. Here, we use a quasi-2D experimental model of a porous medium to characterize biofilm growth dynamics for different pore sizes and flow rates. We present a method to obtain the time-resolved biofilm permeability field from experimental images and use the obtained permeability field to compute the flow field through a numerical model. We observe a biofilm cluster size distribution characterized by a spectrum slope evolving in time between -2 and -1, a fundamental measure that can be used to create spatio-temporal distributions of biofilm clusters for upscaled models. We find a previously undescribed biofilm permeability distribution, which can be used to stochastically generate permeability fields within biofilms. An increase in velocity variance for a decrease in physical heterogeneity shows that the bioclogged porous medium behaves differently than expected from studies on heterogeneity in abiotic porous media.

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.

A Lab on a Chip Experiment for Upscaling Diffusivity of Evolving Porous Media

Reactive transport modelling is a powerful tool to assess subsurface evolution in various energy-related applications. Upscaling, i.e., accounting for pore scale heterogeneities into larger scale analyses, remains one of the biggest challenges of reactive transport modelling. Pore scale simulations capturing the evolutions of the porous media over a wide range of Peclet and Damköhler number in combination with machine learning are foreseen as an efficient methodology for upscaling. However, the accuracy of these pore scale models needs to be tested against experiments. In this work, we developed a lab on a chip experiment with a novel micromodel design combined with operando confocal Raman spectroscopy, to monitor the evolution of porous media undergoing coupled mineral dissolution and precipitation processes due to diffusive reactive fluxes. The 3D-imaging of the porous media combined with pore scale modelling enabled the derivation of upscaled transport parameters. The chemical reaction tested involved the replacement of celestine by strontianite, whereby a net porosity increase is expected because of the smaller molar volume of strontianite. However, under our experimental conditions, the accessible porosity and consequently diffusivity decreased. We propose a transferability of the concepts behind the Verma and Pruess relationship to be applied to also describe changes of diffusivity for evolving porous media. Our results highlight the importance of calibrating pore scale models with quantitative experiments prior to simulations over a wide range of Peclet and Damköhler numbers of which results can be further used for the derivation of upscaled parameters.