Shear viscosity of slightly-elastic concentrated suspensions at low and high shear rates

A quasi-static asymptotic analysis is employed to investigate the elastic effects of fluids on the shear viscosity of highly concentrated suspensions at low and high shear rates. First a brief discussion is presented on the difference between a quasi-static analysis and the periodic-dynamic approach. The critical point is based on the different order-of-contact time between particles. By considering the motions between a particle withN near contact point particles in a two-dimensional “cell” structure and incorporating the concept of shear-dependent maximum packing fraction reveals the structural evolution of the suspension under shear and a newly asymptotic framework is devised. In order to separate the influence of different elastic mechanisms, the second-order Rivlin-Ericksen fluid assumption for describing normal-stress coefficients at low shear rates and Harnoy's constitutive equation for accounting for the stress relaxation mechanism at high shear rates are employed. The derived formulation shows that the relative shear viscosity is characterized by a recoverable shear strain,S _R at low shear rates if the second normal-stress difference can be neglected, and Deborah number,De, at high shear rates. The predicted values of the viscosities increase withS _R , but decrease withDe. The role ofS _R in the matrix is more pronounced than that ofDe. These tendencies are significant when the maximum packing fraction is considered to be shear-dependent. The results are consistent with that of Frankel and Acrivos in the case of a Newtonian suspension, except for when the different divergent threshhold is given as [1 ? (Φ/Φ _m )^1/2] ? 1.     

Influence of temperature on a carbon–fibre epoxy composite subjected to static and fatigue loading under mode-I delamination

In this paper, interlaminar crack initiation and propagation under mode-I with static and fatigue loading of a composite material are experimentally assessed for different test temperatures. The material under study is made of a 3501-6 epoxy matrix reinforced with AS4 unidirectional carbon fibres, with a symmetric laminate configuration [0°]_16/S. In the experimental programme, DCB specimens were tested under static and fatigue loading. Based on the results obtained from static tests, fatigue tests were programmed to analyse the mode-I fatigue behaviour, so the necessary number of cycles was calculated for initiation and propagation of the crack at the different temperatures. G–N curves were determined under fatigue loading, N being the number of cycles at which delamination begins for a given energy release rate. G_ICmax–a, a–N and da/dN–a curves were also determined for different G_cr rates (90%, 85%, 75%, etc.) and different test temperatures: 90 °C, 50 °C, 20 °C, 0 °C, ?30 °C and ?60 °C.     

An experimental study of the first normal stress difference — shear stress relationship in simple shear flow for concentrated shear thickening suspensions

For many polymeric solutions in a simple shear flow a plot of the logarithm of the first normal stress difference N ₁ against the logarithm of the shear stress , for a range of temperatures results in a linear relationship. For such polymer solutions these plots yield a straight line of slope very close to 2 when measured at low shear rates. This relationship is tested against a polymer solution (polyacrylamide in a 50/50 mixture of glycerol/water), a polymer melt (polyethylene), and three concentrated suspensions. These are Superclay (supplied by English China Clays, Cornwall, U.K.) in water, polyvinylchloride in dioctylphthalate and polystyrene latex in water, at volume concentrations of 40% 48%, and 60%, respectively. It was found that the log N ₁ — log relationship is applicable to the polymer solution and melt over a significant range of shear rates and temperatures. In the cases of concentrated suspensions the relationship holds to the point of onset of the shear thickening behaviour. Beyond this point a different relationship exists, however, flow instabilities are apparent. A comment on the contribution of N ₁ and N ₂ to the flow instabilities is made.     

Experimental study on rheological and thermophysical properties of seawater with surfactant additive—part I: rheological properties

The rheological properties of seawater with the addition of surfactant additive (cetyltrimethyl ammonium chloride (CTAC)/sodium salicylate (NaSal)) are measured at different temperatures, including shear viscosity and first normal stress difference (N₁). The effects of the temperature, the salts, and CTAC/NaSal concentration on the rheological properties of test solutions are investigated, and the corresponding influence mechanisms are analyzed. It shows that the addition of salt can decrease the shear viscosities of the solutions, and also decrease N₁ and even eliminate the sharp augment of N₁ above a certain shear rate. The growing elasticity can be characterized by the increase of the initial shear rate for shear-thickening inception. High temperature can also remove the sharp increase of N₁ with salt. Nevertheless, the increase of CTAC/NaSal concentration can withstand the elimination of the sharp augment of N₁.     

