Search

Despite decades of research and development, the optimal efficiency of slurry-packed HPLC columns is still hindered by inherent long-range flow heterogeneity from the wall to the central bulk region of these columns. Here, we show an example of how this issue can be addressed through the straightforward addition of a semidilute amount (500~ppm) of a large, flexible, synthetic polymer (18~MDa partially hydrolyzed polyacrylamide, HPAM) to the mobile phase (1\% NaCl aqueous solution) during operation of a 4.6 mm $\times$ 300 mm column packed with 10~$\mu$m BEH$^{\mathrm{TM}}$ 125~\AA \ Particles. Addition of the polymer imparts elasticity to the mobile phase, causing the flow in the interparticle pore space to become unstable above a threshold flow rate. We verify the development of this elastic flow instability using pressure drop measurements of the friction factor versus Reynolds number. In prior work, we showed that this flow instability is characterized by large spatiotemporal fluctuations in the pore-scale flow velocities that may promote analyte dispersion across the column. Axial dispersion measurements of the quasi non-retained tracer thiourea confirm this possibility: they unequivocally reveal that operating above the onset of the instability improves column efficiency by significantly reducing peak asymmetry. These experiments thereby provide a proof-of-concept demonstration that elastic flow instabilities can be harnessed to mitigate the negative impact of trans-column flow heterogeneities on the efficiency of slurry-packed HPLC columns. While this approach has its own inherent limitations and constraints, our work lays the groundwork for future targeted development of polymers that can impart elasticity when dissolved in commonly used liquid chromatography mobile phases, and can thereby generate elastic flow instabilities to help improve the resolution of HPLC columns.

Diverse chemical, energy, environmental, and industrial processes involve the flow of polymer solutions in porous media. The accumulation and dissipation of elastic stresses as the polymers are transported through the tortuous, confined pore space can lead to the development of an elastic flow instability above a threshold flow rate. This flow instability can generate complex flows with strong spatiotemporal fluctuations, despite the low Reynolds number ($\mathrm{Re} \ll 1$); for example, in 1D ordered arrays of pore constrictions, this unstable flow can be multistable, with distinct pores exhibiting distinct unstable flow states. Here, we examine how this multistability is influenced by fluid rheology. Through experiments using diverse polymer solutions having systematic variations in fluid shear-thinning or elasticity, in pore constriction arrays of varying geometries, we show that the onset of multistability can be described using a single dimensionless parameter. This parameter, the streamwise Deborah number, compares the stress relaxation time of the polymer solution to the time required for the fluid to be advected between pore constrictions. Our work thus helps to deepen understanding of the influence of fluid rheology on elastic instabilities, helping to establish guidelines for the rational design of polymeric fluids with desirable flow behaviors.

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 (EI), 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 $3\times$ reduction in the required mixing length, a $6\times$ increase in solute transverse dispersivity, and can be harnessed to increase the rate at which chemical compounds react by $3\times$ -- enhancements that we rationalize using turbulence-inspired modeling of the underlying transport processes. Our work thereby establishes a simple, robust, versatile, and predictive new way to mix solutes in porous media, with potential applications ranging from large-scale chemical production to environmental remediation.

Inspired by presentations by the authors during a workshop organized at the Princeton Center for Theoretical Science (PCTS) in January 2022, we present a perspective on some of the outstanding questions related to the "physics of the ground beneath our feet." These identified challenges are intrinsically shared with the field of Soft Matter but also have unique aspects when the natural environment is studied., Comment: Perspective Paper, 30 pages, 15 figures

