Dynamic Buckling of a Filament Impacted by a Falling Droplet.

We investigate the buckling dynamics of an elastic filament impacted axially by a falling liquid droplet, and identify the buckling modes through a combination of experimental and theoretical analyses. A phase diagram is constructed on a plane defined by two primary parameters: the falling height and the filament length. Two transition boundaries are observed, with one marking the occurrence of dynamic buckling and the other separating the buckling regime into two distinct modes. Notably, the hydrodynamic viscous force of the liquid dominates during the impact, with the dynamic buckling instability predicted by a single elastoviscous number. The critical load is twice the critical static load, which is, however, lower for the deformable droplet utilized in our study, as compared to a solid object. An additional time-dependent simulation on a longer filament exhibits a higher buckling mode, succeeded by a more distinct coarsening process than our experimental observations.

Symmetry Breaking of Tail-Clamped Filaments in Stokes Flow.

Symmetry breaking (SB) of fluid-structure interaction problems plays an important role in our understanding of animals' locomotive and sensing behaviors. In this Letter, we study the SB of flexible filaments clamped at one end and placed in a spanwise periodic array in Stokes flow. The equilibrium state of the filament along the streamwise direction loses stability and experiences two-dimensional and then three-dimensional SBs as the spanwise distance increases, or as the filament rigidity reduces. For slightly deformed filaments, the viscous and pressure forces are commensurate, while for extremely deformed filaments the viscous force becomes dominant.

Dispersed Fibers Change the Classical Energy Budget of Turbulence via Nonlocal Transfer.

The backreaction of dispersed rigid fibers to turbulence is analyzed by means of a state-of-the-art fully coupled immersed boundary method. The following universal scenario is identified: turbulence at large scales looses a consistent part of its kinetic energy (via a Darcy friction term), which partially reappears at small scales where a new range of energy-containing scales does emerge. Large-scale mixing is thus depleted in favor of a new mixing mechanism arising at the smallest scales. Anchored fibers cause the same backreaction to turbulence as moving fibers of large inertia. Our results thus provide a link between two apparently separated realms: the one of porous media and the one of suspension dynamics.

Numerical simulations of vorticity banding of emulsions in shear flows.

Multiphase shear flows often show banded structures that affect the global behavior of complex fluids e.g. in microdevices. Here we investigate numerically the banding of emulsions, i.e. the formation of regions of high and low volume fractions, alternated in the vorticity direction and aligned with the flow (shear bands). These bands are associated with a decrease of the effective viscosity of the system. To understand the mechanism of experimentally observed banding, we have performed interface-resolved simulations of the two-fluid system. The experiments were performed starting with a random distribution of droplets, which under the applied shear, evolve in time resulting in a phase separation. To numerically reproduce this process, the banded structures are initialized in a narrow channel confined by two walls moving in opposite directions. We find that the initial banded distribution is stable when droplets are free to merge and unstable when coalescence is prevented. In this case, additionally, the effective viscosity of the system increases, resembling the rheological behavior of suspensions of deformable particles. Droplet coalescence, on the other hand, allows emulsions to reduce the total surface of the system and, hence, the energy dissipation associated with the deformation, which in turn reduces the effective viscosity.

The breakdown of Darcy's law in a soft porous material.

We perform direct numerical simulations of the flow through a model of deformable porous medium. Our model is a two-dimensional hexagonal lattice, with defects, of soft elastic cylindrical pillars, with elastic shear modulus G, immersed in a liquid. We use a two-phase approach: the liquid phase is a viscous fluid and the solid phase is modeled as an incompressible viscoelastic material, whose complete nonlinear structural response is considered. We observe that the Darcy flux (q) is a nonlinear function - steeper than linear - of the pressure-difference (ΔP) across the medium. Furthermore, the flux is larger for a softer medium (smaller G). We construct a theory of this super-linear behavior by modelling the channels between the solid cylinders as elastic channels whose walls are made of material with a linear constitutive relation but can undergo large deformation. Our theory further predicts that the flow permeability is an universal function of ΔP/G, which is confirmed by the present simulations.

Flow-assisted droplet assembly in a 3D microfluidic channel.