Measurement of the first and second normal stress differences: Correlation of four experiments on a polyisobutylene/decalin solution “D1”

Pressure distribution measurements for a polyisobutylene/decalin solution D1 in the Truncated Cone-and-Plate (TCP) apparatus are combined with elastic hole pressures obtained for the same solution on the Lodge Stressmeter^® in order to provide two independent estimates of the second normal-stress difference (N ₂). The values ofN ₂ from the TCP apparatus, obtained by numerical differentiation of a function of the center-hole pressure and the pressure gradient, are in good agreement with measurements made on the same sample by Tanner et al. with a direct method, namely the Tilted Trough Experiment, and by Christiansen et al. with a method that requires an extrapolation to the pressure at the free surface of coneand-plate rheogoniometer data obtained with flush-mounted pressure transducers. The viscosities from the modified Stressmeter for low shear rates extend over five decades of shear rate, including a zero-shear-rate region, and agree with the data of Christiansen on a torque-driven flow. The Higashitani-Pritchard-Baird-Lodge (HPBL) equation relatingN ₁–N ₂ to the hole pressure gives good agreement with the data over a certain range of shear stress. The Newtonian hole pressures for several liquids at 20 and 46 °C compare well with a finite-element calculation for a two-dimensional Poiseuille flow. When the elastic hole pressures from the Stressmeter are combined with the extrapolated rim pressures from the TCP Apparatus in order to extract the value ofN ₂, an agreement betweenN ₂ from the center-hole pressure andN ₂ from the rim pressure can only be obtained up to a shear rate of about 40 s^–1, beyond which the value of –N ₂ from the rim pressure diverges abruptly to negative values. It is possible that this constitutes the first quantitative estimate of an edge effect in cone-and-plate rheometry. Alternatively, the elastic hole pressure in cone-and-plate flow is not equivalent to the elastic hole pressure in Poiseuille flow, at least at high shear rates. The data of Christiansen et al. with flush-mounted pressure transducers appear to confirm this second possibility. Finally, a single set of shift factors obeying the Williams, Landel and Ferry equation superposes the viscosity, the first and the second normal-stress difference within experimental scatter, which can be less than 1% for a certain combination of normal-stress differences. The data were recorded at 3, 20, 30, and 46 °C in the shear rate range 1–260 s^–1.     

The evolution of initially uniform shear flow through a nearly two-dimensional 90° curved duct

An experimental study has been made in a nearly two-dimensional 90° curved duct to investigate the effects of interaction between streamline curvature and mean strain on the evolution of turbulence. The initial uniform shear at the entrance to the curved duct was varied by an upstream shear generator to produce five different shear conditions; a uniform flow (UF), a positive weak shear (PW), a positive strong shear (PS), a negative weak shear (NW) and a negative strong shear (NS). The variations of surface pressure and the mean velocity profiles along the downstream direction under different initial shears are carefully measured. The responses of turbulent Reynolds stresses and triple velocity products to the curvature and the mean strain are also investigated. The evolution of turbulence under the curvature with the different shear conditions is described in terms of the turbulent kinetic energy and the various length scales vs the angular distance θ or a curvature parameters S _c which is defined by S _c = (U/R)/(dU/dy- U/R). The results show that the turbulent kinetic energy and the integral length scale are augmented when S _c < 0.054 whereas they are suppressed when S _c > 0.054. It is also observed that the micro-length scales of Taylor and Kolmogoroff are relatively insensitive to the curvature.     

Rheological behavior in the transient state of PP/EPDM blends with carbon nanofillers

The rheological properties in the transient state of PP/EPDM blends with carbon nanofillers had been studied. The carbon nanofillers were incorporated into molten EPDM in an internal mixer at 150 °C. The rheological variables were determined in rotational rheometry at constant temperature of 200 °C. The results suggest that the magnitude of the difference of the normal stress differences (N₁-N₂) of PP/EPDM blends through the time, with and without nanofillers, and has a transition cycle from positive to negative values and vice versa, at constant and at zero shear rate in previously sheared samples. At constant shear rate, the transition cycle is random; meanwhile, it is constant at zero shear rate. This behavior is attributed to the polymeric chain movement, considering that the sheared samples have two molecular reorder processes: an immediate mechanism and another one slower. The fastest reorder process is attributed to the polymeric chains entanglement forming non-stable and stressed molecular structures. In the other hand, the second process is referred to the molecular mobility that takes place inside the stressed entangled polymer, in such a way that its structure tends to molecular stability as the rest time increases.     

