Hydrodynamics control shear-induced pattern formation in attractive suspensions

Structure-property relations in rheology of cellulose nanofibrils-based hydrogels

Hydrogels prepared from self-assembled cellulose nanofibrils (CNFs) are widely used in biomedicine, electronics and environmental technology. Their ability to serve as inks for extrusion-based 3D printing is conventionally evaluated by means of rheological tests. A model is developed that describes the response of CNF gels in small- and large-amplitude oscillatory tests in a unified manner. The model involves a reasonably small number of material parameters, ensures good agreement between results of simulation and observations in oscillatory tests and correctly predicts the stress-strain Lissajous curves, experimental data in hysteresis loop tests, and measurements of the steady-state viscosity. The model is applied to analyze how composition and preparation conditions for CNF gels affect transition from shear thinning to weak strain overshoot in large-amplitude shear oscillatory tests. Based on the model, simple relations are derived for the fractal dimension of CNF clusters and the storage modulus of gels prepared in aqueous solutions of multivalent salts. The validity of these equations is confirmed by comparison of their predictions with observations in independent tests.

Small variations in particle-level interactions lead to large structural heterogeneities in colloidal gels

Colloidal gels typically exhibit mechanical properties akin to a viscoelastic solid, influenced by their underlying particulate network. Hence, the structural and morphological characteristics of the colloidal network have a significant effect on the rigidity of the gel. In this study, we show how seemingly small variations in the particle-level interactions throughout the system result in larger scale structural heterogeneities. While the microscale particle level descriptors of the colloidal network remain largely unaffected by heterogeneous interactions, larger scale properties of a colloidal gel change appreciably. The overall cluster-level mesostructure of a colloidal gel is found to be sensitive to the small variations in the interaction potential at the particle level.

Predicting yielding in attractive colloidal gels

One of the defining characteristics of soft glassy materials is their ability to exhibit a yield stress, which can result in an overall elasto-visco-plastic mechanics. To design soft materials with specific properties, it is essential to gain a comprehensive understanding of the topological and structural failure points that occur during yielding. However, predicting these failure points, which lead to yielding, is challenging due to the dynamic nature of structure development and its cooccurrence with other complicated processes, such as local rearrangements and anisotropy. In this study, we employ a series of tools from network science to investigate colloidal gels as a model for soft glassy materials during yielding. Our findings reveal that edge betweenness centrality can be utilized as a universal predictor for yielding across various state variables, including the volume fraction of solids, the strength, and the range of attraction between colloids.

Effect of Shear History on Solid-Liquid Transition of Particulate Gel Fuels

Investigating the structural evolution of particulate gels is a very challenging task due to their vulnerability and true flow characteristics. In this work, deeper insight into the rheological properties of gel fuels filled with fumed silica (FS) and aluminum microparticles (Al MPs) was gained by changing shear procedures. Firstly, the flow curves were found to no longer follow the monotonic power law and exhibited subtle thixotropic responses. As the shear rate increased, the gel structure underwent a transition from local shear to bulk shear in the nonlinear region after yielding. This finding reveals the prevalence of nonideal local shear in industry. Secondly, the time-dependent rheological responses demonstrated that the strength spectrum of gel fuels depends on the applied shear rate, with stress relaxation more easily observed at lower shear rates. Those results involved the structural disruption, recovery, and equilibrium of particulate gels from two scales of shear rate and shear time.

Nonequilibrium interactions between multi-scale colloids regulate the suspension microstructure and rheology

Understanding nonequilibrium interactions of multi-component colloidal suspensions is critical for many dynamical settings such as self-assembly and material processing. A key question is how the nonequilibrium distributions of individual components influence the effective interparticle interactions and flow behavior. In this work, we develop a first-principle framework to study a bidisperse suspension of colloids and depletants using a Smoluchowski equation and corroborated by Brownian dynamics (BD) simulations. Using nonlinear microrheology as a case study, we demonstrate that effective depletion interactions between driven colloids are sensitive to particle timescales out of equilibrium and cannot be predicted by equilibrium-based pair potentials like Asakura-Oosawa. Furthermore, we show that the interplay between Brownian relaxation timescales of different species plays a critical role in governing the viscosity of multi-component suspensions. Our model highlights the limitations of using equilibrium pair potentials to approximate interparticle interactions in nonequilibrium processes such as hydrodynamic flows and presents a useful framework for studying the transport of driven, interacting suspensions.