Obstructions influence the growth and expansion of bodies in a wide range of settings -- but isolating and understanding their impact can be difficult in complex environments. Here, we study obstructed growth/expansion in a model system accessible to experiments, simulations, and theory: hydrogels swelling around fixed cylindrical obstacles with varying geometries. When the obstacles are large and widely-spaced, hydrogels swell around them and remain intact. In contrast, our experiments reveal that when the obstacles are narrow and closely-spaced, hydrogels fracture as they swell. We use finite element simulations to map the magnitude and spatial distribution of stresses that build up during swelling at equilibrium in a 2D model, providing a route toward predicting when this phenomenon of self-fracturing is likely to arise. Applying lessons from indentation theory, poroelasticity, and nonlinear continuum mechanics, we also develop a theoretical framework for understanding how the maximum principal tensile and compressive stresses that develop during swelling are controlled by obstacle geometry and material parameters. These results thus help to shed light on the mechanical principles underlying growth/expansion in environments with obstructions., Comment: 13 pages, 5 figures; SI: 16 pages, 14 figures

We recast the problem of hydrogel swelling under physical constraints as an energy optimization problem. We apply this approach to compute equilibrium shapes of hydrogel spheres confined within a jammed matrix of rigid beads, and interpret the results to determine how confinement modifies the mechanics of swollen hydrogels. In contrast to the unconfined case, we find a spatial separation of strains within the bulk of the hydrogel as strain becomes localized to an outer region. We also explore the contact mechanics of the gel, finding a transition from Hertzian behavior to non-Hertzian behavior as a function of swelling. Our model, implemented in the Morpho shape optimization environment, can be applied in any dimension, readily adapted to diverse swelling scenarios and extended to use other energies in conjunction.

Drainage, in which a nonwetting fluid displaces a wetting fluid from a porous medium, is well-studied for media with unchanging solid surfaces. However, many media can be eroded by drainage, with eroded material redeposited in pores downstream, altering further flow. Here, we use theory and simulation to examine how these coupled processes both alter the overall fluid displacement pathway and help reshape the solid medium. We find two new drainage behaviors with markedly different characteristics, and quantitatively delineate the conditions under which they arise. Our results thereby help expand current understanding of these rich physics, with implications for applications of drainage in industry and the environment.

Collectives of actively-moving particles can spontaneously separate into dilute and dense phases -- a fascinating phenomenon known as motility-induced phase separation (MIPS). MIPS is well-studied for randomly-moving particles with no directional bias. However, many forms of active matter exhibit collective chemotaxis, directed motion along a chemical gradient that the constituent particles can generate themselves. Here, using theory and simulations, we demonstrate that collective chemotaxis strongly competes with MIPS -- in some cases, arresting or completely suppressing phase separation, or in other cases, generating fundamentally new dynamic instabilities. We establish quantitative principles describing this competition, thereby helping to reveal and clarify the rich physics underlying active matter systems that perform chemotaxis, ranging from cells to robots.

We use a theoretical model to explore how fluid dynamics, in particular, the pressure gradient and wall shear stress in a channel, affect the deposition of particles flowing in a microfluidic network. Experiments on transport of colloidal particles in pressure-driven systems of packed beads have shown that at lower pressure drop, particles deposit locally at the inlet, while at higher pressure drop, they deposit uniformly along the direction of flow. We develop a mathematical model and use agent-based simulations to capture these essential qualitative features observed in experiments. We explore the deposition profile over a two-dimensional phase diagram defined in terms of the pressure and shear stress threshold, and show that two distinct phases exist. We explain this apparent phase transition by drawing an analogy to simple one-dimensional models of aggregation in which the phase transition is calculated analytically., Comment: 10 pages, 8 figures including Supplemental Material

E. coli chemotactic motion in the presence of a chemoattractant field has been extensively studied using wet laboratory experiments, stochastic computational models as well as partial differential equation-based models (PDEs). The most challenging step in bridging these approaches, is establishing a closed form of the so-called chemotactic term, which describes how bacteria bias their motion up chemonutrient concentration gradients, as a result of a cascade of biochemical processes. Data-driven models can be used to learn the entire evolution operator of the chemotactic PDEs (black box models), or, in a more targeted fashion, to learn just the chemotactic term (gray box models). In this work, data-driven Machine Learning approaches for learning the underlying model PDEs are (a) validated through the use of simulation data from established continuum models and (b) used to infer chemotactic PDEs from experimental data. Even when the data at hand are sparse (coarse in space and/or time), noisy (due to inherent stochasticity in measurements) or partial (e.g. lack of measurements of the associated chemoattractant field), we can attempt to learn the right-hand-side of a closed PDE for an evolving bacterial density. In fact we show that data-driven PDEs including a short history of the bacterial density field (e.g. in the form of higher-order in time PDEs in terms of the measurable bacterial density) can be successful in predicting further bacterial density evolution, and even possibly recovering estimates of the unmeasured chemonutrient field. The main tool in this effort is the effective low-dimensionality of the dynamics (in the spirit of the Whitney and Takens embedding theorems). The resulting data-driven PDE can then be simulated to reproduce/predict computational or experimental bacterial density profile data, and estimate the underlying (unmeasured) chemonutrient field evolution., Comment: 44 pages, 20 figures, 3 tables

