On some nearly viscometric flows of viscoelastic liquids

Superposition of oscillatory shear imposed from the boundary and through pressure gradient oscillations and simple shear is investigated. The integral fluid with fading memory shows flow enhancement effects due to the nonlinear structure. Closed-form expressions for the change in the mass transport rate are given at the lowest significant order in the perturbation algorithm. The elasticity of the liquid plays as important a role in determining the enhancement as does the shear dependent viscosity. Coupling of shear thinning and elasticity may produce sharp increases in the flow rate. The interaction of oscillatory shear components may generate a steady flow, either longitudinal or orthogonal, resulting in increases in flow rates akin to resonance, and due to frequency cancellation, even in the absence of a mean gradient. An algorithm to determine the constitutive functions of the integral fluid of order three is outlined.Nomenclature A _n Rivlin-Ericksen tensor of order . - A _k Non-oscillatory component of the first order linear viscoelastic oscillatory velocity field induced by the kth wave in the pressure gradient - d Half the gap between the plates - e _x, e _z Unit vectors in the longitudinal and orthogonal directions, respectively - G(s) Relaxation modulus - G History of the deformation - Stress response functional - I() Enhancement defined as the ratio of the frequency dependent part of the discharge to the frequencyindependent part of it at the third order - I ^*() Enhancement defined as the ratio of the increase in discharge due to oscillations to the total discharge without the oscillations - k Power index in the relaxation modulus G(s) - k _i ^–1 Relaxation times in the Maxwell representation of the quadratic shear relaxation modulus (s ₁, s ₂) - m _i ^–1, n _i ^–1 Relaxation times in the Maxwell representations of the constitutive functions ₁(s ₁,s ₂,s ₃) and ₄ (s ₁, s ₂,s ₃), respectively - P Constant longitudinal pressure gradient - p Pressure field - _mx ,⁽³⁾ _nz ,⁽³⁾ Mean volume transport rates at the third order in the longitudinal and orthogonal directions, respectively - ₀,⁽³⁾, ₁,⁽³⁾ Frequency independent and dependent volume transport rates, respectively, at the third order - s = t- Difference between present and past times t and      

Evaluating Herschel-Bulkley fluids with the back extrusion (annular pumping) technique

A mathematical model was developed to describe the behavior of Herschel-Bulkley fluids in a back extrusion (annular pumping) device. A technique was also developed to determine the rheological properties (yield stress, flow behavior index, and consistency coefficient) of these fluids. Mathematical terms were expressed in four dimensionless terms, and graphical aids and tables were prepared to facilitate the handling of the expressions.Nomenclature a radius of the plunger, m - dv/dr shear rate, s^–1 - F force applied to the plunger, N - F _b buoyancy force, N - F _cb force corrected for buoyancy, N - F _T recorded force just before the plunger is stopped, N - F _Te recorded force after the plunger is stopped, N - g acceleration due to gravity, m/s² - H(t) momentary height between plunger and container bottom, m - K a/R, dimensionless - L length of annular region, m - L(t) depth of plunger penetration, m - n flow behavior index, dimensionless - p static pressure, Pa - P _L pressure in excess of hydrostatic pressure at the plunger base, Pa - p ₀ pressure at entrance to annulus, Pa - P pressure drop per unit of length, Pa/m - Q total volumetric flow rate through the annulus, m³/s - r radial coordinate, measured from common axis of cylinder forming annulus, m - R radius of outer cylinder of annulus, m - s reciprocal of n, dimensionless - t time, s - T dimensionless shear stress, defined in Eq. (3) - T ₀ dimensionless yield stress, defined in Eq. (4) - T _w dimensionless shear stress at the plunger wall - _p velocity of plunger, m/s - velocity, m/s - mass density of fluid, kg/m³ - Newtonian viscosity, Pa s - P p ₀ –p _L , Pa - consistency coefficient, Pa sⁿ - value of where shear stress is zero - _–, ₊ limits of the plug flow region (Fig. 1) - r/R - shear stress, Pa - _y yield stress, Pa - _w shear stress at the plunger wall, Pa - dimensionless flow rate defined in Eq. (24) - dimensionless velocity defined by Eq. (5) - _–, ₊ dimensionless velocity outside the plug flow region - _max dimensionless maximum velocity in the plug flow region - _p dimensionless velocity at the plunger wall     

