skip to main content

Explore Research Products in the NSF-PAR

Title: Transport of complex and active fluids in porous media

Complex and active fluids find broad applications in flows through porous materials. Nontrivial rheology can couple to porous microstructure leading to surprising flow patterns and associated transport properties in geophysical, biological, and industrial systems. Viscoelastic instabilities are highly sensitive to pore geometry and can give rise to chaotic velocity fluctuations. A number of recent studies have begun to untangle how the pore-scale geometry influences the sample-scale flow topology and the resulting dispersive transport properties of these complex systems. Beyond classical rheological properties, active colloids and swimming cells exhibit a range of unique properties, including reduced effective viscosity, collective motion, and random walks, that present novel challenges to understanding their mechanics and transport in porous media flows. This review article aims to provide a brief overview of essential, fundamental concepts followed by an in-depth summary of recent developments in this rapidly evolving field. The chosen topics are motivated by applications, and new opportunities for discovery are highlighted.

  more » « less
Award ID(s):
1554095 1701392
NSF-PAR ID:
10363239
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Society of Rheology
Date Published:
Journal Name:
Journal of Rheology
Volume:
66
Issue:
2
ISSN:
0148-6055
Page Range / eLocation ID:
p. 375-397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Computation)
    null (Ed.)
    In this paper, we consider an important problem for modeling complex coupled phenomena in porous media at multiple scales. In particular, we consider flow and transport in the void space between the pores when the pore space is altered by new solid obstructions formed by microbial growth or reactive transport, and we are mostly interested in pore-coating and pore-filling type obstructions, observed in applications to biofilm in porous media and hydrate crystal formation, respectively. We consider the impact of these obstructions on the macroscopic properties of the porous medium, such as porosity, permeability and tortuosity, for which we build an experimental probability distribution with reduced models, which involves three steps: (1) generation of independent realizations of obstructions, followed by, (2) flow and transport simulations at pore-scale, and (3) upscaling. For the first step, we consider three approaches: (1A) direct numerical simulations (DNS) of the PDE model of the actual physical process called BN which forms the obstructions, and two non-DNS methods, which we call (1B) CLPS and (1C) LP. LP is a lattice Ising-type model, and CLPS is a constrained version of an Allen–Cahn model for phase separation with a localization term. Both LP and CLPS are model approximations of BN, and they seek local minima of some nonconvex energy functional, which provide plausible realizations of the obstructed geometry and are tuned heuristically to deliver either pore-coating or pore-filling obstructions. Our methods work with rock-void geometries obtained by imaging, but bypass the need for imaging in real-time, are fairly inexpensive, and can be tailored to other applications. The reduced models LP and CLPS are less computationally expensive than DNS, and can be tuned to the desired fidelity of the probability distributions of upscaled quantities. 
    more » « less
  2. ; ; ( , Small)
    Abstract

    Polymer solutions are frequently used in enhanced oil recovery and groundwater remediation to improve the recovery of trapped nonaqueous fluids. However, applications are limited by an incomplete understanding of the flow in porous media. The tortuous pore structure imposes both shear and extension, which elongates polymers; moreover, the flow is often at large Weissenberg numbers, Wi, at which polymer elasticity in turn strongly alters the flow. This dynamic elongation can even produce flow instabilities with strong spatial and temporal fluctuations despite the low Reynolds number, Re. Unfortunately, macroscopic approaches are limited in their ability to characterize the pore‐scale flow. Thus, understanding how polymer conformations, flow dynamics, and pore geometry together determine these nontrivial flow patterns and impact macroscopic transport remains an outstanding challenge. This review describes how microfluidic tools can shed light on the physics underlying the flow of polymer solutions in porous media at high Wi and low Re. Specifically, microfluidic studies elucidate how steady and unsteady flow behavior depends on pore geometry and solution properties, and how polymer‐induced effects impact nonaqueous fluid recovery. This work thus provides new insights for polymer dynamics, non‐Newtonian fluid mechanics, and applications such as enhanced oil recovery and groundwater remediation.

     
    more » « less
  3. ; ; ( , Scientific Reports)
    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
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Proceedings of the National Academy of Sciences)
    Multiphase flows in porous media are important in many natural and industrial processes. Pore-scale models for multiphase flows have seen rapid development in recent years and are becoming increasingly useful as predictive tools in both academic and industrial applications. However, quantitative comparisons between different pore-scale models, and between these models and experimental data, are lacking. Here, we perform an objective comparison of a variety of state-of-the-art pore-scale models, including lattice Boltzmann, stochastic rotation dynamics, volume-of-fluid, level-set, phase-field, and pore-network models. As the basis for this comparison, we use a dataset from recent microfluidic experiments with precisely controlled pore geometry and wettability conditions, which offers an unprecedented benchmarking opportunity. We compare the results of the 14 participating teams both qualitatively and quantitatively using several standard metrics, such as fractal dimension, finger width, and displacement efficiency. We find that no single method excels across all conditions and that thin films and corner flow present substantial modeling and computational challenges. 
    more » « less
  5. ; ; ; ; ; ( , Journal of Fluid Mechanics)
    Direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers ( $St = 0.01, 0.1, 0.5, 1, 2$ ) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email