Rheology of shear thickening suspensions and the effects of wall slip in torsional flow

The rheological characterisation of concentrated shear thickening materials suspensions is challenging, as complicated and occasionally discontinuous rheograms are produced. Wall slip is often apparent and when combined with a shear thickening fluid the usual means of calculating rim shear stress in torsional flow is inaccurate due to a more complex flow field. As the flow is no longer “controlled”, a rheological model must be assumed and the wall boundary conditions are redefined to allow for slip. A technique is described where, by examining the angular velocity response in very low torque experiments, it is possible to indirectly measure the wall slip velocity. The suspension is then tested at higher applied torques and different rheometer gaps. The results are integrated numerically to produce shear stress and shear rate values. This enables the measurement of true suspension bulk flow properties and wall slip velocity, with simple rheological models describing the observed complex rheograms.     

A self-similar behavior for the relative viscosity of concentrated suspensions of rigid spheroids

A viscosity model for suspensions of rigid particles with predictive capability over a wide range of particle volume fraction and shear conditions is of interest to quantify the transport of suspensions in fluid flow models. We study the shear viscosity of suspensions and focus on the effect of particle aspect ratio and shear conditions on the rheological behavior of suspensions of rigid bi-axially symmetric ellipsoids (spheroids). We propose a framework that forms the basis to microscopically parameterize the evolution of the suspension microstructures and its effect on the shear viscosity of suspensions. We find that two state variables, the intrinsic viscosity in concentrated limit and the self-crowding factor, control the state of dispersion of the suspension. A combination of these two variables is shown to be invariant with the imposed shear stress (or shear rate) and depends only on the particle aspect ratio. This self-similar behavior, tested against available experimental and numerical data, allows us to derive a predictive model for the relative viscosity of concentrated suspensions of spheroids subjected to low (near zero) strain rates. At higher imposed strain rates, one needs to constrain one of the state variables independently to constrain the state of dispersion of the suspension and its shear dynamic viscosity. Alternatively, the obtained self-similar behavior provides the means to estimate the state variables from the viscosity measurements made in the laboratory, and to relate them to microstructure rearrangements and evolution occurring during deformation.     

Development of a new mixing rheometer for studying rheological behaviour of concentrated energetic suspensions

The overall objective is to present a procedure based on a Couette analogy to quantitatively analyse torque/rotor speed data and extract viscosity/shear-rate curves using a non-conventional geometry. Diphasic flows of energetic concentrated suspensions of melt-cast insensitive explosives exhibit particular rheological properties. The characterization of these complex fluids may be a challenging task when conventional rheometers are used. Placing these dense suspensions in a classic cylindrical geometry may lead to a partial destruction of the internal fluid structure. To prevent that, a “RheoXF” a mixer-type rheometer has been developed: it consists of a mixing device with quite a complex geometry rotating in a cylindrical tank. To evaluate the rheological constants (virtual radius, virtual shear rate and stress constants) of thus mixing rheometer, we used five Newtonian fluids. After this calibration, the rheological characterizations were carried out on five formulations. The unique parameter which changes in these formulations is the batch's origin of a secondary explosive: the 3-nitro-1,2,4-triazole-5-one. These energetic particles differ by their morphology, maximum packing density and may be by their process synthesis. After having determined pseudoplastic parameters, a correlation has been made with the evolution of maximum packing density values calculated with De Larrard model.     

Rheological properties of peanut butter

