Finite Shear Rate Research Articles

Basal boundary conditions for granular surface flows over fragile and brittle erodible beds

Many granular surface flows occur as shear flows of finite thickness, over erodible beds composed of the same granular material. Such beds may be fragile, and offer no more resistance to erosion than to sustained shear. Or they may be brittle, and offer instead an excess resistance to erosion. To take this contrast into account, new basal boundary conditions are proposed. Their implications for parallel flows down infinite slopes are then examined for three different cases: stationary flows; starting; and stopping transients. For all three cases, flow behaviour is altered significantly when beds present an excess resistance to erosion. For stationary flows, non-unique velocity profiles are obtained, implying hysteresis or history-dependence. For starting transients, a power law growth of the flow thickness is predicted, instead of the jump to finite or infinite depth that would otherwise occur. For stopping transients, flows start to decelerate with a finite basal shear rate, even over erodible substrates. Analytical solutions to the corresponding free and moving boundary problems are obtained, and checked against numerical results. Model predictions are then compared with experimental measurements. Overall, good agreement is obtained. In particular, the model describes well the very different erosional responses observed for fragile and brittle beds.

Discontinuous rigidity transition associated with shear jamming in granular simulations.

We investigate the rigidity transition associated with shear jamming in frictionless, as well as frictional, disk packings in the quasi-static regime and at low shear rates. For frictionless disks, the transition under quasi-static shear is discontinuous, with an instantaneous emergence of a system spanning rigid clusters at the jamming transition. For frictional systems, the transition appears continuous for finite shear rates, but becomes sharper for lower shear rates. In the quasi-static limit, it is discontinuous as in the frictionless case. Thus, our results show that the rigidity transition associated with shear jamming is discontinuous, as demonstrated in the past for isotropic jamming of frictionless particles, and therefore a unifying feature of the jamming transition in general.

Dynamic phase diagram of plastically deformed amorphous solids at finite temperature.

The yielding transition that occurs in amorphous solids under athermal quasistatic deformation has been the subject of many theoretical and computational studies. Here, we extend this analysis to include thermal effects at finite shear rate, focusing on how temperature alters avalanches. We derive a nonequilibrium phase diagram capturing how temperature and strain rate effects compete, when avalanches overlap, and whether finite-size effects dominate over temperature effects. The predictions are tested through simulations of an elastoplastic model in two dimensions and in a mean-field approximation. We find a scaling for temperature-dependent softening in the low-strain rate regime when avalanches do not overlap, and a temperature-dependent Herschel-Bulkley exponent in the high-strain rate regime when avalanches do overlap.

Viscuit and the fluctuation theorem investigation of shear viscosity by molecular dynamics simulations: The information and the noise.

The shear viscosity, η, of model liquids and solids is investigated within the framework of the viscuit and Fluctuation Theorem (FT) probability distribution function (PDF) theories, following Heyes et al. [J. Chem. Phys. 152, 194504 (2020)] using equilibrium molecular dynamics (MD) simulations on Lennard-Jones and Weeks-Chandler-Andersen model systems. The viscosity can be obtained in equilibrium MD simulation from the first moment of the viscuit PDF, which is shown for finite simulation lengths to give a less noisy plateau region than the Green-Kubo method. Two other formulas for the shear viscosity in terms of the viscuit and PDF analysis are also derived. A separation of the time-dependent average negative and positive viscuits extrapolated from the noise dominated region to zero time provides another route to η. The third method involves the relative number of positive and negative viscuits and their PDF standard deviations on the two sides for an equilibrium system. For the FT and finite shear rates, accurate analytic expressions for the relative number of positive to negative block average shear stresses is derived assuming a shifted Gaussian PDF, which is shown to agree well with non-equilibrium molecular dynamics simulations. A similar treatment of the positive and negative block average contributions to the viscosity is also shown to match the simulation data very well.

Brittle yielding of amorphous solids at finite shear rates

Amorphous solids display a ductile to brittle transition as the kinetic stability of the quiescent glass is increased, which leads to a material failure controlled by the sudden emergence of a macroscopic shear band in quasi-static protocols. We numerically study how finite deformation rates influence ductile and brittle yielding behaviors using model glasses in two and three spatial dimensions. We find that a finite shear rate systematically enhances the stress overshoot of poorly-annealed systems, without necessarily producing shear bands. For well-annealed systems, the non-equilibrium discontinuous yielding transition is smeared out by finite shear rates and it is accompanied by the emergence of multiple shear bands that have been also reported in metallic glass experiments. We show that the typical size of the bands and the distance between them increases algebraically with the inverse shear rate. We provide a dynamic scaling argument for the corresponding lengthscale, based on the competition between the deformation rate and the propagation time of the shear bands.

