Morphology, repulsion, and ordering of red blood cells in viscoelastic flows under confinement.

Red blood cells (RBC), the primary carriers of oxygen in the body, play a crucial role across several biomedical applications, while also being an essential model system of a deformable object in the microfluidics and soft matter fields. However, RBC behavior in viscoelastic liquids, which holds promise in enhancing microfluidic diagnostic applications, remains poorly studied. We here show that using viscoelastic polymer solutions as a suspending carrier causes changes in the clustering and shape of flowing RBC in microfluidic flows when compared to a standard Newtonian suspending liquid. Additionally, when the local RBC concentration increases to a point where hydrodynamic interactions take place, we observe the formation of equally-spaced RBC structures, resembling the viscoelasticity-driven ordered particles observed previously in the literature, thus providing the first experimental evidence of viscoelasticity-driven cell ordering. The observed RBC ordering, unaffected by polymer molecular architecture, persists as long as the surrounding medium exhibits shear-thinning, viscoelastic properties. Complementary numerical simulations reveal that viscoelasticity-induced repulsion between RBCs leads to equidistant structures, with shear-thinning modulating this effect. Our results open the way for the development of new biomedical technologies based on the use of viscoelastic liquids while also clarifying fundamental aspects related to multibody hydrodynamic interactions in viscoelastic microfluidic flows.

Effect of fluid elasticity on the emergence of oscillations in an active elastic filament.

Many microorganisms propel themselves through complex media by deforming their flagella. The beat is thought to emerge from interactions between forces of the surrounding fluid, the passive elastic response from deformations of the flagellum and active forces from internal molecular motors. The beat varies in response to changes in the fluid rheology, including elasticity, but there are limited data on how systematic changes in elasticity alter the beat. This work analyses a related problem with fixed-strength driving force: the emergence of beating of an elastic planar filament driven by a follower force at the tip of a viscoelastic fluid. This analysis examines how the onset of oscillations depends on the strength of the force and viscoelastic parameters. Compared to a Newtonian fluid, it takes more force to induce the instability in viscoelastic fluids, and the frequency of the oscillation is higher. The linear analysis predicts that the frequency increases with the fluid relaxation time. Using numerical simulations, the model predictions are compared with experimental data on frequency changes in the bi-flagellated alga Chlamydomonas reinhardtii. The model shows the same trends in response to changes in both fluid viscosity and Deborah number and thus provides a possible mechanistic understanding of the experimental observations.

Exploring the relationship between softness and excess entropy in glass-forming systems.

We explore the relationship between a machine-learned structural quantity (softness) and excess entropy in simulations of supercooled liquids. Excess entropy is known to scale well the dynamical properties of liquids, but this quasi-universal scaling is known to breakdown in supercooled and glassy regimes. Using numerical simulations, we test whether a local form of the excess entropy can lead to predictions similar to those made by softness, such as the strong correlation with particles' tendency to rearrange. In addition, we explore leveraging softness to compute excess entropy in the traditional fashion over softness groupings. Our results show that the excess entropy computed over softness-binned groupings is correlated with activation barriers to rearrangement.

Confinement, chaotic transport, and trapping of active swimmers in time-periodic flows.

Microorganisms encounter complex unsteady flows, including algal blooms in marine settings, microbial infections in airways, and bioreactors for vaccine and biofuel production. Here, we study the transport of active swimmers in two-dimensional time-periodic flows using Langevin simulations and experiments with swimming bacteria. We find that long-term swimmer transport is controlled by two parameters, the pathlength of the unsteady flow and the normalized swimmer speed. The pathlength nonmonotonically controls swimmer dispersion dynamics, giving rise to three distinct dispersion regimes. Weak flows hinder swimmer transport by confining cells toward flow manifolds. As pathlength increases, chaotic transport along flow manifolds initiates, maximizing the number of unique flow cells traveled. Last, strong flows trap swimmers at the vortex core, suppressing dispersal. Experiments with Vibrio cholerae showed qualitative agreement with model dispersion patterns. Our results reveal that nontrivial chaotic transport can arise in simple unsteady flows and suggest a potentially optimal dispersal strategy for microswimmers in nature.

To biofilm or not to biofilm.

A new model helps to predict under which conditions a species of bacteria will switch to a static lifestyle.

Sticking Together: Injectable Granular Hydrogels with Increased Functionality via Dynamic Covalent Inter-Particle Crosslinking.