Sur le caractère thermo-extensif de la surface de charge d'un sol non saturéOn the thermo-extensive nature of the yield surface for an unsaturated soil

A study of temperature influence on the yield surface for one unsaturated soil at constant suction is presented. Mechanical consolidation tests are realized at different temperatures on clayey silty sand. A specific triaxial apparatus for unsaturated soils with temperatures included between 30 °C and 60 °C is used. Experimental results show without ambiguity a thermo-extensive nature of the yield surface. The physical interpretation proposed calls for microscopic considerations on the menisci capillary evolution according to temperature and suction. To cite this article: F. Jamin et al., C. R. Mecanique 332 (2004).     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.     

Numerical study of chain conformation on shear banding using diffusive Rolie-Poly model

Shear-banding phenomenon in the entangled polymer systems was investigated in a planar Couette cell with the diffusive Rolie-Poly (ROuse LInear Entangled POLYmers) model, a single-mode constitutive model derived from a tube-based molecular theory. The steady-state shear stress ?? _s was constant in the shear gradient direction while the local shear rate changed abruptly, i.e., split into the bands. We focused on the molecular conformation (also calculated from the Rolie-Poly model) around the band boundary. A band was found also for the conformation, but its boundary was much broader than that for the shear rate. Correspondingly, the first normal stress difference (N ₁) gradually changed in this diffuse boundary of the conformational bands (this change of N ₁ was compensated by a change of the local pressure). For both shear rate and conformation, the boundary widths were quite insensitive to the macroscopic shear rate but changed with various parameters such as the diffusion constant and the relaxation times (the reptation and the Rouse times). The broadness of the conformational banding, associated by the gradual change of N ₁, was attributed to competition between the molecular diffusion (in the shear gradient direction) and the conformational relaxation under a constraint of constant ?? _s.     

A computational study of some rheological influences on the “splashing experiment”

In various attempts to relate the behaviour of highly-elastic liquids in complex flows to their rheometrical behaviour, obvious candidates for study have been the variation of shear viscosity with shear rate, the two normal stress differences N₁ and N₂, especially N₁, the extensional viscosity, and the dynamic moduli G′ and G″. In this paper, we shall confine attention to ‘constant-viscosity’ Boger fluids, and, accordingly, we shall limit attention to N₁, η_E, G′ and G″.We shall concentrate on the “splashing” problem (particularly that which arises when a liquid drop falls onto the free surface of the same liquid). Modern numerical techniques are employed to provide the theoretical predictions. We show that high η_E can certainly reduce the height of the so-called Worthington jet, thus confirming earlier suggestions, but other rheometrical influences (steady and transient) can also have a role to play and the overall picture may not be as clear as it was once envisaged. We argue that this is due in the main to the fact that splashing is a manifestly unsteady flow. To confirm this proposition, we obtain numerical simulations for the linear Jeffreys model.     

Extrudate swell of Boger fluids

Boger fluids are dilute polymer solutions exhibiting high elasticity at low apparent shear rates, which leads to high extrudate swell. Numerical simulations have been undertaken for the flow of three Boger fluids (including benchmark Fluid M1), obeying an integral constitutive equation of the K-BKZ type, capable of describing the behavior of dilute polymer solutions. Their rheology is well captured by the integral model. The flow simulations are performed for planar and axisymmetric geometries without or with gravity. The results provide the extrudate swell and the excess pressure losses (exit correction), as well as the shape and extent of the free surface. All these quantities increase rapidly and monotonically with increasing elasticity level measured by the stress ratio, S_R. It was found that the main reason for the high extrudate swelling is high normal stresses exhibited in shear flow (namely, the first normal-stress difference, N₁). Surprisingly, the elongational parameter of the model or a second normal-stress difference N₂ do not affect the results appreciably. Gravity serves to lower the swelling considerably, and makes the simulations easier and in overall agreement with previous experiments.     

The dependence of viscoelastic flow properties on the structure for styrene-butadiene copolymers