Experimental studies of suspended ferromagnetic fibers in a magnetic field

The analysis of the rotation of a ferromagnetic ellipsoid suspended in a Newtonian fluid and subjected to a uniform magnetic field is extended to include a long, slender cylindrical fiber which is magnetically saturated. Experimental observations of rotating nickel cylinders with aspect ratiosL/D ranging from 5 to 40 agree with the theoretical predictions that: (1) the proper magnetoviscous time constant for the motion is _MV = _s/µ ₀ M _s ² , (2) larger fiber aspect ratios result in considerably longer orientation times; and (3) the strength of the applied external field has only a slight effect on the overall fiber rotation, and has no effect on the maximum angular velocity achieved. Quantitative agreement of theory and experiments is obtained for fibers withL/D 20; for the shorter fibers, the theory tends to overpredict the fiber rotation rate by as much as 30%. D diameter of the cylinder - D ^P (r) position-dependent demagnetization tensor, implicitly defined in eq. (2.5) - D _xx,D _yy,D _zz volume-averaged demagnetizing factors for an ellipsoid equivalent to a uniformly magnetized cylinder, defined in eq. (2.6) - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - k geometric parameter in the hydrodynamic resistance of a body rotating in a Newtonian fluid, eq. (2.2) - L length of the cylinder - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component ofL ^(h) on the cylinder - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.4); thez-component ofL ^(m) on the cylinder - M the magnetization of a magnetic material - M _s the saturation magnitude ofM, approached by all ferromagnetic materials asH _i becomes large - r position vector of a point within a ferromagnetic body - V volume of a magnetic particle - x, y, z rectangular coordinate axes fixed in the cylinder according to figure 1 - angle of inclination of the axis of the cylinder with respect toH ₀ - shear rate - small parameter of slender body theory,=1/ln (2L/D) - _s constant viscosity of the suspending fluid - µ ₀ the magnetic permeability of free space,µ ₀=4 · 10^–7 H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - ₁ the first normal-stress coefficient - ; _z angular velocity of a rotating body; angular velocity of a cylinder about thez-axis, _z =– d/dt     

Measurements of slip at the wall during flow of high-density polyethylene through a rectangular conduit

A hot-film probe has been used to measure slip of high-density polyethylene flowing through a conduit with a rectangular cross section. A transition from no slip to a stick-slip condition has been observed and associated with irregular extrudate shape. Appreciable extrudate roughness was initiated at the same flow rate as that at which the relationship between Nusselt number and Péclet number for heat transfer from the probe departed from the behavior expected for a no-slip condition at the conduit wall. A ₁ constant defined by eq. (A3) - C dimensionless group used in eq. (7) - C _p heat capacity - D constant in eq. (13) - f u _s/u - f _lin defined by eq. (A6) - G storage modulus - G loss modulus - k thermal conductivity - L length of hot film in thex-direction - L _eff effective length of large probe found from eq. (A3) - Nu _L Nusselt number, defined for a lengthL by eq. (2) - (Nu _L)₀ value ofNu _L atPe = 0 (eq. (A 1)) - Pe Péclet number,uL/ - Pe ₀ Péclet number in slip flow, eq. (6) - Pe ₁ Péclet number in shear flow, eq. (4) - q _L average heat flux over hot film of lengthL - R _i resistances defined by figure 8 - R _AB correlation coefficient defined by eq. (14) for signalsA andB - T temperature - T _s temperature of probe surface - T ambient temperature - T T _s –T - u average velocity - u _s slip velocity - V _b voltage indicated in figure 8 - W probe dimension (figure 6) - x distance in flow direction (figure 1) - y distance perpendicular to flow direction (figure 1) - thermal diffusivity,k/C _p - wall shear rate - _5% thickness of lubricating layer during probe calibration for introduction of an error no greater than 5%; see Appendix I - ^* complex viscosity - density - time - _c critical shear stress, eq. (13) - _w wall shear stress - frequency characterizing extrudate distortion (figures 12 and 13), or frequency of oscillation during rheometric characterization (figures 18–20) - * quantity obtained from normalized Nusselt number, eq. (A1), or complex viscosity ^* - A actual (small) probe (see Appendix I) - M model (large) probe (see Appendix I)     

