Alternative Frictional Model for Discontinuous Shear Thickening of Dense Suspensions: Hydrodynamics

Rheological Regimes in Agitated Granular Media under Shear

Agitated granular media have a rich rheology: they exhibit Newtonian behavior at low shear rate and density, develop a yield stress at high density, and cross over to Bagnoldian shear thickening when sheared rapidly-making it challenging to encompass them in one theoretical framework. We measure the rheology of air-fluidized glass particles, spanning 5 orders of magnitude in shear rate. By comparing fluidization-induced to Brownian agitation, we show that all rheological regimes can be delineated by two dimensionless numbers-the Péclet number, Pe, and the ratio of shear-to-fluidization power, Π-and propose a constitutive relation that captures all flow behaviors, qualitatively and quantitatively, in one unified framework.

Suspension physics govern the multiscale dynamics of blood flow in sickle cell disease

In diseases from diabetes to malaria, blood dynamics are significantly altered, resulting in poor clinical outcomes. However, the multiscale mechanisms that determine blood flow in the microcirculation in health and disease are undefined, largely owing to the difficulty in directly linking cell properties to whole-blood rheology. Here, we overcome these difficulties by developing a microfluidic platform to measure red blood cell properties and flow dynamics in the same blood samples from donors. We focus on sickle cell disease (SCD), a genetic disorder that causes red blood cells to stiffen in deoxygenated conditions, with disease pathology driven by oxygen-dependent blood rheology. Our linked cell and whole-blood measurements establish that increases in effective resistances in heterogeneous suspensions are driven by increases in the proportion of stiff cells, similar macroscopically to the behavior of rigid-particle suspensions. Furthermore, by combining simulations with spatially resolved measurements of cell dynamics, we show how the spatio-temporal organization of stiff and deformable cells determines blood rheology and drives disease pathophysiology. In the presence of deformable cells, the stiffened cells marginate towards channel walls, increasing effective wall friction. In fully deoxygenated conditions in which all cells are stiffened, significant heterogeneity in cell volume fraction along the direction of flow causes localized jamming, drastically increasing effective viscous flow resistance. Our work defines the relevant suspension physics required to understand pathological blood rheology in SCD and other diseases affecting red blood cell properties. More broadly, we reveal the multiscale processes that determine emergent rheology in heterogeneous particle suspensions.

Characterizing sliding and rolling contacts between single particles

Contacts between particles in dense, sheared suspensions are believed to underpin much of their rheology. Roughness and adhesion are known to constrain the relative motion of particles, and thus globally affect the shear response, but an experimental description of how they microscopically influence the transmission of forces and relative displacements within contacts is lacking. Here, we show that an innovative colloidal-probe atomic force microscopy technique allows the simultaneous measurement of normal and tangential forces exchanged between tailored surfaces and microparticles while tracking their relative sliding and rolling, unlocking the direct measurement of coefficients of rolling friction, as well as of sliding friction. We demonstrate that, in the presence of sufficient traction, particles spontaneously roll, reducing dissipation and promoting longer-lasting contacts. Conversely, when rolling is prevented, friction is greatly enhanced for rough and adhesive surfaces, while smooth particles coated by polymer brushes maintain well-lubricated contacts. We find that surface roughness induces rolling due to load-dependent asperity interlocking, leading to large off-axis particle rotations. In contrast, smooth, adhesive surfaces promote rolling along the principal axis of motion. Our results offer direct values of friction coefficients for numerical studies and an interpretation of the onset of discontinuous shear thickening based on them, opening up ways to tailor rheology via contact engineering.

Attraction-Enhanced Emergence of Friction in Colloidal Matter

How frictional effects emerge at the microscopic level in particulate materials remains a challenging question, particularly in systems subject to thermal fluctuations due to the transient nature of interparticle contacts. Here, we directly relate particle-level frictional arrest to local coordination in an attractive colloidal model system. We reveal that the orientational dynamics of particles slows down exponentially with increasing coordination number due to the emergence of frictional interactions, the strength of which can be tuned simply by varying the attraction strength. Using a simple computer simulation model, we uncover how the interparticle interactions govern the formation of frictional contacts between particles. Our results establish quantitative relations between friction, coordination, and interparticle interactions. This is a key step toward using interparticle friction to tune the mechanical properties of particulate materials.

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.