Self-assembly of soft matter, such as droplets or colloids, has become a promising scheme to engineer novel materials, model living matter, and explore non-equilibrium statistical mechanics. In this article, we present detailed numerical simulations of few non-Brownian droplets in various flow conditions, specifically, focusing on their self-assembly within a short distance in a three-dimensional (3D) microfluidic channel, cf. [Shen et al., Adv. Sci., 2016, 3(6), 1600012]. Contrary to quasi two-dimensional (q2D) systems, where dipolar interaction is the key mechanism for droplet rearrangement, droplets in 3D confinement produce much less disturbance to the underlying flow, thus experiencing weaker dipolar interactions. Using confined simple shear and Poiseuille flows as reference flows, we show that the droplet dynamics is mostly affected by the shear-induced cross-stream migration, which favors chain structures if the droplets are under an attractive depletion force. For more compact clusters, such as three droplets in a triangular shape, our results suggest that an inhomogeneous cross-sectional inflow profile is further required. Overall, the accelerated self-assembly of a small-size droplet cluster results from the combined effects of strong depletion forces, confinement-mediated shear alignments, and fine-tuned inflow conditions. The deterministic nature of the flow-assisted self-assembly implies the possibility of large throughputs, though calibration of all different effects to directly produce large droplet crystals is generally difficult.

Flexible Fiber Reveals the Two-Point Statistical Properties of Turbulence.

We study the dynamics of a flexible fiber freely moving in a three-dimensional fully developed turbulent field and present a phenomenological theory to describe the interaction between the fiber elasticity and the turbulent flow. This theory leads to the identification of two distinct regimes of flapping, which we validate against direct numerical simulations fully resolving the fiber dynamics. The main result of our analysis is the identification of a flapping regime where the fiber, despite its elasticity, is slaved to the turbulent fluctuations. In this regime the fiber can be used to measure two-point statistical observables of turbulence, including scaling exponents of velocity structure functions, the sign of the energy cascade and the energy flux of turbulence, as well as the characteristic times of the eddies within the inertial range of scales. Our results are expected to have a deep impact on the experimental turbulence research as a new way, accurate and efficient, to measure two-point, and more generally multipoint, statistics of turbulence.

Universal Scaling Laws for Dense Particle Suspensions in Turbulent Wall-Bounded Flows.

The macroscopic behavior of dense suspensions of neutrally buoyant spheres in turbulent plane channel flow is examined. We show that particles larger than the smallest turbulence scales cause the suspension to deviate from the continuum limit in which its dynamics is well described by an effective suspension viscosity. This deviation is caused by the formation of a particle layer close to the wall with significant slip velocity. By assuming two distinct transport mechanisms in the near-wall layer and the turbulence in the bulk, we define an effective wall location such that the flow in the bulk can still be accurately described by an effective suspension viscosity. We thus propose scaling laws for the mean velocity profile of the suspension flow, together with a master equation able to predict the increase in drag as a function of the particle size and volume fraction.

Rheology of Confined Non-Brownian Suspensions.

We study the rheology of confined suspensions of neutrally buoyant rigid monodisperse spheres in plane-Couette flow using direct numerical simulations. We find that if the width of the channel is a (small) integer multiple of the sphere diameter, the spheres self-organize into two-dimensional layers that slide on each other and the effective viscosity of the suspension is significantly reduced. Each two-dimensional layer is found to be structurally liquidlike but its dynamics is frozen in time.

Continuous Growth of Droplet Size Variance due to Condensation in Turbulent Clouds.

We use a stochastic model and direct numerical simulation to study the impact of turbulence on cloud droplet growth by condensation. We show that the variance of the droplet size distribution increases in time as t^{1/2}, with growth rate proportional to the large-to-small turbulent scale separation and to the turbulence integral scales but independent of the mean turbulent dissipation. Direct numerical simulations confirm this result and produce realistically broad droplet size spectra over time intervals of 20 min, comparable with the time of rain formation.

Motion of an elastic capsule in a constricted microchannel.

We study the motion of an elastic capsule through a microchannel characterized by a localized constriction. We consider a capsule with a stress-free spherical shape and impose its steady-state configuration in an infinitely long straight channel as the initial condition for our calculations. We report how the capsule deformation, velocity, retention time, and maximum stress of the membrane are affected by the capillary number, Ca , and the constriction shape. We estimate the deformation by measuring the variation of the three-dimensional surface area and a series of alternative quantities easier to extract from experiments. These are the Taylor parameter, the perimeter and the area of the capsule in the spanwise plane. We find that the perimeter is the quantity that best reproduces the behavior of the three-dimensional surface area. This is maximum at the centre of the constriction and shows a second peak after it, whose location depends on the Ca number. We observe that, in general, area-deformation-correlated quantities grow linearly with Ca , while velocity-correlated quantities saturate for large Ca but display a steeper increase for small Ca . The velocity of the capsule divided by the velocity of the flow displays, surprisingly, two different qualitative behaviors for small and large capillary numbers. Finally, we report that longer constrictions and spanwise wall bounded (versus spanwise periodic) domains cause larger deformations and velocities. If the deformation and velocity in the spanwise wall bounded domains are rescaled by the initial equilibrium deformation and velocity, their behavior is undistinguishable from that in a periodic domain. In contrast, a remarkably different behavior is reported in sinusoidally shaped and smoothed rectangular constrictions indicating that the capsule dynamics is particularly sensitive to abrupt changes in the cross section. In a smoothed rectangular constriction larger deformations and velocities occur over a larger distance.