Many key environmental, industrial, and energy processes rely on controlling fluid transport within subsurface porous media. These media are typically structurally heterogeneous, often with vertically-layered strata of distinct permeabilities -- leading to uneven partitioning of flow across strata, which can be undesirable. Here, using direct in situ visualization, we demonstrate that polymer additives can homogenize this flow by inducing a purely-elastic flow instability that generates random spatiotemporal fluctuations and excess flow resistance in individual strata. In particular, we find that this instability arises at smaller imposed flow rates in higher-permeability strata, diverting flow towards lower-permeability strata and helping to homogenize the flow. Guided by the experiments, we develop a parallel-resistor model that quantitatively predicts the flow rate at which this homogenization is optimized for a given stratified medium. Thus, our work provides a new approach to homogenizing fluid and passive scalar transport in heterogeneous porous media.

Bacteria are ubiquitous in our daily lives, either as motile planktonic cells or as immobilized surface-attached biofilms. These different phenotypic states play key roles in agriculture, environment, industry, and medicine; hence, it is critically important to be able to predict the conditions under which bacteria transition from one state to the other. Unfortunately, these transitions depend on a dizzyingly complex array of factors that are determined by the intrinsic properties of the individual cells as well as those of their surrounding environments, and are thus challenging to describe. To address this issue, here, we develop a generally-applicable biophysical model of the interplay between motility-mediated dispersal and biofilm formation under positive quorum sensing control. Using this model, we establish a universal rule predicting how the onset and extent of biofilm formation depend collectively on cell concentration and motility, nutrient diffusion and consumption, chemotactic sensing, and autoinducer production. Our work thus provides a key step toward quantitatively predicting and controlling biofilm formation in diverse and complex settings.

From pasta to biological tissues to contact lenses, gel and gel-like materials inherently soften as they swell with water. In dry, low-relative-humidity environments, these materials stiffen as they de-swell with water. Here, we use semi-dilute polymer theory to develop a simple power-law relationship between hydrogel elastic modulus and swelling. From this relationship, we predict hydrogel stiffness or swelling at arbitrary relative humidities. Our close predictions of properties of hydrogels across three different polymer mesh families at varying crosslinking densities and relative humidities demonstrate the validity and generality of our understanding. This predictive capability enables more rapid material discovery and selection for hydrogel applications in varying humidity environments., Comment: This is a pre-peer-review manuscript

Viscoelastic fluids are a common subclass of rheologically complex materials that are encountered in diverse fields from biology to polymer processing. Often the flows of viscoelastic fluids are unstable in situations where ordinary Newtonian fluids are stable, owing to the nonlinear coupling of the elastic and viscous stresses. Perhaps more surprisingly, the instabilities produce flows with the hallmarks of turbulence -- even though the effective Reynolds numbers may be $O(1)$ or smaller. We provide perspectives on viscoelastic flow instabilities by integrating the input from speakers at a recent international workshop: historical remarks, characterization of fluids and flows, discussion of experimental and simulation tools, and modern questions and puzzles that motivate further studies of this fascinating subject. The materials here will be useful for researchers and educators alike, especially as the subject continues to evolve in both fundamental understanding and applications in engineering and the sciences., Comment: A perspective article based on a virtual workshop of the Princeton Center for Theoretical Sciences that took place in January 2021