Scaling theory of shear-induced inhomogeneous dilation in granular matter.

Shearing with a finite shear rate a compressed granular system results in a region of grains flowing over a compact, static assembly. Perforce this region is dilated to a degree that depends on the shear rate, the loading pressure, gravity, various material parameters, and the preparation protocol. In spite of numerous studies of granular flows a predictive theory of the amount of dilation is still lacking. Here, we offer a scaling theory that is focused on such a prediction as a function of shear rate and the dissipative parameters of the granular assembly. The resulting scaling laws are universal with respect to changing the interparticle force laws.

Rheology of Inelastic Hard Spheres at Finite Density and Shear Rate.

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 σ to the shear rate γ[over ˙] and predictions of the yield stress complete our discussion of granular rheology derived from first principles.

Beyond linear elasticity: jammed solids at finite shear strain and rate.

The shear response of soft solids can be modeled with linear elasticity, provided the forcing is slow and weak. Both of these approximations must break down when the material loses rigidity, such as in foams and emulsions at their (un)jamming point - suggesting that the window of linear elastic response near jamming is exceedingly narrow. Yet precisely when and how this breakdown occurs remains unclear. To answer these questions, we perform computer simulations of stress relaxation and shear start-up tests in athermal soft sphere packings, the canonical model for jamming. By systematically varying the strain amplitude, strain rate, distance to jamming, and system size, we identify characteristic strain and time scales that quantify how and when the window of linear elasticity closes, and relate these scales to changes in the microscopic contact network.

Shearrate diffusion and constitutive relations during transients in simple shear

Granular matter, consisting of hard, frictional, cohesionless spheres, sheared in a simple shear geometry with smooth walls undergoes a velocity driven transition from a jammed or creeping state (low wall velocity) to a flow state with a finite shear rate in the bulk (high wall velocity). In the flow state, the state variables volume fraction \(\nu \), inertial number I and the macroscopic friction \(\mu \) of the bulk follow an exponential transient. The characteristic time of this progression grows with the wall velocity and the system size and is typically large compared to the inverse shear rate. It is shown that I, first being stationary in the shear zones, spreads diffusively into the bulk. The other state variables follow according to the constitutive laws, well known from the steady state.

A mesoscopic model for the rheology of soft amorphous solids, with application to microchannel flows

We have studied a mesoscopic model for the flow of amorphous solids. The model is based on key features identified at the microscopic level, namely periods of elastic deformation interspersed with localised rearrangements of particles that induce long-range elastic deformations. These long-range deformations are derived following a continuum mechanics approach, in the presence of solid boundaries, and are included in full in the model. Indeed, they mediate spatial cooperativity in the flow, whereby a localised rearrangement may lead a distant region to yield. In particular, we have simulated a channel flow and found manifestations of spatial cooperativity that are consistent with published experimental observations for concentrated emulsions in microchannels. Two categories of effects are distinguished. On the one hand, the coupling of regions subject to different shear rates, for instance, leads to finite shear rate fluctuations in the seemingly unsheared "plug" in the centre of the channel. On the other hand, there is convincing experimental evidence of a specific rheology near rough walls. We discuss the diverse possible physical origins for this effect, and we suggest that it may be associated with the bumps of particles into surface asperities as they slide along the wall.

Non-linear oscillatory rheological properties of a generic continuum foam model: Comparison with experiments and shear-banding predictions

The occurrence of shear bands in a complex fluid is generally understood as resulting from a structural evolution of the material under shear, which leads (from a theoretical perspective) to a non-monotonic stationary flow curve related to the coexistence of different states of the material under shear. In this paper we present a scenario for shear-banding in a particular class of complex fluids, namely foams and concentrated emulsions, which differs from other scenarios in two important ways. First, the appearance of shear bands is shown to be possible both without any intrinsic physical evolution of the material (e.g. via a parameter coupled to the flow such as concentration or entanglements) and without any finite critical shear rate below which the flow does not remain stationary and homogeneous. Secondly, the appearance of shear bands depends on the initial conditions, i.e. the preparation of the material. In other words, it is history dependent. This behaviour relies on the tensorial character of the underlying model (2D or 3D) and is triggered by an initially inhomogeneous strain distribution in the material. The shear rate displays a discontinuity at the band boundary whose amplitude is history dependent and thus depends on the sample preparation.