Laminar, turbulent, and inertial shear-thickening regimes in channel flow of neutrally buoyant particle suspensions.

The aim of this Letter is to characterize the flow regimes of suspensions of finite-size rigid particles in a viscous fluid at finite inertia. We explore the system behavior as a function of the particle volume fraction and the Reynolds number (the ratio of flow and particle inertia to viscous forces). Unlike single-phase flows, where a clear distinction exists between the laminar and the turbulent states, three different regimes can be identified in the presence of a particulate phase, with smooth transitions between them. At low volume fractions, the flow becomes turbulent when increasing the Reynolds number, transitioning from the laminar regime dominated by viscous forces to the turbulent regime characterized by enhanced momentum transport by turbulent eddies. At larger volume fractions, we identify a new regime characterized by an even larger increase of the wall friction. The wall friction increases with the Reynolds number (inertial effects) while the turbulent transport is weakly affected, as in a state of intense inertial shear thickening. This state may prevent the transition to a fully turbulent regime at arbitrary high speed of the flow.

A microfluidic device to sort capsules by deformability: a numerical study.

Guided by extensive numerical simulations, we propose a microfluidic device that can sort elastic capsules by their deformability. The device consists of a duct embedded with a semi-cylindrical obstacle, and a diffuser which further enhances the sorting capability. We demonstrate that the device can operate reasonably well under changes in the initial position of the capsule. The efficiency of the device remains essentially unaltered under small changes of the obstacle shape (from semi-circular to semi-elliptic cross-section). Confinement along the direction perpendicular to the plane of the device increases its efficiency. This work is the first numerical study of cell sorting by a realistic microfluidic device.

Rotational propulsion enabled by inertia.

The fluid mechanics of small-scale locomotion has recently attracted considerable attention, due to its importance in cell motility and the design of artificial micro-swimmers for biomedical applications. Most studies on the topic consider the ideal limit of zero Reynolds number. In this paper, we investigate a simple propulsion mechanism --an up-down asymmetric dumbbell rotating about its axis of symmetry-- unable to propel in the absence of inertia in a Newtonian fluid. Inertial forces lead to continuous propulsion for all finite values of the Reynolds number. We study computationally its propulsive characteristics as well as analytically in the small-Reynolds-number limit. We also derive the optimal dumbbell geometry. The direction of propulsion enabled by inertia is opposite to that induced by viscoelasticity.

Elasto-inertial focusing and particle migration in high aspect ratio microchannels for high-throughput separation.

The combination of flow elasticity and inertia has emerged as a viable tool for focusing and manipulating particles using microfluidics. Although there is considerable interest in the field of elasto-inertial microfluidics owing to its potential applications, research on particle focusing has been mostly limited to low Reynolds numbers (Re<1), and particle migration toward equilibrium positions has not been extensively examined. In this work, we thoroughly studied particle focusing on the dynamic range of flow rates and particle migration using straight microchannels with a single inlet high aspect ratio. We initially explored several parameters that had an impact on particle focusing, such as the particle size, channel dimensions, concentration of viscoelastic fluid, and flow rate. Our experimental work covered a wide range of dimensionless numbers (0.05 < Reynolds number < 85, 1.5 < Weissenberg number < 3800, 5 < Elasticity number < 470) using 3, 5, 7, and 10 µm particles. Our results showed that the particle size played a dominant role, and by tuning the parameters, particle focusing could be achieved at Reynolds numbers ranging from 0.2 (1 µL/min) to 85 (250 µL/min). Furthermore, we numerically and experimentally studied particle migration and reported differential particle migration for high-resolution separations of 5 µm, 7 µm and 10 µm particles in a sheathless flow at a throughput of 150 µL/min. Our work elucidates the complex particle transport in elasto-inertial flows and has great potential for the development of high-throughput and high-resolution particle separation for biomedical and environmental applications.