Granular hydrogels are an exciting class of microporous and injectable biomaterials that are being explored for many biomedical applications, including regenerative medicine, 3D printing, and drug delivery. Granular hydrogels often possess low mechanical moduli and lack structural integrity due to weak physical interactions between microgels. This has been addressed through covalent inter-particle crosslinking; however, covalent crosslinking often occurs through temporal enzymatic methods or photoinitiated reactions, which may limit injectability and material processing. To address this, a hyaluronic acid (HA) granular hydrogel is developed with dynamic covalent (hydrazone) inter-particle crosslinks. Extrusion fragmentation is used to fabricate microgels from photocrosslinkable norbornene-modified HA, additionally modified with either aldehyde or hydrazide groups. Aldehyde and hydrazide-containing microgels are mixed and jammed to form adhesive granular hydrogels. These granular hydrogels possess enhanced mechanical integrity and shape stability over controls due to the covalent inter-particle bonds, while maintaining injectability due to the dynamic hydrazone bonds. The adhesive granular hydrogels are applied to 3D printing, which allows the printing of structures that are stable without any further post-processing. Additionally, the authors demonstrate that adhesive granular hydrogels allow for cell invasion in vitro. Overall, this work demonstrates the use of dynamic covalent inter-particle crosslinking to enhance injectable granular hydrogels.

Bacteria hinder large-scale transport and enhance small-scale mixing in time-periodic flows.

Understanding mixing and transport of passive scalars in active fluids is important to many natural (e.g., algal blooms) and industrial (e.g., biofuel, vaccine production) processes. Here, we study the mixing of a passive scalar (dye) in dilute suspensions of swimming Escherichia coli in experiments using a two-dimensional (2D) time-periodic flow and in a simple simulation. Results show that the presence of bacteria hinders large-scale transport and reduces overall mixing rate. Stretching fields, calculated from experimentally measured velocity fields, show that bacterial activity attenuates fluid stretching and lowers flow chaoticity. Simulations suggest that this attenuation may be attributed to a transient accumulation of bacteria along regions of high stretching. Spatial power spectra and correlation functions of dye-concentration fields show that the transport of scalar variance across scales is also hindered by bacterial activity, resulting in an increase in average size and lifetime of structures. On the other hand, at small scales, activity seems to enhance local mixing. One piece of evidence is that the probability distribution of the spatial concentration gradients is nearly symmetric with a vanishing skewness. Overall, our results show that the coupling between activity and flow can lead to nontrivial effects on mixing and transport.

Aerodynamics-assisted, efficient and scalable kirigami fog collectors.

To address the global water shortage crisis, one of the promising solutions is to collect freshwater from the environmental resources such as fog. However, the efficiency of conventional fog collectors remains low due to the viscous drag of fog-laden wind deflected around the collecting surface. Here, we show that the three-dimensional and centimetric kirigami structures can control the wind flow, forming quasi-stable counter-rotating vortices. The vortices regulate the trajectories of incoming fog clusters and eject extensive droplets to the substrate. As the characteristic structural length is increased to the size of vortices, we greatly reduce the dependence of fog collection on the structural delicacy. Together with gravity-directed gathering by the folds, the kirigami fog collector yields a collection efficiency of 16.1% at a low wind speed of 0.8 m/s and is robust against surface characteristics. The collection efficiency is maintained even on a 1 m2 collector in an outdoor setting.

The perpetual fragility of creeping hillslopes.

Soil creeps imperceptibly but relentlessly downhill, shaping landscapes and the human and ecological communities that live within them. What causes this granular material to 'flow' at angles well below repose? The unchallenged dogma is churning of soil by (bio)physical disturbances. Here we experimentally render slow creep dynamics down to micron scale, in a laboratory hillslope where disturbances can be tuned. Surprisingly, we find that even an undisturbed sandpile creeps indefinitely, with rates and styles comparable to natural hillslopes. Creep progressively slows as the initially fragile pile relaxes into a lower energy state. This slowing can be enhanced or reversed with different imposed disturbances. Our observations suggest a new model for soil as a creeping glass, wherein environmental disturbances maintain soil in a perpetually fragile state.

Crystalline shielding mitigates structural rearrangement and localizes memory in jammed systems under oscillatory shear.

The nature of yield in amorphous materials under stress has yet to be fully elucidated. In particular, understanding how microscopic rearrangement gives rise to macroscopic structural and rheological signatures in disordered systems is vital for the prediction and characterization of yield and the study of how memory is stored in disordered materials. Here, we investigate the evolution of local structural homogeneity on an individual particle level in amorphous jammed two-dimensional (athermal) systems under oscillatory shear and relate this evolution to rearrangement, memory, and macroscale rheological measurements. We define the structural metric crystalline shielding, and show that it is predictive of rearrangement propensity and structural volatility of individual particles under shear. We use this metric to identify localized regions of the system in which the material's memory of its preparation is preserved. Our results contribute to a growing understanding of how local structure relates to dynamic response and memory in disordered systems.

Bacterial activity hinders particle sedimentation.

