Strain stiffening of non-colloidal hard sphere suspensions dispersed in Newtonian fluid near liquid-and-crystal coexistence region

Concentrated hard sphere suspensions often show an interesting nonlinear behavior, called strain stiffening, in which the viscosity or modulus starts to increase at critical strain amplitude. Sudden increase of rheological properties is similar to shear thickening; however, the particle dynamics in the strain stiffening under oscillatory shear flow does not necessarily coincide with the mechanism of shear thickening under step shear flow. In this study, we have systematically investigated the nonlinear rheology of non-colloidal (>1???m) hard sphere suspensions dispersed in Newtonian fluid near liquid-and-crystal coexistence region in order to better understand the strain stiffening behavior. The suspensions near liquid-and-crystal coexistence region are known to locally form the closed packing structure. The critical strain amplitude which is the onset of strain stiffening was different for the storage and loss modulus. But they converged to each other as the suspension forms a more crystalline structure. The critical strain amplitude was independent of medium viscosity, imposed angular frequency, and particle size, but was strongly dependent upon particle volume fraction. The onset of strain stiffening was explained in terms of shear-induced collision due to particle motion in the closed packing structure. Nonlinear stress wave-forms, which reflect the micro-structural change, were observed with the onset of strain stiffening. During the strain stiffening, enhanced elastic stress before and after flow reversal was observed which originates from changes in the suspension microstructure. Nonlinearity of the shear stress in terms of Fourier intensity was extremely increased up to 0.55. Beyond the strain stiffening, the suspension responded liquid-like and the nonlinearity decreased but the elastic shear stress was still indicating the microstructure rearrangement within a cycle.     

Dynamic properties of shear thickening colloidal suspensions

The transient shear rheology (i.e., frequency and strain dependence) is compared to the steady rheology for a model colloidal dispersion through the shear thickening transition. Reversible shear thickening is observed and the transition stress compares well to theoretical predictions. Steady and transient shear thickening are observed to occur at the same value of the average stress. The critical strain for shear thickening is found to depend inversely on the frequency at fixed applied stress for low frequencies (high strains), but is limited to an apparent minimum critical strain at higher frequencies. This minimum critical strain is shown to be an artifact of slip. Lissajous plots illustrate the transition in material properties through the shear thickening transition, and the energy dissipated by a shear thickening suspension is analyzed as a function of strain amplitude.     

Microstructure of shear-thickening concentrated suspensions determined by flow-USANS

Reversible shear thickening in colloidal suspensions is a consequence of the formation of hydroclusters due to the dominance of short-ranged lubrication hydrodynamic interactions at relatively high shear rates. Here, we develop and demonstrate a new method of flow-ultra small angle neutron scattering to probe the colloidal microstructure under steady flow conditions on length scales suitable to characterize the formation of hydroclusters. Results are presented for a model near hard-sphere colloidal suspension of 260 nm radius (10% polydisperse) sterically stabilized silica particles in poly(ethylene glycol) at shear rates in the shear thinning and shear thickening regime for dilute, moderately concentrated, and concentrated (ordered) suspensions. Hydrocluster formation is observed as correlated, broadly distributed density fluctuations in the suspension with a characteristic length scale of a few particle diameters. An order–disorder transition is observed to be coincident with shear thickening for the most concentrated sample, but the shear-thickened state shows hydrocluster formation. These structural observations are correlated to the behavior of the shear viscosity and discussed within the framework of theory, simulation, and prior experiments.     

Investigating the transient response of a shear thickening fluid using the split Hopkinson pressure bar technique

The split Hopkinson pressure bar (SHPB) technique is implemented to evaluate the transient response of a colloidal suspension exhibiting shear thickening at strain rates and timescales never before explored in a laboratory instrument. These suspensions are shown to exhibit a discontinuous transition from fluid-like (shear thinning) to solid-like (shear thickening) behavior when evaluated using rotational rheometry. The effect of loading rate on this transition time is studied for a particle volume fraction of 0.54 using the SHPB technique. It is shown that the time required for transition to occur decreases logarithmically with loading rate. From these results, we conclude that transition is not triggered by a characteristic shear rate, but rather a critical shear strain is required. Results from SHPB experiments performed up to Peclet numbers of order 10⁷ are presented and discussed for 0.50, 0.52, and 0.54 particle volume fraction suspensions.     

