Search

This study presents an experimental framework for large amplitude oscillatory extension (LAOE) to investigate nonlinear material properties of complex fluids. Using a microfluidic optimized shape cross-slot extensional rheometer, we generate approximately homogeneous planar extensional flows driven by programmable syringe pumps operating in oscillatory or pulsatile sinusoidal modes. Micro-particle image velocimetry and simultaneous pressure drop measurements are employed to analyze the time-dependent flow field and elastic stress response. For Newtonian fluids, a linear relationship between the applied strain rate and pressure drop is observed across a wide range of oscillation amplitudes and frequencies. In contrast, dilute polymer solutions exhibit significant deviations, with excess pressure drops and divergence between average strain rates along extension and compression axes during the LAOE cycle. By spanning a broad range of Weissenberg and Deborah numbers, we identify unique Lissajous curves and critical conditions for the onset of nonlinearities under oscillatory extension. Numerical simulations, assuming homogeneous flow, underpin the experimental findings, validating the robustness of our microfluidic approach. This study demonstrates the utility of oscillatory extensional flows for probing the nonlinear rheological behavior of soft materials, offering quantitative insights into their extensional properties under nonlinear flow conditions.

Elastic turbulence can lead to to increased flow resistance, mixing and heat transfer. Its control -- either suppression or promotion -- has significant potential, and there is a concerted ongoing effort by the community to improve our understanding. Here we explore the dynamics of uncertainty in elastic turbulence, inspired by an approach recently applied to inertial turbulence in Ge et al. (2023) \textit{J. Fluid Mech.} 977:A17. We derive equations for the evolution of uncertainty measures, yielding insight on uncertainty growth mechanisms. Through numerical experiments, we identify four regimes of uncertainty evolution, characterised by I) rapid transfer to large scales, with large scale growth rates of $\tau^{6}$ (where $\tau$ represents time), II) a dissipative reduction of uncertainty, III) exponential growth at all scales, and IV) saturation. These regimes are governed by the interplay between advective and polymeric contributions (which tend to amplify uncertainty), viscous, relaxation and dissipation effects (which reduce uncertainty), and inertial contributions. In elastic turbulence, reducing Reynolds number increases uncertainty at short times, but does not significantly influence the growth of uncertainty at later times. At late times, the growth of uncertainty increases with Weissenberg number, with decreasing polymeric diffusivity, and with the logarithm of the maximum length scale, as large flow features adjust the balance of advective and relaxation effects. These findings provide insight into the dynamics of elastic turbulence, offering a new approach for the analysis of viscoelastic flow instabilities., Comment: submitted to JFM

The FENE-P and sPTT viscoelastic models are widely used for modelling of complex fluids. Although they are derived from distinct micro-structural theories, these models can become mathematically identical in steady and homogeneous flows with a particular choice of the value of the model parameters. However, even with this choice of parameter values, the model responses are known to differ from each other in transient flows. In this work, we investigate the responses of the FENE-P and sPTT constitutive models in Large Amplitude Oscillatory Shear (LAOS). In steady-shear, the shear stress scales with the non-dimensional group $Wi/(aL) \ (Wi\sqrt{\epsilon})$ for the FENE-P (sPTT) model, where $Wi$ is the Weissenberg number, $L^2$ is the limit of extensibility in the FENE-P model ($a$ being $L^2/(L^2-3)$) and $\epsilon$ is the extensibility parameter in the sPTT model. Our numerical and analytical results show that, in LAOS, the FENE-P model only shows this universality for large values of $L^2$ whereas the sPTT model shows this universality for all values of$\epsilon$. In the strongly non-linear region, there is a drastic difference between the responses of the two models, with the FENE-P model exhibiting strong shear stress overshoots which manifest as self-intersecting secondary loops in the viscous Lissajous curves. We quantify the non-linearity exhibited by each constitutive model using the Sequence of Physical Processes framework. Despite the high degree of non-linearity exhibited by the FENE-P model, we also show using fully non-linear 1D simulations that it does not shear band in LAOS within the range of conditions studied., Comment: Submitted to Journal of Fluid Mechanics

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

In this work we consider theoretically the problem of a Newtonian droplet moving in an otherwise quiescent infinite viscoelastic fluid under the influence of an externally applied temperature gradient. The outer fluid is modelled by the Oldroyd-B equation, and the problem is solved for small Weissenberg and Capillary numbers in terms of a double perturbation expansion. We assume microgravity conditions and neglect the convective transport of energy and momentum. We derive expressions for the droplet migration speed and its shape in terms of the properties of both fluids. In the absence of shape deformation, the droplet speed decreases monotonically for sufficiently viscous inner fluids, while for fluids with a smaller inner-to-outer viscosity ratio, the droplet speed first increases and then decreases as a function of the Weissenberg number. For small but finite values of the Capillary number, the droplet speed behaves monotonically as a function of the applied temperature gradient for a fixed ratio of the Capillary and Weissenberg numbers. We demonstrate that this behaviour is related to the polymeric stresses deforming the droplet in the direction of its migration, while the associated changes in its speed are Newtonian in nature, being related to a change in the droplet's hydrodynamic resistance and its internal temperature distribution. When compared to the results of numerical simulations, our theory exhibits a good predictive power for sufficiently small values of the Capillary and Weissenberg numbers., Comment: 28 pages, 7 figures, accepted for publication in Journal of Non-Newtonian Fluid Mechanics

