Purely-elastic flow instabilities and elastic turbulence in microfluidic cross-slot devices

Comparison of viscoelastic flows in two- and three-dimensional serpentine channels

Polymer solutions in the dilute regime play a significant role in industrial applications. Due to the intricate rheological properties of these highly viscoelastic fluids, especially in complex flow geometries, a thorough numerical analysis of their flow dynamics is imperative. In this research, we present a numerical investigation of purely elastic instability occurring in two- and three-dimensional serpentine channels under conditions where fluid inertia is negligible and across a broad spectrum of polymer relaxation times. Our findings reveal a strong qualitative agreement between the existing experimental results obtained from dilute solutions of flexible polymers in microfluidic devices and the numerical simulations conducted in two and three dimensions using the Oldroyd-B model. Spatial flow observations and statistical analysis of temporal flow features indicate that this purely elastic turbulent flow exhibits nonhomogeneous, non-Gaussian, and anisotropic characteristics across all scales. Additionally, our comparison of two- and three-dimensional simulations demonstrates that the elastic instability is primarily driven by the curvature of the streamlines induced by the flow geometry, rather than the weak secondary flow in the azimuthal direction. Therefore, our two-dimensional numerical simulations successfully replicate, at least qualitatively, the features observed in three-dimensional experiments. Furthermore, spectral analysis suggests that, in comparison to elastic turbulence in the dilute regime, the range of scales for the excited fluctuations is narrower.

Rheological effects on purely-elastic flow asymmetries in the cross-slot geometry

Viscoelastic flows in the cross-slot geometry can undergo a transition from a steady symmetric to a steady asymmetric flow state, ostensibly due to purely-elastic effects arising beyond a critical flow rate, or Weissenberg number Wi. However, some reports suggest that shear thinning of the fluid's viscosity may also play an important role in this transition. We employ a series of polymer solutions of varying rheological properties to investigate in detail how the interplay between fluid elasticity and shear thinning affects the onset and development of asymmetric flows in the cross-slot. Flow velocimetry is performed on each of the polymer solutions, and is used to assess the degree of flow asymmetry I in the cross-slot as a function of both Wi and a dimensionless parameter S quantifying the flow-rate-dependent extent of shear thinning. Typically, the flow field breaks symmetry as Wi is increased beyond a critical value, but the magnitude of I is found to also be dependent on S. For a few specific polymer solutions, the flow field recovers symmetry above a second, higher critical Wi as S becomes small. The experimental results are summarized in a flow state diagram in Wi-S space, showing the relationship between flow asymmetry and fluid rheology. Finally, to gain a deeper understanding of the effects of shear thinning, numerical simulations are performed using the linear simplified Phan-Thien-Tanner model. We demonstrate that the degree of both shear thinning and elasticity of the fluid, and their interplay, are important factors controlling elastic instabilities in the cross-slot geometry.

Nonlinear Phenomena in Microfluidics

This review focuses on experimental work on nonlinear phenomena in microfluidics, which for the most part are phenomena for which the velocity of a fluid flowing through a microfluidic channel does not scale proportionately with the pressure drop. Examples include oscillations, flow-switching behaviors, and bifurcations. These phenomena are qualitatively distinct from laminar, diffusion-limited flows that are often associated with microfluidics. We explore the nonlinear behaviors of bubbles or droplets when they travel alone or in trains through a microfluidic network or when they assemble into either one- or two-dimensional crystals. We consider the nonlinearities that can be induced by the geometry of channels, such as their curvature or the bas-relief patterning of their base. By casting posts, barriers, or membranes─situated inside channels─from stimuli-responsive or flexible materials, the shape, size, or configuration of these elements can be altered by flowing fluids, which may enable autonomous flow control. We also highlight some of the nonlinearities that arise from operating devices at intermediate Reynolds numbers or from using non-Newtonian fluids or liquid metals. We include a brief discussion of relevant practical applications, including flow gating, mixing, and particle separations.

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.

Statistical properties of two-dimensional elastic turbulence

We numerically investigate the spatial and temporal statistical properties of a dilute polymer solution in the elastic turbulence regime, i.e., in the chaotic flow state occurring at vanishing Reynolds and high Weissenberg numbers. We aim at elucidating the relations between measurements of flow properties performed in the spatial domain with the ones taken in the temporal domain, which is a key point for the interpretation of experimental results on elastic turbulence and to discuss the validity of Taylor's hypothesis. To this end, we carry out extensive direct numerical simulations of the two-dimensional Kolmogorov flow of an Oldroyd-B viscoelastic fluid. Static pointlike numerical probes are placed at different locations in the flow, particularly at the extrema of mean flow amplitude. The results in the fully developed elastic turbulence regime reveal large velocity fluctuations, as compared to the mean flow, leading to a partial breakdown of Taylor's frozen-field hypothesis. While second-order statistics, probed by spectra and structure functions, display consistent scaling behaviors in the spatial and temporal domains, the third-order statistics highlight robust differences. In particular the temporal analysis fails to capture the skewness of streamwise longitudinal velocity increments. Finally, we assess both the degree of statistical inhomogeneity and isotropy of the flow turbulent fluctuations as a function of scale. While the system is only weakly nonhomogenous in the cross-stream direction, it is found to be highly anisotropic at all scales.