The ability of many living systems to actively self-propel underlies critical biomedical, environmental, and industrial processes. While such active transport is well-studied in uniform settings, environmental complexities such as geometric constraints, mechanical cues, and external stimuli such as chemical gradients and fluid flow can strongly influence transport. In this chapter, we describe recent progress in the study of active transport in such complex environments, focusing on two prominent biological systems -- bacteria and eukaryotic cells -- as archetypes of active matter. We review research findings highlighting how environmental factors can fundamentally alter cellular motility, hindering or promoting active transport in unexpected ways, and giving rise to fascinating behaviors such as directed migration and large-scale clustering. In parallel, we describe specific open questions and promising avenues for future research. Furthermore, given the diverse forms of active matter -- ranging from enzymes and driven biopolymer assemblies, to microorganisms and synthetic microswimmers, to larger animals and even robots -- we also describe connections to other active systems as well as more general theoretical/computational models of transport processes in complex environments., Comment: Preprint of a chapter in the book \textit{Out-of-Equilibrium Soft Matter: Active Fluids}, to be published by the Royal Society of Chemistry; revised version has added references and corrected typos

The spreading of bacterial populations is central to processes in agriculture, the environment, and medicine. However, existing models of spreading typically focus on cells in unconfined settings--despite the fact that many bacteria inhabit complex and crowded environments, such as soils, sediments, and biological tissues/gels, in which solid obstacles confine the cells and thereby strongly regulate population spreading. Here, we develop an extended version of the classic Keller-Segel model of bacterial spreading that incorporates the influence of confinement in promoting both cell-solid and cell-cell collisions. Numerical simulations of this extended model demonstrate how confinement fundamentally alters the dynamics and morphology of spreading bacterial populations, in good agreement with recent experimental results. In particular, with increasing confinement, we find that cell-cell collisions increasingly hinder the initial formation and the long-time propagation speed of chemotactic pulses. Moreover, also with increasing confinement, we find that cellular growth and division plays an increasingly dominant role in driving population spreading--eventually leading to a transition from chemotactic spreading to growth-driven spreading via a slower, jammed front. This work thus provides a theoretical foundation for further investigations of the influence of confinement on bacterial spreading. More broadly, these results help to provide a framework to predict and control the dynamics of bacterial populations in complex and crowded environments.

Imbibition, the displacement of a nonwetting fluid by a wetting fluid, plays a central role in diverse energy, environmental, and industrial processes. While this process is typically studied in homogeneous porous media with uniform permeabilities, in many cases, the media have multiple parallel strata of different permeabilities. How such stratification impacts the fluid dynamics of imbibition, as well as the fluid saturation after the wetting fluid breaks through to the end of a given medium, is poorly understood. We address this gap in knowledge by developing an analytical model of imbibition in a porous medium with two parallel strata, combined with a pore network model that explicitly describes fluid crossflow between the strata. By numerically solving these models, we examine the fluid dynamics and fluid saturation left after breakthrough. We find that the breakthrough saturation of nonwetting fluid is minimized when the imposed capillary number Ca is tuned to a value Ca$^*$ that depends on both the structure of the medium and the viscosity ratio between the two fluids. Our results thus provide quantitative guidelines for predicting and controlling flow in stratified porous media, with implications for water remediation, oil/gas recovery, and applications requiring moisture management in diverse materials.

In contexts ranging from embryonic development to bacterial ecology, cell populations migrate chemotactically along self-generated chemical gradients, often forming a propagating front. Here, we theoretically show that the stability of such chemotactic fronts to morphological perturbations is determined by limitations in the ability of individual cells to sense and thereby respond to the chemical gradient. Specifically, cells at bulging parts of a front are exposed to a smaller gradient, which slows them down and promotes stability, but they also respond more strongly to the gradient, which speeds them up and promotes instability. We predict that this competition leads to chemotactic fingering when sensing is limited at too low chemical concentrations. Guided by this finding and by experimental data on E. coli chemotaxis, we suggest that the cells' sensory machinery might have evolved to avoid these limitations and ensure stable front propagation. Finally, as sensing of any stimuli is necessarily limited in living and active matter in general, the principle of sensing-induced stability may operate in other types of directed migration such as durotaxis, electrotaxis, and phototaxis.