The rotation of a suspended axisymmetric ellipsoid in a magnetic field

Under the influence of a uniform and parallel magnetic field, a ferromagnetic fiber suspended in a Newtonian fluid rotates to align with the field direction. This study examines the field-induced rotation process for an individual non-Brownian axisymmetric ellipsoid suspended in a stagnant Newtonian fluid. Theoretical predictions are derived by a perturbation analysis for the limiting case where the strength of the applied magnetic field far exceeds the saturation magnetization of the ellipsoid. Numerical calculations are performed for the more general problem of an ellipsoid with known isotropic, non-hysteretic magnetic properties, using nickel and a stainless steel as examples. The analysis encompasses materials with field-induced, nonlinear magnetic properties, distinguishing these results from the simpler cases where the particle magnetization is either independent of, or linearly dependent on, the strength of the applied external field. In this study, predictions indicate that when the ellipsoid is magnetically saturated, the particle rotation is governed by the magnetoviscous time constant, _MV = _s/₀ M _s ² . It is found that the rotation rate depends strongly on the aspect ratio,a/b, of the ellipsoid, but only weakly on the dimensionless magnetization,M _s/H ₀. A geometric parameter for an ellipsoid, defined in eq. (2.5) - a, b major, minor semi-axes of an axisymmetric ellipsoid - D demagnetization tensor for an ellipsoid - D ^M magnetometric demagnetization tensor, the volume-average ofD ^P (r) - D ^P (r) position dependent demagnetization tensor, implicitly defined in eq. (2.12) - D _xx,D _yy,D _zz demagnetization factors, the diagonal elements ofD. Values for ellipsoids are defined in eq. (2.15) - F ^(m) magnetic force exerted on a body in a magnetic field - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - h _i ;h _ix,h _iy Cartesian components of dimensionless internal magnetic field,h _i =H _i /H ₀ - I moment of inertia tensor - k geometric parameter for hydrodynamic resistance of a body rotating in a Newtonian fluid given in eq. (2.3) - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component of the hydrodynamic torque - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.10); thez-component of the magnetic torque - M;M the magnetization, or dipole moment density, of a magnetic material; the magnitude ofM - M _s the saturation value ofM, approached by all ferromagnetic materials asH _i becomes large (figure 3) - m _s the dimensionless saturation magnetization,M _s/H ₀ - r position vector of a point within a ferromagnetic body - s dummy integration variable in eq. (2.5) - t time - U magnetoquasistatic potential energy of a magnetic body in a magnetic field, given in eq. (2.8) - u curve-fitting variable in eq. (4.1);u = logH _i - V volume of a magnetic particle; for an axisymmetric ellipsoid,V = (4/3) ab ² - x, y, z rectangular coordinate axes fixed in the ellipsoid (figure 1) - angle of inclination of the major axis of the ellipsoid with respect toH ₀ - _s viscosity of the Newtonian suspending medium - µ ₀ the magnetic permeability of free space,µ ₀ =4 · 10^–7H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - the magnetic susceptibility of a magnetic material, = M/H _i - ; _z angular velocity of a rotating body; angular velocity about thez-axis of an ellipsoid, _z=–d/dt     

Effect of solids concentration in dense suspensions of silica-iron(III) oxide and water

The rheological properties of dense suspensions, of silica, iron (III) oxide and water, were studied over a range of solids concentrations using a viscometer, which was modified so as to prevent settling of the solid components. Over the conditions studied, the material behaved according to power—law flow relationships. As the concentrations of silica and iron(III) oxide were increased, an entropy term in the flow equation was identified which had a silica dependent and an iron (III) oxide dependent component. This was attributed to a tendency to order into some form of structural regularity. A, A, B, C pre-exponential functions (K Pa^–n s^–1) - C _ox volume fraction iron (III) oxide - Q activation energy (kJ mol^–1) - R gas constant (kJ mol^–1 K^–1) - R _v silica/water volume ratio - T temperature (K) - n power-law index - H enthalpy (kJ mol^–1) - S entropy change (kJ mol^–1 K^–1) - shear strain rate (s^–1) - shear stress (Pa)     