Rheological properties of concentrated latex suspensions of poly(styrene-butadiene)

Concentrated suspensions of charged latex particles of poly(styrene-butadiene) have been used as model systems to investigate the influence of surface charges on the rheology of colloidal suspensions. The suspensions were found to behave as elastic solids at small strains and to require a finite stress to flow. This was related to an ordered structure of the suspensions at rest, resulting from electrostatic and van der Waals forces. Important shear-thinning effects were observed as a consequence of structure rearrangements under shear. At a fixed shear rate, the steady-shear viscosity as a function of the ionic strength exhibits a minimum. Under oscillatory shear flow, the behavior of the concentrated suspensions was found to be non-linear above a very small strain amplitude. The non-linear output signal from dynamic experiments was analyzed using a fast Fourier transform algorithm. A maximum in the third harmonic intensity as a function of the strain amplitude was observed and the intensity of higher harmonics decreased with increasing ionic strength. The behavior of the suspensions could be adequately described using the structural model of Yziquel et al. (Yziquel F, Carreau PJ, Moan M, Tanguy PA (1999) Rheological modeling of concentrated colloidal suspensions. J Non-Newtonian Fluid Mech 86:133–155).     

Strain-thickening transition in ferric-oxide suspensions under oscillatory shear

Shear-strain-thickening transition under oscillatory flow was observed in flocculated ferric-oxide suspensions in mineral oil. The value of the dynamic modulus of the suspensions that was measured at small strain amplitude after cessation of shear also became higher when the strain amplitude of the applied shear had been within or above the transition region.The ferric-oxide powders used were an acicular submicron maghemite (magnetic) and the hematite (non-magnetic) that was converted from the maghemite by heat treatment. The powders were treated with a dispersing agent and the suspensions were prepared in 33% by particle weight. The strain-thickening transition was observed in both the magnetic and the non-magnetic suspensions. However, the onset of the strain-thickening in the magnetic suspension was found at about one decade larger strain amplitude than that in the non-magnetic analog suspension, indicating particle interactions affect to the appearance of the phenomenon.A qualitative interpretation was made in view of site percolation for the enhancement of modulus at rest after the application of the large-amplitude oscillatory shear, where the process of the strain-thickening transition under shear and the development of the modulus after stopping the shear was described with a floc model in which the flocculation phase dilates as a result of the reduction of the particle linkages under higher shear.     

Shear-thickening flow of suspensions of carbon nanofibers in aqueous PVA solutions

The suspensions of carbon nanofibers in aqueous poly(vinyl alcohol) solutions were prepared in the presence of spherical carbon black particles, and the steady-shear viscosity and dynamic viscoelasticity were measured for complex suspensions. Although the single suspensions of carbon black are highly stable, the flocculation of carbon nanofibers is promoted by the addition of carbon black particles. The complex suspensions show remarkable shear thickening in the steady-flow and strain hardening in oscillatory shear with large amplitude. The nonlinear responses strongly depend on the carbon black concentration, whereas the dynamic viscoelasticity at low strains in the linear ranges is not significantly influenced. As the highly elastic effects arise from the long-range motion of particles, the possible mechanism may be the orientation of nanofibers in strong shear fields. The suspensions show the time-dependent behavior of viscosity when the time-scale of measurements is shorter than that of structural recovery to the isotropic states.     

Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Modeling adiabatic shear failure from energy considerations