We perform Brownian dynamics simulations of active stiff polymers undergoing run-reverse dynamics, and so mimic bacterial swimming, in porous media. In accord with recent experiments of \emph{Escherichia coli}, the polymer dynamics are characterized by trapping phases interrupted by directed hopping motion through the pores. We find that the effective translational diffusivities of run-reverse agents can be enhanced up to two orders in magnitude, compared to their non-reversing counterparts, and exhibit a non-monotonic behavior as a function of the reversal rate, which we rationalize using a coarse-grained model. Furthermore, we discover a geometric criterion for the optimal spreading, which emerges when their run lengths are comparable to the longest straight path available in the porous medium. More significantly, our criterion unifies results for porous media with disparate pore sizes and shapes and thus provides a fundamental principle for optimal transport of microorganisms and cargo-carriers in densely-packed biological and environmental settings.

Diverse processes -- e.g., environmental pollution, groundwater remediation, oil recovery, filtration, and drug delivery -- involve the transport of colloidal particles in porous media. Using confocal microscopy, we directly visualize this process in situ and thereby identify the fundamental mechanisms by which particles are distributed throughout a medium. At high injection pressures, hydrodynamic stresses cause particles to be continually deposited on and eroded from the solid matrix -- notably, forcing them to be distributed throughout the entire medium. By contrast, at low injection pressures, the relative influence of erosion is suppressed, causing particles to localize near the inlet of the medium. Unexpectedly, these macroscopic distribution behaviors depend on imposed pressure in similar ways for particles of different charges, although the pore-scale distribution of deposition is sensitive to particle charge. These results reveal how the multiscale interactions between fluid, particles, and the solid matrix control how colloids are distributed in a porous medium., Comment: This article is published in Science Advances and supplementary material is available at https://advances.sciencemag.org/content/6/46/eabc2530

Hydrogels hold promise in agriculture as reservoirs of water in dry soil, potentially alleviating the burden of irrigation. However, confinement in soil can markedly reduce the ability of hydrogels to absorb water and swell, limiting their widespread adoption. Unfortunately, the underlying reason remains unknown. By directly visualizing the swelling of hydrogels confined in three-dimensional granular media, we demonstrate that the extent of hydrogel swelling is determined by the competition between the force exerted by the hydrogel due to osmotic swelling and the confining force transmitted by the surrounding grains. Furthermore, the medium can itself be restructured by hydrogel swelling, as set by the balance between the osmotic swelling force, the confining force, and intergrain friction. Together, our results provide quantitative principles to predict how hydrogels behave in confinement, potentially improving their use in agriculture as well as informing other applications such as oil recovery, construction, mechanobiology, and filtration., Comment: This article is published in Science Advances and supplementary material is available at https://advances.sciencemag.org/content/7/7/eabd2711

Models of disease spreading are critical for predicting infection growth in a population and evaluating public health policies. However, standard models typically represent the dynamics of disease transmission between individuals using macroscopic parameters that do not accurately represent person-to-person variability. To address this issue, we present a dynamic network model that provides a straightforward way to incorporate both disease transmission dynamics at the individual scale as well as the full spatiotemporal history of infection at the population scale. We find that disease spreads through a social network as a traveling wave of infection, followed by a traveling wave of recovery, with the onset and dynamics of spreading determined by the interplay between disease transmission and recovery. We use these insights to develop a scaling theory that predicts the dynamics of infection for diverse diseases and populations. Furthermore, we show how spatial heterogeneities in susceptibility to infection can either exacerbate or quell the spread of disease, depending on its infectivity. Ultimately, our dynamic network approach provides a simple way to model disease spreading that unifies previous findings and can be generalized to diverse diseases, containment strategies, seasonal conditions, and community structures., Comment: In press, Frontiers in Physics (2021)