Flow Direction-Dependent Elastic Instability in a Symmetry-Breaking Microchannel

This paper reports flow direction-dependent elastic instability in a symmetry-breaking microchannel. The microchannel consisted of a square chamber and a nozzle structure. A viscoelastic polyacrylamide solution was used for the instability demonstration. The instability was realized as the viscoelastic flow became asymmetric and unsteady in the microchannel when the flow exceeded a critical Weissenberg number. The critical Weissenberg number was found to be different for the forward-directed flow and the backward-directed flow in the microchannel.

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.

Nonlinear microfluidics: device physics, functions, and applications

The microfluidic flow is typically laminar due to the dominant viscous effects. At Reynolds numbers far below 1 (Re ≪ 1), the fluid inertia can be neglected. For the steady flow of incompressible Newtonian fluids, it approaches linear Stokes flow. At intermediate Re, there exists a weak-inertia flow regime where secondary flows such as Dean vortices are accessible for microfluidic manipulations. Apart from the fluid inertia, other nonlinear factors such as the non-Newtonian fluid properties, concurrent flow of dissimilar fluids, compliant fluidic structures and stimuli-responsive materials can also cause intriguing flow behaviours. Through proper designs, they can be applied for a variety of microfluidic components including mixers, valves, oscillators, stabilizers and auto-regulators etc., greatly enriching the microfluidic flow control and manipulation strategies. Due to its unique working characteristics and advantages, nonlinear microfluidics has increasingly attracted extensive attention. This review presents a systematic survey on this subject. The designs of typical nonlinear microfluidic devices, their working mechanisms, key applications, and the perspective of their future developments will be discussed. The nonlinear microfluidic techniques are believed to play an essential role in the next generation of highly-integrated, automated, and intelligent microfluidics.

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.

Three-dimensional structures and symmetry breaking in viscoelastic cross-channel flow

Using holographic particle tracking, we report the three-dimensional flow structure organizing the viscoelastic instability in cross-channel flow. Beyond a critical Wi, the advective core flow undergoes an out-of-plane instability marked by the emergence of tertiary flow, resembling that of the toroidal vortices in Taylor-Couette geometry. The out-of-plane flow component distorts the separatrix between the impinging inflow streams, triggering symmetry breaking normal to the extension plane. As extensional rate increases, progressively higher order modes of the separatrix are observed, akin to Euler buckling of a rigid column. The disturbances propagate upstream via stress fluctuations despite viscous dissipation. These complex flow structures may be generic to elastic turbulence in mixed flows.

Signature of elastic turbulence of viscoelastic fluid flow in a single pore throat

When a viscoelastic fluid, such as an aqueous polymer solution, flows through a porous medium, the fluid undergoes a repetitive expansion and contraction as it passes from one pore to the next. Above a critical flow rate, the interaction between the viscoelastic nature of the polymer and the pore configuration results in spatial and temporal flow instabilities reminiscent of turbulentlike behavior, even though the Reynolds number Re≪1. To investigate whether this is caused by many repeated pore body-pore throat sequences, or simply a consequence of the converging (diverging) nature present in a single pore throat, we performed experiments using anionic hydrolyzed polyacrylamide (HPAM) in a microfluidic flow geometry representing a single pore throat. This allows the viscoelastic fluid to be characterized at increasing flow rates using microparticle image velocimetry in combination with pressure drop measurements. The key finding is that the effect, popularly known as "elastic turbulence," occurs already in a single pore throat geometry. The critical Deborah number at which the transition in rheological flow behavior from pseudoplastic (shear thinning) to dilatant (shear thickening) strongly depends on the ionic strength, the type of cation in the anionic HPAM solution, and the nature of pore configuration. The transition towards the elastic turbulence regime was found to directly correlate with an increase in normal stresses. The topology parameter, Q_{f}, computed from the velocity distribution, suggests that the "shear thickening" regime, where much of the elastic turbulence occurs in a single pore throat, is a consequence of viscoelastic normal stresses that cause a complex flow field. This flow field consists of extensional, shear, and rotational features around the constriction, as well as upstream and downstream of the constriction. Furthermore, this elastic turbulence regime, has high-pressure fluctuations, with a power-law decay exponent of up to |-2.1| which is higher than the Kolmogorov value for turbulence of |-5/3|.

Oscillations of a cantilevered micro beam driven by a viscoelastic flow instability

