No Arabic abstract
We find in complementary experiments and event driven simulations of sheared inelastic hard spheres that the velocity autocorrelation function $psi(t)$ decays much faster than $t^{-3/2}$ obtained for a fluid of elastic spheres at equilibrium. Particle displacements are measured in experiments inside a gravity driven flow sheared by a rough wall. The average packing fraction obtained in the experiments is 0.59, and the packing fraction in the simulations is varied between 0.5 and 0.59. The motion is observed to be diffusive over long times except in experiments where there is layering of particles parallel to boundaries, and diffusion is inhibited between layers. Regardless, a rapid decay of $psi(t)$ is observed, indicating that this is a feature of the sheared dissipative fluid, and is independent of the details of the relative particle arrangements. An important implication of our study is that the non-analytic contribution to the shear stress may not be present in a sheared inelastic fluid, leading to a wider range of applicability of kinetic theory approaches to dense granular matter.
Considering a granular fluid of inelastic smooth hard spheres we discuss the conditions delineating the rheological regimes comprising Newtonian, Bagnoldian, shear thinning, and shear thickening behavior. Developing a kinetic theory, valid at finite shear rates and densities around the glass transition density, we predict the viscosity and Bagnold coefficient at practically relevant values of the control parameters. The determination of full flow curves relating the shear stress $sigma$ to the shear rate $dotgamma$, and predictions of the yield stress complete our discussion of granular rheology derived from first principles.
We present molecular dynamics simulations of pseudo hard sphere fluid (generalized WCA potential with exponents (50, 49) proposed by Jover et al. J. Chem. Phys 137, (2012)) using GROMACS package. The equation of state and radial distribution functions at contact are obtained from simulations and compared to the available theory of true hard spheres (HS) and available data on pseudo hard spheres. The comparison shows agreements with data by Jover et al. and the Carnahan-Starling equation of HS. The shear viscosity is obtained from the simulations and compared to the Enskog expression and previous HS simulations. It is demonstrated that the PHS potential reproduces the HS shear viscosity accurately.
Non-Newtonian transport properties of an inertial suspension of inelastic rough hard spheres under simple shear flow are determined from the Boltzmann kinetic equation. The influence of the interstitial gas on rough hard spheres is modeled via a Fokker-Planck generalized equation for rotating spheres accounting for the coupling of both the translational and rotational degrees of freedom of grains with the background viscous gas. The generalized Fokker-Planck term is the sum of two ordinary Fokker-Planck differential operators in linear $mathbf{v}$ and angular $boldsymbol{omega}$ velocity space. As usual, each Fokker-Planck operator is constituted by a drag force term (proportional to $mathbf{v}$ and/or $boldsymbol{omega}$) plus a stochastic Langevin term defined in terms of the background temperature $T_text{ex}$. The Boltzmann equation is solved by two different but complementary approaches: (i) by means of Grads moment method, and (ii) by using a Bhatnagar-Gross-Krook (BGK)-type kinetic model adapted to inelastic rough hard spheres. As occurs in the case of emph{smooth} inelastic hard spheres, our results show that both the temperature and the non-Newtonian viscosity increase drastically with increasing the shear rate (discontinuous shear thickening effect) while the fourth-degree velocity moments also exhibit an $S$-shape. In particular, while high levels of roughness may slightly attenuate the jump of the viscosity in comparison to the smooth case, the opposite happens for the rotational temperature. As an application of these results, a linear stability analysis of the steady simple shear flow solution is also carried out showing that there are regions of the parameter space where the steady solution becomes linearly unstable.
Dense suspensions of model hard-sphere-like colloids, with different particle sizes, are examined experimentally in the glass state, under shear and extensional rheology. Under steady shear flow we detect Discontinuous Shear Thickening (DST) above a critical shear rate. Start-up shear experiments show stress overshoots in the vicinity of the onset of DST related with a change in microscopic morphology, as the sample shows dilatancy effects. The analysis of the normal stress together with direct sample observation by high speed camera, indicates the appearance of positive N1 and dilation behavior at the shear thickening onset. Dilatancy effects are detected also under extensional flow. The latter was studied through capillary breakup and filament stretching experimental setups, where liquid-like response is seen for strain rate lower than a critical strain rate and solid like-behavior for higher strain rates. Monitoring the filament thinning processes under different conditions (volume fractions and strain rates) we have created a state diagram where all responses of a hard-sphere suspension (Newtonian, shear thinning, shear thickening, dilatant) are shown. We finally compare the shear thickening response of these hard-sphere-like suspensions and glasses in shear with that in extensional flow.
The Boltzmann equation for inelastic Maxwell models is considered to determine the velocity moments through fourth degree in the simple shear flow state. First, the rheological properties (which are related to the second-degree velocity moments) are {em exactly} evaluated in terms of the coefficient of restitution $alpha$ and the (reduced) shear rate $a^*$. For a given value of $alpha$, the above transport properties decrease with increasing shear rate. Moreover, as expected, the third-degree and the asymmetric fourth-degree moments vanish in the long time limit when they are scaled with the thermal speed. On the other hand, as in the case of elastic collisions, our results show that, for a given value of $alpha$, the scaled symmetric fourth-degree moments diverge in time for shear rates larger than a certain critical value $a_c^*(alpha)$ which decreases with increasing dissipation. The explicit shear-rate dependence of the fourth-degree moments below this critical value is also obtained.