Collective migration -- the directed, coordinated motion of many self-propelled agents -- is a fascinating emergent behavior exhibited by active matter that has key functional implications for biological systems. Extensive studies have elucidated the different ways in which this phenomenon may arise. Nevertheless, how collective migration can persist when a population is confronted with perturbations, which inevitably arise in complex settings, is poorly understood. Here, by combining experiments and simulations, we describe a mechanism by which collectively migrating populations smooth out large-scale perturbations in their overall morphology, enabling their constituents to continue to migrate together. We focus on the canonical example of chemotactic migration of Escherichia coli, in which fronts of cells move via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We identify two distinct modes in which chemotaxis influences the morphology of the population: cells in different locations along a front migrate at different velocities due to spatial variations in (i) the local nutrient gradient and in (ii) the ability of cells to sense and respond to the local nutrient gradient. While the first mode is destabilizing, the second mode is stabilizing and dominates, ultimately driving smoothing of the overall population and enabling continued collective migration. This process is autonomous, arising without any external intervention; instead, it is a population-scale consequence of the manner in which individual cells transduce external signals. Our findings thus provide insights to predict, and potentially control, the collective migration and morphology of cell populations and diverse other forms of active matter.

Colloidal particles hold promise for mobilizing and removing trapped immiscible fluids from porous media, with implications for key energy and water applications. Most studies focus on accomplishing this goal using particles that can localize at the immiscible fluid interface. Therefore, researchers typically seek to optimize the surface activity of particles, as well as their ability to freely move through a pore space with minimal deposition onto the surrounding solid matrix. Here, we demonstrate that deposition can, surprisingly, promote mobilization of a trapped fluid from a porous medium without requiring any surface activity. Using confocal microscopy, we directly visualize both colloidal particles and trapped immiscible fluid within a transparent, three-dimensional (3D) porous medium. We find that as non-surface active particles deposit on the solid matrix, increasing amounts of trapped fluid become mobilized. We unravel the underlying physics by analyzing the extent of deposition, as well as the geometry of trapped fluid droplets, at the pore scale: deposition increases the viscous stresses on trapped droplets, overcoming the influence of capillarity that keeps them trapped. Given an initial distribution of trapped fluid, this analysis enables us to predict the extent of fluid mobilized through colloidal deposition. Taken together, our work reveals a new way by which colloids can be harnessed to mobilize trapped fluid from a porous medium.

Diverse processes rely on the viscous flow of polymer solutions through porous media. In many cases, the macroscopic flow resistance abruptly increases above a threshold flow rate in a porous medium---but not in bulk solution. The reason why has been a puzzle for over half a century. Here, by directly visualizing the flow in a transparent 3D porous medium, we demonstrate that this anomalous increase is due to the onset of an elastic instability. We establish that the energy dissipated by the unstable flow fluctuations, which vary across pores, generates the anomalous increase in flow resistance through the entire medium. Thus, by linking the pore-scale onset of unstable flow to macroscopic transport, our work provides generally-applicable guidelines for predicting and controlling polymer solution flows.

The seemingly simple process of inhalation relies on a complex interplay between muscular contraction in the thorax, elasto-capillary interactions in individual lung branches, propagation of air between different connected branches, and overall air flow into the lungs. These processes occur over considerably different length and time scales; consequently, linking them to the biomechanical properties of the lungs, and quantifying how they together control the spatiotemporal features of inhalation, remains a challenge. We address this challenge by developing a computational model of the lungs as a hierarchical, branched network of connected liquid-lined flexible cylinders coupled to a viscoelastic thoracic cavity. Each branch opens at a rate and a pressure that is determined by input biomechanical parameters, enabling us to test the influence of changes in the mechanical properties of lung tissues and secretions on inhalation dynamics. By summing the dynamics of all the branches, we quantify the evolution of overall lung pressure and volume during inhalation, reproducing the shape of measured breathing curves. Using this model, we demonstrate how changes in lung muscle contraction, mucus viscosity and surface tension, and airway wall stiffness---characteristic of many respiratory diseases, including those arising from COVID-19, cystic fibrosis, chronic obstructive pulmonary disease, asthma, and emphysema---drastically alter inhaled lung capacity and breathing duration. Our work therefore helps to identify the key factors that control breathing dynamics, and provides a way to quantify how disease-induced changes in these factors lead to respiratory distress., Comment: In press