Shear thickening in presence of adhesive contact forces: The singularity of cornstarch

HYPOTHESIS

A number of dense particle suspensions experience a dramatic increase in viscosity with the shear stress, up to a solid-like response. This shear-thickening process is understood as a transition under flow of the nature of the contacts - from lubricated to frictional - between initially repellent particles. Most systems are now assumed to fit in with this scenario, which is questionable.

EXPERIMENT

Using an in-house pressure sensor array, we provide a spatio-temporal map of the normal stresses in the flows of two shear-thickening fluids: a stabilized calcium carbonate suspension, known to fit in with the standard scenario, and a cornstarch suspension, which spectacular thickening behavior remains poorly understood.

FINDINGS

We evidence in cornstarch a unique, stable heterogeneous structure, which moves in the velocity direction and does not appear in calcium carbonate. Its nature changes from a stress wave to a rolling solid jammed aggregate at high solid fraction and small gap width. The modeling of these heterogenities points to an adhesive force between cornstarch particles at high stress, also evidenced in microscopic measurements. Cornstarch being also attractive at low stress, it stands out of the classical shear-thickening frame, and might be part of a larger family of adhesive and attractive shear-thickening fluids.

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.

Surface instabilities generated by a slider pulled across a granular bed

We report an instability of a slider slowly dragged at the surface of a granular bed in a quasistatic regime. The boat-shaped slider sits on the granular medium under its own weight and is free to translate vertically and to rotate around the pitch axis while a constant horizontal speed is imposed. For a wide range of parameters (mass, length, shape, velocity) a regular pattern of peaks and troughs spontaneously emerges as the slider travels forward. This instability is studied through experiments using a conveyor belt and by means of two-dimensional discrete elements method simulations. We show that the wavelength and amplitude of the pattern scale as the length of the slider. We also observe that the ripples disappear for low and high masses, indicating an optimal confining pressure. The effect of the shape, more specifically the inclination of the front spatula, is studied and found to drastically influence both the wavelength and the amplitude. Finally, we show that the mechanical details (friction, cohesion) of the contact point between the slider and the pulling device is critical and remains to be fully understood.

Topological origins of yielding in short-ranged weakly attractive colloidal gels

Yielding of the particulate network in colloidal gels under applied deformation is accompanied by various microstructural changes, including rearrangement, bond rupture, anisotropy, and reformation of secondary structures. While much work has been done to understand the physical underpinnings of yielding in colloidal gels, its topological origins remain poorly understood. Here, employing a series of tools from network science, we characterize the bonds using their orientation and network centrality. We find that bonds with higher centralities in the network are ruptured the most at all applied deformation rates. This suggests that a network analysis of the particulate structure can be used to predict the failure points in colloidal gels a priori.

Shear induced tuning and memory effects in colloidal gels of rods and spheres

Shear history plays an important role in determining the linear and nonlinear rheological response of colloidal gels and can be used for tuning their structure and flow properties. Increasing the colloidal particle aspect ratio lowers the critical volume fraction for gelation due to an increase in the particle excluded volume. Using a combination of rheology and confocal microscopy, we investigate the effect of steady and oscillatory preshear history on the structure and rheology of colloidal gels formed by silica spheres and rods of length L and diameter D (L/D = 10) dispersed in 11 M CsCl solution. We use a non-dimensional Mason number, Mn (=F_visc./F_attr.), to compare the effect of steady and oscillatory preshear on gel viscoelasticity. We show that after preshearing at intermediate Mn, attractive sphere gel exhibits strengthening, whereas attractive rod gel exhibits weakening. Rheo-imaging of gels of attractive rods shows that at intermediate Mn, oscillatory preshear induces large compact rod clusters in the gel microstructure, compared to steady preshear. Our study highlights the impact of particle shape on gel structuring under flow and viscoelasticity after shear cessation.