The numerical simulation of dynamic structural failure by localized shear is quite complex in terms of constitutive models and choice of adequate failure criteria, along with a pronounced mesh-sensitivity. As a result, the existing numerical procedures are usually quite sophisticated, so that their application for design purposes is still limited. This study is based on the implementation of a simple energy-based criterion, which was developed on experimental considerations (Rittel et al., 2006), and uses a minimal number of adjustable parameters. According to this criterion, a material point starts to fail when the total strain energy density reaches a critical value. Thereafter, the strength of the element decreases gradually to zero to mimic the actual structural behavior. The criterion was embedded into commercial finite element software and tested by simulating numerically four typical high-rate experiments. The first is the dynamic torsion test of a tubular specimen. The second concerns the failure mode transition in mode II fracture of an edge crack in plain strain. The last two involve dynamic shear localization under high rate compression of a cylindrical and a shear compression specimen. A very good adequation was found both qualitatively and quantitatively. Qualitatively, in terms of failure path selection, and quantitatively, in terms of local strains, temperatures and critical impact velocity. The proposed approach is enticing from an engineering perspective aimed at predicting the onset and propagation of dynamic shear localization in actual structures.     

The impact of non-DLVO forces on the onset of shear thickening of concentrated electrically stabilized suspensions

This paper exposes an extension of an activation model previously published by the authors. When particles arranged along the compression axis of a sheared suspension, they may overcome the electrostatic repulsion and form force chains associated with shear thickening. A percolation-based consideration allows an estimation of the impact of the force chains on a flowing suspension. It suggests that similar to mode coupling models, the suspension becomes unstable before the critical stress evaluated from the activation model is reached. The percolated force chains lead to discontinuous shear thickening. The model predictions are compared with results from two experimental studies on aqueous suspensions of inorganic oxides; in one of them, hydration repulsion and in the other hydrophobic attraction can be expected. It is shown that the incorporation of non-Derjaguin–Landau–Verwey–Overbeek forces greatly improve predictions of the shear thickening instability.     

Elastic properties of polyethylene melts at high shear rates with respect to extrusion

At high shear rates a steady state of shear flow with constant shear rate, constant shear stress, and constant recoverable shear strain is observed in the short-time sandwich rheometer after some few shear units already. The melt exhibits rather high elastic shear deformations and the recovery occurs at much higher speed than it is observed in the newtonian range. The ratio of first normal stress difference and twice the shear stress, being equal to the recoverable strain in the second-order fluid limit, significantly underestimates the true elastic shear strains at high shear rates. The observed shear rate dependence of shear stress and first normal stress difference as well as of the (constrained) elastic shear strain is correctly described on the basis of a discrete relaxation time spectrum. In simple shear a stick-slip transition at the metal walls is found. Necessary for the onset of slip is a critical value of shear stress and a certain amount of elastic shear deformation or orientation of the melt.     

Non-linear viscoelastic response of magnetic fiber suspensions in oscillatory shear

This paper reports the first study on the large amplitude oscillatory shear flow for magnetic fiber suspensions subject to a magnetic field perpendicular to the flow. The suspensions used in our experiments consisted of cobalt microfibers of the average length of 37 μm and diameter of 4.9 μm, dispersed in a silicon oil. Rheological measurements have been carried out at imposed stress using a controlled stress magnetorheometer. The stress dependence of the shear moduli presented a staircase-like decrease with, at least, two viscoelastic quasi-plateaus corresponding to the onset of microscopic and macroscopic scale rearrangement of the suspension structure, respectively. The frequency behavior of the shear moduli followed a power-law trend at low frequencies and the storage modulus showed a high-frequency plateau, typical for Maxwell behavior. Our simple single relaxation time model fitted reasonably well the rheological data. To explain a relatively high viscous response of the fiber suspension, we supposed a coexistence of percolating and pivoting aggregates. Our simulations revealed that the former became unstable beyond some critical stress and broke in their middle part. At high stresses, the free aggregates were progressively destroyed by shear forces that contributed to a drastic decrease of the moduli. We have also measured and predicted the output strain waveforms and stress–strain hysteresis loops. With the growing stress, the shape of the stress–strain loops changed progressively from near-ellipsoidal one to the rounded-end rectangular one due to a progressive transition from a linear viscoelastic to a viscoplastic Bingham-like behavior.     