Imbibition plays a central role in diverse energy, environmental, and industrial processes. In many cases, the medium has multiple parallel strata of different permeabilities; however, how this stratification impacts imbibition is poorly understood. We address this gap in knowledge by directly visualizing forced imbibition in three-dimensional (3D) porous media with two parallel strata. We find that imbibition is spatially heterogeneous: for small capillary number Ca, the wetting fluid preferentially invades the fine stratum, while for Ca above a threshold value, the fluid instead preferentially invades the coarse stratum. This threshold value depends on the medium geometry, the fluid properties, and the presence of residual wetting films in the pore space. These findings are well described by a linear stability analysis that incorporates crossflow between the strata. Thus, our work provides quantitative guidelines for predicting and controlling flow in stratified porous media., Comment: In press

Polymer solutions are often injected in porous media for applications such as oil recovery and groundwater remediation. As the fluid navigates the tortuous pore space, elastic stresses build up, causing the flow to become unstable at sufficiently large injection rates. However, it is poorly understood how the spatial and temporal characteristics of this unstable flow depend on pore space geometry. We elucidate this dependence by systematically varying the spacing between pore constrictions in a model porous medium. We find that when the pore spacing is large, unstable eddies form upstream of each pore, similar to observations of an isolated pore. By contrast, when the pore spacing is sufficiently small, the flow exhibits a surprising bistability, stochastically switching between two distinct unstable flow states. We hypothesize that this unusual behavior arises from the interplay between flow-induced polymer elongation and relaxation of polymers as they are advected through the porous medium. Consistent with this idea, we find that the flow state in a given pore persists for long times. Moreover, we find that the flow state is correlated between neighboring pores; however, these correlations do not persist long-range. Our results thus help to elucidate the rich array of flow behaviors that can arise in polymer solution flow through porous media., Comment: In press; this version is pre-review. Updated version will be posted in six months

Understanding how bacteria move in porous media is critical to applications in healthcare, agriculture, environmental remediation, and chemical sensing. Recent work has demonstrated that E. coli, which moves by run-and-tumble dynamics in a homogeneous medium, exhibits a new form of motility when confined in a disordered porous medium: hopping-and-trapping motility, in which cells perform rapid, directed hops punctuated by intervals of slow, undirected trapping. Here, we use direct visualization to shed light on how these processes depend on pore-scale confinement and cellular activity. We find that hopping is determined by pore-scale confinement, and is independent of cellular activity; by contrast, trapping is determined by the competition between pore-scale confinement and cellular activity, as predicted by an entropic trapping model. These results thus help to elucidate the factors that regulate bacterial motion in porous media, and could help aid the development of new models of motility in heterogeneous environments.

Hydrated granular packings often crack into discrete clusters of grains when dried. Despite its ubiquity, accurate prediction of cracking remains elusive. Here, we elucidate the previously overlooked role of individual grain shrinkage---a feature common to many materials---in determining crack patterning using both experiments and simulations. By extending the classical Griffith crack theory, we obtain a scaling law that quantifies how cluster size depends on the interplay between grain shrinkage, stiffness, and size---applicable to a diverse array of shrinkable, granular packings.

Polymer solutions are frequently used in enhanced oil recovery and groundwater remediation to improve the recovery of trapped non-aqueous 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 non-trivial flow patterns and impact macroscopic transport remains an outstanding challenge. Here, we describe 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 non-aqueous fluid recovery. This work provides new insights for polymer dynamics, non-Newtonian fluid mechanics, and applications such as enhanced oil recovery and groundwater remediation., Comment: 45 pages, 12 figures

