Dramatic effect of fluid chemistry on cornstarch suspensions: Linking particle interactions to macroscopic rheology

Stress and flow inhomogeneity in shear-thickening suspensions

HYPOTHESIS

The viscosity of dense suspensions surges when the applied stress surpasses a material-specific critical threshold. There is growing evidence that the thickening transition involves non-uniform flow and stress with considerable spatiotemporal complexity. Nevertheless, it is anticipated that dense suspensions of calcium carbonate particles with purely repulsive interactions may not conform to this scenario, as indicated by local pressure measurements with millimeter spatial resolution.

EXPERIMENT

Here we utilize Boundary Stress Microscopy (BSM), a technique capable of resolving stresses down to the micron scale, to search for evidence of stress heterogeneity. In addition, we measure the flow field at the lower boundary of the suspension where the boundary stress is measured.

FINDINGS

We find localized regions of high-stresses that are extended in the vorticity direction and propagate in the flow direction at a speed approximately half that of the rheometer's top plate. These high-stress regions proliferate with the applied stress accounting for the increased viscosity. Furthermore, the velocity of particles at the lower boundary of the suspension shows a significant and complex nonaffine flow that accompanies regions of high-stresses. Hence, our findings demonstrate that stress and flow inhomogeneity are intrinsic characteristics of shear-thickening suspensions, regardless of the nature of interparticle interactions.

Rheological Properties of Dense Particle Suspensions of Starches: Shear Thickening, Shear Jamming, and Shock Absorption Properties

Concentrated suspensions of Brownian and non-Brownian particles display distinctive rheological behavior highly dependent on shear rate and shear stress. Cornstarch suspensions, composed of starch particles from corn plants, served as a model for concentrated non-Brownian suspensions, demonstrating discontinuous shear thickening (DST) and dynamic shear jamming (SJ). However, starch particles from other plant sources have not yet been investigated, despite their different sizes and shapes. This study is focused on the evaluation of the effects of the structural parameters of starch particles by preparing concentrated suspensions of starch particles from 13 different plants at particle fractions of 25-50% and their rheological behavior through steady shear, pull-out, and ball-drop tests. Starch particles can be roughly classified as polygonal and ellipsoidal. The DST and SJ behavior typically reported for concentrated cornstarch suspensions were confirmed for other starch particles in both particle groups. The ball-drop test demonstrated excellent shock absorption properties for 11 concentrated suspensions of starch particles, except for sago palms. In the case of concentrated suspensions of starch particles, the particle fraction and shear applied were the dominant factors that significantly affected the rheological behavior, whereas the particle shape was not a primary contributor. The findings of this study drive further investigation on the effect of liquid and particle surface properties in concentrated particle suspensions on DST and SJ behaviors.

Effect of Liquid Properties on the Non-Newtonian Rheology of Concentrated Silica Suspensions: Discontinuous Shear Thickening, Shear Jamming, and Shock Absorbance

Concentrated particle suspensions exhibit rheological behavior, such as discontinuous shear thickening (DST) and dynamic shear jamming (SJ), which affect applications such as soft armors. Although the origin of this behavior in shear-activated particle-particle interactions has been identified, the effect of chemical factors, especially the role of liquids, on this behavior remains unexplored. Hydrogen bonding in suspensions has been proposed to be essential for frictional contacts between particles, and therefore, most studies on DST and SJ have focused on aqueous and protic organic media with a definite hydrogen bonding ability. To identify an alternative molecular mechanism, this study explored the effects of liquid polarity and an aprotic nature on the rheological behavior of concentrated suspensions of silica microparticles. Owing to their excellent particle dispersion, the DST behavior of polar liquids was observed, independent of protic and aprotic liquids. In contrast, nonpolar liquids formed particle agglomerates because of the particle-particle attraction and became a paste at a high particle fraction. The SJ behavior was confirmed for three aprotic organic liquids (propylene carbonate, 1,3-dimethyl-2-imidazolidinone, and 1,3-dimethylpropyleneurea), suggesting the hydrogen bonding ability of these aprotic liquids. The diverse mechanisms of shear-activated interactions between particles present material design possibilities for the non-Newtonian rheology of concentrated particle suspensions.

Dense bidisperse suspensions under non-homogeneous shear