In situ monitoring the sol–gel transition for polyacrylamide gel

The effects of heating rate and shear stress on the sol–gel transition for a polyacrylamide gel during in situ preparation were investigated by dynamic rheology. The gelation evolution was also studied through monitoring of the static gelation process by means of temperature, pH, and conductivity variation. The gels were prepared by polymerization of acrylamide crosslinked with N,N′-methylenebisacrylamide using a redox initiation based on potassium persulfate/ascorbic acid. In situ gelation process was studied using oscillatory deformation tests at constant frequency of 0.1 Hz for different heating rate (from 0.5 to 5 °C/min) and shear stress (in the range of 0.1 to 10 Pa). The paper was presented at the third annual rheology conference, AERC 2006, Crete, Greece, 27–29 April 2006.     

Rheology and flow-birefringence from viscoelastic polymer-clay solutions

 The shear orientation of viscoelastic clay-polymer solutions was investigated by means of rheology and flow birefringence (Δn). The polymer chains are in dynamic adsorption/desorption equilibrium with the clay particles to form a “network”. The elastic behavior of the network was characterized by constant stress, oscillatory shear, and stress relaxation experiments. Constant stress experiments indicated a yield stress upon which shear flow started and no strain recovery could be observed. Oscillatory shear experiments showed a broad elastic region followed by flow when a critical strain was reached. Stress relaxation experiments showed several relaxation times when the same critical strain was reached. Experiments under steady flow characterized the transient behavior of the network. With increasing steady shear rate a pronounced minimum in birefringence was observed at a critical shear rate. The shear rate dependent viscosity showed near power law behavior and no corresponding critical feature. While birefringence detects orientational effects on a microscopic length scale, rheology averages over macroscopic changes in the sample. The same degree of orientation could be achieved under constant shear rate or constant stress conditions. Received: 25 January 2001 Accepted: 22 May 2001     

A microstructural investigation of the nonlinear response of electrorheological suspensions

In the preceding paper, Part 1, the transition from linear to nonlinear behavior for electrorheological (ER) suspensions under start-up of steady shear flow was found to first arise from the slight rearrangement of unstable structures. In this paper, we investigate the transition to nonlinear behavior for ER suspensions under oscillatory shear flow, focusing on the role of the rearrangement of unstable structures, and employing experimental and simulation results. Again, we find that nonlinear deformation first arises from these rearrangements, as opposed to the gross rearrangement or rupture of particulate chains. The Fourier transform of the simulated time-dependent shear stress is employed to quantify the dependence of the critical strain on the deformation frequency and electric field strength. The predicted behavior is consistent with experimental trends. Methods for verifying the predictions are discussed, as well as possible avenues for exploiting this information in improved operating strategies and improved ER fluids.     

Shear-thickening behavior of high molecular weight poly(ethylene oxide) solutions

Simple shear rheological properties of solutions of a high molecular weight (8 × 10⁶ g/mol) poly(ethylene oxide) (PEO) and its mixtures with sodium dodecyl sulfate (SDS) have been studied. Shear-thickening effects set in at a critical shear rate for PEO solutions. This particular behavior has not been reported for aqueous solutions of PEO, to our knowledge. The effect is attributed to PEO flow-induced self-aggregation. The experiments were performed in different operation modes (strain rate and stress controlled) and with different geometries (double wall Couette and Couette) and identical viscosities were obtained, which rules out flow instabilities as possible cause for the shear-thickening effect. Shear thickening was observed in the temperature range 15–50°C. Flow-induced PEO degradation occurs for shear rates in the shear-thickening regime, which indicates substantial chain deformation and accumulated stresses in the molecule when shear thickening occurs. Addition of SDS to the PEO solutions induces the formation of surfactant polymer complexes that preserve the characteristic shear-thickening effect.     

Shear thickening in electrically-stabilized colloidal suspensions

