non-Brownian Suspensions Research Articles

Unifying suspension and granular shear-induced self-diffusion

Shear-induced self-diffusion is a fundamental mode of transport in granular flows. Yet its critical behaviour and dependence on the particle solid fraction are still unclear. Here, we rationalize these dependencies by performing two-dimensional pressure-imposed numerical simulations of dense non-Brownian frictional suspensions. Our results, combined with existing numerical data on inertial granular flows, show that the shear-induced diffusion coefficients of both systems can be captured by a single function of the distance to jamming. They further show that the grain diffusive behaviour is underpinned by a specific random walk process, having a constant elementary step length driven at a frequency that increases with the solid fraction. The proposed scaling laws pave the way for a better understanding of mixing processes in granular media.

Dispersion and particle pressure in sedimenting suspensions with hydrodynamic interactions and particle inertia

Particle inertia is a feature of suspension dynamics that is often neglected. However, its neglect changes the order of the governing equations of particle motion. In this work, we examine the impact of the inertia of particles over their velocity fluctuations and the resulting stresses. In order to observe effects of particle inertia, we consider dusty-gas suspensions of spherical weakly Brownian microparticles sedimenting in a Newtonian fluid at low Reynolds numbers. We first investigate the motion of a single spherical particle. The simulations for this simple case allow us to define the appropriate physical parameters of the flow and to validate the numerical calculation of velocity statistics over a range of Péclet and Stokes numbers. Next, we perform Langevin dynamics simulations with periodic boundary conditions to integrate the equations governing the translational and rotational motions of N hydrodynamically interacting particles suspended in the viscous fluid. The first aim of this paper is to examine the behavior of the velocity variance of inertial particles, with small fluctuations produced by Brownian motion at the moderate Péclet numbers. The interplay between mild Brownian forces and hydrodynamic interactions decreases the anisotropy of velocity fluctuations observed in non-Brownian suspensions (Pe≫1) of sedimenting particles. The simulation results suggest that particle inertia attenuates the magnitude of velocity fluctuations produced by both Brownian motion and viscous hydrodynamic interactions. Additionally, we study the long-time behavior of the velocity fluctuations by calculating their autocorrelation functions and the particle diffusivities. The numerical simulations show clear evidence of particle pressure arising from the flow disturbance produced by hydrodynamic interactions in a suspension. From the particle velocity variance obtained in the present numerical simulations, we propose a simple model for particle-phase pressure as a linear function of the particle volume fraction ϕ⩽5%, which is usually a closure quantity required in models of more complex particulate systems such as fluidized beds.

Mass and momentum balance during particle migration in the pressure-driven flow of frictional non-Brownian suspensions

The transient shear-induced particle migration of frictional non-Brownian suspensions is studied using particle-resolved simulations. The numerical method – the fictitious domain method – is well suited to heterogeneous flows thanks to a frame-invariant formulation of the subgrid (lubrication) corrections that does not involve the ambient flow (Orsi et al., J. Comput. Phys., vol. 474, 2023, 111823). The paper aims to give an accurate quantitative picture of the mass and momentum balance during the flow. The various assumptions and local constitutive laws that together form the suspension balance model (SBM) are thoroughly examined. To this purpose, the various quantities of interest are locally averaged in space and time, and their profile across the channel is extensively studied, with specific attention to the time evolution of the different contributions, either hydrodynamic in nature or from contact interactions, to the shear and normal stresses. The latter, together with the velocity gradient in the wall-normal direction and the volume fraction profile, yield the local constitutive laws, which are compared with their counterpart obtained in homogeneous shear flow. A fair agreement is observed except in a layering area at the boundaries and at the very centre of the channel. In addition, the main assumption of the SBM, i.e. the local relation between the hydrodynamic force on the particles and the particle flux, is meticulously investigated. The hydrodynamic force is found to be mainly a drag, except in the lower range of the probed volume fractions, where a non-drag contribution is observed.

Shear-induced migration of rigid spheres in a Couette flow

Concentration inhomogeneities occur in many flows of non-Brownian suspensions. Their modeling necessitates the description of the relative motion of the particle phase and of the fluid phase, as well as the accounting for their interaction, which is the object of the suspension balance model (SBM). We systematically investigate the dynamics and the steady state of shear-induced migration in a wide-gap Couette flow for a wide range of particle volume fraction, and we test the ability of the SBM to account for the observations. We use a model suspension for which macroscopic particle stresses are known. Surprisingly, the observed magnitude of migration is much lower than that predicted by the SBM when the particle stress in the SBM is equated to the macroscopic particle stress. Another noteworthy observation is the quasi-absence of migration for semidilute suspensions. From the steady-state volume fraction profiles, we derive the local particle normal stress responsible for shear-induced migration according to the SBM. However, the observed dynamics of migration is much faster than that predicted by the SBM when using this stress in the model. More generally, we show that it is not possible to build a local friction law consistent with both the magnitude and the dynamics of migration within the standard SBM framework. This suggests that there is a missing term in the usual macroscopic constitutive law for the particle normal stress driving migration. The SBM is indeed capable of accurately predicting both the magnitude and the dynamics of migration when a tentative phenomenological term involving a concentration gradient is added to the particle normal stresses determined in macroscopic experiments.