Sedimentation in active fluids has come into focus due to the ubiquity of swimming micro-organisms in natural and industrial processes. Here, we investigate sedimentation dynamics of passive particles in a fluid as a function of bacteria E. coli concentration. Results show that the presence of swimming bacteria significantly reduces the speed of the sedimentation front even in the dilute regime, in which the sedimentation speed is expected to be independent of particle concentration. Furthermore, bacteria increase the dispersion of the passive particles, which determines the width of the sedimentation front. For short times, particle sedimentation speed has a linear dependence on bacterial concentration. Mean square displacement data shows, however, that bacterial activity decays over long experimental (sedimentation) times. An advection-diffusion equation coupled to bacteria population dynamics seems to capture concentration profiles relatively well. A single parameter, the ratio of single particle speed to the bacteria flow speed can be used to predict front sedimentation speed.

A perfusion phantom for ASL MRI based on impinging jets.

PURPOSE: We present a novel perfusion phantom for validation of arterial spin labeled (ASL) perfusion MRI methods and protocols. METHODS: Impinging jets, driven by a peristaltic pump, were used to achieve perfusion-like mixing of magnetically labeled inflowing fluid within a perfusion compartment. The phantom was validated by varying pump rates and obtaining ASL-MRI data at multiple postlabeling delays using a pseudo-continuous ASL sequence with a 3D stack-of-spirals readout. An additional data set was acquired using a pseudo-continuous ASL sequence with a 2D EPI readout. Phantom sensitivity to pseudo-continuous ASL labeling efficiency was also tested. RESULTS: Fluid dynamics simulations predicted that maximum mixing would occur near the central axis of the perfusion compartment. Experimentally observed signal changes within this region were reproducible and well fit by the standard Buxton general kinetic model. Simulations and experimental data showed no label outflow from the perfusion chamber and calculated perfusion rates, averaged over the entire phantom volume, agreed with the expected volumetric flow rates provided by the flow pump. Phantom sensitivity to pseudo-continuous ASL labeling parameters was also demonstrated. CONCLUSION: Perfusion-like signal can be simulated using impinging jets to create a well-mixed compartment. Observed perfusion and transit time values were reproducible and within the physiological range for brain perfusion. This phantom design has a broad range of potential applications in both basic and clinical research involving ASL MRI.

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.

Scaling of relaxation and excess entropy in plastically deformed amorphous solids.

When stressed sufficiently, solid materials yield and deform plastically via reorganization of microscopic constituents. Indeed, it is possible to alter the microstructure of materials by judicious application of stress, an empirical process utilized in practice to enhance the mechanical properties of metals. Understanding the interdependence of plastic flow and microscopic structure in these nonequilibrium states, however, remains a major challenge. Here, we experimentally investigate this relationship, between the relaxation dynamics and microscopic structure of disordered colloidal solids during plastic deformation. We apply oscillatory shear to solid colloidal monolayers and study their particle trajectories as a function of shear rate in the plastic regime. Under these circumstances, the strain rate, the relaxation rate associated with plastic flow, and the sample microscopic structure oscillate together, but with different phases. Interestingly, the experiments reveal that the relaxation rate associated with plastic flow at time t is correlated with the strain rate and sample microscopic structure measured at earlier and later times, respectively. The relaxation rate, in this nonstationary condition, exhibits power-law, shear-thinning behavior and scales exponentially with sample excess entropy. Thus, measurement of sample static structure (excess entropy) provides insight about both strain rate and constituent rearrangement dynamics in the sample at earlier times.

Quantification of plasticity via particle dynamics above and below yield in a 2D jammed suspension.

