Disorder Suppresses Chaos in Viscoelastic Flows

Harnessing elastic instabilities for enhanced mixing and reaction kinetics in porous media

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 [Formula: see text] reduction in the required mixing length, a [Formula: see text] increase in solute transverse dispersivity, and can be harnessed to increase the rate at which chemical compounds react by [Formula: see text]-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.

Geometry-Dependent Elastic Flow Dynamics in Micropillar Arrays

Regular device-scale DNA waves for high DNA concentrations and flow velocities have been shown to emerge in quadratic micropillar arrays with potentially strong relevance for a wide range of microfluidic applications. Hexagonal arrays constitute another geometry that is especially relevant for the microfluidic pulsed-field separation of DNA. Here, we report on the differences at the micro and macroscopic scales between the resulting wave patterns for these two regular array geometries and one disordered array geometry. In contrast to the large-scale regular waves visible in the quadratic array, in the hexagonal arrays, waves occur in a device-scale disordered zig-zag pattern with fluctuations on a much smaller scale. We connect the large-scale pattern to the microscopic flow and observe flow synchronization that switches between two directions for both the quadratic and hexagonal arrays. We show the importance of order using the disordered array, where steady-state stationary and highly fluctuating flow states persist in seemingly random locations across the array. We compare the flow dynamics of the arrays to that in a device with sparsely distributed pillars. Here, we observe similar vortex shedding, which is clearly observable in the quadratic and disordered arrays. However, the shedding of these vortices couples only in the flow direction and not laterally as in the dense, ordered arrays. We believe that our findings will contribute to the understanding of elastic flow dynamics in pillar arrays, helping us elucidate the fundamental principles of non-Newtonian fluid flow in complex environments as well as supporting applications in engineering involving e.g., transport, sorting, and mixing of complex fluids.

Particle dispersion through porous media with heterogeneous attractions

Porous media used in many practical applications contain natural spatial variations in composition and surface charge that lead to heterogeneous physicochemical attractions between the media and transported particles. We performed Stokesian dynamics (SD) simulations to examine the effects of heterogeneous attractions on quiescent diffusion and hydrodynamic dispersion of particles within geometrically ordered arrays of nanoposts. We find that transport under quiescent conditions occurs by two mechanisms, diffusion through the void space and intermittent hopping between the attractive wells of different nanoposts. As the attraction heterogeneity increases, the latter mechanism becomes dominant, resulting in an increase in the particle trajectory tortuosity, deviations from Gaussian behavior in the particle displacement distributions, and a decrease in the long-time particle diffusivity. Similarly, under flow conditions corresponding to low Péclet number (Pe), increased attraction heterogeneity leads to transient localization near the nanoposts, resulting in a broadening of the particle distribution and enhanced longitudinal dispersion in the direction of flow. At high Pe where advection strongly dominates, however, the longitudinal dispersion coefficient is insensitive to attraction heterogeneity and exhibits Taylor-Aris dispersion behavior. Our findings provide insight into how heterogeneous interactions may influence particle transport in complex 3-D porous media.

Using symmetry to control viscoelastic waves in pillar arrays

Solutions of macromolecules exhibit viscoelastic properties and unlike Newtonian fluids, they may break time-reversal symmetry at low Reynolds numbers resulting in elastic turbulence. Furthermore, under some conditions, instead of the chaotic turbulence, the result is large-scale waves in the form of cyclic spatial and temporal concentration variations, as has been shown for macromolecular DNA flowing in microfluidic pillar arrays. We here demonstrate how altering the symmetry of the individual pillars can be used to influence the symmetry of these waves. We control the extent of instabilities in viscoelastic flow by leveraging the effects of the symmetry of the pillars on the waves, demonstrating suppressed viscoelastic fluctuations with relevance for transport and sorting applications, or conversely opening up for enhanced viscoelasticity-mediated mixing. The onset of waves, which changes flow resistance, occurs at different Deborah numbers for flow in different directions through the array of triangular pillars, thus breaking the symmetry of the flow resistance along the device, opening up for using the occurrence of the waves to construct a fluidic diode.

Anomalous crystalline ordering of particles in a viscoelastic fluid under high shear