The rheological properties of two types of commercial peanut butter have been studied. Both products are concentrated suspensions, and differ by the presence of additives. The first type, referred to as “100% peanuts,” is an unstabilized suspension consisting of solid peanut particles in peanut oil which is a Newtonian fluid. The second type, referred to as “smooth,” consists of the same suspension stabilized with a vegetable oil and contains other ingredients such as salt and sugar in very small quantities. A mean volume particle diameter of 6.6 μm has been determined, the particle diameter distribution was found to be narrow, and the solids volume fraction was estimated to be 0.6. Slip encountered in rheometry was greatly reduced by gluing sandpaper to the parallel plates of the rheometer. Both samples behaved like plastic materials and apparent yield stresses of 24 Pa and 370 Pa have been determined for the unstabilized and the stabilized suspensions, respectively. No linear domain was found for both suspensions and the non-linearity was confirmed by deformed Lissajous curves and higher odd harmonics in the output signal of small amplitude oscillatory shear experiments. The stabilized suspension behaved more like a solid, the elastic modulus being larger than the loss modulus and almost independent of the frequency. This solid-like behavior is supposedly caused by strong repulsive (steric) forces induced by the stabilizing agent. Received: 29 September 1999 Accepted: 9 August 2000     

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.     

Shear rheometry and visualization of glass fiber suspensions

The nonlinear rheological behavior of short glass fiber suspensions has been investigated in this work by rotational rheometry and flow visualization. A Newtonian and a Boger fluid (BF) were used as suspending media. The suspensions exhibited shear thinning in the semidilute regime and weaker shear thinning in the transition to the concentrated one. Normal stresses and relative viscosity were higher for the BF suspensions than for the Newtonian ones presumably due to enhanced hydrodynamic interactions resulting from BF elasticity. In addition, relative viscosity of the suspensions increased rapidly with fiber content, suggesting that the rheological behavior in the concentrated regime is dominated by mechanical contacts between fibers. Visualization of individual fibers and their interactions under flow allowed the detection of aggregates, which arise from adhesive contacts. The orientation states of the fibers were quantified by a second order tensor and fast Fourier transforms of the flow field images. Fully oriented states occurred for shear rates around 20 s^− 1. Finally, the energy required to orient the fibers was higher in step forward than in reversal flow experiments due to a change in the spatial distribution of fibers, from isotropic to planar oriented, during the forward experiments.     

Experimental and conceptual problems in the rheological characterization of wheat flour doughs

Wheat flour dough is an industrially important material and a better understanding of its rheological behavior could have long ranging impact on the agricultural and the food processing industries. However, rheological characterization of dough is proving to be difficult due to a range of testing issues and anomalies in flow behavior. In a cone-and-plate rheometer wheat flour doughs “roll-out” of the gap before steady state viscosities can be established, as discussed by Bloksma and Nieman (1975). However, the mirror image of the transient viscosity-time plot obtained using a cone-and-plate viscometer has been used to obtain an estimate of steady shear viscosity behavior (Gleissle, 1975). To check this transient methodology for doughs, a second method, in addition to cone-and-plate transient flow, for determination of the shear viscosity, was needed. For this, capillary extrusion was chosen. Both a piston-driven and pressure driven capillary rheometer were employed. End corrections were determined to provide information on both the shear viscosity and, following Binding (1988), the extensional viscosity of the doughs. There are few data available on end corrections for doughs, though published data by Kieffer indicate that the corrections are unexpectedly very high. In this present work it was found that the end correction experiments were very difficult and imprecise in part due to the time-dependent nature of the doughs and difficulties in preparing replicate batches required to compare dies of differing L/R values. Further it was unexpectedly found that the samples, though prepared by normal mixing procedures to the “optimum” level, were so heterogeneous that large fluctuations in the pressure at constant output rate (in the piston-driven rheometer) and in output rate at constant pressure (in the pressure-driven instrument) were observed. These fluctuations could be eliminated by overmixing of the doughs, but overmixed doughs are of little practical interest. Although the problems encountered in this work were significant, it was encouraging that even these preliminary studies indicate that rheological measurements are effective in differentiating between spring and winter wheats. Defining a constitutive model for dough rheology still remains a major challenge, as results from one type of testing do not corroborate the findings from a different type of testing. Received: 19 May 1998 Accepted: 27 July 1998     

An experimental study of multimodal glass suspension rheology to test and validate a polydisperse suspension viscosity model

