Turbulent flows have been used for millennia to mix solutes; a familiar example is stirring cream into coffee. However, many energy, environmental, and industrial processes rely on the mixing of solutes in porous media where confinement suppresses inertial turbulence. As a result, mixing is drastically hindered, requiring fluid to permeate long distances for appreciable mixing and introducing additional steps to drive mixing that can be expensive and environmentally harmful. Here, we demonstrate that this limitation can be overcome just by adding dilute amounts of flexible polymers to the fluid. Flow-driven stretching of the polymers generates an elastic instability, driving turbulent-like chaotic flow fluctuations, despite the pore-scale confinement that prohibits typical inertial turbulence. Using in situ imaging, we show that these fluctuations stretch and fold the fluid within the pores along thin layers (“lamellae”) characterized by sharp solute concentration gradients, driving mixing by diffusion in the pores. This process results in a reduction in the required mixing length, a increase in solute transverse dispersivity, and can be harnessed to increase the rate at which chemical compounds react by —enhancements that we rationalize using turbulence-inspired modeling of the underlying transport processes. Our work thereby establishes a simple, robust, versatile, and predictive way to mix solutes in porous media, with potential applications ranging from large-scale chemical production to environmental remediation. 
                        more » 
                        « less   
                    This content will become publicly available on December 10, 2025
                            
                            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. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2046015
- PAR ID:
- 10587196
- Publisher / Repository:
- United States National Academy of Sciences (NAS)
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 121
- Issue:
- 50
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract Large discrepancies between well-mixed reaction rates and effective reactions rates estimated under fluid flow conditions have been a major issue for predicting reactive transport in porous media systems. In this study, we introduce a framework that accurately predicts effective reaction rates directly from pore structural features by combining 3D pore-scale numerical simulations with machine learning (ML). We first perform pore-scale reactive transport simulations with fluid–solid reactions in hundreds of porous media and calculate effective reaction rates from pore-scale concentration fields. We then train a Random Forests model with 11 pore structural features and effective reaction rates to quantify the importance of structural features in determining effective reaction rates. Based on the importance information, we train artificial neural networks with varying number of features and demonstrate that effective reaction rates can be accurately predicted with only three pore structural features, which are specific surface, pore sphericity, and coordination number. Finally, global sensitivity analyses using the ML model elucidates how the three structural features affect effective reaction rates. The proposed framework enables accurate predictions of effective reaction rates directly from a few measurable pore structural features, and the framework is readily applicable to a wide range of applications involving porous media flows.more » « less
- 
            Having a basic understanding of non-Newtonian fluid flow through porous media, which usually consist of series of expansions and contractions, is of importance for enhanced oil recovery, groundwater remediation, microfluidic particle manipulation, etc. The flow in contraction and/or expansion microchannel is unbounded in the primary direction and has been widely studied before. In contrast, there has been very little work on the understanding of such flow in an expansion–contraction microchannel with a confined cavity. We investigate the flow of five types of non-Newtonian fluids with distinct rheological properties and water through a planar single-cavity microchannel. All fluids are tested in a similarly wide range of flow rates, from which the observed flow regimes and vortex development are summarized in the same dimensionless parameter spaces for a unified understanding of the effects of fluid inertia, shear thinning, and elasticity as well as confinement. Our results indicate that fluid inertia is responsible for developing vortices in the expansion flow, which is trivially affected by the confinement. Fluid shear thinning causes flow separations on the contraction walls, and the interplay between the effects of shear thinning and inertia is dictated by the confinement. Fluid elasticity introduces instability and asymmetry to the contraction flow of polymers with long chains while suppressing the fluid inertia-induced expansion flow vortices. However, the formation and fluctuation of such elasto-inertial fluid vortices exhibit strong digressions from the unconfined flow pattern in a contraction–expansion microchannel of similar dimensions.more » « less
- 
            A fundamental understanding of the flow of polymer solutions through the pore spaces of porous media is relevant and significant to enhanced oil recovery and groundwater remediation. We present in this work an experimental study of the fluid rheological effects on non-Newtonian flows in a simple laboratory model of the real-world pores—a rectangular sudden contraction–expansion microchannel. We test four different polymer solutions with varying rheological properties, including xanthan gum (XG), polyvinylpyrrolidone (PVP), polyethylene oxide (PEO), and polyacrylamide (PAA). We compare their flows against that of pure water at the Reynolds ( R e ) and Weissenburg ( W i ) numbers that each span several orders of magnitude. We use particle streakline imaging to visualize the flow at the contraction–expansion region for a comprehensive investigation of both the sole and the combined effects of fluid shear thinning, elasticity and inertia. The observed flow regimes and vortex development in each of the tested fluids are summarized in the dimensionless W i − R e and χ L − R e parameter spaces, respectively, where χ L is the normalized vortex length. We find that fluid inertia draws symmetric vortices downstream at the expansion part of the microchannel. Fluid shear thinning causes symmetric vortices upstream at the contraction part. The effect of fluid elasticity is, however, complicated to analyze because of perhaps the strong impact of polymer chemistry such as rigidity and length. Interestingly, we find that the downstream vortices in the flow of Newtonian water, shear-thinning XG and elastic PVP solutions collapse into one curve in the χ L − R e space.more » « less
- 
            Blunt, MJ (Ed.)Multiphase fluid flow in porous media is important to a wide variety of processes of fundamental scientific and practical importance. Developing a model for the pore space of porous media represents the first step for simulating such flows. With rapid increase in the computation power and advances in instrumentation and imaging processes, it has become feasible to carry out simulation of multiphase flow in two- and three-dimensional images of porous media, hence dispensing with development of models of pore space that are based on approximating their morphology. Image-based simulations are, however, very time consuming. We describe an approach for speeding-up image-based simulation of multiphase flow in porous media based on curvelet transformations, which are specifically designed for processing of images that contain complex curved surfaces. Most porous media contain correlations in their morphology and, therefore, their images carry redundant information that, in the curvelet transform space, can be removed efficiently and accurately in order to obtain a coarser image with which the computations are far less intensive. We utilize the methodology to simulate two-phase flow of oil and water in two-dimensional digital images of sandstone and carbonate samples, and demonstrate that while the results with the curvelet-processed images are as accurate as those with the original ones, the computations are speeded up by a factor of 110–150. Thus, the methodology opens the way toward achieving the ultimate goal of simulation of multiphase flow in porous media, namely, making image-based computations a standard practice.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