Addition of particles to a viscoelastic suspension dramatically alters the properties of the mixture, particularly when it is sheared or otherwise processed. Shear-induced stretching of the polymers results in elastic stress that causes a substantial increase in measured viscosity with increasing shear, and an attractive interaction between particles, leading to their chaining. At even higher shear rates, the flow becomes unstable, even in the absence of particles. This instability makes it very difficult to determine the properties of a particle suspension. Here, we use a fully immersed parallel plate geometry to measure the high-shear-rate behavior of a suspension of particles in a viscoelastic fluid. We find an unexpected separation of the particles within the suspension resulting in the formation of a layer of particles in the center of the cell. Remarkably, monodisperse particles form a crystalline layer which dramatically alters the shear instability. By combining measurements of the velocity field and torque fluctuations, we show that this solid layer disrupts the flow instability and introduces a single-frequency component to the torque fluctuations that reflects a dominant velocity pattern in the flow. These results highlight the interplay between particles and a suspending viscoelastic fluid at very high shear rates.

Stress and stretching regulate dispersion in viscoelastic porous media flows

In this work, we study the role of viscoelastic instability in the mechanical dispersion of fluid flow through porous media at high Péclet numbers. Using microfluidic experiments and numerical simulations, we show that viscoelastic instability in flow through a hexagonally ordered (staggered) medium strongly enhances dispersion transverse to the mean flow direction with increasing Weissenberg number (Wi). In contrast, preferential flow paths can quench the elastic instability in disordered media, which has two important consequences for transport: first, the lack of chaotic velocity fluctuations reduces transverse dispersion relative to unstable flows. Second, the amplification of flow along preferential paths with increasing Wi causes strongly-correlated stream-wise flow that enhances longitudinal dispersion. Finally, we illustrate how the observed dispersion phenomena can be understood through the lens of Lagrangian stretching manifolds, which act as advective transport barriers and coincide with high stress regions in these viscoelastic porous media flows.

Short and long-range cyclic patterns in flows of DNA solutions in microfluidic obstacle arrays

We observe regular patterns emerging across multiple length scales with high-concentration DNA solutions in microfluidic pillar arrays at low Reynolds numbers and high Deborah numbers. Interacting vortices between pillars lead to long-range order in the form of large travelling waves consisting of DNA at high concentration and extension. Waves are formed in quadratic arrays of pillars, while randomizing the position of the pillar in each unit cell of a quadratic array leads to suppression of the long-range patterns. We find that concentrations exceeding the overlap concentration of the DNA enables the waves, and exploring the behavior of the waves as a function of flow rate, buffer composition, concentration and molecular length, we identify elastic effects as central to the origin of the waves. Our work may not only help increase the low throughput that often limits sample processing in microfluidics, it may also provide a platform for further studies of the underlying viscoelastic mechanisms.

Experimental insights into elasto-inertial transitions in Taylor-Couette flows

Since the seminal work of Taylor in 1923, Taylor-Couette (TC) flow has served as a paradigm to study hydrodynamic instabilities and bifurcation phenomena. Transitions of Newtonian TC flows to inertial turbulence have been extensively studied and are well understood, while in the past few years, there has been an increasing interest in TC flows of complex, viscoelastic fluids. The transitions to elastic turbulence (ET) or elasto-inertial turbulence (EIT) have revealed fascinating dynamics and flow states; depending on the rheological properties of the fluids, a broad spectrum of transitions has been reported, including rotating standing waves, flame patterns (FP), and diwhirls (DW). The nature of these transitions and the relationship between ET and EIT are not fully understood. In this review, we discuss experimental efforts on TC flows of viscoelastic fluids. We outline the experimental methods employed and the non-dimensional parameters of interest, followed by an overview of inertia, elasticity and elasto-inertia-driven transitions to turbulence and their modulation through shear thinning or particle suspensions. The published experimental data are collated, and a map of flow transitions to EIT as a function of the key fluid parameters is provided, alongside perspectives for the future work. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (part 1)'.

Upstream wall vortices in viscoelastic flow past a cylinder

We report a novel inertia-less, elastic flow instability for a viscoelastic, shear-thinning wormlike micellar solution flowing past a microcylinder in a channel with blockage ratio B_R = 2R/W = 0.5 and aspect ratio α = H/W ≈ 5, where R ≈ 100 μm is the cylinder radius, W is the channel width, and H is the channel height. The instability manifests upstream of the cylinder and changes form with increasing Weissenberg number over the range 0.5 ≲ Wi = Uλ/R ≲ 900, where U is the average flow velocity and λ is the terminal relaxation time of the fluid. Beyond a first critical Wi, the instability begins as a bending of the streamlines near the upstream pole of the cylinder that breaks the symmetry of the flow. Beyond a second critical Wi, small, time-steady, and approximately symmetric wall-attached vortices form upstream of the cylinder. Beyond a third critical Wi, the flow becomes time dependent and pulses with a characteristic frequency commensurate with the breakage timescale of the wormlike micelles. This is accompanied by a breaking of the symmetry of the wall-attached vortices, where one vortex becomes considerably larger than the other. Finally, beyond a fourth critical Wi, a vortex forms attached to the upstream pole of the cylinder whose length fluctuates in time. The flow is highly time dependent, and the cylinder-attached vortex and wall-attached vortices compete dynamically for space and time in the channel. Our results add to the rapidly growing understanding of viscoelastic flow instabilities in microfluidic geometries.