Shear thickening in suspensions of particles with dynamic brush layers

Control of frictional interactions among liquid-suspended particles has led to tunable, strikingly non-Newtonian rheology via the formation of strong flow constraints as particles come into close proximity under shear. Typically, these frictional interactions have been in the form of physical contact, controllable via particle shape and surface roughness. We investigate a different route, where molecular bridging between nearby particle surfaces generates a controllable constraint to relative particle movement. This is achieved with surface-functionalized colloidal particles capable of forming dynamic covalent bonds with telechelic polymers that comprise the suspending fluid. At low shear stress this results in particles coated with a uniform polymer brush layer. Beyond an onset stress σ* the telechelic polymers become capable of bridging and generate shear thickening. Over the size range investigated, we find that the dynamic brush layer leads to dependence of σ* on particle diameter that closely follows a power law with exponent -1.76. In the shear thickening regime, we observe an enhanced dilation in measurements of the first normal stress difference N1 and reduction in the extrapolated volume fraction required for jamming, both consistent with an effective particle friction that increases with decreasing particle diameter. These results are discussed in light of predictions for suspensions of hard spheres and of polymer-grafted particles.

Decoupling of rotation and translation at the colloidal glass transition

In dense particle systems, the coupling of rotation and translation motion becomes intricate. Here, we report the results of confocal fluorescence microscopy where simultaneous recording of translational and rotational particle trajectories from a bidisperse colloidal dispersion is achieved by spiking the samples with rotational probe particles. The latter consist of colloidal particles containing two fluorescently labeled cores suited for tracking the particle's orientation. A comparison of the experimental data with event driven Brownian simulations gives insights into the system's structure and dynamics close to the glass transition and sheds new light onto the translation-rotation coupling. The data show that with increasing volume fractions, translational dynamics slows down drastically, whereas rotational dynamics changes very little. We find convincing agreement between simulation and experiments, even though the simulations neglect far-field hydrodynamic interactions. An additional analysis of the glass transition following mode coupling theory works well for the structural dynamics but indicates a decoupling of the diffusion of the smaller particle species. Shear stress correlations do not decorrelate in the simulated glass states and are not affected by rotational motion.

Order-disorder transition during shear thickening in bidisperse dense suspensions

Shear thickening of multimodal suspensions has proven difficult to understand because the rheology depends largely on the microscopic details of stress-induced frictional contacts at different particle size distributions (PSDs). Our discrete particle simulations below a critical volume fraction ϕc over a broad range of shear rates and PSDs elucidate the basic mechanism of order-disorder transition. Around the theoretical optimal PSD (relative content of small particles ζ1= 0.26), particles order into a layered structure in the Newtonian regime. At the onset of shear thickening, this layered structure transforms to a disordered one, accompanied by an abrupt viscosity jump. Minor increase in large-large particle contacts after the order-disorder transition causes apparent increase in radial force along the compressional axis. Bidisperse suspensions with less regular but stable layered structure at ζ1= 0.50 show good fluidity in the shear thickening regime. This work shows that in inertial flows where particle collisions dominate, order-disorder transition could play an essential role in shear thickening for bidisperse suspensions.

Rheology of Highly Filled Polymer Compositions-Limits of Filling, Structure, and Transport Phenomena

The current state of the rheology of various polymeric and other materials containing a high concentration of spherical solid filler is considered. The physics of the critical points on the concentration scale are discussed in detail. These points determine the features of the rheological behavior of the highly filled materials corresponding to transitions from a liquid to a yielding medium, elastic-plastic state, and finally to an elastic solid-like state of suspensions. Theoretical and experimental data are summarized, showing the limits of the most dense packing of solid particles, which is of key importance for applications and obtaining high-quality products. The results of model and fine structural studies of physical phenomena that occur when approaching the point of filling the volume, including the occurrence of instabilities, are considered. The occurrence of heterogeneity in the form of individual clusters is also described. These heterogeneous objects begin to move as a whole that leads to the appearance of discontinuities in the suspension volume or wall slip. Understanding these phenomena is a key for particle technology and multiphase processing.