Role of the constriction angle on the clogging by bridging of suspensions of particles

Confined flows of particles can lead to clogging, and therefore failure, of various fluidic systems across many applications. As a result, design guidelines need to be developed to ensure that clogging is prevented or at least delayed. In this Letter, we investigate the influence of the angle of reduction in the cross section of the channel on the bridging of semidilute and dense non-Brownian suspensions of spherical particles. We observe a decrease of the clogging probability with the reduction of the constriction angle. This effect is more pronounced for dense suspensions close to the maximum packing fraction where particles are in contact in contrast to semidilute suspensions. We rationalize this difference in terms of arch selection. We describe the role of the constriction angle and the flow profile, providing insights into the distinct behavior of semidilute and dense suspensions. Published by the American Physical Society 2024

Rheology of bidisperse non-Brownian suspensions.

We study the rheology of bidisperse non-Brownian suspensions using particle-based simulation, mapping the viscosity as a function of the size ratio of the species, their relative abundance, and the overall solid content. The variation of the viscosity with applied stress exhibits shear-thickening phenomenology irrespective of composition, though the stress-dependent limiting solids fraction governing the viscosity and its divergence point are nonmonotonic in the mixing ratio. Contact force data demonstrate an asymmetric exchange in the dominant stress contribution from large-large to small-small particle contacts as the mixing ratio of the species evolves. Combining a prior model for shear thickening with one for composition-dependent jamming, we obtain a full description of the rheology of bidisperse non-Brownian suspensions capable of predicting effects such as the viscosity reduction observed upon adding small particles to a suspension of large particles.

Heterogeneous solvent dissipation coupled with particle rearrangement in shear-thinning non-Brownian suspensions.

Dense non-Brownian suspensions exhibit significant shear thinning, although a comprehensive understanding of the full scope of this phenomenon remains elusive. This study numerically reveals intimate heterogenous coupled dynamics between many-body particle motions and solvent hydrodynamics in shear-thinning non-Brownian suspensions. In our simulation systems, we do not account for frictional contact forces, reflecting experimental conditions under low shear rates where shear thinning occurs, while hydrodynamic interactions are directly incorporated using the Smoothed Profile Method. We demonstrate the spatially correlated viscous dissipation and particle motions; they share the same characteristic length, which decreases with increasing shear rate. We further show that, at lower shear rates, significant particle density changes are induced against the incompressibility of the solvent, suggesting the cooperative creation and annihilation of gaps and flow channels. We discuss that hydrodynamic interactions may substantially restrict particle rearrangements even in highly dense suspensions, influencing the quantitative aspects of macroscopic rheology.

Dynamic driving eliminates volume fraction inhomogeneity and apparent yield stress in flowing dense non-Brownian suspensions

Dynamic driving eliminates volume fraction inhomogeneity and apparent yield stress in flowing dense non-Brownian suspensions

Porous structures impact on particle dynamics of non-Brownian and noncolloidal suspensions

This study aims to provide valuable insights into the impact of porous structures on particle dynamics in non-Brownian, non-colloidal suspension flows at very low Reynolds numbers. Two experimental approaches, Particle Image Velocimetry (PIV) with refractive index matching and Optical Flow Tracking Velocimetry (OFTV) were employed to analyze very dilute suspensions over various porous media models. The study considered three different porous structures with permeabilities ranging from 0.7 to 0.9 and three different thicknesses ranging from 0.2 cm to 0.5 cm, while the suspension bulk volume fraction was maintained at 3 %. In the PIV analysis, we observed that decreasing the porous permeability resulted in the maximum velocity location within the free flow region moving closer towards the interface between the flow and the porous media. We further quantified the effect of the porous structure on the suspension by characterizing interface properties, such as dimensionless slip velocity, shear rate, and slip length. These interface properties were found to be influenced by both the thickness and permeability of the porous media. Next, we analyzed particle migration due to the presence of porous structures using OFTV for very dilute suspensions of 1 %, 2 %, and 3 %, considering a porous medium with known physical properties and thickness. The study revealed two local concentration maxima: one within the free flow region on top of the rod arrays used to create the porous structure and a second along the rods' centerline inside the porous media model.

Effects of confinement-induced non-Newtonian lubrication forces on the rheology of a dense suspension