Secondary flows are ubiquitous in channel flows, where small velocity components perpendicular to the main velocity appear due to the complexity of the channel geometry and/or that of the flow itself such as from inertial or non-Newtonian effects, etc. We investigate here the inertialess secondary flow of viscoelastic fluids in curved microchannels of rectangular cross-section and constant but alternating curvature: the so-called "serpentine channel" geometry. Numerical calculations (Poole et al, 2013) have shown that in this geometry, in the absence of elastic instabilities, a steady secondary flow develops that takes the shape of two counter-rotating vortices in the plane of the channel cross-section. We present the first experimental visualization evidence and characterization of these steady secondary flows, using a complementarity of microPIV in the plane of the channel, and confocal visualisation of dye-stream transport in the cross-sectional plane. We show that the measured streamlines and the relative velocity magnitude of the secondary flows are in qualitative agreement with the numerical results. In addition to our techniques being broadly applicable to the characterisation of three-dimensional flow structures in microchannels, our results are important for understanding the onset of instability in serpentine viscoelastic flows., Comment: This is a post-peer-review, pre-copyedit version of an article published in Microfluidics and Nanofluidics. The final authenticated version is available online at: http://dx.doi.org/10.1007/s10404-019-2195-0

Peritrichous bacteria swim in viscous fluids by rotating multiple helical flagellar filaments. As the bacterium swims forward, all its flagella rotate in synchrony behind the cell in a helical bundle. When the bacterium changes its direction, the flagellar filaments unbundle and reorient the cell for a short period of time before returning to their bundled state and resuming swimming. This rapid bundling and unbundling is, at its heart, a mechanical process whereby hydrodynamic interactions balance with elasticity to determine the time-varying deformation of the filaments. Inspired by this biophysical problem, we present in this paper what is perhaps the simplest model of bundling whereby two, or more, straight elastic filaments immersed in a viscous fluid rotate about their centreline, inducing rotational flows which tend to bend the filaments around each other. We derive an integro-differential equation governing the shape of the filaments resulting from mechanical balance in a viscous fluid at low Reynolds number. We show that such equation may be evaluated asymptotically analytically in the long-wavelength limit, leading to a local partial differential equation governed by a single dimensionless Bundling number. A numerical study of the dynamics predicted by the model reveals the presence of two configuration instabilities with increasing Bundling numbers: first to a crossing state where filaments touch at one point and then to a bundled state where filaments wrap along each other in a helical fashion. We next provide an intuitive physical model for the crossing instability and show that it may be used to predict analytically its threshold, and adapted to address the transition to a bundling state. We then use a macro-scale experimental implementation of the two-filament configuration in order to validate our theoretical predictions and obtain excellent agreement.

Experimental measurements and numerical simulations are made on fluid flow through cross-slot devices with a range of aspect (depth:width) ratios, 0.4 < alpha < 3.87. For low Reynolds numbers Re, the flow is symmetric and a sharp boundary exists between fluid streams entering the cross-slot from opposite directions. Above an alpha-dependent critical value Re_c, the flow undergoes a symmetry-breaking bifurcation (though remains steady and laminar) and a spiral vortex structure develops about the central axis of the outflow channel. An order parameter characterizing the instability grows according to a sixth-order Landau potential, and shows a progression from second order to first order transitions as alpha increases. A tricritical point occurs for alpha ~ 0.55. The spiral vortex acts as a mixing region in the flow field and this phenomenon can be used to drive enhanced mixing in microfluidic devices.

We propose a novel device capable of measuring the relaxation time of viscoelastic fluids as small as 1\,ms. In contrast to most rheometers, which by their very nature are concerned with producing viscometric or nearly-viscometric flows, here we make use of an elastic instability which occurs in the flow of viscoelastic fluids with curved streamlines. To calibrate the rheometer we combine simple scaling arguments with relaxation times obtained from first normal-stress difference data measured in a classical shear rheometer. As an additional check we also compare these relaxation times to those obtained from Zimm theory and good agreement is observed. Once calibrated, we show how the serpentine rheometer can be used to access smaller polymer concentrations and lower solvent viscosities where classical measurements become difficult or impossible to use due to inertial and/or resolution limitations. In the absence of calibration the serpentine channel can still be a very useful comparative or index device., Comment: accepted for for publication in Lab on a chip