Stress-activated friction in sheared suspensions probed with piezoelectric nanoparticles

A hallmark of concentrated suspensions is non-Newtonian behavior, whereby the viscosity increases dramatically once a characteristic shear rate or stress is exceeded. Such strong shear thickening is thought to originate from a network of frictional particle-particle contact forces, which forms under sufficiently large stress, evolves dynamically, and adapts to changing loads. While there is much evidence from simulations for the emergence of this network during shear thickening, experimental confirmation has been difficult. Here, we use suspensions of piezoelectric nanoparticles and exploit the strong local stress focusing within the network to activate charge generation. This charging can then be detected in the measured ac conductance and serve as a signature of frictional contact formation. The direct link between stress-activated frictional particle interactions and piezoelectric suspension response is further demonstrated by tracking the emergence of structural memory in the contact network under oscillatory shear and by showing how stress-activated friction can drive mechano-transduction of chemical reactions with nonlinear reaction kinetics. Taken together, this makes the ac conductance of piezoelectric suspensions a sensitive in-situ reporter of the micromechanics associated with frictional interactions.

Dynamic-bond-induced sticky friction tailors non-Newtonian rheology

Frictional network formation has become a new paradigm for understanding the non-Newtonian shear-thickening behavior of dense suspensions. Recent studies have exclusively focused on interparticle friction that instantaneously vanishes when applied shear is ceased. Herein, we investigate a friction that emerges from dynamic chemical bridging of functionalized particle surfaces sheared into close proximity. This enables tailoring of both friction magnitude and the time release of the frictional coupling. The experiments use dense suspensions of thiol-functionalized particles suspended in ditopic polymers endcapped with benzalcyanoacetamide Michael-acceptors. The subsequent room temperature, catalyst-free dynamic thia-Michael reactions can form bridging interactions between the particles with dynamic covalent bonds that linger after formation and release in the absence of shear. This chemical friction mimics physical friction but is stickier, leading to tunable rheopexy. The effect of sticky friction on dense suspension rheology is explored by varying the electronic nature of the benzalcyanoacetamide moiety, the molecular weight of the ditopic polymers, the amount of a competitive bonding compound, and temperature. These results demonstrate how dynamic-bond-induced sticky friction can be used to systematically control the time dependence of the non-Newtonian suspension rheology.

Scaling relationships between viscosity and diffusivity in shear-thickening suspensions

Dense suspensions often exhibit a dramatic response to large external deformation. The recent body of work has related this behavior to transition from an unconstrained lubricated state to a constrained frictional state. Here, we use numerical simulations to study the flow behavior and shear-induced diffusion of frictional non-Brownian spheres in two dimensions under simple shear flow. We first show that both viscosity η and diffusivity D/ of the particles increase under characteristic shear stress, which is associated with lubrication to frictional transition. Subsequently, we propose a one-to-one relationship between viscosity and diffusivity using the length scale ξ associated with the size of collective motions (rigid clusters) of the particles. We demonstrate that η and D/ are controlled by ξ in two distinct flow regimes, i.e. in the frictionless and frictional states, where the one-to-one relationship is described as a crossover from D/ ∼ η (frictionless) to η1/3 (frictional). We also confirm that the proposed power laws are insensitive to the interparticle friction and system size.

Rheology of 3D printable ceramic suspensions: effects of non-adsorbing polymer on discontinuous shear thickening