Plasticity and dynamical heterogeneity in driven glassy materials

Many amorphous glassy materials exhibit complex spatio-temporal mechanical response and rheology, characterized by an intermittent stress strain response and a fluctuating velocity profile. Under quasistatic and athermal deformation protocols this heterogeneous plastic flow was shown to be composed of plastic events of various sizes, ranging from local quadrupolar plastic rearrangements to system spanning shear bands. In this paper, through numerical study of a 2D Lennard-Jones amorphous solid, we generalize the study of the heterogeneous dynamics of glassy materials to the finite shear rate (gamma not equal to 0) and temperature case (T not equal to 0). In practice, we choose an effectively athermal limit (T approximately 0) and focus on the influence of shear rate on the rheology of the glass. In line with previous works we find that the model Lennard-Jones glass follows the rheological behavior of a yield stress fluid with a Herschel-Bulkley response of the form, sigma = sigmaY + c1gamma(beta). The global mechanical response obtained through the use of Molecular Dynamics is shown to converge in the limit gamma --> 0 to the quasistatic limit obtained with an energy minimization protocol. The detailed analysis of the plastic deformation at different shear rates shows that the glass follows different flow regimes. At sufficiently low shear rates the mechanical response reaches a shear-rate-independent regime that exhibits all the characteristics of the quasistatic response (finite-size effects, cascades of plastic rearrangements, yield stress, ...). At intermediate shear rates the rheological properties are determined by the externally applied shear rate and the response deviates from the quasistatic limit. Finally at higher shear the system reaches a shear-rate-independent homogeneous regime. The existence of these three regimes is also confirmed by the detailed analysis of the atomic motion. The computation of the four-point correlation function shows that the transition from the shear-rate-dominated to the quasistatic regime is accompanied by the growth of a dynamical cooperativity length scale xi that is shown to diverge with shear rate as xi is proportional to gamma(-nu), with nu approximately 0.2 -0.3. This scaling is compared with the prediction of a simple model that assumes the diffusive propagation of plastic events.

Dynamics of Fluid Vesicles in Oscillatory Shear Flow

The dynamics of fluid vesicles in oscillatory shear flow was studied using differential equations of two variables: the Taylor deformation parameter and inclination angle $\theta$. In a steady shear flow with a low viscosity $\eta_{\rm {in}}$ of internal fluid, the vesicles exhibit steady tank-treading motion with a constant inclination angle $\theta_0$. In the oscillatory flow with a low shear frequency, $\theta$ oscillates between $\pm \theta_0$ or around $\theta_0$ for zero or finite mean shear rate $\dot\gamma_{\rm m}$, respectively. As shear frequency $f_{\gamma}$ increases, the vesicle oscillation becomes delayed with respect to the shear oscillation, and the oscillation amplitude decreases. At high $f_{\gamma}$ with $\dot\gamma_{\rm m}=0$, another limit-cycle oscillation between $\theta_0-\pi$ and $-\theta_0$ is found to appear. In the steady flow, $\theta$ periodically rotates (tumbling) at high $\eta_{\rm {in}}$, and $\theta$ and the vesicle shape oscillate (swinging) at middle $\eta_{\rm {in}}$ and high shear rate. In the oscillatory flow, the coexistence of two or more limit-cycle oscillations can occur for low $f_{\gamma}$ in these phases. For the vesicle with a fixed shape, the angle $\theta$ rotates back to the original position after an oscillation period. However, it is found that a preferred angle can be induced by small thermal fluctuations.

Shear stresses of colloidal dispersions at the glass transition in equilibrium and in flow

We consider a model dense colloidal dispersion at the glass transition, and investigate the connection between equilibrium stress fluctuations, seen in linear shear moduli, and the shear stresses under strong flow conditions far from equilibrium, viz., flow curves for finite shear rates. To this purpose, thermosensitive core-shell particles consisting of a polystyrene core and a cross-linked poly(N-isopropylacrylamide) shell were synthesized. Data over an extended range in shear rates and frequencies are compared to theoretical results from integrations through transients and mode coupling approaches. The connection between nonlinear rheology and glass transition is clarified. While the theoretical models semiquantitatively fit the data taken in fluid states and the predominant elastic response of glass, a yet unaccounted dissipative mechanism is identified in glassy states.