The yield transition of amorphous materials is characterized by a swift increase of energy dissipation. The connection between particle dynamics, dissipation, and overall material rheology, however, has still not been elucidated. Here, we take a new approach relating trajectories to yielding, using a custom built interfacial stress rheometer, which allows for measurement of shear moduli (G',G'') of a dense athermal suspension's microstructure while simultaneously tracking particle trajectories undergoing cyclic shear. We find an increase in total area traced by particle trajectories as the system is stressed well below to well above yield. Trajectories may be placed into three categories: reversibly elastic paths; reversibly plastic paths, associated with smooth limit cycles; and irreversibly plastic paths, in which particles do not return to their original position. We find that above yield, reversibly plastic trajectories are predominantly found near to the shearing surface, whereas reversibly elastic paths are more prominent near the stationary wall. This spatial transition between particles acting as liquids to those acting as solids is characteristic of a 'melting front', which is observed to shift closer to the wall with increasing strain. We introduce a non-dimensional measure of plastic dissipation based on particle trajectories that scales linearly with strain amplitude both above and below yield, and that is unity at the rheological yield point. Surprisingly, this relation collapses for three systems of varying degrees of disorder.

Active Reversible Swimming of Magnetically Assembled "Microscallops" in Non-Newtonian Fluids.

Miniaturized devices capable of active swimming at low Reynolds numbers are of fundamental importance and possess potential biomedical utility. The design of colloidal microswimmers requires not only miniaturizing reconfigurable structures but also understanding their interactions with media at low Reynolds numbers. We investigate the dynamics of "microscallops" made of asymmetric magnetic cubes, which are assembled and actuated using magnetic fields. One approach to achieving directional propulsion is to break the symmetry of the viscous forces by coupling the reciprocal motions of such microswimmers with the nonlinear rheology inherent in non-Newtonian fluids. When placed in shear-thinning fluids, the local viscosity gradient resulting from nonuniform shear stresses exerted by time-asymmetric strokes of the microscallops generates propulsive thrust through an effect we term "self-viscophoresis". Surprisingly, we found that the direction of propulsion changes with the size and structure of these assemblies. We analyze the origins of their directional propulsion and explain the variable propulsion direction in terms of multiple counterbalancing domains of shear dissipation around the microscale structures. The principles governing the locomotion of these microswimmers may be extended to other reconfigurable microbots assembled from colloidal-scale units.

Formation of stable aggregates by fluid-assembled solid bridges.

When a colloidal suspension is dried, capillary pressure may overwhelm repulsive electrostatic forces, assembling aggregates that are out of thermal equilibrium. This poorly understood process confers cohesive strength to many geological and industrial materials. Here we observe evaporation-driven aggregation of natural and synthesized particulates, probe their stability under rewetting, and measure bonding strength using an atomic force microscope. Cohesion arises at a common length scale (∼5 µm), where interparticle attractive forces exceed particle weight. In polydisperse mixtures, smaller particles condense within shrinking capillary bridges to build stabilizing "solid bridges" among larger grains. This dynamic repeats across scales, forming remarkably strong, hierarchical clusters, whose cohesion derives from grain size rather than mineralogy. These results may help toward understanding the strength and erodibility of natural soils, and other polydisperse particulates that experience transient hydrodynamic forces.

Flow Resistance and Structures in Viscoelastic Channel Flows at Low Re.

The flow of viscoelastic fluids in channels and pipes remains poorly understood, particularly at low Reynolds numbers. Here, we investigate the flow of polymeric solutions in straight channels using pressure measurements and particle tracking. The flow friction factor f_{η} versus flow rate exhibits two regimes: a transitional regime marked by rapid increase in drag, and a turbulentlike regime characterized by a sudden decrease in drag and a weak dependence on flow rate. Lagrangian trajectories show finite transverse modulations not seen in Newtonian fluids. These curvature perturbations far downstream can generate sufficient hoop stresses to sustain the flow instabilities in the parallel shear flow.

Quenching active swarms: effects of light exposure on collective motility in swarming Serratia marcescens.

Swarming colonies of the light-responsive bacteria Serratia marcescens grown on agar exhibit robust fluctuating large-scale flows that include arrayed vortices, jets and sinuous streamers. We study the immobilization and quenching of these collective flows when the moving swarm is exposed to intense wide-spectrum light with a substantial ultraviolet component. We map the emergent response of the swarm to light in terms of two parameters-light intensity and duration of exposure-and identify the conditions under which collective motility is impacted. For small exposure times and/or low intensities, we find collective motility to be negligibly affected. Increasing exposure times and/or intensity to higher values suppresses collective motility but only temporarily. Terminating exposure allows bacteria to recover and eventually reestablish collective flows similar to that seen in unexposed swarms. For long exposure times or at high intensities, exposed bacteria become paralysed and form aligned, jammed regions where macroscopic speeds reduce to zero. The effective size of the quenched region increases with time and saturates to approximately the extent of the illuminated region. Post-exposure, active bacteria dislodge immotile bacteria; initial dissolution rates are strongly dependent on duration of exposure. Based on our experimental observations, we propose a minimal Brownian dynamics model to examine the escape of exposed bacteria from the region of exposure. Our results complement studies on planktonic bacteria, inform models of patterning in gradated illumination and provide a starting point for the study of specific wavelengths on swarming bacteria.

Upstream vortex and elastic wave in the viscoelastic flow around a confined cylinder.

The viscoelastic flow past a cylinder is a classic benchmark problem that is not completely understood. Using novel 3D holographic particle velocimetry, we report three main discoveries of the elastic instability upstream of a single cylinder in viscoelastic channel flow. First, we observe that upstream vortices initiate at the corner between the cylinder and the wall and grow with increasing flow rate. Second, beyond a critical Weissenberg, the flow upstream becomes unsteady and switches between two bistable configurations, leading to symmetry breaking in the cylinder axis direction that is highly three-dimensional in nature. Lastly, we find that the disturbance of the elastic instability propagates relatively far upstream via an elastic wave, and is weakly correlated with that in the cylinder wake. The wave speed and the extent of the instability increase with Weissenberg number, indicating an absolute instability in viscoelastic fluids.