We study the rheological behaviour of bidisperse suspensions in three dimensions under a non-uniform shear flow, made by the superimposition of a linear shear and a sinusoidal disturbance. Our results show that (i) only a streamwise disturbance in the shear-plane alters the suspension dynamics by substantially reducing the relative viscosity, (ii) with the amplitude of the disturbance determining a threshold value for the effect to kick-in and its wavenumber controlling the amount of reduction and which of the two phases is affected. We show that, (iii) the rheological changes are caused by the effective separation of the two phases, with the large or small particles layering in separate regions. We provide a physical explanation of the phase separation process and of the conditions necessary to trigger it. We test the results in the whole flow curve, and we show that the mechanism remains substantially unaltered, with the only difference being the nature of the interactions between particles modified by the phase separation.

Smart Polymers for Soft Materials: From Solution Processing to Organic Solids

Polymeric materials are ubiquitous in our everyday life, where they find a broad range of uses-spanning across common household items to advanced materials for modern technologies. In the context of the latter, so called "smart polymers" have received a lot of attention. These systems are soluble in water below their lower critical solution temperature Tℓ and often exhibit counterintuitive solvation behavior in mixed solvents. A polymer is known as smart-responsive when a slight change in external stimuli can significantly change its structure, functionm and stability. The interplay of different interactions, especially hydrogen bonds, can also be used for the design of lightweight high-performance organic solids with tunable properties. Here, a general scheme for establishing a structure-property relationship is a challenge using the conventional simulation techniques and also in standard experiments. From the theoretical side, a broad range of all-atom, multiscale, generic, and analytical techniques have been developed linking monomer level interaction details with macroscopic material properties. In this review, we briefly summarize the recent developments in the field of smart polymers, together with complementary experiments. For this purpose, we will specifically discuss the following: (1) the solution processing of responsive polymers and (2) their use in organic solids, with a goal to provide a microscopic understanding that may be used as a guiding tool for future experiments and/or simulations regarding designing advanced functional materials.

Solvents govern rheology and jamming of polymeric bead suspensions

Despite their apparent simplicity, suspensions of hard spheres in a Newtonian fluid show complex non-Newtonian behaviors and remain poorly understood. Recent works have pointed out the crucial role of interparticle contact forces in these behaviors. Here, we show that the same (polystyrene) particles, when immersed in different Newtonian solvents, show different behaviors at both the microscopic and macroscopic scales. Thanks to interparticle force measurements in each solvent together with rheological measurements, we show how the fine details of the pairwise particle interactions impact the macroscopic behavior. The rheological properties (shear thinning, shear thickening, jamming solid fraction value) of the suspensions, made up of same particles, are shown to depend on the nature of the solvent. Here, we highlight several mechanisms at the particle scale: the swelling of polymeric particles in an organic solvent, the role of colloidal repulsive forces and inertia for particles in a water solution, and the variation of the friction coefficient as a function of the load for particles immersed in silicone oils. Our study provides new quantitative data to test micromechanical models and simulations. It questions the interpretation of previous experimental works. Finally, it shows the need to systematically characterize the interparticle normal and tangential forces when studying a given suspension of hard spheres in a Newtonian fluid.

Aging of cornstarch particles suspended in aqueous solvents at room temperature

Starch suspensions are often used as model systems to demonstrate extreme shear-thickening effects. We study the aging of cornstarch particles in aqueous suspensions at room temperature by granulometry and rheological measurements. When starch is diluted in glycerol, no long-term changes are observed. The situation differs when water is used as solvent. For volume fractions up to 20 vol %, when the cornstarch suspensions in water are stored under continual agitation, we observe an increase in viscosity. When the cornstarch suspension is aged under quiescent conditions, no evolution of the particle size is observed. In the concentrated situation, the rheological properties vary independent of the storage condition. We show that the increase in viscosity is related to air trapped in the pore space and to the swelling of the granules and leakage of the amylopectin component of the starch into the surrounding water. The relative importance of the two processes depends upon the particle concentration and upon the energy brought to the system.

Density waves in shear-thickening suspensions

Shear thickening corresponds to an increase of the viscosity as a function of the shear rate. It is observed in many concentrated suspensions in nature and industry: water or oil saturated sediments, crystal-bearing magma, fresh concrete, silica suspensions, and cornstarch mixtures. Here, we reveal how shear-thickening suspensions flow, shedding light onto as yet non-understood complex dynamics reported in the literature. When shear thickening is important, we show the existence of density fluctuations that appear as periodic waves moving in the direction of flow and breaking azimuthal symmetry. They come with strong normal stress fluctuations of the same periodicity. The flow includes small areas of normal stresses of the order of tens of kilopascals and areas of normal stresses of the order of hundreds of pascals. These stress inhomogeneities could play an important role in the damage caused by thickening fluids in the industry.

Testing the Wyart-Cates model for non-Brownian shear thickening using bidisperse suspensions