Shear jamming and fragility in fractal suspensions under confinement

Under applied stress, the viscosity of many dense particulate suspensions increases drastically, a response known as discontinuous shear-thickening (DST). In some cases, the applied stress can even transform the suspension into a solid-like shear jammed state. Although shear jamming (SJ) has been probed for dense suspensions with particles having well-defined shapes, such a phenomenon for fractal objects has not been explored. Here, using rheology and in situ optical imaging, we study the flow behaviour of ultra-dilute fractal suspensions of multi-walled carbon nanotubes (MWCNT) under confinement. We show a direct transition from flowing to SJ state without a precursory DST in fractal suspensions at an onset volume fraction, ϕ ∼ 0.5%, significantly lower than that of conventional dense suspensions (ϕ ∼ 55%). The ultra-low concentration enables us to demonstrate the fragility and associated contact dynamics of the SJ state, which remain experimentally unexplored in suspensions. Furthermore, using a generalized Wyart-Cates model, we propose a generic phase diagram for fractal suspensions that captures the possibility of SJ without prior DST over a wide range of shear stress and volume fractions.

Flow-induced "waltzing" red blood cells: Microstructural reorganization and the corresponding rheological response

We investigate flow-induced structural organization in a dilute suspension of tumbling red blood cells (RBCs) under confined shear flow. For small Reynolds (Re = 0.1) and capillary numbers (Ca), with fully coupled hydrodynamic interaction (HI) and without interparticle adhesion, we find that HI between the biconcave discoid particles prompts the formation of layered RBC chains and synchronized rotating RBC pairs, referred here as "waltzing doublets." As the volume fraction ϕ increases, more waltzing doublets appear in RBC files. Stronger shear stress disrupts structural arrangements at higher Ca. We find that the flow-induced organization of waltzing doublets changes how the suspension viscosity varies with ϕ qualitatively. The intrinsic viscosity is particularly sensitive to microstructural rearrangement, increasing (decreasing) with ϕ at low (high) Ca that correlates with the change in the fraction of doublets. We verified flow-induced collective motion with comparison to two-cell simulations in which the cell volume fraction is controlled by varying the domain volume.

Brownian dynamics simulations of shear-induced aggregation of charged colloidal particles in the presence of hydrodynamic interactions

HYPOTHESIS

In spite of the abundant literature on Brownian simulations of the aggregation behavior of colloidal suspensions both under quiescent conditions and in the presence of shear, few works performed simulations including the effect of hydrodynamic interactions. Even fewer works have investigated the effects of shear on the aggregation of electrostatically-stabilized colloidal suspensions.

SIMULATIONS

In this work, we employed Brownian dynamics simulations implementing the Rotne-Prager-Yamakawa approximation to account for hydrodynamic interactions and investigated the aggregation kinetics of electrostatically-stabilized colloidal suspensions exposed to simple shear, for various Péclet number values, particle volume fractions and surface potential values.

RESULTS

The increase in Péclet number (i.e., in the shear rate), leads to an overall increase in the aggregation rate and the formation of large aggregates that, for sufficiently high volume fractions, rapidly grow, leading to either breakup and restructuring phenomena or percolation of the system. In some cases, a bimodal distribution of the cluster population was observed. Our simulations further indicate that at the highest Péclet, the aggregation dynamics is independent of the energy barrier and entirely controlled by shear. A comparison with a simple BD method reveals that neglecting long-range hydrodynamic interactions leads to a substantial underestimation of the aggregation rate.

Flow-Switched Bistability in a Colloidal Gel with Non-Brownian Grains

