Thixotropy, yielding and ultrasonic Doppler velocimetry in pulp fibre suspensions

This paper reports thixotropy in concentrated pulp fibre suspensions and studies their transient flow behaviour using conventional rheometry coupled with a velocimetry technique. Specifically, an ultrasonic Doppler velocimeter is used in conjunction with a rate-controlled rheometer to deduce the local velocity profiles of pulp fibre suspensions. Pulp suspensions are found to exhibit a plateau in their flow curves where a slight increase in the shear stress generates a jump in the corresponding shear rate, implying the occurrence of shear banding. The velocity profiles were found to be discontinuous in the vicinity of the yielding radius where the Herschel–Bulkley model failed to predict the flow. Shear history and the time of rest prior to the measurement were found to play a significant role on the rheology and the local velocity profiles of pulp suspensions.     

Rheological characterization and constitutive modeling of bread dough

Bread dough (a flour–water system) has been rheologically characterized using a parallel-plate, an extensional, and a capillary rheometer at room temperature. Based on the linear and nonlinear viscoelastic and viscoplastic data, two constitutive equations have been applied, namely a viscoplastic Herschel–Bulkley model and a viscoelastoplastic K–BKZ model with a yield stress. For cases where time effects are unimportant, the viscoplastic Herschel–Bulkley model can be used. For cases where transient effects are important, it is more appropriate to use the K-BKZ model with the addition of a yield stress. Finally, the wall slip behavior of dough was studied in capillary flow, and an appropriate slip law was formulated. These models characterize the rheological behavior of bread dough and constitute the basic ingredients for flow simulation of dough processing, such as extrusion, calendering, and rolling.     

Application of Rheological Models in Prediction of Turbulent Slurry Flow

The paper deals with fully developed steady turbulent flow of slurry in a circular straight and smooth pipe. The Kaolin slurry consists of very fine solid particles, so the solid particles concentration, and density, and viscosity are assumed to be constant across the pipe. The mathematical model is based on the time averaged momentum equation. The problem of closure was solved by the Launder and Sharma k-ε turbulence model (Launder and Sharma, Lett Heat Mass Transf 1:131–138, 1974) but with a different turbulence damping function. The turbulence damping function, used in the mathematical model in the present paper, is that proposed by Bartosik (1997). The mathematical model uses the apparent viscosity concept and the apparent viscosity was calculated using two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley. The main aim of the paper is to compare measurements and predictions of the frictional head loss and velocity distribution, taking into account two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley, if the Kaolin slurry possesses low, moderate, and high yield stress. Predictions compared with measurements show an observable advantage of the Herschel–Bulkley rheological model over the Bingham model particularly if the bulk velocity decreases.     

Flow-induced crystallization of high-density polyethylene: the effects of shear and uniaxial extension

The effects of shear, uniaxial extension and temperature on the flow-induced crystallization of two different types of high-density polyethylene (a metallocene and a ZN-HDPE) are examined using rheometry. Shear and uniaxial extension experiments were performed at temperatures below and well above the peak melting point of the polyethylenes in order to characterize their flow-induced crystallization behavior at rates relevant to processing (elongational rates up to 30 s^− 1 and shear rates 1 to 1,000 s^− 1 depending on the application). Generally, strain and strain rate found to enhance crystallization in both shear and elongation. In particular, extensional flow was found to be a much stronger stimulus for polymer crystallization compared to shear. At temperatures well above the melting peak point (up to 25°C), polymer crystallized under elongational flow, while there was no sign of crystallization under simple shear. A modified Kolmogorov crystallization model (Kolmogorov, Bull Akad Sci USSR, Class Sci, Math Nat 1:355–359, 1937) proposed by Tanner and Qi (Chem Eng Sci 64:4576–4579, 2009) was used to describe the crystallization kinetics under both shear and elongational flow at different temperatures.     

Deformation and orientation during shear and elongation of a polycarbonate/carbon nanotubes composite in the melt

In this study, we focused on the elongational rheology and the morphology of an electrically conductive polycarbonate/multiwalled carbon nanotubes (2 wt%) composite in the melt. In shear and melt elongation, the influence of the carbon nanotubes was large when the externally applied stress was small. Consequently, the elastic interactions resulting from the carbon nanotubes dominated in the low frequency range of the shear oscillations. The elongational viscosity of the composite was only moderately influenced by the addition of 2 wt% carbon nanotubes. Transmission electron microscopy investigations of the stretched composite showed that isolated carbon nanotubes were oriented in elongation. In recovery after melt elongation, the recovered stretch of the composite was much smaller than the recovered stretch of pure polycarbonate. This effect is caused by the carbon nanotubes network, which prohibited large extensions of the macromolecules and led to a yield stress of the composite.     