Experimental measurements of non-colloidal multimodal suspension viscosities are performed over a wide range of mixing ratios and used to test the robustness and predictive capability of a recent viscosity model (Mwasame et al. in Phys Fluids 28:061701, 2016b), subsequently referred to as the MWB model. Three unimodally distributed particle suspensions with narrow size distributions are blended to make the bimodal and trimodal suspensions used in the rheological experiments. We demonstrate how predictions for mixture viscosities can be made using the MWB model only requiring the volume-weighted average particle sizes and viscosity correlations of the individual unimodal suspensions comprising the blend. The resultant model predictions are found to be in good agreement with measured bimodal and trimodal viscosity data to within expected experimental uncertainty. The datasets provided here can be used to validate future modeling efforts, and the MWB model can be used to optimize the viscosity of multimodal suspension mixtures for specific performance criteria.     

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.     

Effects of polymer–fiber interactions on rheology and flow behavior of suspensions of semi-flexible fibers in polymeric liquids

Rheological properties of suspensions of fibers in polymeric fluids are influenced by fiber–polymer interactions. In this paper, we investigate this influence from both experimental and modeling standpoints. In the experimental part of this investigation, we have changed the fiber–polymer interactions by treating the surface of the fibers. The resulting effects are observed using scanning electron microscopy and dynamic mechanical analysis techniques and quantified from the measurements of the viscosity in the start-up of shear flows and dynamic tests in the linear viscoelastic range region. The results are interpreted with the help of a mesoscopic rheological model developed for suspensions of fibers in viscoelastic fluids.     

Bimodal model of slurry viscosity with application to coal-slurries. Part 1. Theory and experiment

The rheological behavior of stable slurries is shown to be characterized by a bimodal model that represents a slurry as made up of a coarse fraction and a colloidal size fine fraction. According to the model, the two fractions behave independently of each other, and the non-Newtonian behavior of the viscosity is solely caused by the colloidal fraction, while the coarse fraction increases the viscosity level through hydrodynamic interactions. Data from experiments run with colloidal coal particles of about 2–3 µm average size dispersed in water show the viscosity of these colloidal suspensions to exhibit a highly shearrate-dependent behavior and, in the high shear limit, to match very closely the viscosity of suspensions of uniform size rigid spheres although the coal volume fraction must be determined semi-empirically. Different amounts of coarse coal particles are added to the colloidal suspension and the viscosity of the truly bimodal slurries measured as a function of shear rate. In agreement with the bimodal model, the measured shear viscosities show the coarse fraction to behave independently of the colloidal fraction and its contribution to the viscosity rise to be independent of the shear rate. It is shown that the shear rate exerted on the colloidal fraction is higher than that applied by the viscometer as a result of hydrodynamic interactions between the coarse particles, and that it is this effective higher shear rate which is necessary to apply in the correlations. For determining the coal volume fraction a relatively simple and quite accurate measurement technique is developed for determining the density and void fraction of coarse porous particles; the technique directly relates volume fraction to mass fraction.     

Rheological characterization of a solder paste for surface mount applications

Solder pastes used in surface mount soldering techniques (SMT) are very complex suspensions containing high volumes of metallic powder in a carrier fluid. The rheological complexity results largely from the carrier fluid itself, which is a suspension of colloidal particles. In this work, we have characterized the rheological properties of a typical carrier fluid and its solder paste containing 64 vol.% metallic powder. A six-blade vane geometry was used to avoid wall slip and sample fracture. All measurements were carried out following pre-shearing and rest time in order to obtain reproducible results. Steady shear experiments showed that the solder paste was highly shear-thinning and thixotropic. In oscillatory shear, the linear viscoelastic domain was found to be very narrow for both the suspending fluid and the paste. Frequency sweep tests in the linear domain revealed a gel-like structure with a nearly constant G′ for the suspending fluid and a slightly increasing G′ for the solder paste. From creep experiments, a yield stress of about 40 Pa was determined for the suspending fluid at temperatures between 25 and 40°C, and of 100 Pa at 4°C. A much larger yield stress, 480 Pa, was determined for the solder paste at 25°C.     