We show that mixing a colloidal gel with larger, non-Brownian grains generates novel flow-switched bistability. Using a combination of confocal microscopy and rheology, we find that prolonged moderate shear results in liquefaction by collapsing the gel into disjoint globules, whereas fast shear gives rise to a yield-stress gel with granular inclusions upon flow cessation. We map out the state diagram of this new "mechanorheological material" with varying granular content and demonstrate that its behavior is also found in separate mixture using different particles and solvents.

X-ray computerized tomography observation of Lycopodium paste incorporating memory of shaking

In a uniform layer consisting of a mixture of granular material and liquid, it is known that desiccation cracks exhibit various anisotropic patterns that depend on the nature of the shaking that the layer experienced before drying. The existence of this effect implies that information regarding the direction of shaking is retained as a kind of memory in the arrangements of granular particles. In this work we make measurements in paste composed of Lycopodium powder using microfocus x-ray computerized tomography (μCT) in order to investigate the three-dimensional arrangements of particles. We find shaking-induced anisotropic arrangements of neighboring particles and density fluctuations forming interstices mainly in the lower part of the layer. We compare the observed properties of these arrangements with numerical results obtained in the study of a model of non-Brownian particles under shear deformation. In the experimental system, we also observe crack tips in the μCT images and confirm that these cracks grow along interstices in the direction perpendicular to the initial shaking.

Tuning local microstructure of colloidal gels by ultrasound-activated deformable inclusions

Colloidal gels possess a memory of previous shear events, both steady and oscillatory. This memory, embedded in the microstructure, affects the mechanical response of the gel, and therefore enables precise tuning of the material properties under careful preparation. Here we demonstrate how the dynamics of a deformable inclusion, namely a bubble, can be used to locally tune the microstructure of a colloidal gel. We examine two different phenomena of bubble dynamics that apply a local strain to the surrounding material: dissolution due to gas diffusion, with a characteristic strain rate of ∼10^-3 s^-1; and volumetric oscillations driven by ultrasound, with a characteristic frequency of ∼10⁴ s^-1. We characterise experimentally the microstructure of a model colloidal gel around bubbles in a Hele-Shaw geometry using confocal microscopy and particle tracking. In bubble dissolution experiments, we observe the formation of a pocket of solvent next to the bubble surface, but marginal changes to the microstructure. In experiments with ultrasound-induced bubble oscillations, we observe a striking rearrangement of the gel particles into a microstructure with increased local ordering. High-speed bright-field microscopy reveals the occurrence of both high-frequency bubble oscillations and steady microstreaming flow; both are expected to contribute to the emergence of the local order in the microstructure. These observations open the way to local tuning of colloidal gels based on deformable inclusions controlled by external pressure fields.

Nanoparticle Assembly as a Materials Development Tool

Nanoparticle assembly is a complex and versatile method of generating new materials, capable of using thousands of different combinations of particle size, shape, composition, and ligand chemistry to generate a library of unique structures. Here, a history of particle self-assembly as a strategy for materials discovery is presented, focusing on key advances in both synthesis and measurement of emergent properties to describe the current state of the field. Several key challenges for further advancement of nanoparticle assembly are also outlined, establishing a roadmap of critical research areas to enable the next generation of nanoparticle-based materials synthesis.

Shear-induced migration of confined flexible fibers

We report an experimental study of the shear-induced migration of flexible fibers in suspensions confined between two parallel plates. Non-Brownian fiber suspensions are imaged in a rheo-microscopy setup, where the top and the bottom plates counter-rotate and create a Couette flow. Initially, the fibers are near the bottom plate due to sedimentation. Under shear, the fibers move with the flow and migrate towards the center plane between the two walls. Statistical properties of the fibers, such as the mean values of the positions, orientations, and end-to-end lengths of the fibers, are used to characterize the behaviors of the fibers. A dimensionless parameter Λ_eff, which compares the hydrodynamic shear stress and the fiber stiffness, is used to analyze the effective flexibility of the fibers. The observations show that the fibers that are more likely to bend exhibit faster migration. As Λ_eff increases (softer fibers and stronger shear stresses), the fibers tend to align in the flow direction and the motions of the fibers transition from tumbling and rolling to bending. The bending fibers drift away from the walls to the center plane. Further increasing Λ_eff leads to more coiled fiber shapes, and the bending is more frequent and with larger magnitudes, which leads to more rapid migration towards the center. Different behaviors of the fibers are quantified with Λ_eff, and the structures and the dynamics of the fibers are correlated with the migration.