The flow properties for 300 kg/m³ solutions of four-arm, star-branched random and block styrene-butadiene-styrene copolymers in n-butylbenzene are presented. The viscosity, η, first and second normal stress differences, N₁ and N₂, and the steady sher compliance, J, were determined as a function of the shear rate from cone-plate shearing data obtained with a stiffened Model R17 Weissenberg Rheogoniometer. The normal stress differences were determined from total normal force and plate pressure distribution data. Four sensitive, miniature, variable-capacitance pressure transducers mounted in the 7.4-cm plate with their approximately 2.4 mm in diameter pressure-sensing diaphragms flush with the plate surface provided data for the pressure distribution on the plate. In general, the data extend from the zero shear-rate viscosity region somewhat into the shear thinning region. Based on limited data, the zero shear viscosities for the random copolymers increased with (M₃)^3.2, whereas those for the 28% block and 38% block polymers increased with (M_w)^?;6, respectively. The latter high exponents are believed to be a consequence of a network with junctions formed by dispersed phase polystyrene block domains. The sign for N₂ was opposite that of N₁ and the ratio N₂/M₁ for all of the star copolymers averaged –0.214 with a standard deviation of 1.5%. This value is within ±1% of the ratios for tetrachain polybutadienes and polystyrenes and is significantly lower than the –0.29 for linear polybutadiene and polystyrene solutions in normal butylbenzene. The N₂/N₁ ratio did not vary significantly over the shear-rate range investigated.     

A comparison of three different methods for measuring both normal stress differences of viscoelastic liquids in torsional rheometers

A novel pressure sensor plate (normal stress sensor (NSS) from RheoSense, Inc.) was adapted to an Advanced Rheometrics Expansion System rheometer in order to measure the radial pressure profile for a standard viscoelastic fluid, a poly(isobutylene) solution, during cone–plate and parallel-plate shearing flows at room temperature. We observed in our previous experimental work that use of the NSS in cone-and-plate shearing flow is suitable for determining the first and second normal stress differences N ₁ and N ₂ of various complex fluids. This is true, in part, because the uniformity of the shear rate at small cone angles ensures the existence of a simple linear relationship between the pressure [i.e., the vertical diagonal component of the total stress tensor (Π₂₂)] and the logarithm of the radial position r (Christiansen and coworkers, Magda et al.). However, both normal stress differences can also be calculated from the radial pressure distribution measured in parallel-plate torsional flows. This approach has rarely been attempted, perhaps because of the additional complication that the shear rate value increases linearly with radial position. In this work, three different methods are used to investigate N ₁ and N ₂ as a function of shear rate in steady shear flow. These methods are: (1) pressure distribution cone–plate (PDCP) method, (2) pressure distribution parallel-plate (PDPP) method, and (3) total force cone–plate parallel-plate (TFCPPP) method. Good agreement was obtained between N ₁ and N ₂ values obtained from the PDCP and PDPP methods. However, the measured N ₁ values were 10–15% below the certified values for the standard poly(isobutylene) solution at higher shear rates. The TFCPPP method yielded N ₁ values that were in better agreement with the certified values but gave positive N ₂ values at most shear rates, in striking disagreement with published results for the standard poly(isobutylene) solution.

Pressure hole errors in the annular flow measurement of the second normal stress difference

Summary For viscometric, axial, annular flow, the second normal stress differenceN ₂ is related to the difference in normal thrust across the annular space,T _rr . Past attempts using this method have yielded values ofN ₂ for polymer solutions which are different in magnitudeand opposite in sign from those obtained in other experiments. This inconsistency is attributed to errors resulting from the use of pressure holes in the measurement ofT _rr , and is supported by a second-order fluid analysis.The present work focuses on the measurement of the effect of pressure hole errors on the determination ofN ₂ with aqueous polymer solutions. In the measurement ofT _rr , simultaneous use is made of both pressure holes and miniature pressure transducers to measure and account for pressure hole errors. Results indicate that hole errors are sufficiently large to reverse the sign of the computedN ₂. This technique is therefore suggested as a preferred method for determiningN ₂, especially at high shear rates.With 6 figures     

Activation energy of viscoelastic flow emulsion inks as a function of shear rate

The plastic viscosities of five production batches of a nonthixotropic pseudoplasticGestetner 217 type emulsion duplicating ink were determined at 40 RPM (700 sec^–1) using aFerranti-Shirley cone-plate viscometer. The activation energiesE were determined from the slope of anArrhenius plot over the range of 1/2 °C to 40 °C. The value ofE was not affected by batch-to-batch variations.Working with one ink only, the same procedure was then applied to shear rates from 5 to 220 RPM (87.5 to 3850 sec^–1) over the same temperature range.E was found to be constant below about 60 RPM (1000 sec^–1) but dropped rapidly with increasing shear rate above that threshold. A plot ofE vs. In shear rate gave a straight line parallel to the shear rate axis, intersected by a second straight line with a negative slope. The value ofE corresponding to the point of intersection is denotedE _L , the limiting activation energy of viscoelastic flow in a rotational shear field.It is suggested that the rate of shear corresponding toE _L is that necessary to eliminate the last traces of slippage flow in the system and to convert it to 100% laminar flow. The theoretical implications of this suggestion are discussed.Presented at the Joint Meeting of the British Society of Rheology and Research Association of British Paint, Colour and Varnish Manufacturers at Teddington, April 29, 1964.     