Concentrated suspensions of particles at volume fractions (ϕ) ≥ 0.5 often exhibit complex rheological behavior, transitioning from shear thinning to shear thickening as the shear stress or shear rate is increased. These suspensions can be extruded to form 3D structures, with non-adsorbing polymers often added as rheology modifiers to improve printability. Understanding how non-adsorbing polymers affect the suspension rheology, particularly the onset of shear thickening, is critical to the design of particle inks that will extrude uniformly. In this work, we examine the rheology of concentrated aqueous suspensions of colloidal alumina particles and the effects of adding non-adsorbing polyvinylpyrrolidone (PVP). First, we show that suspensions with ϕ_alumina = 0.560-0.575 exhibited discontinuous shear thickening (DST), where the viscosity increased by up to two orders of magnitude above an onset stress (τ_min). Increasing ϕ_alumina from 0.550 to 0.575 increased the viscosity and yield stress in the shear thinning regime and decreased τ_min. Next, PVP was added at concentrations within the dilute and semi-dilute non-entangled regimes of polymer conformation (ϕ_PVP = 0.005-0.050) to suspensions with constant ϕ_alumina = 0.550. DST was observed in all cases and increasing ϕ_PVP increased the viscosity and yield stress. Interestingly, increasing ϕ_PVP also increased τ_min. We posit that the free PVP chains act as lubricants between alumina particles, increasing the stress needed to induce thickening. Finally, we demonstrate through direct comparisons of suspensions with and without PVP how non-adsorbing polymer addition can extend the extrusion processing window due to the increase in τ_min.

Microscopic mechanism study of the rheological behavior of non-Newtonian fluids based on dissipative particle dynamics

Non-Newtonian fluid rheological properties are a hot research topic for realizing intelligent applications. In order to investigate the microscopic mechanism and structural evolution process of the nonlinear rheological behavior of non-Newtonian fluids, this paper systematically investigates two continuous nonlinear rheological behaviors of non-Newtonian fluids, namely shear-thickening and shear-thinning rheological properties, using a non-Newtonian fluid system composed of polyethylene glycol (PEG) mixed with nano-silica (Nano-SiO₂) by a dissipative particle dynamics (DPD) method. It is shown that at low shear rates, the molecular chains of PEG in the fluid are stretched due to shear flow and the molecular structure is transformed into an ordered state; and the effective hydrodynamic radius of Nano-SiO₂ beads decreases, which makes the translational friction coefficient of the beads decrease and the system mobility increases, exhibiting shear-thinning behavior. When the shear rate exceeds the critical value, the contact and collision probability between Nano-SiO₂ beads in the non-Newtonian fluid increases; a large number of silicon hydroxyl groups exist on the surface of Nano-SiO₂, which form a large number of hydrogen bonds when they are close to each other and constrain the particle separation, resulting in a large aggregation of Nano-SiO₂ beads, leading to an increase in the effective kinetic radius of Nano-SiO₂ beads and an increase in the coefficient of translational friction, forming a blockage of the fluid system and exhibiting a shear-thickening behavior. Our study provides insights for understanding the rheological behavior of non-Newtonian fluids from a microscopic perspective, and contributes to the intelligent application of non-Newtonian fluids.

Structure and Dynamics of Force Clusters and Networks in Shear Thickening Suspensions

Dense suspensions can exhibit shear thickening in response to large deformation. A consensus has emerged over the past few years on the formation of force networks, that span the entire system size, that lead to increased resistance to motion. Nonetheless, the characteristics of these networks are to a large extent poorly understood. Here, force networks formed in continuous and discontinuous shear thickening dense suspensions (CST and DST, respectively) are studied. We first show the evolution of the network formation and its topological heterogeneities as the applied stress increases. Subsequently, we identify force communities and coarse grain the suspension into a cluster network, and show that cluster-level dynamics are responsible for stark differences between the CST and DST behavior. Our results suggest that the force clusters formed in the DST regime are considerably more constrained in their motion, while CST clusters are loosely connected to their surrounding clusters.

Origin of Two Distinct Stress Relaxation Regimes in Shear Jammed Dense Suspensions

Many dense particulate suspensions show a stress induced transformation from a liquidlike state to a solidlike shear jammed (SJ) state. However, the underlying particle-scale dynamics leading to such striking, reversible transition of the bulk remains unknown. Here, we study transient stress relaxation behaviour of SJ states formed by a well-characterized dense suspension under a step strain perturbation. We observe a strongly nonexponential relaxation that develops a sharp discontinuous stress drop at short time for high enough peak-stress values. High resolution boundary imaging and normal stress measurements confirm that such stress discontinuity originates from the localized plastic events, whereas system spanning dilation controls the slower relaxation process. We also find an intriguing correlation between the nature of transient relaxation and the steady-state shear jamming phase diagram obtained from the Wyart-Cates model.