Observations on the rheological response of alkali activated fly ash suspensions: the role of activator type and concentration

This paper reports the influence of activator type and concentration on the rheological properties of alkali-activated fly ash suspensions. A thorough investigation of the rheological influences (yield stress and plastic viscosity) of several activator parameters, including: (i) the cation type and concentration of alkali hydroxide and (ii) the alkali-to-binder ratio (n) and silica modulus (M_s), and (iii) the volume of the activation solution, on the suspension rheology is presented. The results indicate a strong dependence on the cation and its concentration in the activation solution. The viscosity of the activation solution and the volumetric solution-to-powder ratio are shown to most strongly influence the plastic viscosity of the suspension. The suspension yield stress is predominantly influenced by the changes in fly ash particle surface charge and the ionic species in the activator. A shift from non-Newtonian to Newtonian flow behavior is noted in the case of silicate-based suspensions for M_s?≤?1.5. This behavior, which is not observed at higher M_S values, or when the fly ash is dispersed in hydroxide solutions or pure water, is hypothesized to be caused by colloidal siliceous species present in this system, or surface charge effects on the fly ash particles. Comparisons of the rheological response of alkali-activated suspensions to that of portland cement-water suspensions are also reported.     

Rheological properties of suspensions of TiO2 particles in polymer solutions. 1. Shear viscosity

This paper presents results on the rheological behaviour of suspensions of two kinds of TiO₂ particles in two different polymer solutions. The particles differ in their hydrophilic or hydrophobic properties. The dispersing media are a solution of high molecular weight polyisobutylene in decalin and a solution of a low molecular weight polybutene in decalin. The concentrations of polymer are adjusted in order to get the same zero shear viscosity. The shear viscosity measurements display an apparent yield stress in some cases. The existence and the values of the yield stress depend on the volume fraction of solid particles and on the type of particles. The evolution of the intrinsic viscosity and of the maximum packing fraction vs the shear rate is interpreted in terms of evolution of the size and of the shape of aggregates of particles under shear. The effect of temperature on the development of the yield stress is also discussed. The results are completed by microscopic observations.     

Breakup of a capillary bridge of suspensions

The breakup of liquid bridges under the action of capillary forces is used for studying the rheology of suspensions under stretching. The experiments were performed with suspensions of finegrained (3–30 μm) sand in glycerin for sand volume fractions up to 0.465. The bridge thinning process was registered using an electro-optical measuring device and videofilming. The results were analyzed on the basis of a theory developed earlier for the thinning of a liquid bridge under the action of capillary forces. It is found that, for fairly slow stretching realized in the initial stage of the thinning, the rheological behavior of the suspensions considered agrees with the model of a Newtonian viscous fluid. Along with this, the measured effective viscosity of the suspension turned out to be approximately two-fold greater than the suspension viscosity under shear. The origin of this discrepancy is analyzed. With increase in the stretching rate, in the final stage of the thinning, the weakening of the suspension occurs, which is manifested in the formation of a local rapidly thinning neck in the bridge, similar to that observed in the breakup of plastic materials.     

Dynamics and rheology of a dilute suspension of dipolar nonspherical particles in an external field: Part 1. Steady shear flows

Suspensions of small nonspherical particles having dipolar moments exhibit non-Newtonian behavior under the influence of shear and external fields. Numerical methods are presented for calculating the rheological and rheo-optical properties of dilute suspensions of Brownian particles having permanent dipoles subject to time-dependent shear and external fields. The numerical methods employ the Galerkin method of weighted residuals to solve the differential equation for the particle orientation distribution function. The steady-state shear flow intrinsic viscosity of suspensions of particles with sufficiently extreme aspect ratio is predicted to exhibit a maximum value attained for intermediate shear rates at selected field orientations. These numerical results provide valuable insight into the coupling which occurs between the effects of rotary Brownian motion, the hydrodynamic resistance of nonspherical particles, and the external torque exerted on dipolar particles. The results are applicable to both suspensions of magnetic particles and electrically dipolar particles.     