Nanoparticle transport within non-Newtonian fluid flow in porous media

Control over dispersion of nanoparticles in polymer solutions through porous media is important for subsurface applications such as soil remediation and enhanced oil recovery. Dispersion is affected by the spatial heterogeneity of porous media, the non-Newtonian behavior of polymer solutions, and the Brownian motion of nanoparticles. Here, we use the Euler-Lagrangian method to simulate the flow of nanoparticles and inelastic non-Newtonian fluids (described by Meter model) in a range of porous media samples and injection rates. In one case, we use a fine mesh of more than 3 million mesh points to model nanoparticles transport in a sandstone sample. The results show that the velocity distribution of nanoparticles in the porous medium is non-Gaussian, which leads to the non-Fickian behavior of nanoparticles dispersion. Due to pore-space confinement, the long-time mean-square displacement of nanoparticles depends nonlinearly on time. Additionally, the gradient of shear stress in the pore space of the porous medium dictates the transport behavior of nanoparticles in the porous medium. Furthermore, the Brownian motion of nanoparticles increases the dispersion of nanoparticles along the longitudinal and transverse direction.

Towards Predicting the Onset of Elastic Turbulence in Complex Geometries

Abstract

Flow of complex fluids in porous structures is pertinent in many biological and industrial processes. For these applications, elastic turbulence, a viscoelastic instability occurring at low Re—arising from a non-trivial coupling of fluid rheology and flow geometry—is a common and relevant effect because of significant over-proportional increase in pressure drop and spatio-temporal distortion of the flow field. Therefore, significant efforts have been made to predict the onset of elastic turbulence in flow geometries with constrictions. The onset of flow perturbations to fluid streamlines is not adequately captured by Deborah and Weissenberg numbers. The introduction of more complex dimensionless numbers such as the M-criterion, which was meant as a simple and pragmatic method to predict the onset of elastic instabilities as an order-of-magnitude estimate, has been successful for simpler geometries. However, for more complex geometries which are encountered in many relevant applications, sometimes discrepancies between experimental observation and M-criteria prediction have been encountered. So far these discrepancies have been mainly attributed to the emergence from disorder. In this experimental study, we employ a single channel with multiple constrictions at varying distance and aspect ratios. We show that adjacent constrictions can interact via non-laminar flow field instabilities caused by a combination of individual geometry and viscoelastic rheology depending (besides other factors) explicitly on the distance between adjacent constrictions. This provides intuitive insight on a more conceptual level why the M-criteria predictions are not more precise. Our findings suggest that coupling of rheological effects and fluid geometry is more complex and implicit and controlled by more length scales than are currently employed. For translating bulk fluid, rheology determined by classical rheometry into the effective behaviour in complex porous geometries requires consideration of more than only one repeat element. Our findings open the path towards more accurate prediction of the onset of elastic turbulence, which many applications will benefit.

Article Highlights

Oscillative Trapping of a Droplet in a Converging Channel Induced by Elastic Instability

Permanent trapping of an oscillating, nonwetting droplet is observed in a converging-diverging microchannel when aqueous, viscoelastic fluids are injected. Classical theories based on the balance between capillary and viscous forces suggest that the droplet size should decrease with increasing flow rates of a displacing Newtonian fluid, and the droplet should be completely displaced at high injection rates. However, droplets in viscoelastic fluids cannot be removed by increasing flow rates due to the oscillation. The oscillation amplitude linearly increases with the Deborah number (De), which further inhibits the droplet's passing through the constriction, "permanently." Our microfluidic experiments show that the onset of oscillation is determined by a critical De, which is near 1. We derive a linear relationship for the trapped droplet length with Ec^{1/3}, where Ec is the elastocapillary number, by introducing the elastic force into the force balance in addition to the capillary and viscous forces.

Elastic turbulence generates anomalous flow resistance in porous media

Many energy, environmental, industrial, and microfluidic processes rely on the flow of polymer solutions through porous media. Unexpectedly, the macroscopic flow resistance often 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 flow in a transparent 3D porous medium, we demonstrate that this anomalous increase is due to the onset of an elastic instability in which the flow exhibits strong spatiotemporal fluctuations reminiscent of inertial turbulence, despite the small Reynolds number. Our measurements enable us to quantitatively establish that the energy dissipated by pore-scale fluctuations generates the anomalous increase in the overall flow resistance. Because the macroscopic resistance is one of the most fundamental descriptors of fluid flow, our results both help deepen understanding of complex fluid flows and provide guidelines to inform a broad range of applications.