Microstructure of continuous shear thickening colloidal suspensions determined by rheo-VSANS and rheo-USANS

Research on shear thickening colloidal suspensions demonstrates that measurements of the microstructure can elucidate the source of the rheological material properties in the shear thickened state as well as critically test simulations and theory based on a variety of mechanisms such as enhanced lubrication hydrodynamics, elastohydrodynamics, and contact friction. Prior work on continuous shear thickening dispersions with a well-defined shear thickened state identified the formation of hydroclusters as characteristic of this state, determined the anisotropy in the nearest neighbor distribution, and used this information to test prevailing theories and simulations. However, important questions remain about the mesoscale (i.e., particle cluster scale) microstructure of the shear thickened state. Here we employ neutron scattering methods applied to shearing colloidal dispersions of spherical particles with two extremes of friction and lubrication surface properties to resolve the longer-length scale microstructure in the shear thickened state. Hydroclusters are shown to be highly localized, in agreement with prior neutron scattering and direct optical measurements, but in disagreement with the most recent simulations that predict a longer-range structure formation. These results combined with prior measurements provide experimental evidence about the length scale of microstructure formation in continuous shear thickening suspensions necessary to improve our understanding of the phenomenon as well as guide theoretical investigations that quantitatively link nanoscale forces to macroscopic properties in the shear thickened state.

Jamming Distance Dictates Colloidal Shear Thickening

We report experimental and computational observations of dynamic contact networks for colloidal suspensions undergoing shear thickening. The dense suspensions are comprised of sterically stabilized poly(methyl methacrylate) colloids that are spherically symmetric and have varied surface roughness. Confocal rheometry and dissipative particle dynamics simulations show that the shear thickening strength β scales exponentially with the scaled deficit contact number and the scaled jamming distance. Rough colloids, which experience additional rotational constraints, require an average of 1.5-2 fewer particle contacts as compared to smooth colloids, in order to generate the same β. This is because the surface roughness enhances geometric friction in such a way that the rough colloids do not experience a large change in the free volume near the jamming point. The available free volume for colloids of different roughness is related to the deficiency from the maximum number of nearest neighbors at jamming under shear. Our results further suggest that the force per contact is different for particles with different morphologies.

Supramolecular assembly inspired molecular engineering to dynamically tune non-Newtonian fluid:from quasi-static flowability-free to shear thickening

Shear thickening fluids (STFs) have been the research focus for decades because of the prospect as a damping ingredient. However, their inherent liquid character confines their practical applications. In this work, inspired by the assembly engineering, novel gelatinous shear thickening fluids (GSTFs) are fabricated by integrating low molecular weight gelators (LMWGs) into STFs and investigated by rheological experiments. The results show that the apparent performances of GSTFs are determined by the LMWGs content. LMWGs inside GSTFs can assemble into three-dimensional network that can constraint the flowability of liquid molecular and their content dominate the density and strength of assembly network. At a moderate content, GSTFs exhibit desired properties with restricted quasi-static flowability and almost undamaged dynamic shear thickening character. While a higher content will disappear shear thickening and a lower content cannot gelate STFs. Besides, three different LMWGs are employed to gelate STFs and all they can gelate STFs in spite of the distinct minimum gelation concentration, indicating the universality for GSTFs preparation and the superiority of a reasonable molecular structure of LMWGs. Further, the temperature sweep experiments suggest that GSTFs can endure higher temperature without flowing due to its higher gel-sol transition temperature. Basing on these advanced mechanical properties, we believe that the GSTFs with more expected characters have significance for the study of non-Newtonian fluids and will broaden the special application field of STFs.

Shear thickening behavior in dense repulsive and attractive suspensions of hard spheres