New Model to Predict the Concentration-Dependent Viscosity of Monoclonal Antibody Solutions

Concentrated monoclonal antibody solutions exhibit high solution viscosity, which is experimentally measured to be ∼1-2 orders of magnitude higher than the viscosity of water. However, physical processes responsible for the high antibody viscosity are not fully understood. We show that fluid occlusion due to the trapped solvent molecules within the boundaries formed by the aggregated antibodies is responsible for the elevated solution viscosity. We develop a theory to predict the viscosity of monoclonal antibodies based on the geometry of the antibody molecule and the aggregate morphology. We validate our theory with experiments and highlight useful insights obtained from the viscosity equation which could help in controlling the drug viscosity at the molecular design stage itself.

Migration and Morphology of Colloidal Gel Clusters in Cylindrical Channel Flow

We report the cluster-level structural parameters of colloidal thermogelling nanoemulsions in channel flow as a function of attractive interactions and local shear stress. The spatiotemporal evolution of the gel microstructure is obtained by directly visualizing the dispersed phase near the edge of a cylindrical channel. We observe the flow of the nanoemulsion gels in a range of radial positions (r) and shear stresses between 70 and 220 Pa, finding that the r-dependent cluster sizes are due to a balance between shear forces that yield bonds and attractive interactions that rebuild the inter-colloid bonds. In addition, the largest clusters appear to be affected by confinement and accumulate toward the central axis of the channel, resulting in a volume fraction gradient. Cluster size and volume fraction variabilities are most prominent when the attractive interactions are the strongest. Specifically, a distinct transition from sparse, fluidized clusters near the walls to concentrated, large clusters toward the center is observed. These two structural states coincide with a velocity-based transition from higher shear rates near the walls to lower shear rates toward the center of the channel. We find a compounding effect where larger gel clusters, formed under strong attractions and low shear stresses, are susceptible to shear-induced migration that intensifies r-dependent heterogeneity and deviations in the flow behavior from predictive models.

Life and death of colloidal bonds control the rate-dependent rheology of gels

Colloidal gels exhibit rich rheological responses under flowing conditions. A clear understanding of the coupling between the kinetics of the formation/rupture of colloidal bonds and the rheological response of attractive gels is lacking. In particular, for gels under different flow regimes, the correlation between the complex rheological response, the bond kinetics, microscopic forces, and an overall micromechanistic view is missing in previous works. Here, we report the bond dynamics in short-range attractive particles, microscopically measured stresses on individual particles and the spatiotemporal evolution of the colloidal structures in different flow regimes. The interplay between interparticle attraction and hydrodynamic stresses is found to be the key to unraveling the physical underpinnings of colloidal gel rheology. Attractive stresses, mostly originating from older bonds dominate the response at low Mason number (the ratio of shearing to attractive forces) while hydrodynamic stresses tend to control the rheology at higher Mason numbers, mostly arising from short-lived bonds. Finally, we present visual mapping of particle bond numbers, their life times and their borne stresses under different flow regimes.

Shear driven vorticity aligned flocs in a suspension of attractive rigid rods