An empirical model predicting the viscosity of highly concentrated, bimodal dispersions with colloidal interactions

The relationship between particle size distribution and viscosity of concentrated dispersions is of great industrial importance, since it is the key to get high solids dispersions or suspensions. The problem is treated here experimentally as well as theoretically for the special case of strongly interacting colloidal particles. An empirical model based on a generalized Quemada equation is used to describe η as a function of volume fraction for mono- as well as multimodal dispersions. The pre-factor η˜ accounts for the shear rate dependence of η and does not affect the shape of the η vs φ curves. It is shown here for the first time that colloidal interactions do not show up in the maximum packing parameter and φ_max can be calculated from the particle size distribution without further knowledge of the interactions among the suspended particles. On the other hand, the exponent ɛ is controlled by the interactions among the particles. Starting from a limiting value of 2 for non-interacting either colloidal or non-colloidal particles, ɛ generally increases strongly with decreasing particle size. For a given particle system it then can be expressed as a function of the number average particle diameter. As a consequence, the viscosity of bimodal dispersions varies not only with the size ratio of large to small particles, but also depends on the absolute particle size going through a minimum as the size ratio increases. Furthermore, the well-known viscosity minimum for bimodal dispersions with volumetric mixing ratios of around 30/70 of small to large particles is shown to vanish if colloidal interactions contribute significantly. Received: 7 June 2000/Accepted: 12 February 2001     

Wall slip of mayonnaises in viscometers

The shear flow of mayonnaise is generally characterized by an apparent yield stress, shear thinning in steady flow, stress overshoots upon inception of flow and other time-dependent effects. These observations are usually understood to be the result of structural rearrangement within the material. Additionally and separately, the possibility that emulsions may exhibit apparent wall slip on a microscopic scale at a solid-liquid boundary has been reported by some researchers. Thus, observed rheological behavior is likely to be the result of the interplay between these two phenomena. In the present work, it is demonstrated that when measurements are sought to be made on mayonnaise using rotational viscometers visible wall slip occurs, rendering such instruments ineffective for the purpose of making viscosity measurements even at shear rates as low as 10^–3s^–1. The factors that influence the onset and extent of slip are investigated with the help of parallel plate viscometers, and it is concluded that the observed “yielding” of mayonnaise is actually an artifact of the onset of macroscopic slip. Slip effects are also found in capillary flow but are ameliorated with increasing shear rate. To circumvent these problems, it is proposed that extensional viscometry be employed for determining the flow behavior of mayonnaises. Received: 18 August 1997 Accepted: 1 April 1998     

Rheological characterization of bimodal colloidal dispersions

In order to investigate the effect of the particle size distribution on the rheological properties of concentrated colloidal dispersions both steady-state shear and oscillatory measurements have been performed on well-characterized bimodal dispersions of sterically stabilized PMMA particles. Replacing a minor amount of large particles by small ones in a concentrated dispersion, keeping the total effective volume fraction constant, decreases the viscosity quite drastically. On the other hand, replacing a small amount of small particles by big ones hardly effects the viscosity at all. This behavior can be attributed to the deformability of the stabilizing polymer layer. A procedure is proposed to calculate the limiting viscosities in a bimodal colloidal dispersion starting from the characteristics of the monodisperse systems. A good agreement has been obtained between the calculated values and the experimental results. The linear viscoelastic properties of the concentrated dispersions have been investigated by means of oscillatory measurements. The plateau values of the storage modulus for the bimodal dispersions decrease with an increasing fraction of the coarse particles. By substituting the bimodal dispersion by an equivalent monodisperse system the storage modulus can be superimposed on the values for the monodisperse suspensions when plotted as a function of the mean interparticle distance.