A theory is presented for the onset of shear thickening in colloidal suspensions of particles, stabilized by an electrostatic repulsion. Based on an activation model, a critical shear stress can be derived for the onset of shear thickening in dense suspensions for a constant potential and a constant charge approach of the spheres. Unlike previous models, the total interaction potential is taken into account (sum of attraction and repulsion). The critical shear stress is related to the maximum of the total interaction potential scaled by the free volume per particle. A comparison with experimental investigations shows the applicability of the theory.     

Extensional rheology of a shear-thickening cornstarch and water suspension

A filament-stretching rheometer is used to measure the extensional viscosity of a shear-thickening suspension of cornstarch in water. The experiments are performed at a concentration of 55 wt.%. The shear rheology of these suspensions demonstrates a strong shear-thickening behavior. The extensional rheology of the suspensions demonstrates a Newtonian response at low extension rates. At moderate strain rates, the fluid strain hardens. The speed of the strain hardening and the extensional viscosity achieved increase quickly with increasing extension rate. Above a critical extension rate, the extensional viscosity goes through a maximum and the fluid filaments fail through a brittle fracture at a constant tensile stress. The glassy response of the suspension is likely the result of jamming of particles or clusters of particles at these high extension rates. This same mechanism is responsible for the shear thickening of these suspensions. In capillary breakup extensional rheometry, measurement of these suspensions demonstrates a divergence in the extensional viscosity as the fluid stops draining after a modest strain is accumulated.     

Reological behaviour of concentrated monodisperse suspensions as a function of preshear conditions and temperature: an experimental study

 The influence of preshearing on the rheological behaviour of model suspensions was investigated with a stress-controlled cone-and-plate rheometer. The used matrix fluids showed Newtonian behaviour over the whole range of applied shear stresses. Highly monodisperse spherical glass spheres with various particle diameters were used as fillers. By applying steady preshearing at a low preshear stress, where a diffusion of particles can be expected, it was found for all model suspensions investigated at volume fractions ranging from 0.20 to 0.35 that the time-temperature superposition in the steady shear and in the dynamic mode holds within the chosen temperature range. Furthermore, all presheared model suspensions displayed a high and a low frequency range which are either separated by a shoulder or by a plateau value of G′ at intermediate frequencies. It could clearly be demonstrated that the low frequency range strongly depends on the preshear conditions. Hence, the features observed in the low frequency range can be attributed to a structure formation of a particulate network. In the high frequency range a frequency-dependent behaviour was observed which obeys the classical behaviour of Newtonian fluids (G′∝ω², G′′∝ω). The resulting temperature shift factors from the dynamic and the steady shear mode are identical and independent of the volume fraction and the particle size of the filler. Received: 29 November 2000 Accepted: 12 July 2001     

A Study of Large Plastic Deformations in Dual Phase Steel Using Digital Image Correlation and FE Analysis

Large plastic deformation in sheets made of dual phase steel DP800 is studied experimentally and numerically. Shear testing is applied to obtain large plastic strains in sheet metals without strain localisation. In the experiments, full-field displacement measurements are carried out by means of digital image correlation, and based on these measurements the strain field of the deformed specimen is calculated. In the numerical analyses, an elastoplastic constitutive model with isotropic hardening and the Cockcroft–Latham fracture criterion is adopted to predict the observed behaviour. The strain hardening parameters are obtained from a standard uniaxial tensile test for small and moderate strains, while the shear test is used to determine the strain hardening for large strains and to calibrate the fracture criterion. Finite Element (FE) calculations with shell and brick elements are performed using the non-linear FE code LS–DYNA. The local strains in the shear zone and the nominal shear stress-elongation characteristics obtained by experiments and FE simulations are compared, and, in general, good agreement is obtained. It is demonstrated how the strain hardening at large strains and the Cockcroft–Latham fracture criterion can be calibrated from the in-plane shear test with the aid of non-linear FE analyses. An erratum to this article can be found at