Shear thickening in stable dense colloidal suspensions is a reversible phenomenon and no hysteresis is observed in the flow curve measurements. However, a reduction in the stability of colloids promotes particle aggregation and introduces a time dependent rheological response. In this work, by using a model colloidal system of hard spherical silica particles (average diameter of 415 nm) with varying particle volume fractions 0.2 ≤ ϕ ≤ 0.56, we study the effect of particle stability on the hysteresis of the shear thickening behavior of these suspensions. The particle stability is manipulated by adding a simple monovalent salt (sodium chloride) in the silica suspension with varying concentrations α ∈ [0,0.5] M. For repulsive and weakly attractive suspensions, the flow behavior is history independent and the shear thickening behavior does not exhibit hysteresis. However, significant hysteresis is observed in rheological measurements for strongly attractive suspensions, with shear history playing a critical role due to the dynamic nature of particle clusters, resulting in time dependent hysteresis behavior. By performing numerical simulations, we find that this hysteresis behavior arises due to the competition among shear, electrostatic repulsive, van der Waals attractive, and frictional contact forces. The critical shear stress (i.e., the onset of shear thickening) decreases with increasing salt concentrations, which can be captured by a scaling relationship based on the force balance between particle-particle contact force and electrostatic repulsive force. Our combined experimental and simulation results imply the formation of particle contacts in our sheared suspensions.

Unifying disparate rate-dependent rheological regimes in non-Brownian suspensions

A typical dense non-Brownian particulate suspension exhibits shear thinning (decreasing viscosity) at low shear rate or stress followed by a Newtonian plateau (constant viscosity) at intermediate shear rate or stress values which transitions to shear thickening (increasing viscosity) beyond a critical shear rate or stress value and finally undergoes a second shear thinning transition at extremely high shear rate or stress values. In this study, we unify and quantitatively reproduce all the disparate rate-dependent regimes and the corresponding transitions for a dense non-Brownian suspension with increasing shear rate or stress. We employ discrete particle dynamics simulations based on the proposed mechanism to elucidate its accuracy. We find that a competition between interparticle interactions of hydrodynamic and nonhydrodynamic origins and the switching in the dominant stress scale with increasing the shear rate or stress lead to each of the above transitions. Inclusion of traditional hydrodynamic interactions, attractive or repulsive Derjaguin-Landau-Verwey-Overbeek (DLVO) interactions the interparticle contact interactions, and a constant friction (or other constraint mechanism) reproduces the initial thinning as well as the shear thickening transition. However, to quantitatively capture the intermediate Newtonian plateau and the second shear thinning, an additional nonhydrodynamic interaction of non-DLVO origin and a decreasing coefficient of friction, respectively, are essential, thus providing an explanation for the presence of the intermediate Newtonian plateau along with reproducing the second shear thinning in a single framework. Expressions utilized for various interactions and friction are determined from experimental measurements and hence result in excellent quantitative agreement between the simulations and previous experiments.

Machine learning of lubrication correction based on GPR for the coupled DPD-DEM simulation of colloidal suspensions

Hydrodynamic interactions have a major impact on the suspension properties, but they are absent in atomic and molecular fluids due to a lack of intervening medium at close range. To reproduce the correct hydrodynamic interactions, lubrication correction is essential to compensate the missing short-range hydrodynamics from the fluids. However, lubrication correction requires many simulations in particle-based simulations of colloidal suspensions. To address the problem, we employ an active learning strategy based on Gaussian process regression (GPR) for normal and tangential lubrication corrections to significantly reduce the number of necessary simulations and apply the correction to the coupled multiscale simulation of monodisperse hard-sphere colloidal suspensions. In particular, a single-particle dissipative particle dynamics (DPD) model with parameter correction is used to describe the solvent-solvent and colloid-solvent interactions, and a discrete element method (DEM) model to depict the colloid-colloid frictional contacts. The lubrication correction results demonstrate that only six and four independent simulations (observation points for GPR training) are required to achieve accurate normal and tangential lubrication corrections, respectively. To validate the machine learning of lubrication correction based on GPR, we investigate the self-diffusion coefficients of colloids, suspension rheology and microstructure using the coupled DPD-DEM model with GPR lubrication correction. Our simulation results show that the machine learning of lubrication correction based on GPR is effective and the lubrication corrected DPD-DEM model is indeed capable of accurately capturing hydrodynamic interactions and correctly reproducing dynamical and rheological properties of colloidal suspensions. Moreover, the machine learning of lubrication correction based on GPR is not limited to the coupled DPD-DEM simulation of colloidal suspensions presented here, but can be easily applied to other particle-based simulations of particulate suspensions.