The influence of temperature on rheological properties of blood mixtures with different volume expanders—implications in numerical arterial hemodynamics simulations

During the complicated cardiac surgery on a non-beating heart with cardiopulmonary bypass, protection of the heart is accomplished by injecting cold cardioplegic solutions. In most forms of circulatory shock, it is necessary to immediately restore the circulating volume. Intravenous solutions of volume expanders, such as hydroxyethyl starch and dextrans, are used to increase the volume of fluid in the circulating blood. In this work, blood samples of six donors were obtained and used to prepare mixtures with different volume expanders in concentrations ranging from 10 to 50 vol./vol.%. The flow curves of all mixtures in the temperature range from 4°C to 37°C were constructed and fitted to the Herschel–Bulkley model, in order to extract the shear thinning and yield stress parameters. To assess the influence of the observed changes in the rheological properties of blood on the hemodynamics in arterial vasculature, a realistic three-dimensional rigid-wall computational model was constructed from MRI images of the right carotid bifurcation obtained in vivo from a healthy male volunteer. The time-varying flow field was numerically computed using the Newtonian model as well as the Herschel–Bulkley model with the Papanastasiou regularization. The numerical simulations indicate only moderate changes in the time-averaged flow field that become accentuated when the instantaneous flow field is considered. We also found that although the influence of temperature, hematocrit, and volume expanders on hemodynamics is significant, this can primarily be attributed to the changes in the nominal viscosity of the flow medium.     

Comparison of viscous and elastic properties of polyolefin melts in shear and elongation

Viscous and elastic properties of a linear polypropylene (PP) and a long-chain branched low-density polyethylene (LDPE) have been investigated by creep and creep–recovery experiments in shear and elongation. The data obtained verify the ratios between the linear values of the viscosities and the steady-state elastic compliances in shear and elongation predicted by the theory of linear viscoelasticity. In the nonlinear range, no simple correlation between the viscous behaviour in shear and elongation exists. The elongational viscosity of the PP decreases with increasing stress analogously to the shear thinning observed; the linear range extends to higher stresses in elongation than in shear, however. The LDPE shows thinning in shear and strain hardening in elongational flow. For the LDPE, a linear steady-state elastic tensile compliance corresponding to one third of the linear steady-state elastic compliance in shear was determined. For the PP, this theoretically predicted value is approximately reached. Analogous to the viscous behaviour, the linear range extends to higher stresses in elongation than in shear. For both materials, the steady-state elastic compliances in the nonlinear range decrease with increasing stress in shear as well as in elongation. However, the decrease in elongation is more pronounced.     

The influence of flow confinement on the rheological properties of complex fluids

The rheological properties of plastic and viscoplastic complex fluids have been investigated using constant velocity squeeze flow rheometry at small gaps in order to examine the influence of the resulting flow confinement. The fluids investigated were aqueous carbopol suspensions and three commercial products (Tesco English mustard, Tesco value lemon curd and Heinz tomato ketchup (HTK)). The bulk rheological properties were measured using parallel plate rheometry. Rough plates were used to eliminate wall slip in both rheometric configurations. The commercial products are Herschel–Bulkley fluids and, for gaps less than a critical value of ~200 μm, the yield stresses tended to large values. There was a corresponding trend to small values for the flow consistencies except for HTK, which has a negligibly small bulk value. A possible explanation is that the micro-structural correlation lengths of these fluids are of the order of that of the critical separation. The yield stresses of the carbopol suspensions also increased sharply at the critical gap but then decreased gradually with decreasing gap separation. The transition may correspond to the formation of a jammed state with the subsequent plastic flow at smaller gap separations arising from localised shearing.     

Shear localisation with 2D viscous froth and its relation to the continuum model

Simulations of monodisperse and polydisperse (μ ₂(A) = 0.13±0.002) 2D foam samples undergoing simple shear are performed using the 2D viscous froth (VF) model. These simulations clearly demonstrate shear localisation. The dependence of localisation length on the product λV (shearing velocity V times the wall drag coefficient λ) is examined and is shown to agree qualitatively with published experimental data. A wide range of localisation lengths is found at low λV, an effect which is attributed to the existence of distinct yield and limit stresses. The general continuum model is extended to incorporate such an effect and its parameters are subsequently related to those of the VF model. A Herschel–Bulkley exponent of a = 0.3 is shown to accurately describe the observed behaviour. The localisation length is found to be independent of λV for monodisperse foam samples.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

Linear viscoelastic model for elongational viscosity by control theory