Profile blunting and flow blockage in a yield-stress fluid: A molecular dynamics study

The flow of a simple glass forming system (a 80:20 binary Lennard-Jones mixture) through a planar channel is studied via molecular dynamics simulations. The flow is driven by an external body force similar to gravity. Previous studies show that the model exhibits both a static [F. Varnik, J. Chem. Phys. 120, 2788 (2004)] and a dynamic [F. Varnik and O. Henrich, Phys. Rev. B 73, 174209 (2006)] yield stress in the glassy phase. These observations are corroborated by the present work, where we investigate how the presence of a yield stress may affect the system behavior in a Poiseuille-type flow geometry. In particular, we observe a blunted velocity profile across the channel: A relatively wide region in the channel center flows with a constant velocity (zero shear rate) followed by a nonlinear change of the shear rate as the walls are approached. The observed velocity gradients are compared to those obtained from the knowledge of the shear stress across the channel and the flow curves (stress versus shear rate), the latter being determined in our previous simulations of homogeneous shear flow. Furthermore, using the value of the (dynamic) yield stress known from previous simulations, we estimate the threshold body force for a complete arrest of the flow. Indeed, a blockage is observed as the imposed force falls below this threshold value. Small but finite shear rates are observed at stresses above the dynamic but below the static yield stress. We discuss the possible role of the stick-slip-like motion for this observation.

On the orientation of lamellar block copolymer phases under shear

Self-assembled lamellar structures composed of block copolymers are simulated by molecular dynamics. The response of a bulk system to external shear is investigated, in particular, the average energy, the entropy production, and the stability of the lamellae's orientation. We distinguish two orientations, a parallel orientation in which the normal to the lamellae sheets lies in the direction of the shear gradient, and a perpendicular orientation in which the normal lies perpendicular to the shear gradient and shear direction. The perpendicular phase is stable throughout all shear rates. The parallel phase has higher internal energy and larger entropy production than the perpendicular phase and moreover becomes unstable at relatively small shear rates. The perpendicular orientation should therefore be more stable at any finite shear rate. Surface effects are probably responsible for the stability of the parallel phase observed experimentally at small shear rates.

The flow-phase diagram of Doi-Hess theory for sheared nematic polymers II: finite shear rates

The purpose of this paper is to extend the rheological predictions of the Doi-Hess kinetic theory for sheared nematic polymers from the anomalous weak shear regime (Forest et al. 2004a) to arbitrary shear rates, and to associate salient rheological and optical properties with the solution space of kinetic theory. Using numerical bifurcation software (AUTO), we provide the phase diagram of all stable monodomain orientational probability distribution functions (PDFs) and their phase transitions (bifurcations) vs nematic concentration (N) and normalized shear rate (Peclet number, Pe) for Pe≥1. Shear stresses, normal stress differences, the peak direction of the orientational distribution, and birefringence order parameters are calculated and illustrated for each type of PDF attractor: steady flow-aligning, both in and out of the flow deformation plane and along the vorticity axis; unsteady limit cycles, where the peak orientation direction rotates in-plane or around the vorticity axis or in bi-stable orbits tilted between them; and chaotic attractors first observed in kinetic simulations by Grosso et al. (2001). We pay particular attention to correlations between rheological features and the variety of monodomain phase transitions. Together with the weak flow regime, these results provide a nearly complete picture of the rheological consequences of the Doi-Hess kinetic theory for sheared monodomains of rigid, extreme aspect ratio, nematic rods or platelets.

Wall slip and yielding in pasty materials