There is a growing consensus that shear thickening of concentrated dispersions is driven by the formation of stress-induced frictional contacts. The Wyart-Cates (WC) model of this phenomenon, in which the microphysics of the contacts enters solely via the fraction f of contacts that are frictional, can successfully fit flow curves for suspensions of weakly polydisperse spheres. However, its validity for "real-life", polydisperse suspensions has yet to be seriously tested. By performing systematic simulations on bidisperse mixtures of spheres, we show that the WC model applies only in the monodisperse limit and fails when substantial bidispersity is introduced. We trace the failure of the model to its inability to distinguish large-large, large-small and small-small frictional contacts. By fitting our data using a polydisperse analogue of f that depends separately on the fraction of each of these contact types, we show that the WC picture of shear thickening is incomplete. Systematic experiments on model shear-thickening suspensions corroborate our findings, but highlight important challenges in rigorously testing the WC model with real systems. Our results prompt new questions about the microphysics of thickening for both monodisperse and polydisperse systems.

Competing Timescales Lead to Oscillations in Shear-Thickening Suspensions

Competing timescales generate novelty. Here, we show that a coupling between the timescales imposed by instrument inertia and the formation of interparticle frictional contacts in shear-thickening suspensions leads to highly asymmetric shear-rate oscillations. Experiments tuning the presence of oscillations by varying the two timescales support our model. The observed oscillations give access to a shear-jamming portion of the flow curve that is forbidden in conventional rheometry. Moreover, the oscillation frequency allows us to quantify an intrinsic relaxation time for particle contacts. The coupling of fast contact network dynamics to a slower system variable should be generic to many other areas of dense suspension flow, with instrument inertia providing a paradigmatic example.

Discontinuous rate-stiffening in a granular composite modeled after cornstarch and water

Cornstarch in water exhibits impact-activated solidification (IAS) and strong discontinuous shear thickening, with "shear jamming". However, these phenomena are absent in cornstarch in ethanol. Here we show that cornstarch granules swell under ambient conditions. We postulate that this granule swelling is linked to an interparticle force scale that introduces a discontinuous rate-dependence to the generation of stable contacts between granules. We studied this force scale by coating sand with ~ 2 μm-thick polydimethysiloxane, creating a material that exhibits a similar IAS and discontinuous deformation rate-stiffening despite being a granular composite, not a suspension. This result suggests rate-dependence can be tuned by coating granular materials, introducing an interparticle force scale from rate-dependent properties present in the coating material. Our work provides insights into the unique behavior of cornstarch in water, bridges our understanding of suspensions and dry granular materials, and introduces a method to make discontinuous rate-dependent materials without suspending particles.

From Yielding to Shear Jamming in a Cohesive Frictional Suspension

Simulations are used to study the steady shear rheology of dense suspensions of frictional particles exhibiting discontinuous shear thickening and shear jamming, in which finite-range cohesive interactions result in a yield stress. We develop a constitutive model that combines yielding behavior and shear thinning at low stress with the frictional shear thickening at high stresses, in good agreement with the simulation results. This work shows that there is a distinct difference between solids below the yield stress and in the shear-jammed state, as the two occur at widely separated stress levels, with an intermediate region of stress in which the material is flowable.

Interparticle hydrogen bonding can elicit shear jamming in dense suspensions

Dense suspensions of hard particles in a liquid can exhibit strikingly counter-intuitive behaviour, such as discontinuous shear thickening (DST)^1-7 and reversible shear jamming (SJ) into a state where flow is arrested and the suspension is solid-like^8-12. A stress-activated crossover from hydrodynamic interactions to frictional particle contacts is key for these behaviours^2-4,6,7,9,13. However, in experiments, many suspensions show only DST, not SJ. Here we show that particle surface chemistry plays a central role in creating conditions that make SJ readily observable. We find the system's ability to form interparticle hydrogen bonds when sheared into contact elicits SJ. We demonstrate this with charge-stabilized polymer microspheres and non-spherical cornstarch particles, controlling hydrogen bond formation with solvents. The propensity for SJ is quantified by tensile tests¹² and linked to an enhanced friction by atomic force microscopy. Our results extend the fundamental understanding of the SJ mechanism and open avenues for designing strongly non-Newtonian fluids.

Constraint-Based Approach to Granular Dispersion Rheology

We present a phenomenological model for granular suspension rheology in which particle interactions enter as constraints to relative particle motion. By considering constraints that are formed and released by stress respectively, we derive a range of experimental flow curves in a single treatment and predict singularities in viscosity and yield stress consistent with literature data. Fundamentally, we offer a generic description of suspension flow that is independent of bespoke microphysics.

Testing constitutive relations by running and walking on cornstarch and water suspensions