Flows involving different types of chain branches have been modelled as functions of the uniaxial elongation using the recently generated constitutive model and molecular dynamics for linear viscoelasticity of polymers. Previously control theory was applied to model the relationship between the relaxation modulus, dynamic and shear viscosity, transient flow effects, power law and Cox–Merz rule related to the molecular weight distribution (MWD) by melt calibration. Temperature dependences and dimensions of statistical chain tubes were also modelled. The present study investigated the elongational viscosity. We introduced earlier the rheologically effective distribution (RED), which relates very accurately and linearly to the viscoelastic properties. The newly introduced effective strain-hardening distribution (RED_H) is related to long-chain branching. This RED_H is converted to real long-chain branching distribution by melt calibration and a simple relation formula. The presented procedure is very effective at characterizing long-chain branches, and also provides information on their structure and distribution. Accurate simulations of the elongational viscosities of low-density polyethylene, linear low-density polyethylene and polypropylene, and new types of MWDs are presented. Models are presented for strain-hardening that includes the monotonic increase and overshoot effects. Since the correct behaviour at large Hencky strains is still unclear, these theoretical models may aid further research and measurements.     

Numerical simulation of entry flow of the IUPAC-LDPE melt

Numerical simulations have been undertaken for the creeping entry flow of a well-characterized polymer melt (IUPAC-LDPE) in a 4:1 axisymmetric and a 14:1 planar contraction. The fluid has been modeled using an integral constitutive equation of the K-BKZ type with a spectrum of relaxation times (Papanastasiou–Scriven–Macosko or PSM model). Numerical values for the constants appearing in the equation have been obtained from fitting shear viscosity and normal stress data as measured in shear and elongational data from uniaxial elongation experiments. The numerical solutions show that in the axisymmetric contraction the vortex in the reservoir first increases with increasing flow rate (or apparent shear rate), goes through a maximum and then decreases following the behavior of the uniaxial elongational viscosity. For the planar contraction, the vortex diminishes monotonically with increasing flow rate following the planar extensional viscosity. This kinematic behavior is not in agreement with recent experiments. The PSM strain-memory function of the model is then modified to account for strain-hardening in planar extension. Then the vortex pattern shows an increase in both axisymmetric and planar flows. The results for planar flow are compared with recent experiments showing the correct trend.     

Investigation of the rheological properties of short glass fiber-filled polypropylene in extensional flow

The behavior of short glass fiber–polypropylene suspensions in extensional flow was investigated using three different commercial instruments: the SER wind-up drums geometry (Extensional Rheology System) with a strain-controlled rotational rheometer, a Meissner-type rheometer (RME), and the Rheotens. Results from uniaxial tensile testing have been compared with data previously obtained using a planar slit die with a hyperbolic entrance. The effect of three initial fiber orientations was investigated: planar random, fully aligned in the stretching flow direction and perpendicular to it. The elongational viscosity increased with fiber content and was larger for fibers initially oriented in the stretching direction. The behavior at low elongational rates showed differences among the various experimental setups, which are partly explained by preshearing history and nonhomogenous strain rates. However, at moderate and high rates, the results are comparable, and the behavior is strain thinning. Finally, a new constitutive equation for fibers suspended into a fluid obeying the Carreau model is used to predict the elongational viscosity, and the predictions are in good agreement with the experimental data.     

Rheology of pulp suspensions using ultrasonic Doppler velocimetry

Conventional rheometry coupled with local velocity measurements (ultrasonic Doppler velocimetry) are used to study the flow behaviour of various commercial pulp fibre suspensions at fibre mass concentrations ranging from 1 to 5 wt.%. Experimental data obtained using a stress-controlled rheometer by implementing a vane in large cup geometry exhibits apparent yield stress values which are lower than those predicted before mainly due to existence of apparent slip. Pulp suspensions exhibit shear-thinning behaviour up to a high shear rate value after which Newtonian behaviour prevails. Local velocity measurements prove the existence of significant wall slippage at the vane surface. The velocimetry technique is also used to study the influence of pH and lignin content on the flow behaviour of pulp suspensions. The Herschel–Bulkley constitutive equation is used to fit the local steady-state velocity profiles and to predict the steady-state flow curves obtained by conventional rheometry. Consistency between the various sets of data is found for all suspensions studied, including apparent yield stress, apparent wall slip and complete flow curves.     

Constitutive equations for the Doi–Edwards model without independent alignment

We present two representations of the Doi–Edwards model without Independent Alignment explicitly expressed in terms of the Finger strain tensor, its inverse and its invariants. The two representations provide explicit expressions for the stress prior to and after Rouse relaxation of chain stretch, respectively. The maximum deviations from the exact representations in simple shear, biaxial extension and uniaxial extension are of order 2%. Based on these two representations, we propose a framework for Doi–Edwards models including chain stretch in the memory integral form.     