We carried out systematic rheometrical tests under controlled stress with smooth and rough parallel disks, along with magnetic resonance imaging (MRI) tests in coaxial cylinder geometry, with foam and a model concentrated emulsion. At low shear stress wall slip appears to occur but the bulk fluid remains static, as proved by the fact that in this regime the apparent shear rate obtained for a given shear stress is inversely proportional to the gap between the disks. At high shear stress data with different surface types and gaps coincide, suggesting that wall slip is negligible in this regime. In parallel, MRI results show that, in contrast with the apparent, simple, yielding behavior observed in usual rheometry, there is an abrupt transition from a finite shear rate to a static one at critical stress. This critical shear rate precisely corresponds to the transition between the two regimes of slip. This suggests that different flow regimes occur with these materials: (1) at low stress, with smooth surfaces a layer of a different material is sheared along the solid surfaces whereas the rest of the fluid does not flow; with rough surfaces there is no flow; (2) beyond a critical stress, for both surface types, the bulk fluid starts to flow but the shear is localized in a thin layer; then the thickness of this layer increases when stress is applied; (3) for both surface types homogeneous flow is obtained only beyond slightly larger stress, which is associated to a critical, apparent shear rate.We carried out systematic rheometrical tests under controlled stress with smooth and rough parallel disks, along with magnetic resonance imaging (MRI) tests in coaxial cylinder geometry, with foam and a model concentrated emulsion. At low shear stress wall slip appears to occur but the bulk fluid remains static, as proved by the fact that in this regime the apparent shear rate obtained for a given shear stress is inversely proportional to the gap between the disks. At high shear stress data with different surface types and gaps coincide, suggesting that wall slip is negligible in this regime. In parallel, MRI results show that, in contrast with the apparent, simple, yielding behavior observed in usual rheometry, there is an abrupt transition from a finite shear rate to a static one at critical stress. This critical shear rate precisely corresponds to the transition between the two regimes of slip. This suggests that different flow regimes occur with these materials: (1) at low stress, with smooth surfaces...

Creep instability of the weak layer and natural slab avalanche triggerings

The weak layer connecting the snow slab with the older snow substrate is modelled as an open cell ice foam made of an array of ice bonds. The bond rupture and rewelding rates under load are derived using simple analytical equations. We show that the weak layer experiences stable flow for a wide range of loads, bond brittleness and bond rewelding rates. However, a critical point may be reached, corresponding to a finite shear rate of the weak layer, beyond which the shear rate increases catastrophically, leading to a natural avalanche release. This critical situation is equivalent to a ductile to britte transition, whose activation energy is one-half of that of diffusion in ice. A measurement of the exponent in the shear strain vs. stress relation may give a hint on the imminence of the avalanche release. A similar calculation performed at the tip of a pre-existing basal crack gives a simple criterion for crack stability and an analytical expression for snow fracture toughness.

Brownian dynamics simulations of the stress and molecular configuration of polymers in exponential and linearly-ramped shear flow

The rheological and optical properties of dilute polymer solutions during the startup and subsequent relaxation in exponential shear flow are studied using Brownian dynamics simulations of freely draining, flexible bead-rod chains. We find that exponential shear flow is capable of effecting large molecular deformation for finite shear rates and strains, and furthermore, the polymer undergoes an initial stress and index of refraction relaxation that is much faster than a single exponential decay upon cessation of flow. When we compare the first normal stress difference ( τ 11– τ 22) to the corresponding components of the index of refraction tensor, n 11– n 22, we find that there is hysteresis when plotted over the course of a full startup and relaxation experiment. This is very similar to the hysteresis originally discovered by Doyle et al. [J. Non-Newtonian Fluid Mech. 76 (1998) 79–110] for uniaxial extensional flow. Moreover, the shear components of the stress and index of refraction tensors ( τ 12 and n 12) also exhibit hysteretic behavior. However, whereas the hysteresis loop for the 11–22 component traverses up the left branch during startup and relaxes on the right branch, that for the 12 component is in the completely opposite direction, i.e. the right branch corresponds to startup and the left branch relaxation. The presence of stress-index of refraction hysteresis clearly has important implications for the stress-optic rule. We calculate the stress-optic coefficient for the startup process, and find that the normal and shear components give the same value, but that it is strongly dependent on Wi. This result is consistent with that found by Doyle et al. during the startup of steady shear flow [Doyle and Shaqfeh, Dynamic simulation of freely-draining, flexible bead-rod chains: startup of extensional and shear flow, J. Non-Newtonian Fluid Mech. 76 (1998) 43–78]. We also simulated shear flows with a linearly ramped shear rate to compare to our results for exponential shear. We find that for comparable final shear rate and time, both exponential and linear ramping shear produce very similar hysteresis effects in both the shear and normal components suggesting that there may be many unsteady shear flow capable of creating large polymer stretch and stress-birefringence hysteresis. Finally, we compare our bead-rod chain simulations to the FENE dumbbell model, and we find that while the latter can qualitatively capture the normal stress hysteresis, it predicts negligible shear stress hysteresis at equivalent shear rates and shear strains.