The interaction of flexible structures with viscoelastic flows can result in very rich dynamics. In this paper, we present the results of the interactions between the flow of a viscoelastic polymer solution and a cantilevered beam in a confined microfluidic geometry. Cantilevered beams with varying length and flexibility were studied. With increasing flow rate and Weissenberg number, the flow transitioned from a fore-aft symmetric flow to a stable detached vortex upstream of the beam, to a time-dependent unstable vortex shedding. The shedding of the unstable vortex upstream of the beam imposed a time-dependent drag force on the cantilevered beam resulting in flow-induced beam oscillations. The oscillations of the flexible beam were classified into two distinct regimes: a regime with a clear single vortex shedding from upstream of the beam resulting in a sinusoidal beam oscillation pattern with the frequency of oscillation increasing monotonically with Weissenberg number, and a regime at high Weissenberg numbers characterized by 3D viscoelastic instabilities where the frequency of oscillations plateaued. The critical onset of the flow transitions, the mechanism of vortex shedding and the dynamics of the cantilevered beam response are presented in detail here as a function of beam flexibility and flow viscoelasticity.

Numerical Study on the Characteristics of Boger Type Viscoelastic Fluid Flow in a Micro Cross-Slot under Sinusoidal Stimulation

The cross-slot geometry plays an important role in the study of nonlinear effects of viscoelastic fluids. The flow of viscoelastic fluid in a micro cross-slot with a high channel aspect ratio (AR, the ratio of channel depth to width) can be divided into three types, which are symmetric flow, steady-state asymmetric flow and time-dependent flow under the inlet condition with a constant velocity. However, the flow pattern of a viscoelastic fluid in the cross-slot when a stimulation is applied at inlets has been rarely reported. In this paper, the response of cross-slot flow under an external sinusoidal stimulation is studied by numerical simulations of a two-dimensional model representing the geometry with a maximum limit of AR. For the cases under constant inlet velocity conditions, three different flow patterns occur successively with the increase of Weissenberg number (Wi). For the cases under sinusoidal varying inlet velocity conditions, when the stimulation frequency is far away from the natural frequency of a viscoelastic fluid, the frequency spectrum of velocity fluctuation field shows the characteristics of a fundamental frequency and several harmonics. However, the harmonic frequency disappears when the stimulation frequency is close to the natural frequency of the viscoelastic fluid. Besides, the flow pattern shows spatial symmetry and changes with time. In conclusion, the external stimulation has an effect on the flow pattern of viscoelastic fluid in the 2D micro cross-slot channel, and a resonance occurs when the stimulation frequency is close to the natural frequency of the fluid.

Flow of wormlike micellar solutions around microfluidic cylinders with high aspect ratio and low blockage ratio

We employ time-resolved flow velocimetry and birefringence imaging methods to study the flow of a well-characterized shear-banding wormlike micellar solution around a novel glass-fabricated microfluidic circular cylinder. In contrast with typical microfluidic cylinders, our geometry is characterized by a high aspect ratio α = H/W = 5 and a low blockage ratio β = 2r/W = 0.1, where H and W are the channel height and width, and the cylinder radius r = 20 μm. The small cylinder radius allows access up to very high Weissenberg numbers 1.9 ≤ Wi = λMU/r ≤ 3750 (where λM is the Maxwell relaxation time) while inertial effects remain entirely negligible (Reynolds number, Re < 10-4). At low Wi values, the flow remains steady and symmetric and a birefringent region (indicating micellar alignment and tensile stress) develops downstream of the cylinder. Above a critical value Wic ≈ 60 the flow transitions to a steady asymmetric state, characterized as a supercritical pitchfork bifurcation, in which the fluid takes a preferential path around one side of the cylinder. At a second critical value Wic2 ≈ 130, the flow becomes time-dependent, with a characteristic frequency f0 ≈ 1/λM. This initial transition to time dependence has characteristics of a subcritical Hopf bifurcation. Power spectra of the measured fluctuations become complex as Wi is increased further, showing a gradual slowing down of the dynamics and emergence of harmonics. A final transition at very high Wic3 corresponds to the re-emergence of a single peak in the power spectrum but at much higher frequency. We discuss this in terms of possible flow-induced breakage of micelles into shorter species with a faster relaxation time.

Particle-laden two-dimensional elastic turbulence

The aggregation properties of heavy inertial particles in the elastic turbulence regime of an Oldroyd-B fluid with periodic Kolmogorov mean flow are investigated by means of extensive numerical simulations in two dimensions. Both the small- and large-scale features of the resulting inhomogeneous particle distribution are examined, focusing on their connection with the properties of the advecting viscoelastic flow. We find that particles preferentially accumulate on thin highly elastic propagating structures and that this effect is the largest for intermediate values of particle inertia. We provide a quantitative characterization of this phenomenon that allows to relate it to the accumulation of particles in filamentary highly strained flow regions producing clusters of correlation dimension close to 1. At larger scales, particles are found to undergo turbophoretic-like segregation. Indeed, our results indicate a close relationship between the profiles of particle density and fluid velocity fluctuations. The large-scale inhomogeneity of the particle distribution is interpreted in the framework of a model derived in the limit of small, but finite, particle inertia. The qualitative characteristics of different observables are, to a good extent, independent of the flow elasticity. When increased, the latter is found, however, to slightly reduce the globally averaged degree of turbophoretic unmixing.