Capillary fingering is a displacement process that can occur when a non-wetting fluid displaces a wetting fluid from a homogeneous disordered porous medium. Here, we investigate how this process is influenced by a pore size gradient. Using microfluidic experiments and computational pore-network models, we show that the non-wetting fluid displacement behavior depends sensitively on the direction and the magnitude of the gradient. The fluid displacement depends on the competition between a pore size gradient and pore-scale disorder; indeed, a sufficiently large gradient can completely suppress capillary fingering. By analyzing capillary forces at the pore scale, we identify a non-dimensional parameter that describes the physics underlying these diverse flow behaviors. Our results thus expand the understanding of flow in complex porous media, and suggest a new way to control flow behavior via the introduction of pore size gradients., Comment: In press, Physical Review Fluids (2019)

Many clays, soils, biological tissues, foods, and coatings are shrinkable, granular materials: they are composed of packed, hydrated grains that shrink when dried. In many cases, these packings crack during drying, critically hindering applications. However, while cracking has been widely studied for bulk gels and packings of non-shrinkable grains, little is known about how packings of shrinkable grains crack. Here, we elucidate how grain shrinkage alters cracking during drying. Using experiments with model shrinkable hydrogel beads, we show that differential shrinkage can dramatically alter crack evolution during drying---in some cases, even causing cracks to spontaneously "self-close". In other cases, packings shrink without cracking or crack irreversibly. We developed both granular and continuum models to quantify the interplay between grain shrinkage, poromechanics, packing size, drying rate, capillarity, and substrate friction on cracking. Guided by the theory, we also found that cracking can be completely altered by varying the spatial profile of drying. Our work elucidates the rich physics underlying cracking in shrinkable, granular packings, and yields new strategies for controlling crack evolution., Comment: In press

Diverse processes--e.g. bioremediation, biofertilization, and microbial drug delivery--rely on bacterial migration in disordered, three-dimensional (3D) porous media. However, how pore-scale confinement alters bacterial motility is unknown due to the opacity of typical 3D media. As a result, models of migration are limited and often employ ad hoc assumptions. Here we reveal that the paradigm of run-and-tumble motility is dramatically altered in a porous medium. By directly visualizing individual Escherichia coli, we find that the cells are intermittently and transiently trapped as they navigate the pore space, exhibiting diffusive behavior at long time scales. The trapping durations and the lengths of "hops" between traps are broadly distributed, reminiscent of transport in diverse other disordered systems; nevertheless, we show that these quantities can together predict the long-time bacterial translational diffusivity. Our work thus provides a revised picture of bacterial motility in complex media and yields principles for predicting cellular migration.

We use confocal microscopy to directly visualize the simultaneous flow of both a wetting and a non-wetting fluid through a model three-dimensional (3D) porous medium. We find that, for small flow rates, both fluids flow through unchanging, distinct, connected 3D pathways; in stark contrast, at sufficiently large flow rates, the non-wetting fluid is broken up into discrete ganglia. By performing experiments over a range of flow rates, using fluids of different viscosities, and with porous media having different geometries, we show that this transition can be characterized by a state diagram that depends on the capillary numbers of both fluids, suggesting that it is controlled by the competition between the viscous forces exerted on the flowing oil and the capillary forces at the pore scale. Our results thus help elucidate the diverse range of behaviors that arise in two-phase flow through a 3D porous medium., Comment: Physics of Fluids (2014)

We use confocal microscopy to directly visualize the formation and complex morphologies of trapped non-wetting fluid ganglia within a model 3D porous medium. The wetting fluid continues to flow around the ganglia after they form; this flow is characterized by a capillary number, Ca. We find that the ganglia configurations do not vary for small Ca; by contrast, as Ca is increased above a threshold value, the largest ganglia start to become mobilized and are ultimately removed from the medium. By combining our 3D visualization with measurements of the bulk transport, we show that this behavior can be quantitatively understood by balancing the viscous forces exerted on the ganglia with the pore-scale capillary forces that keep them trapped within the medium. Our work thus helps elucidate the fluid dynamics underlying the mobilization of a trapped non-wetting fluid from a 3D porous medium.