In this work, we propose a functionalised bi-viscous lubrication model to study the material properties of concentrated non-Brownian suspensions and explore the possible confinement-induced non-Newtonian effects of the lubricant in the rheological response of this type of suspensions. From tribological studies, it is well-known that even macroscopically Newtonian liquids under strong confinement might exhibit properties which deviate significantly from their bulk behaviour. When two surfaces separated by an extremely small gap (still large compared to the molecular size) are sheared, strong shear-thinning of the lubricant viscosity at low shear-rates is observed, in spite of its Newtonian-like bulk response. This is connected to a significant increase of the zero-shear-rate viscosity under extreme confinement. We start from an effective lubrication algorithm recently proposed and develop a new gap-size-dependent interparticle bi-viscous lubrication model, able to capture qualitatively the main phenomenology of confined lubricants. We solve the lubrication interaction between particles iteratively via a semi-implicit splitting scheme. Since the handling of lubrication is made implicitly here, the method copes efficiently with large increases of the inter-particle effective viscosities, which would otherwise lead to simulation blow-up or the use of vanishing time-steps in standard explicit schemes. We analyse the rheological response of the suspension systematically in terms of model parameters. In contrast to pure Newtonian lubrication interactions, distinct shear-thinning and shear-thickening regimes in the relative suspension viscosity are observed, which are discussed in terms of particle microstructure coupled with the complex rheology of the confined lubricant. In addition, normal-stress response is negative in both N1 and N2, which is difficult to achieve with standard contact frictional models.

Fresh cement as a frictional non-Brownian suspension

Cement is an essential construction material due to its ability to flow before later setting, however the rheological properties must be tightly controlled. Despite this, much understanding remains empirical. Using a combination of continuous and oscillatory shear flow, we compare fresh Portland cement suspensions to previous measurements on model non-Brownian suspensions to gain a micro-physical understanding. Comparing steady and small-amplitude oscillatory shear, we reveal two distinct jamming concentrations, ϕμ and ϕrcp, where the respective yield stresses diverge. As in model suspensions, the steady-shear jamming point is notably below the oscillatory jamming point, ϕμ<ϕrcp, suggesting that it is tied to frictional particle contacts. These results indicate that recently established models for the rheology of frictional, adhesive non-Brownian suspensions can be applied to fresh cement pastes, offering a new framework to understand the role of additives and fillers. Such micro-physical understanding can guide formulation changes to improve performance and reduce environmental impact.

Shear thickening and hysteresis in dense suspensions: The effect of particle shape

The flow of dense suspension of non-Brownian particles has been considered by various studies affected by their significance in a variety of industries and natural phenomena. In this study, we investigate the effect of polyhedron morphology on shear thickening, shear jamming, and hysteresis characteristics of non-Brownian suspension of acrylate particles. Particles with the same chemical nature and three different shapes of spherical (aspect ratio Γ = 1), elliptical paraboloid (Γ ≈ 1), and boat-shaped (Γ ≈ 3) are fabricated via photopolymerization-based methods. Studied suspensions show the shear-thinning behavior at low shear stresses and shear thickening behavior at the higher range of shear stress. Also, the strength of observed shear thickening is enhanced for the suspensions of polyhedron particles, which can be attributed to the heightened degree of interparticle frictional contacts. Furthermore, it is found that angularity not only shifts the predicted frictionless and frictional jamming packing fractions to lower values but also expands the shear jamming packing fraction range. Finally, a history-dependent hysteresis is observed in all samples due to the different particle spatial structures forming in ascending and descending flow modes. The observed hysteresis loops strongly depend on the volume fraction and diminish near the jamming packing fraction due to the restricted mobility space of particles. In addition, the tumbling of elongated particles also can decrease the hysteresis loop by enhancing viscosity in the ascending flow mode, where the structures are not fully developed.

Simulating the rheology of dense suspensions using pairwise formulation of contact, lubrication and Brownian forces

Dense suspensions of solid particles in viscous liquid are ubiquitous in both industry and nature, and there is a clear need for efficient numerical routines to simulate their rheology and microstructure. Particles of micron size present a particular challenge: at low shear rates, colloidal interactions control their dynamics while at high rates, granular-like contacts dominate. While there are established particle-based simulation schemes for large-scale non-Brownian suspensions using only pairwise lubrication and contact forces, common schemes for colloidal suspensions generally are more computationally costly and thus restricted to relatively small system sizes. Here, we present a minimal particle-based numerical model for dense colloidal suspensions that incorporates Brownian forces in pairwise form alongside contact and lubrication forces. We show that this scheme reproduces key features of dense suspension rheology near the colloidal-to-granular transition, including both shear thinning due to entropic forces at low rates and shear thickening at high rates due to contact formation. This scheme is implemented in LAMMPS, a widely used open source code for parallelised particle-based simulations, with a runtime that scales linearly with the number of particles, making it amenable for large-scale simulations.