Ermittlung des dehnrheologischen Verhaltens thermoplastischer Kautschuke

Zusammenfassung Mit Hilfe eines einfachen Dehnrheometers wurden Dehnfließversuche mit thermoplastischen Kautschukblends (Ibuflex-SEH und Santoprene) durchgeführt. Die auf diese Weise erhaltenen Dehnviskositätsdaten wurden mit den von kautschukmodifizierten Thermoplasten (PVC, PE und PS) verglichen. Nach der Anlaufphase zeigt Ibuflex-SEH ein stationäres Fließverhalten, während die Dehnviskosität von Santoprene als Funktion der Zeit monoton ansteigt, ohne einen stationären Zustand zu erreichen. Die Ursachen des Unterschieds im Fließverhalten der untersuchten Polymermaterialien wurden diskutiert.

---

Extensional flow experiments were carried out on a simple extensional rheometer with the thermoplastic elastomer blends: Ibuflex-SEH and Santoprene. The extensional viscosity data obtained were comparable with those of elastomer modified PVC, PE and PS. After the transient stage Ibuflex-SEH shows a steady-state flow behaviour, whereas the extensional viscosity of Santoprene as a function of time increases monotonically without reaching a stationary state. The causes for these differences were discussed.

Erweiterte Fassung eines Vortrags anläßlich der Jahrestagung der Deutschen Rheologischen Gesellschaft vom 13.–15. Mai 1985 in Berlin     

Der elektroviskose Effekt als Folge elektrostatischer Kraft

Zusammenfassung Die Viskosität vieler Flüssigkeiten ändert sich als Folge elektrischer Felder. Die größten elektroviskosen Effekte findet man unter äußerem Feld bei Suspensionen aus dielektrischen Flüssigkeiten und feinen Festkörperteilchen.Die Flüssigkeiten zeigen ausgeprägtes Bingham-Verhalten. Solange die auf die Wand ausgeübte Scherspannung kleiner als die von der Feldstärke abhängige Haftspannung ist, verhalten sich die Materialien wie Festkörper. Man beobachtet, daß Elektroviskosität proportional zur Feldstärke im Quadrat und umgekehrt proportional zur Schergeschwindigkeit ist, und daß sich die Teilchen der Flüssigkeiten in Richtung der Feldstärke orientieren und faserige Strukturen bilden. In dieser Arbeit erklären wir diese Eigenschaften durch ein einfaches Modell. In diesem Modell resultiert die Scherkraftzunahme aus der elektrostatischen Kraft, die sich als Gradient der im Feld abgespeicherten Energie ergibt.The electroviscosity effect following electrostatic force

---

Apparent viscosity of many fluids changes under an electric field. The largest effects are found with dispersions of fine particles in dielectric liquids subjected to an externally applied field. The large effect is apparently due to induced fibration when the particles align themselves as fibers in field direction. The dispersions exhibit Bingham-behavior, i.e., as long as the shear stress is smaller than the field-dependent yield stress, the material behaves as a solid, a behavior which can be made plausible as the electrodes are mechanically linked together by the fibers (Winslow 1949). In this paper arguments are put forward that the induced fibration is the consequence of electrostatic forces acting on the particles and pulling them into regions of higher field strength. The resulting configuration is an equilibrium configuration and it is postulated that even under shear load, the electric field tends to maintain this configuration such that only the particles close to the electrodes are sheared off, which then form a suspension near the electrode surfaces. A phenomenological theory based on the assumption that the fibers are columns of rectangular particles gives quantitative results for the electroviscosity and the yield stress which agree with experimental observations.

     

The steady shear viscosity of filled polymeric liquids described by a linear superposition of two relaxation mechanisms