A combination of rheology, optical microscopy and computer simulations was used to investigate the microstructural changes of a semi-dilute suspension of attractive rigid rods in an imposed shear flow. The aim is to understand the relation of the microstructure with the viscoelastic response, and the yielding and flow behaviour in different shear regimes of gels built from rodlike colloids. A semi-dilute suspension of micron sized, rodlike silica particles suspended in 11 M CsCl salt solution was used as a model system for attractive rods' gel. Upon application of steady shear the gel microstructure rearranges in different states and exhibits flow instabilities depending on shear rate, attraction strength, volume fraction and geometrical confinement. At low rod volume fractions, the suspension forms large, vorticity aligned, particle rich flocs that roll in the flow-vorticity plane, an effect that is due to an interplay between hydrodynamic interactions and geometrical confinement as suggested by computer simulations. Experimental data allow the creation of a state diagram, as a function of volume fraction and shear rates, identifying regimes of stable (or unstable) floc formation and of homogeneous gel or broken clusters. The transition is related to dimensionless Mason number, defined as the ratio of shear forces to interparticle attractive force.

Effect of Particle Specific Surface Area on the Rheology of Non-Brownian Silica Suspensions

Industrial formulations very often involve particles with a broad range of surface characteristics and size distributions. Particle surface asperities (roughness) and porosity increase particle specific surface area and significantly alter suspension rheology, which can be detrimental to the quality of the end product. We examine the rheological properties of two types of non-Brownian, commercial precipitated silicas, with varying specific surface area, namely PS52 and PS226, suspended in a non-aqueous solvent, glycerol, and compare them against those of glass sphere suspensions (GS2) with a similar size distribution. A non-monotonic effect of the specific surface area (S) on suspension rheology is observed, whereby PS52 particles in glycerol are found to exhibit strong shear thinning response, whereas such response is suppressed for glass sphere and PS226 particle suspensions. This behaviour is attributed to the competing mechanisms of particle-particle and particle-solvent interactions. In particular, increasing the specific surface area beyond a certain value results in the repulsive interparticle hydration forces (solvation forces) induced by glycerol overcoming particle frictional contacts and suppressing shear thinning; this is evidenced by the response of the highest specific surface area particles PS226. The study demonstrates the potential of using particle specific surface area as a means to tune the rheology of non-Brownian silica particle suspensions.

Temperature-Induced Aggregation in Portlandite Suspensions

Temperature is well known to affect the aggregation behavior of colloidal suspensions. This paper elucidates the temperature dependence of the rheology of portlandite (calcium hydroxide: Ca(OH)₂) suspensions that feature a high ionic strength and a pH close to the particle's isoelectric point. In contrast to the viscosity of the suspending medium (saturated solution of Ca(OH)₂ in water), the viscosity of Ca(OH)₂ suspensions is found to increase with elevating temperature. This behavior is shown to arise from the temperature-induced aggregation of polydisperse Ca(OH)₂ particulates because of the diminution of electrostatic repulsive forces with increasing temperature. The temperature dependence of the suspension viscosity is further shown to diminish with increasing particle volume fraction as a result of volumetric crowding and the formation of denser fractal structures in the suspension. Significantly, the temperature-dependent rheological response of suspensions is shown to be strongly affected by the suspending medium's properties, including ionic strength and ion valence, which affect aggregation kinetics. These outcomes provide new insights into aggregation processes that affect the temperature-dependent rheology of portlandite-based and similar suspensions that feature strong charge screening behavior.

Time-dependent shear rate inhomogeneities and shear bands in a thixotropic yield-stress fluid under transient shear

We study the rheological responses and shear-rate inhomogeneities and shear banding behaviors of a thixotropic fumed silica suspension in shear startup tests and flow reversal tests. We find that this suspension under transient shear exhibits not only viscoelasticity, yielding, kinematic hardening, and thixotropy, but also time-dependent shear inhomogeneities including bands when the apparent shear rate is below a critical value between 0.1 and 0.25 s-1. Through multiple shear startup tests and flow reversal tests, we find that thixotropy promotes flow heterogeneity while kinematic hardening suppresses it. We propose a simple thixo-plastic constitutive equation that can qualitatively predict the important features of the rheological response and banding dynamics in shear startup tests and flow reversal tests.