Bifurcations in flows of complex fluids around microfluidic cylinders

Flow around a cylinder is a classical problem in fluid dynamics and also one of the benchmarks for testing viscoelastic flows. The problem is of wide relevance to understanding many microscale industrial and biological processes and applications, such as porous media and mucociliary flows. In recent years, we have developed model microfluidic geometries consisting of very slender cylinders fabricated in glass by selective laser-induced etching. The cylinder radius is small compared with the channel width, which allows the effects of the stagnation points in the flow to dominate over the effects of squeezing between the cylinder and the channel walls. Furthermore, the cylinders are contained in high aspect ratio microchannels that render the flow field approximately two-dimensional (2D) and therefore conveniently permit comparison between experiments and 2D numerical simulations. A number of different viscoelastic fluids including wormlike micellar and various polymer solutions have been tested in our devices. Of particular interest to us has been the occurrence of a striking, steady-in-time, flow asymmetry that occurs for certain non-Newtonian fluids when the dimensionless Weissenberg number (quantifying the importance of elastic over viscous forces in the flow) increases above a critical value. In this perspective review, we present a summary of our key findings related to this novel flow instability and present our current understanding of the mechanism for its onset and growth. We believe that the same fundamental mechanism may also underlie some important non-Newtonian phenomena observed in viscoelastic flows around particles, drops, and bubbles, or through geometries composed of multiple bifurcation points such as cylinder arrays and other porous media. Knowledge of the instability we discuss will be important to consider in the design of optimally functional lab-on-a-chip devices in which viscoelastic fluids are to be used.

Stagnation points control chaotic fluctuations in viscoelastic porous media flow

Viscoelastic flows through porous media become unstable and chaotic beyond critical flow conditions, impacting widespread industrial and biological processes such as enhanced oil recovery and drug delivery. Understanding the influence of the pore structure or geometry on the onset of flow instability can lead to fundamental insights into these processes and, potentially, to their optimization. Recently, for viscoelastic flows through porous media modeled by arrays of microscopic posts, Walkama et al. [D. M. Walkama, N. Waisbord, J. S. Guasto, Phys. Rev. Lett 124, 164501 (2020)] demonstrated that geometric disorder greatly suppressed the strength of the chaotic fluctuations that arose as the flow rate was increased. However, in that work, disorder was only applied to one originally ordered configuration of posts. Here, we demonstrate experimentally that, given a slightly modified ordered array of posts, introducing disorder can also promote chaotic fluctuations. We provide a unifying explanation for these contrasting results by considering the effect of disorder on the occurrence of stagnation points exposed to the flow field, which depends on the nature of the originally ordered post array. This work provides a general understanding of how pore geometry affects the stability of viscoelastic porous media flows.

Flow of Non-Newtonian Fluids in a Single-Cavity Microchannel

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.

Tristability in Viscoelastic Flow Past Side-by-Side Microcylinders

Viscoelastic flows through microscale porous arrays exhibit complex path selection and switching phenomena. However, understanding this process is limited by a lack of studies linking between a single object and large arrays. Here, we report experiments on viscoelastic flow past side-by-side microcylinders with variable intercylinder gap. With increasing flow rate, a sequence of two imperfect symmetry-breaking bifurcations forces selection of either one or two of the three possible flow paths around the cylinders. Tuning the gap length through the value where the first bifurcation becomes perfect reveals regions of bistability and tristability in a dimensionless flow rate-gap length phase diagram.

Active open-loop control of elastic turbulence

We demonstrate through numerical solutions of the Oldroyd-B model in a two-dimensional Taylor-Couette geometry that the onset of elastic turbulence in a viscoelastic fluid can be controlled by imposed shear-rate modulations, one form of active open-loop control. Slow modulations display rich and complex behavior where elastic turbulence is still present, while it vanishes for fast modulations and a laminar response with the Taylor-Couette base flow is recovered. We find that the transition from the laminar to the turbulent state is supercritical and occurs at a critical Deborah number. In the state diagram of both control parameters, Weissenberg versus Deborah number, we identify the region of elastic turbulence. We also quantify the transition by the flow resistance, for which we derive an analytic expression in the laminar regime within the linear Oldroyd-B model. Finally, we provide an approximation for the transition line in the state diagram introducing an effective critical Weissenberg number in comparison to constant shear. Deviations from the numerical result indicate that the physics behind the observed laminar-to-turbulent transition is more complex under time-modulated shear flow.