Emergence of a hexagonal pattern in shear-thickening suspensions under orbital oscillations

A dense particle suspension under shear may lose its uniform state to large local density and stress fluctuations, which challenge the mean-field description of the system. Here, we explore the novel dynamics of a non-Brownian suspension under orbital oscillations, where localized density waves along the flow direction appear beyond an excitation frequency threshold and self-organize into a hexagonal pattern across the system. The spontaneous occurrence of the inhomogeneity pattern arises from a coupling between particle advection and the shear-thickening nature of the suspension. Through linear stability analysis, we show that they overcome the stabilizing effects of particle pressure at sufficient particle volume fraction and oscillation frequency. In addition, the long-standing density waves degenerate into random fluctuations when replacing the free surface with rigid confinement. It indicates that the shear-thickened state is intrinsically heterogeneous, and the boundary conditions are crucial for developing local disturbance.

Interparticle friction in sheared dense suspensions: Comparison of the viscous and frictional rheology descriptions

In the literature, two different frameworks exist for describing the rheology of solid/liquid suspensions: (1) the “viscous” framework in terms of the relative suspension viscosity, ηr, as a function of the reduced solid volume fraction, ϕ/ϕm, with ϕm the maximum flowable packing fraction, and (2) the “frictional” framework in terms of a macroscopic friction coefficient, μ, as a function of the viscous number, Iv, defined as the ratio of the viscous shear to the wall-normal particle stress. Our goal is to compare the two different frameworks, focusing on the effect of friction between particles. We have conducted a particle-resolved direct numerical simulation study of a dense non-Brownian suspension of neutrally buoyant spheres in slow plane Couette flow. We varied the bulk solid volume fraction from ϕb=0.1 to 0.6 and considered three different Coulomb friction coefficients: μc=0, 0.2, and 0.39. We find that ηr scales well with ϕ/ϕm, with ϕm obtained from fitting the Maron–Pierce correlation. We also find that μ scales well with Iv. Furthermore, we find a monotonic relation between ϕ/ϕm and Iv, which depends only weakly on μc. Since ηr=μ/Iv, we thus find that the two frameworks are largely equivalent and that both account implicitly for Coulomb friction. However, we find that the normal particle stress differences, N1 and N2, when normalized with the total shear stress and plotted against either ϕ/ϕm or Iv, remain explicitly dependent on μc in a manner that is not yet fully understood.

Hydrodynamic irreversibility of non-Brownian suspensions in highly confined duct flow

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

Two-Dimensional Crystals far from Equilibrium.

When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behavior with no equilibrium counterpart. We study a two-dimensional particle model initially proposed to describe driven non-Brownian suspensions undergoing nonequilibrium absorbing phase transitions. We show that when the transition occurs at large density, the dynamics produces long-range crystalline order. In the ordered phase, long-range translational order is observed because equipartition of energy is lacking, phonons are suppressed, and density fluctuations are hyperuniform. Our study offers an explicit microscopic model where nonequilibrium violations of the Mermin-Wagner theorem stabilize crystalline order in two dimensions.

Rheology of non-Brownian suspensions: a rough contact story

Rheology of non-Brownian suspensions: a rough contact story

Effect of DLVO interactions on the rheology and microstructure of non-Brownian suspensions

Effect of DLVO interactions on the rheology and microstructure of non-Brownian suspensions

Physics-informed neural networks for understanding shear migration of particles in viscous flow

We harness the physics-informed neural network (PINN) approach to extend the utility of phenomenological models for particle migration in shear flow. Specifically, we propose to constrain the neural network training via a model for the physics of shear-induced particle migration in suspensions. Then, we train the PINN against experimental data from the literature, showing that this approach provides both better fidelity to the experiments, and a novel understanding of the relative roles of the hypothesized migration fluxes. We first verify the PINN approach for solving the inverse problem of radial particle migration in a non-Brownian suspension in an annular Couette flow. In this classical case, the PINN yields the same value (as reported in the literature) for the ratio of the two parameters of the empirical model. Next, we apply the PINN approach to analyze experiments on particle migration in both non-Brownian and Brownian suspensions in Poiseuille slot flow, for which a definitive calibration of the phenomenological migration model has been lacking. Using the PINN approach, we identify the unknown/empirical parameters in the physical model through the inverse solver capability of PINNs. Specifically, the values are significantly different from those for the Couette cell, highlighting an inconsistency in the literature that uses the latter value for Poiseuille flow. Importantly, the PINN results also show that the inferred values of the empirical model’s parameters vary with the shear Péclet number and the particle bulk volume fraction of the suspension, instead of being constant as assumed in some previous literature.