Evaluatin of the higashitani and pritchard analysis of the hole pressure using flow birefingence

Higashitani and Pritchard (H-P) carried out an analysis of the hole pressure (P_H) for viscoelastic fluids which leads to expressions relating P_H to the shear stress (σ), the wall shear stress (σ_w), and the normal stress differneces (N₁ and N₂). Although very good agreement their theory and experimental results has been obtained for several polymer solutions and three polymer melts, it is known that at least two of the key assumptions in the theory are violated. In this study flow birefingence has been used to determine the stress field (i.e. σ and the normal stress difference, σ₁₁ — σ₂₂) in the region of a slot placed perpendicular to the flow direction for a polystyrene melt. Values of σ and σ₁₁ — σ₂₂ were then used to evaluate the integrand in the expression relating P_e1 where P_e1 is the hole pressure measured at the base of the rectangular slot, σ and N₁. Values of P_e1 evaluated using flow birefringence data agreed well with those obtained using the same integral expression and cone-and-plate values of N₁ and σ, and with directly measured values of P_e1. This agreement occurred even though the stress field was found to be asymmetric around the centerline of the slot and with secondary flow in the slot. A detailed evaluation of the values of N₁/2σ, which constitute the integral in the H-P theory, along the centerline of the slot revealed most of the contributions to the integral canceled in integrating from the base of the slot where σ is zero at the centerline of the slit-die. The main contributions to the integral occurred from the integration taken from the centerline of the slit-die to the upperdie     

Similarity solutions in free convection boundary-layer flows adjacent to vertical permeable surfaces in porous media: II prescribed surface heat flux

The equation which governs the similarity solution for free convection boundary-layer flow along a vertical permeable surface with prescribed surface heating and mass transfer rate is discussed. The solution is seen to depend on two non-dimensional parameters;m, the power-law exponent, and γ, the mass transfer parameter. It is shown that solutions exist for allm>?1 for γ>0 (fluid injection) whereas for γ<0 (fluid withdrawal), solution exist form>m ₀(γ), wherem ₀ is determined as a function of γ. Solutions for large mass transfer rates are obtained, for both γ>0 and γ<0. For γ>0 the form of the asymptotic solution for γ large is seen to depend on the value ofm. Solutions form large are derived, these are seen to be different depending on whether γ is positive or negative.     

Measuring and assessing first and second normal stress differences of polymeric fluids with a modular cone-partitioned plate geometry

We propose a simple, robust method to measure both the first and second normal stress differences of polymers, hence obtaining the full set of viscometric material functions in nonlinear shear flow. The method is based on the use of a modular cone-partitioned plate (CPP) setup with two different diameters of the inner plate, mounted on a rotational strain-controlled rheometer. The use of CPP allows extending the measured range of shear rates without edge fracture problems. The main advantage of such a protocol is that it overcomes limitations of previous approaches based on CPP (moderate temperatures not exceeding 120 °C, multiple measurements of samples with different volume) and yields data over a wide temperature range by performing a two-step measurement on two different samples with the same volume. The method was tested with two entangled polystyrene solutions at elevated temperatures, and the results were favorably compared with both the limited literature data on the second normal stress difference and the predictions obtained with a recent tube-based model of entangled polymers accounting for shear flow-induced molecular tumbling. Limitations and possible improvements of the proposed simple experimental protocol are also discussed.

Open image in new window

Graphical abstract The effects of edge fracture in start-up shear experiments can be circumvented with the use of a cone-partitioned plate (CPP) geometry. Such a device consists of an inner measuring plate surrounded by an outer nonmeasuring corona. The radius of the sample exceeds that of the measuring plate so that the measured volume is not affected by edge instability. However, the measured first normal stress difference is an apparent one (N_app,1), owing to the contribution of the nonmeasured part of the sample. The figure depicts a schematic design of a modular CPP geometry. Such a fixture is built in a way that the inner tool and the outer partition can be easily replaced, in order to have different measuring diameters (i.e., 6 and 10 mm). From the corresponding signals of the N_app,1, the effective first and second normal stress differences can be calculated.