Filled polymeric liquids often exhibit apparent yielding and shear thinning in steady shear flow. Yielding results from non-hydrodynamic particle—particle interactions, while shear thinning results from the non-Newtonian behavior of the polymer melt. A simple equation, based on the linear superposition of two relaxation mechanisms, is proposed to describe the viscosity of filled polymer melts over a wide range of shear rates and filler volume fraction.The viscosity is written as the sum of two generalized Newtonian liquid models. The resulting equation can describe a wide range of shear-thinning viscosity curves, and a hierarchy of equations is obtained by simplifying the general case. Some of the parameters in the equation can be related to the properties of the unfilled liquid and the solid volume fraction. One adjustable parameter, a yield stress, is necessary to describe the viscosity at low rates where non-hydrodynamic particle—particle interaction dominate. At high shear rates, where particle—particle interactions are dominated by interparticle hydrodynamics, no adjustable parameters are necessary. A single equation describes both the high and low shear rate regimes. Predictions of the equation closely fit published viscosity data of filled polymer melts. n power-law index - n ₁,n ₂ power-law index of first (second) term - shear rate - steady shear viscosity - ₀ zero-shear rate viscosity - _{0, 1}, _{0, 2} zero-shear rate viscosity of first (second) term - time constant - ₁, ₂ time constant of first (second) term - µ _r relative viscosity of filled Newtonian liquid - ₀ yield stress - ø solid volume fraction - ø _m maximum solid volume fraction     

Viscous heating correction for thermally developing flows in slit die viscometry