The ability of a person to run on the surface of a suspension of cornstarch and water has fascinated scientists and the public alike. However, the constitutive relation obtained from traditional steady-state rheology of cornstarch and water suspensions has failed to explain this behavior. In another paper we presented an averaged constitutive relation for impact rheology consisting of an effective compressive modulus of a system-spanning dynamically jammed structure [R. Maharjan et al., this issue, Phys. Rev. E 97, 052602 (2018)10.1103/PhysRevE.97.052602]. Here we show that this constitutive model can be used to quantitatively predict, for example, the trajectory and penetration depth of the foot of a person walking or running on cornstarch and water. The ability of the constitutive relation to predict the material behavior in a case with different forcing conditions and flow geometry than it was obtained from suggests that the constitutive relation could be applied more generally. We also present a detailed calculation of the added mass effect to show that while it may be able to explain some cases of people running or walking on the surface of cornstarch and water for pool depths H>1.2 m and foot impact velocities V_{I}>1.7 m/s, it cannot explain observations of people walking or running on the surface of cornstarch and water for smaller H or V_{I}.

Constitutive relation for the system-spanning dynamically jammed region in response to impact of cornstarch and water suspensions

We experimentally characterize the impact response of concentrated suspensions consisting of cornstarch and water. We observe that the suspensions support a large normal stress-on the order of MPa-with a delay after the impactor hits the suspension surface. We show that neither the delay nor the magnitude of the stress can yet be explained by either standard rheological models of shear thickening in terms of steady-state viscosities, or impact models based on added mass or other inertial effects. The stress increase occurs when a dynamically jammed region of the suspension in front of the impactor propagates to the opposite boundary of the container, which can support large stresses when it spans between solid boundaries. We present a constitutive relation for impact rheology to relate the force on the impactor to its displacement. This can be described in terms of an effective modulus but only after the delay required for the dynamically jammed region to span between solid boundaries. Both the modulus and the delay are reported as a function of impact velocity, fluid height, and weight fraction. We report in a companion paper the structure of the dynamically jammed region when it spans between the impactor and the opposite boundary [Allen et al., Phys. Rev. E 97, 052603 (2018)10.1103/PhysRevE.97.052603]. In a direct follow-up paper, we show that this constitutive model can be used to quantitatively predict, for example, the trajectory and penetration depth of the foot of a person walking or running on cornstarch and water [Mukhopadhyay et al., Phys. Rev. E 97, 052604 (2018)10.1103/PhysRevE.97.052604].

Depleted depletion drives polymer swelling in poor solvent mixtures

Establishing a link between macromolecular conformation and microscopic interaction is a key to understand properties of polymer solutions and for designing technologically relevant "smart" polymers. Here, polymer solvation in solvent mixtures strike as paradoxical phenomena. For example, when adding polymers to a solvent, such that all particle interactions are repulsive, polymer chains can collapse due to increased monomer-solvent repulsion. This depletion induced monomer-monomer attraction is well known from colloidal stability. A typical example is poly(methyl methacrylate) (PMMA) in water or small alcohols. While polymer collapse in a single poor solvent is well understood, the observed polymer swelling in mixtures of two repulsive solvents is surprising. By combining simulations and theoretical concepts known from polymer physics and colloidal science, we unveil the microscopic, generic origin of this collapse-swelling-collapse behavior. We show that this phenomenon naturally emerges at constant pressure when an appropriate balance of entropically driven depletion interactions is achieved.

Oscillatory rheology of dense, athermal suspensions of nearly hard spheres below the jamming point

The viscosity of a dense suspension has contributions from hydrodynamics and particle interactions, both of which depend upon the flow-induced arrangement of particles into fragile structures. Here, we study the response of nearly hard sphere suspensions to oscillatory shear using simulations and experiments in the athermal, non-inertial limit. Three distinct regimes are observed as a function of the strain amplitude γ₀. For γ₀ < 10^-1, initially non-contacting particles remain separated and the suspension behaves similarly to a Newtonian fluid, with the shear stress proportional to the strain rate, and the normal stresses close to zero. For γ₀ > 10¹, the microstructure becomes well-established at the beginning of each shear cycle and the rheology is quasi-Newtonian: the shear stress varies with the rate, but flow-induced structures lead to non-zero normal stresses. At intermediate γ₀, particle-particle contacts break and reform across entire oscillatory cycles, and we probe a non-linear regime that reveals the fragility of the material. Guided by these features, we further show that oscillatory shear may serve as a diagnostic tool to isolate specific stress contributions in dense suspensions, and more generally in those materials whose rheology has contributions with both hydrodynamic and non-hydrodynamic origin.