Shear thickening in non-Brownian suspensions: an excluded volume effect.

Shear thickening appears as an increase of the viscosity of a dense suspension with the shear rate, sometimes sudden and violent at high volume fraction. Its origin for noncolloidal suspension with non-negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate that fluid inertia causes a strong microstructure anisotropy that results in the formation of a shadow region with no relative flux of particles. We show that shear thickening at finite inertia can be explained as an increase of the effective volume fraction when considering the dynamically excluded volume due to these shadow regions.

A dual resolution phase-field solver for wetting of viscoelastic droplets.

We present a new and efficient phase-field solver for viscoelastic fluids with moving contact line based on a dual-resolution strategy. The interface between two immiscible fluids is tracked by using the Cahn-Hilliard phase-field model, and the viscoelasticity incorporated into the phase-field framework. The main challenge of this approach is to have enough resolution at the interface to approach the sharp-interface methods. The method presented here addresses this problem by solving the phase field variable on a mesh twice as fine as that used for the velocities, pressure, and polymer-stress constitutive equations. The method is based on second-order finite differences for the discretization of the fully coupled Navier-Stokes, polymeric constitutive, and Cahn-Hilliard equations, and it is implemented in a 2D pencil-like domain decomposition to benefit from existing highly scalable parallel algorithms. An FFT-based solver is used for the Helmholtz and Poisson equations with different global sizes. A splitting method is used to impose the dynamic contact angle boundary conditions in the case of large density and viscosity ratios. The implementation is validated against experimental data and previous numerical studies in 2D and 3D. The results indicate that the dual-resolution approach produces nearly identical results while saving computational time for both Newtonian and viscoelastic flows in 3D.

Nonequilibrium thermodynamics and the optimal path to turbulence in shear flows.

We determine the initial condition on the laminar-turbulent boundary closest to the laminar state using nonlinear optimization for plane Couette flow. Resorting to the general evolution criterion of nonequilibrium systems we optimize the route to the statistically steady turbulent state, i.e., the state characterized by the largest entropy production. This is the first time information from the fully turbulent state is included in the optimization procedure. We demonstrate that the optimal initial condition is localized in space for realistic flow domains.

Yield-stress fluids in porous media: a comparison of viscoplastic and elastoviscoplastic flows.

A numerical and theoretical study of yield-stress fluid flows in two types of model porous media is presented. We focus on viscoplastic and elastoviscoplastic flows to reveal some differences and similarities between these two classes of flows. Small elastic effects increase the pressure drop and also the size of unyielded regions in the flow which is the consequence of different stress solutions compare to viscoplastic flows. Yet, the velocity fields in the viscoplastic and elastoviscoplastic flows are comparable for small elastic effects. By increasing the yield stress, the difference in the pressure drops between the two classes of flows becomes smaller and smaller for both considered geometries. When the elastic effects increase, the elastoviscoplastic flow becomes time-dependent and some oscillations in the flow can be observed. Focusing on the regime of very large yield stress effects in the viscoplastic flow, we address in detail the interesting limit of 'flow/no flow': yield-stress fluids can resist small imposed pressure gradients and remain quiescent. The critical pressure gradient which should be exceeded to guarantee a continuous flow in the porous media will be reported. Finally, we propose a theoretical framework for studying the 'yield limit' in the porous media.

Integral representation of channel flow with interacting particles.

We construct a boundary integral representation for the low-Reynolds-number flow in a channel in the presence of freely suspended particles (or droplets) of arbitrary size and shape. We demonstrate that lubrication theory holds away from the particles at horizontal distances exceeding the channel height and derive a multipole expansion of the flow which is dipolar to the leading approximation. We show that the dipole moment of an arbitrary particle is a weighted integral of the stress and the flow at the particle surface, which can be determined numerically. We introduce the equation of motion that describes hydrodynamic interactions between arbitrary, possibly different, distant particles, with interactions determined by the product of the mobility matrix and the dipole moment. Further, the problem of three identical interacting spheres initially aligned in the streamwise direction is considered and the experimentally observed "pair exchange" phenomenon is derived analytically and confirmed numerically. For nonaligned particles, we demonstrate the formation of a configuration with one particle separating from a stable pair. Our results suggest that in a dilute initially homogenous particulate suspension flowing in a channel the particles will eventually separate into singlets and pairs.