In the thermally developing region, d _yy /dx| _y=h varies along the flow direction x, where _yy denotes the component of stress normal to the y-plane; y = ±h at the die walls. A finite element method for two-dimensional Newtonian flow in a parallel slit was used to obtain an equation relating d _yy /dx/ _y=h and the wall shear stress ₀ at the inlet; isothermal slit walls were used for the calculation and the inlet liquid temperature T₀ was assumed to be equal to the wall temperature. For a temperature-viscosity relation /₀ = [1+(T–T₀]^–1, a simple expression [(hd _yy /dx/ _y=h )/ _w0] = 1–[1-F _c(Na)] [M()+P(Pr) ·Q(Gz ^–1)] was found to hold over the practical range of parameters involved, where Na, Gz, and Pr denote the Nahme-Griffith number, Graetz number, and Prandtl number; is a dimensionless variable which depends on Na and Gz. An order-of-magnitude analysis for momentum and energy equations supports the validity of this expression. The function F _c(Na) was obtained from an analytical solution for thermally developed flow; F _c(Na) = 1 for isothermal flow. M(), P(Pr), and Q(Gz) were obtained by fitting numerical results with simple equations. The wall shear rate at the inlet can be calculated from the flow rate Q using the isothermal equation.Notation x,y Cartesian coordinates (Fig. 2) - , dimensionless spatial variables [Eq. (16)] - dimensionless variable, : = Gz(x)^–1 - dimensionless variable [Eq. (28)] - t,t ^* time, dimensionless time [Eq. (16)] - , velocity vector, dimensionless velocity vector - _x , velocity in x-direction, dimensionless velocity - _y , velocity in y-direction, dimensionless velocity - V average velocity in x-direction - _yy , ^* normal stress on y-planes, dimensionless normal stress - shear stress on y-planes acting in x-direction - _w , _w ^* value of shear stress stress at the wall, dimensionless wall shear stress - _w0, _w0 ^* wall shear stress at the inlet, dimensionless variable - , ^* rate-of-strain tensor, dimensionless tensor - wall shear rate, wall shear rate at the inlet - Q flow rate - T, T ₀, temperature, temperature at the wall and at the inlet, dimensionless temperature - h, w half the die height, width of the die - l,L the distance between the inlet and the slot region, total die length - T ₂, T ₃, T ₄ pressure transducers in the High Shear Rate Viscometer (HSRV) (Fig. 1) - P, P2, P3 pressure, liquid pressures applied to T ₂ and T ₃ - , ₀, ^* viscosity, viscosity at T = T ₀, dimensionless viscosity - viscosity-temperature coefficient [Eq. (8)] - k thermal conductivity - C _p specific heat at constant pressure - Re Reynolds number - Na Nahme-Griffith number - Gz Graetz number - Pr Prandtl number     

Rheological characteristics of glass-fibre-filled polypropylene melts

The rheological properties of glass fibre-filled polypropylene melts have been investigated. A high pressure capillary rheometer has been used for the experimental study. The effect of shear rate, temperature, and fibre concentration on the melt viscosity and viscoelastic properties have been studied. An equation has been proposed to correlate the melt viscosity with shear rate, temperature and fibre content. A master curve relation on this basis has been brought out using the shift factora _T . a _T shift factor (=/ _r ) - A _i coefficients of the polynomical of eq. (1) (i = 0, 1, 2, ,n) - B constant in the AFE equation (eq. (2)) (Pa s) - B constant in eq. (3) - D extrudate diameter - d capillary diameter - activation energy at constant shear rate (kcal/mole) - E activation energy at constant shear stress (kcal/mole) - T melt temperature (K) - X fraction glass fibre by weight - shear rate (s^–1) - shear viscosity (Pa s) - normal stress coefficient (Pa s²) - ₁ – ₂ first normal-stress difference (Pa) - shear stress (Pa) - r at reference temperature     

Newtonian stratified flow through an abrupt expansion

The present paper is devoted to the numerical simulations of stratified generalized Newtonian flow. The main feature of the algorithm is to include the coordinates of the interface as primitive unknowns of the flow problem together with the nodal velocity components and pressures calculated on a deformable finite element mesh. Newton-Raphson's iterative method is used for solving the non-linear problem. Special attention is given to the location of the intersection between the interface and a solid boundary, when it is unknown at the outset. The method is applied to the creeping stratified flow of two Newtonian fluids through a one-to-four abrupt expansion; the results remain valid for the contraction problem. It is found that, under appropriate conditions, the small Newtonian corner vortex can increase by a large factor, both in size and intensity; the degrees of freedom of the simulation are the ratios of flow-rate and of viscosity of the fluids.     

Generalized characteristic of polymer blend melt viscosity

A very simple reduction procedure is suggested for the blend viscosities of different polymer pairs. This procedure is based on the comparison of the blend viscosity, normalized either to the matrix or to the disperse phase viscosity, with the viscosities ratio of the initial polymers ( _m / _d ). We have obtained, for 13 different pairs containing 30% of the second component, the universal linear dependencies, mutual analysis of which allows connection of their special points with the stream morphology. The fibrillous morphology takes place in the range of _m / _d = 0, 1–5. Simultaneous, the thin skin consisting of the disperse phase polymers is formed. These results confirm the predominant role of the viscosities ratio in fibrillar composite material formation in comparison with the interphase tension phenomena.     

Experimental investigations of the stability limit of the helical flow of pseudoplastic liquids

The stability of the laminar helical flow of pseudoplastic liquids has been investigated with an indirect method consisting in the measurement of the rate of mass transfer at the surface of the inner rotating cylinder. The experiments have been carried out for different values of the geometric parameter = R ₁/R ₂ (the radius ratio) in the range of small values of the Reynolds number,Re < 200. Water solutions of CMC and MC have been used as pseudoplastic liquids obeying the power law model. The results have been correlated with the Taylor and Reynolds numbers defined with the aid of the mean viscosity value. The stability limit of the Couette flow is described by a functional dependence of the modified critical Taylor number (including geometric factor) on the flow indexn. This dependence, general for pseudoplastic liquids obeying the power law model, is close to the previous theoretical predictions and displays destabilizing influence of pseudoplasticity on the rotational motion. Beyond the initial range of the Reynolds numbers values (Re>20), the stability of the helical flow is not affected considerably by the pseudoplastic properties of liquids. In the range of the monotonic stabilization of the helical flow the stability limit is described by a general dependence of the modified Taylor number on the Reynolds number. The dependence is general for pseudoplastic as well as Newtonian liquids.Nomenclature C _i concentration of reaction ions, kmol/m³ - d = R ₂ –R ₁ gap width, m - F _M () Meksyn's geometric factor (Eq. (1)) - F ₀ Faraday constant, C/kmol - i _l density of limit current, A/m³ - k _c mass transfer coefficient, m/s - n flow index - R ₁,R ₂ inner, outer radius of the gap, m - Re = V _m ·2d·/µ _m Reynolds number - Ta _c = _c ·d^3/2·R ₁ ^1/2 ·/µ _m Taylor number - Z _i number of electrons involved in electrochemical reaction - = R ₁/R ₂ radius ratio - µ apparent viscosity (local), Ns/m² - µ _m mean apparent viscosity value (Eq. (3)), Ns/m² - µ _i apparent viscosity value at a surface of the inner cylinder, Ns/m² - density, kg/m³ - _c angular velocity of the inner cylinder (critical value), 1/s     

The squeezing flow of a model suspension fluid

The paper is concerned with the squeezing flow of a model suspension fluid. The numerical solution obtained by a time-dependent Boundary Element Method is compared to an asymptotic solution at large radius. It is found that the kinematics are Newtonian in character, and the fibres quickly align themselves radially. Consequently, the squeezing force is only weakly dependent on the initial orientations of the fibres and the device can be used for measuring the effective viscosity of the suspension. The effective viscosity found from the squeezing flow agrees surprisingly well with experimental data and numerical data derived from the falling sphere geometry at low volume fractions ( < 0.1).     

Pulsatile blood flow through arterioles

An analytical study was made to examine the effect of vascular deformability on the pulsatile blood flow in arterioles through the use of a suitable mathematical model. The blood in arterioles is assumed to consist of two layers — both Newtonian but with differing coefficients of viscosity. The flow characteristics of blood as well as the resistance to flow have been determined using the numerical computations of the resulting expressions. The applicability of the model is illustrated using numerical results based on the existing experimental data. r, z coordinate system - u, axial/longitudinal velocity component of blood - p pressure exerted by blood - _b density of blood - µ viscosity of blood - t time - , displacement components of the vessel wall - T _t0,T ₀ known initial stresses - density of the wall material - h thickness of the vessel wall - T _t,T stress components of the vessel - K _l,K _r components of the spring coefficient - C _l,C _r components of the friction coefficient - M _a additional mass of the mechanical model - r ₁ outer radius of the vessel - thickness of the plasma layer - r ₁ inner radius of the vessel - circular frequency of the forced oscillation - k wave number - E ₀,E _t, , _t material parameters for the arterial segment - µ _p viscosity of the plasma layer - Q total flux - Q _p flux across the plasma zone - Q _h flux across the core region - Q mean flow rate - resistance to flow - P pressure difference - l length of the segment of the vessel     

The evolution of linear viscoelasticity during the vulcanization of polyethylene

The evolution of linear viscoelasticity during the vulcanization of polyethylene is studied through the gel point. The material in the vicinity of the gel point is described by two scaling laws: one characterizes the viscoelasticity at the critical point and a second characterizes the evolution of viscoelasticity near the gel point. Time Resolved Mechanical Spectroscopy is used to observe both scaling phenomena. The material at the gel point displays power law relaxation over five decades of time with a power-law relaxation exponent equal to 0.32. This study conforms with previous findings that this exponent is composition-dependent. The longest relaxation time diverges in the vicinity of the gel point as _max |p _c - p| ^–1/, and we find = 0.2. This result conforms with previous reports that this exponent may be independent of composition. The Arrhenius flow activation energy for this material undergoes three-fold changes during crosslinking up to the gel point. A single-adjustable-parameter stretched-exponential-power law relaxation function is an inadequate representation of crosslinked materials over any significant range of extent of the reaction up to the gel point.     

Shear creep and recovery of a technical polystyrene

The torsional creep and recoverable bahaviour of a technical polystyrene is reported over seven orders of magnitude of the value of the compliance from 10^–8 to 10^–1 Pa^–1 and over more than seven decades in time. The results for the recoverable compliance J_R (t) reveal a dispersion region seen between the glass transition and the steady-state recoverable compliance J_e. The limiting value of the final dispersion J_e = 4.7 · 10^–4 Pa^–1 indicates a broad molecular-weight distribution. The steady-state recoverable compliance J_e is independent of the temperature. The temperature dependence of the final dispersion was found to be indistinguishable from that of viscous flow. However, this temperature dependence differs significantly from that of the glass-rubber transition. A proposal has been made for the construction of creep compliance and recoverable compliance over an extended time scale.     

The hole pressure problem: On the Higashitani-Pritchard theory for transverse and axial slots

Analysis of the Higashitani-Pritchard (HP) theory for a transverse slot provides insight into the nature of the errors involved. The effects of geometry, elasticity and asymmetry on the errors in the HP theory are also discussed. Inertial effects and the effect of asymmetry are included in the HP theory for a transverse slot, through modification of the pressure term in the momentum equation; the extended HP theory proposed in this work provides evidence as to the reasons for the failure of the Newtonian hole pressure as a correction term at higher Reynolds numbers. Finally, the HP theory is used to analyse a control volume of fluid in the axial slot geometry and deduce a rigorous sufficient condition for secondary flows to exist. This analysis illustrates the nature of the role played by the second normal-stress difference (N ₂) in causing secondary flows, and suggests the possibility of the direction of the secondary flows being related to the sign ofN ₂. A computational study of all three aspects of the HP theory investigated in this work seems worth-while. x, y, z Cartesian coordinates - u, Cartesian velocity components inx andy directions - velocity vector - q ₁,q ₂,q ₃ orthogonal curvilinear coordinates - h ₁,h ₂,h ₃ scale factors of curvilinear coordinate system - P isotropic pressure - N ₁ first normal-stress difference - N ₂ second normal-stress difference - R _e Reynolds number - R _L hole-based Reynolds number (eqs. (5.3) and (5.4)) - W _e Weissenberg number - H channel height - W hole width - L hole depth - S path of zero slope of streamlines (pathS) - P _H hole pressure (eq. (2.1)) - P _I,P _E inertial and elastic pressures (eqs. (5.7–9)) - P _R pressure term used in correction for asymmetry (eqs. (5.19–21)) - P _IH,P _EH inertial and elastic hole pressures (eqs. (5.17) and (5.12)) - P _RH,P _SH,P _TH correction terms used in extended HP theory (eqs. (5.23–25)) - total stress tensor - extra-stress tensor - del operator - unit tensor - _ij, _ij components of total and extra stress tensors in curvilinear coordinate system - _w wall shear stress at pointB of figure 1 - unperturbed wall shear rate - density of fluid - viscosity of fluid - ₂ second normal-stress coefficient - _ij radius of curvature of surface of constantq _i in theq _j direction     

Viscoelastic swirling flow with free surface in cylindrical chambers

The flow of a viscoelastic liquid driven by the steadily rotating bottom cover of a cylindrical cup is investigated. The flow field and the shape of the free surface are determined at the lowest significant orders of the regular domain perturbation in terms of the angular velocity of the bottom cap. The meridional field superposed on a primary azimuthal field shows a structure of multiple cells. The velocity field and the shape of the free surface are strongly effected by the cylinder aspect ratio and the elasticity of the liquid. The use of this flow configuration as a free surface rheometer to determine the first two Rivlin-Ericksen constants is shown to be promising.Nomenclature R, ,Z Coordinates in the physical domain D - , , Coordinates in the rest stateD ₀ - r, ,z Dimensionless coordinates in the rest stateD ₀ - Angular velocity - Zero shear viscosity - Surface tension coefficient - Density - Dimensionless surface tension parameter - ₁, ₂ The first two Rivlin-Ericksen constants - Stream function - Dimensionless second order meridional stream function - ^* Dimensionless second normal stress function - 2 Dimensionless sum of the first and second normal stress functions - N ₁,N ₂ The first and second normal stress functions - n Unit normal vector - D Stretching tensor - A _n nth order Rivlin-Ericksen tensor - S Extra-stress - u Velocity field - U Dimensionless second order meridional velocity field - V Dimensionless first order azimuthal velocity field - p Pressure - Modified pressure field - P Dimensionless second order pressure field - J Mean curvature - a Cylinder radius - d Liquid depth at rest - D Dimensionless liquid depth at rest - h Free surface height - H Dimensionless free surface height at the second order