Dynamic simulation of non-spherical particulate suspensions

Particle-level simulation has been employed to investigate rheology and microstructure of non-spherical particulate suspensions in a simple shear flow. Non-spherical particles in Newtonian fluids are modeled as three-dimensional clusters of neutrally buoyant, non-Brownian spheres linked together by Hookean-type constraint force. Rotne–Prager correction to velocity disturbance has been employed to account for far-field hydrodynamic interactions. An isolated rod-like particle in simple shear flow exhibits a periodic orientation distribution, commonly referred to as Jeffery orbit. Lubrication-like repulsive potential between clusters have been included in simulation of rod-like suspensions at various aspect ratios over dilute to semi-dilute volume fractions. Shear viscosity evaluated by orientation distribution qualitatively agrees with one obtained by direct computation of shear stress.     

On the rheology of a dilute suspension of rigid spheres in a weakly viscoelastic matrix fluid

The complete solution for the pressure and the velocity field up to O(De) of a dilute suspension of neutrally buoyant, non-Brownian rigid spheres suspended in an unbounded, weakly viscoelastic matrix fluid, where is the solid volume fraction and De is the Deborah number of the matrix fluid, is presented. The spheres are subjected to an arbitrary linear velocity profile at infinity. The analytical solution is used for the prediction of the bulk stress, and specifically for the calculation of the first and the second normal stress differences in simple shear and uniaxial elongational flows. A comparison of the results with available values reported in the literature is also offered. The final expressions for the bulk normal stress differences in shear and uniaxial elongational flow fully agree with those reported earlier by Greco et al., J. Non-Newton. Fluid Mech., 147 (2007) 1–10.     

Shear rheology of carbon nanotube suspensions

The shear rheology of carbon nanotube suspensions is reviewed from the perspective of colloid and polymer science. In the semi-dilute to concentrated regimes, the nature of the equilibrium or quiescent state is often dominated by nanotube entanglement and strong attractive inter-particle interactions that favor the formation of a disordered network or gel. The strength of this network with respect to the applied stress dictates the development of mesoscale structural anisotropy, first through a global yield stress and then through a critical stress for homogenization. For concentrated suspensions, the nematic liquid-crystalline order anticipated for high-aspect-ratio rigid rods has been observed in a few select scenarios. The opportunity for deeper theoretical insight is emphasized and intuitive physical arguments are offered that might serve as a foundation for future study.     

Elongational flows of some non-colloidal suspensions

The rheology of a system must be explored not only in viscometric flows, but also in other flow classes, and so, we present some results for the axisymmetric elongational flow of non-colloidal suspensions of spheres. We compare our results with data from shear flows using the same matrices and spheres. We have experimented with non-colloidal suspensions of 40-μm diameter polystyrene spheres with volume fractions (?) varying from 0.3 to 0.5. Two matrix fluids were used—one was a near-Newtonian polydimethyl siloxane of 12 Pa-s viscosity and the other was a variant of the M1 Boger fluid sample of Sridhar which we call M1*. We did not find that the Trouton ratio for either of these fluids was 3; generally, the ratio was larger. We investigated the role of sphere roughness using spheres roughened to 5.3 % of the radius in a 50 % suspension in silicone oil and found an increase of elongational viscosity of about 65 % which is comparable with the 60 % increase in shear viscosity with roughness noted previously. For the silicone oil matrix, we found no rate effect, with very little strain-hardening. By contrast, the M1-type matrix suspensions showed strain-hardening and an increase of elongational viscosity with elongation rate.     

Rheological properties and percolation in suspensions of multiwalled carbon nanotubes in polycarbonate

This paper is concerned with several issues related to the rheological behavior of polycarbonate/multiwalled carbon nanotube nanocomposites. The composites were prepared by diluting a masterbatch of 15 wt.% nanotubes using melt-mixing method, and the dispersion was analyzed by SEM, TEM, and AFM techniques. To understand the percolated structure, the nanocomposites were characterized via a set of rheological, electrical, and thermal conductivity measurements. The rheological measurements revealed that the structure and properties were temperature dependent; the percolation threshold was significantly lower at higher temperature suggesting stronger nanotube interactions. The nanotube networks were also sensitive to the steady shear deformation particularly at high temperature. Following preshearing, the elastic modulus decreased markedly suggesting that the nanotubes became more rigid. These results were analyzed using simple models for suspensions of rod-like particles. Finally, the rheological, electrical, and thermal conductivity percolation thresholds were compared. As expected, the rheological threshold was smaller than the thermal and electrical threshold.