In silico biophysics and hemorheology of blood hyperviscosity syndrome

Hyperviscosity syndrome (HVS) is characterized by an increase of the blood viscosity by up to seven times the normal blood viscosity, resulting in disturbances to the circulation in the vasculature system. HVS is commonly associated with an increase of large plasma proteins and abnormalities in the properties of red blood cells, such as cell interactions, cell stiffness, and increased hematocrit. Here, we perform a systematic study of the effect of each biophysical factor on the viscosity of blood by employing the dissipative particle dynamic method. Our in silico platform enables manipulation of each parameter in isolation, providing a unique scheme to quantify and accurately investigate the role of each factor in increasing the blood viscosity. To study the effect of these four factors independently, each factor was elevated more than its values for a healthy blood while the other factors remained constant, and viscosity measurement was performed for different hematocrits and flow rates. Although all four factors were found to increase the overall blood viscosity, these increases were highly dependent on the hematocrit and the flow rates imposed. The effect of cell aggregation and cell concentration on blood viscosity were predominantly observed at low shear rates, in contrast to the more magnified role of cell rigidity and plasma viscosity at high shear rates. Additionally, cell-related factors increase the whole blood viscosity at high hematocrits compared with the relative role of plasma-related factors at lower hematocrits. Our results, mapped onto the flow rates and hematocrits along the circulatory system, provide a correlation to underpinning mechanisms for HVS findings in different blood vessels.

Signature of Transition between Granular Solid and Fluid: Rate-Dependent Stick Slips in Steady Shearing

Despite extensive studies on either smooth granular-fluid flow or the solidlike deformation at the slow limit, the change between these two extremes remains largely unexplored. By systematically investigating the fluctuations of tightly packed grains under steady shearing, we identify a transition zone with prominent stick-slip avalanches. We establish a state diagram, and propose a new dimensionless shear rate based on the speed dependence of interparticle friction and particle size. With fluid-immersed particles confined in a fixed volume and forced to "flow" at viscous numbers J decades below reported values, we answer how a granular system can transition to the regime sustained by solid-to-solid friction that goes beyond existing paradigms based on suspension rheology.

Consititutive modeling of time-dependent flows in frictional suspensions

This paper summarizes recent joint work towards a constitutive modelling framework for dense granular suspensions. The aim is to create a time-dependent, tensorial theory that can implement the physics described in steady state by the Wyart-Cates model. This model of shear thickening suspensions supposes that lubrication films break above a characteristic normal force so that frictional contact forces come into play: the resulting non-sliding constraints can be enough to rigidify a system that would flow freely at lower stresses [1]. Implementing this idea for time-dependent flows requires the introduction of new concepts including a configuration-dependent ‘jamming coordinate’, alongside a decomposition of the velocity gradient tensor into compressive and extensional components which then enter the evolution equation for particle contacts in distinct ways. The resulting approach [2, 3] is qualitatively successful in addressing (i) the collapse of stress during flow reversal in shear flow, and (ii) the ability of transverse oscillatory flows to unjam the system. However there is much work required to refine this approach towards quantitative accuracy, by incorporating more of the physics of contact evolution under flow as determined by close interrogation of particle-based simulations.

Transition from steady shear to oscillatory shear rheology of dense suspensions

Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.

Shear Thickening and Jamming of Dense Suspensions: The "Roll" of Friction

Particle-based simulations of discontinuous shear thickening (DST) and shear jamming (SJ) suspensions are used to study the role of stress-activated constraints, with an emphasis on resistance to gearlike rolling. Rolling friction decreases the volume fraction required for DST and SJ, in quantitative agreement with real-life suspensions with adhesive surface chemistries and "rough" particle shapes. It sets a distinct structure of the frictional force network compared to only sliding friction, and from a dynamical perspective leads to an increase in the velocity correlation length, in part responsible for the increased viscosity. The physics of rolling friction is thus a key element in achieving a comprehensive understanding of strongly shear-thickening materials.