Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

Rheology of molten polymers

Summary TheCross equation describes the flow of pseudoplastic liquids in terms of an upper and a lower Newtonian viscosity corresponding to infinite and zero shear, and ₀, and of a third material constant related to the mechanism of rupture of linkages between particles in the intermediate, non-Newtonian flow regime, Calculation of of bulk polymers is important, since it cannot be determined experimentally. The equation was applied to the melt flow data of two low density polyethylenes at three temperatures.Using data in the non-Newtonian region covering 3 decades of shear rate to extrapolate to the zero-shear viscosity resulted in errors amounting to about onethird of the measured ₀ values. The extrapolated upper Newtonian viscosity was found to be independent of temperature within the precision of the data, indicating that it has a small activation energy.The ₀ values were from 100 to 1,400 times larger than the values at the corresponding temperatures.The values of were large compared to the values found for colloidal dispersions and polymer solutions, but decreased with increasing temperature. This shows that shear is the main factor in reducing chain entanglements, but that the contribution of Brownian motion becomes greater at higher temperatures.

---

Zusammenfassung Die Gleichung vonCross beschreibt das Fließverhalten von pseudoplastischen Flüssigkeiten durch drei Konstante: Die obereNewtonsche Viskosität (bei sehr hohen Schergeschwindigkeiten), die untereNewtonsche Viskosität ₀ (bei Scherspannung Null), und eine Materialkonstante, die vom Brechen der Bindungen zwischen Partikeln im nicht-Newtonschen Fließbereich abhängt. Die Berechnung von ist wichtig für unverdünnte Polymere, wo man sie nicht messen kann.Die Gleichung wurde auf das Fließverhalten der Schmelzen von zwei handelsüblichen Hochdruckpolyäthylenen bei drei Temperaturen angewandt. Die Werte von ₀, durch Extrapolation von gemessenen scheinbaren Viskositäten im Schergeschwindigkeitsbereich von 10 bis 4000 sec^–1 errechnet, wichen bis 30% von den gemessenen ₀-Werten ab. Die Aktivierungsenergie der war so klein, daß die-Werte bei den drei Temperaturen innerhalb der Genauigkeit der Extrapolation anscheinend gleich waren. Die ₀-Werte waren 100 bis 1400 mal größer als die-Werte.Im Verhältnis zu kolloidalen Dispersionen und verdünnten Polymerlösungen war das der Schmelzen groß, nahm aber mit steigender Temperatur ab. Deshalb wird die Verhakung der Molekülketten hauptsächlich durch Scherbeanspruchung vermindert, aber der Beitrag derBrownschen Bewegung nimmt mit steigender Temperatur zu.

     

Extensional flow of HDPE/LDPE blends

Elongational viscosity data, obtained through the converging flow analysis by Cogswell, are presented for two types of HDPE/LDPE blends at various compositions and different temperatures. The results relative to the homopolymer parents compare favourably with literature data obtained also with different and more sophisticated techniques. Those relative to blends show peculiar features for the two cases: when the newtonian viscosity of the LDPE is higher all the blends show a behaviour typical of the LDPE with a maximum in _el / ₀ enhanced at small percentage of HDPE; when the newtonian viscosity of the LDPE is similar to that of the HDPE there is only a gradual change in the properties.     

Melt viscosity of linear and branched poly(butyleneisophthalate)

Linear and branched poly(butyleneisophthalate) samples were synthesized and characterized in terms of the intrinsic viscosity, the molecular weight and the melt viscosity over a wide range of shear rates at 200 °C. An exponent of about 4.6 in the equation relating ₀ to was found for linear samples; this high value is probably due to the high content of cyclic oligomers in low molecular weight samples. Both linear and branched samples exhibited Newtonian behaviour over a rather wide range of shear rates, but for any given melt-viscosity, the branched samples became shear thinning at lower shear rates than the linear ones. Our experimental data were found to fit a previously proposed correlation between the melt viscosity ratio ( _{0, b} / _{0, 1} ) and a branching index quite well.     

Die temperatur- und druckinvariante Darstellung der Scherviskosität von geschmolzenem Polymethylmethacrylat

Zusammenfassung An Polymethylmethacrylaten mit zahlenmittleren Molekulargewichten zwischen 8000 und 145000 g/Mol wurde die Scherviskosität bei Temperaturen im Bereich von 130–190 °C und bei Drücken bis zu 1000 kp/cm² in Abhängigkeit vom Schergefälle gemessen. Trägt man die an ein und derselben Probe gemessenen relativen Viskositäten/ ₀ (₀ untere Newtonsche Viskosität) über ₀· doppeltlogarithmisch auf, so erhält man eine einzige temperatur- und druckinvariante Kurve (master curve). Mit sinkendem Molekulargewicht machen sich jedoch bei großen Abszissenwerten zunehmend Abweichungen von der Temperatur- und Druckinvarianz bemerkbar. Die an Proben mit verschiedenem Molekulargewicht ermittelten invarianten Kurven weisen Wendepunkte auf, die sich mit wachsendem Molekulargewicht immer mehr zu niedrigeren Werten von/ ₀ verlagern. Durch Verschieben in Abszissenrichtung kann man die Kurven bis zu ihren Wendepunkten miteinander zur Deckung bringen und erhält auf diese Weise einen begrenzten Bereich der Molekulargewichtsinvarianz.Die Temperatur- und Druckinvarianz läßt sich sowohl mit den vonPrandtl undEyring als auch mit den vonF. Bueche entwickelten Modellen für das Fließen von Hochpolymeren begründen. Aus der Temperatur- und Druckinvarianz folgt, daß der Temperatur- und der Druckkoeffizient der Viskosität mit zunehmendem Schergefälle in gleicher Weise abnehmen und ein Minimum durchlaufen, wie auch experimentell bestätigt wird.

---

Summary The shear viscosity of polymethylmethacrylates (number average molecular weights from 8000 to 145000 g/mole) was measured in dependence of shear rate in the temperature range from 130–190 °C and with pressures up to 1000 kp/cm². If for one specimen ₀ ( ₀ lower Newtonian viscosity) is plotted over ₀· , both in logarithmic scales, one gets a curve, which is independent of temperature and pressure (master curve). With decreasing molecular weights increasing deviations from the master curve are found. The master curves found with specimens of different molecular weights have inflection points, which shift to lower values of ₀ with increasing molecular weight. By shifting of the whole curves in direction of the abscissa the parts until to the inflection point might be reduced to a single plot which is independent of molecular weight.The independence of temperature and pressure might be derived from the models developed by eitherPrandtl andEyring orF. Bueche for the flow in high polymers. From this follows that the temperature- and the pressure-coefficient of viscosity will similarly decrease with increasing shear rate until to a minimum, which was experimentally verified.

---

---

Erweiterte Fassung eines auf der Jahrestagung der Deutschen Rheologen am 21. Mai 1968 in Berlin gehaltenen Vortrags.

---

Herrn Dr.H. Kausch und Herrn Dipl.-Ing.H. Schönewald sind wir für wertvolle Diskussionen zu Dank verpflichtet. Der Arbeitsgemeinschaft industrieller Forschungsvereinigungen e.V. (AIP) danken wir für die finanzielle Unterstützung dieser Arbeit.     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration     

The secondary flow induced around a sphere in an oscillating stream of elastico-viscous liquid

The present paper is devoted to the theoretical study of the secondary flow induced around a sphere in an oscillating stream of an elastico-viscous liquid. The boundary layer equations are derived following Wang's method and solved by the method of successive approximations. The effect of elasticity of the liquid is to produce a reverse flow in the region close to the surface of the sphere and to shift the entire flow pattern towards the main flow. The resistance on the surface of the sphere and the steady secondary inflow increase with the elasticity of the liquid.Nomenclature a radius of the sphere - b ^ik contravariant components of a tensor - e contravariant components of the rate of strain tensor - F() see (47) - G total nondimensional resistance on the surface of the sphere - g _ik covariant components of the metric tensor - f, g, h secondary flow components introduced in (34) - k ₀ measure of relaxation time minus retardation time (elastico-viscous parameter) - K =k ₀ ²/V ₀ ² , nondimensional parameter characterizing the elasticity of the liquid - n measure of the ratio of the boundary layer thickness and the oscillation amplitude - N, T defined in (44) - p arbitrary isotropic pressure - p _ik covariant components of the stress tensor - p ^ik contravariant components of the stress tensor associated with the change of shape of the material - R =V ₀ a/v, the Reynolds number - S =a/V ₀, the Strouhall number - r, , spherical polar coordinates - u, v, w r, , component of velocity - t time - V(, t) potential velocity distribution around the sphere - V ₀ characteristic velocity - u, v, t, y, P nondimensional quantities defined in (15) - reciprocal of s - density - defined in (32) - defined in (42) - ₀ limiting viscosity for very small changes in deformation velocity - complex conjugate of - oscillation frequency - = ₀/, the kinematic coefficient of viscosity - , defined in (52) - (, y) stream function defined in (45) - =(NT/2n)^1/2 y - /t convective time derivative (1) _ik      

A three-parameter model describing the experimental relation between shear stress and shear rate for laminar flow of aqueous polymer solutions

Summary A three-parameter model is introduced to describe the shear rate — shear stress relation for dilute aqueous solutions of polyacrylamide (Separan AP-30) or polyethylenoxide (Polyox WSR-301) in the concentration range 50 wppm – 10,000 wppm. Solutions of both polymers show for a similar rheological behaviour. This behaviour can be described by an equation having three parameters i.e. zero-shear viscosity ₀, infinite-shear viscosity , and yield stress ₀, each depending on the polymer concentration. A good agreement is found between the values calculated with this three-parameter model and the experimental results obtained with a cone-and-plate rheogoniometer and those determined with a capillary-tube rheometer.

---

Zusammenfassung Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit von strukturviskosen Flüssigkeiten wird durch ein Modell mit drei Parametern beschrieben. Mit verdünnten wäßrigen Polyacrylamid-(Separan AP-30) sowie Polyäthylenoxidlösungen (Polyox WSR-301) wird das Modell experimentell geprüft. Beide Polymerlösungen zeigen im untersuchten Schergeschwindigkeitsbereich von ein ähnliches rheologisches Verhalten. Dieses Verhalten kann mit drei konzentrationsabhängigen Größen, nämlich einer Null-Viskosität ₀, einer Grenz-Viskosität und einer Fließgrenze ₀ beschrieben werden. Die Ergebnisse von Experimenten mit einem Kegel-Platte-Rheogoniometer sowie einem Kapillarviskosimeter sind in guter Übereinstimmung mit den Werten, die mit dem Drei-Parameter-Modell berechnet worden sind.

---

a Pa^–1 physical quantity defined by:a = {1 – ( / ₀)}/ ₀ - c l concentration (wppm) - D m capillary diameter - L m length of capillary tube - P Pa pressure drop - R m radius of capillary tube - u m s^–1 average velocity - v _r m s^–1 local axial velocity at a distancer from the axis of the tube - shear rate (–dv _r /dr) - local shear rate in capillary flow - s^–1 wall shear rate in capillary flow - Pa s dynamic viscosity - _a Pa s apparent viscosity defined by eq. [2] - ( _a ) Pa s apparent viscosity in capillary tube at a distanceR from the axis - ₀ Pa s zero-shear viscosity defined by eq. [4] - Pa s infinite-shear viscosity defined by eq. [5] - l ratior/R - kg m^– density - Pa shear stress - ₀ Pa yield stress - _r Pa local shear stress in capillary flow - _R Pa wall shear stress in capillary flow _R = (PR/2L) - _v m³ s^–1 volume rate of flow With 8 figures and 1 table     

A study on the melt viscosity of linear and branched poly(butyleneterephthalate)

Linear and branched PBTP samples were synthesized and characterized in terms of the intrinsic viscosity, the melt-flow-index and, for some, the melt viscosity over a range of shear rates at 250 °C.An exponent of 3.2 in the equation relating to was found for linear samples. Both linear and branched samples exhibited Newtonian behaviour over a wide range of shear rates, but for any given melt-viscosity the branched samples became shear thinning at lower shear rates than the linear ones. Correlation between a branching index,, and melt-visocity ratio (0,b/0,l) was in agreement with a previous theoretical study.     

Application of melt flow index for rheological characterization of branched and linear polymers

A procedure for evaluating rheological characteristics, such as the master curves log/ ₀ vs. log % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xd9Gqpe0x% c9q8qqaqFn0dXdir-xcvk9pIe9q8qqaq-xir-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaaieGaceWFZo% Gbaiaaaaa!3B59!\[\dot \gamma \] ₀ and flow curves, using the melt flow index is described for branched and linear polymers. Experimental data on the melt flow index and branching degree are needed for this purpose, as well as some polymer constants, i.e. coefficients of the ₀ vs. MFI relation and coefficients of fluidity dependence on molecular characteristics. An example is given for bisphenol A polycarbonate.     

The very slow flow of a Powell-Eyring fluid around a sphere

Summary The very slow flow of a Powell-Eyring type non-Newtonian fluid around a sphere is investigated by a variational technique. The result, a correction factor that is applied to the Stokes' equation, is given as a plot and as an equation which is empirically fit to the plot. Also, a comparison of the very slow flows of a simplified viscoelastic Oldroyd fluid and the Powell-Eyring fluid is made which indicates that in a certain restricted region of the very slow flows, both models give essentially the same results. The Oldroyd and Powell-Eyring model parameters are interrelated by forcing both models to fit the same tube flow viscosity data.Nomenclature B dimensionless quantity, v /R - C dimensionless second invariant - c ₁ constant determined by variational method - D dimensionless variational integral - D _2j , D _j+k position-independent variables used in specification of trial functions - E _2j , E _j+k position-independent variables used in specification of trial functions - f friction factor - f _corr friction factor correction - F _drag drag force on sphere - g, g ₀, g ₁ general trial function; first and second terms in the general trial function - G, H terms in the expression for C - j index - J variational integral - k index - K term in the expression for C - p, q integers - r integer, radial coordinate - R radius of sphere - Re Reynolds number - Re ₀ Reynolds number at point of zero shear rate - Re Reynolds number at infinite distance from sphere - Re _NN Reynolds number based on variable part of viscosity - u, v dimensionless position coordinates - V volume considered - v _i ith velocity component - v _r , v , v _z velocity components in the r, , and z-directions - v approach velocity of the fluid - x/ parameter in Powell-Eyring model - x _i i-position coordinate - parameter in Powell-Eyring model - rate of deformation - , _c , _N , ₀ coefficient of viscosity; cross viscosity; parameter in Powell-Eyring model; viscosity in limit of zero shear rate - spherical coordinate - , _ij rate of deformation tensor; ij-component of rate of deformation tensor - ₁, ₂ parameters in Oldroyd model - Newtonian viscosity - ₁, ₂ parameters in Oldroyd model - dimensionless radial coordinate, r/R - second invariant - fluid density - spherical coordinate - stream function     

Measurement of the relaxation time in Darcy flow

The inertia of a liquid flowing through a porous medium is normally ignored, but if the acceleration is great, it may be important. The relaxation time, defined so that it alone accounts for the inertia, has been determined experimentally with a simple oscillator. A U-Tube is provided with a porous plug and filled with a liquid. During pendulation of the liquid, the frequency and the damping define the relaxation time. The measured value of the relaxation time is about 10 times the theoretical estimate derived from Navier-Stokes equation.Symbols E modulus of elasticity - E _D dissipated energy - E _k kinetic energy - g acceleration of gravity - G pressure gradient - h height - K ₀ permeability - L length of porous plug - n porosity - P dissipated power - pressure - R half the tube length - R _c radius of the tube bend - r radial coordinate - r _o radius of the tube - s coordinate along a streamline in the tube - t time - v flux per unit area - it relaxation time - , auxiliary variables - , v dynamic and kinematic viscosity - , velocity potential for inviscid flow and gravity potential - dissipation function - displacement of the liquid - , _o frequency of damped and undamped oscillations     

General viscosity-composition relationship for dispersions,solutions and binary liquid systems

An equation for the viscosity of a mixture of two imaginary Newtonian liquids is derived. In the derivation the mathematical assumption is used that the effective activation energy for viscous flow of a binary liquid mixture is a linear combination of the reciprocals of the activation energy of the components. It contains two dependent fitting constants and has the same structure as the Mooney equation for dispersions of spherical solid particles, the Huggins equation for polymer solutions and is identical to an equation by Hoffmann and Rother, when written in the variables that the last authors used.As a consequence it can be shown that the viscosity of binary liquid mixtures, liquid resion solutions, dispersions of solid spherical particles and polymer solutions can be described very well by one and the same equation, up to the highest concentrations.It has further been found that the viscosity of dispersions of non-spherical particles, solutions of solids in organic solvents and solutions of electrolytes and non-electrolytes in water can also be described by this formula. The equation permits the construction of a straight line on which all liquids can be plotted.An algebraic analysis of the equation shows that each series of viscosity composition data can be placed in one of three rheological groups independent of the type of fraction that is used to characterize the composition.Seventy-four binary systems, covering a wide range of liquids have been used to show the applicability of the developed equation.It has been found that in most cases the data are best described by splitting them into two regions, each with its own set of dependent constants. General symbol for the fraction or concentration of the component with the higher viscosity determining the composition of a binary mixture [—] - _v Volume fraction of the component with the higher viscosity [—] - _w Weight fraction of the component with the higher viscosity [—] - _mw Molecular weight fraction of the component with the higher viscosity [—] - c Concentration of the component with the higher viscosity [g/cm³] - E ₂,E ₁,E Activation energy for viscous flow referring to the component with the higher viscosity, the lower viscosity and the viscosity of the binary mixtures, respectively [J] - ₂, ₁, Experimental parameter (with the dimension of energy) referring to the component with the higher viscosity, the component with the lower viscosity and to the binary mixtures, respectively [J] - ₁, ₂ Viscosity of the component with the lower and the higher viscosity, respectively [Pa · s] - Viscosity of a binary mixture [Pa · s] - [] The usual intrinsic viscosity of the component with the highest viscosity [cm³/g] - _r / ₁ [—] - _{sp r} – 1 [—] - [] -intrinsic viscosity [—] - [] _v Volume intrinsic viscosity [—] - [] _w Weight intrinsic viscosity [—] - [] _c Concentration intrinsic viscosity, identical to [] [cm³/g] - T _e Temperature at which the two liquids have the same viscosity [K] - _e Viscosity at temperatureT _e [Pa · s] - P ₁,P ₂ Density of the component with the lower and the higher viscosity, respectively - R Gas constant [J · Mol^–1 · K^–1]     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

In-line measurement of rheological properties of polymer melts

Shear viscosity (), first normal stress difference (N ₁), and extensional viscosity ( _E ) of polymer melts measured under processing conditions are important in process modeling, quality control, and process control. A slit rheometer that could simultaneously measure , N ₁, and the planar extensional viscosity ( _p ) was designed and tested by attaching it in-line to a laboratory model single-screw extruder. A tube (circular cross-section) rheometer to measure and the uniaxial extensional viscosity ( _u ) simultaneously was also designed and tested. Two commercial grades of LDPE (low density polyethylene) with melt index values of 6 and 12 were used as test materials for the study. Exit and hole pressure methods were used to estimate N ₁, and the entrance pressure drop method using the analyses of Cogswell, Binding, and Gibson (the last analysis used with the axisymmetric case only) was used to estimate _E .The hole pressure method was considered better than the exit pressure method to estimate N ₁ (due to the greater susceptibility of the latter to experimental errors). From the hole pressure method N ₁ was obtained from 100 kPa to 500 kPa over a range of shear rates from 40 s^–1 to 700 s^–1. Among the analyses used to estimate the extensional viscosity, Cogswell's is recommended due to its simpler equations without loss of much information compared to the other analyses. The range of extension rates achieved was 1 to 30 s^–1. The combination of the hole pressure and entrance pressure drop methods in a slit rheometer is a feasible design for a process rheometer, allowing the simultaneous measurement of the shear viscosity, first normal stress difference and planar extensional viscosity under processing conditions. Similarly, combining the entrance pressure drop measurements with a tube rheometer is also feasible and convenient.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Determination of the steady-state shear viscosity from measurements of the apparent viscosity for some common types of viscometers

Considering a number of model fluids, the relation between the (measurable) apparent viscosity _a and the (true) shear viscosity is studied for some commonly used viscometers, like capillary, slit, plate-plate and concentric cylinders (including the influence of the bottom of the cylinder), as well as for one laboratory type of viscometer. As long as is a purely monotonic function, a shift factor < 1 allows one to deduce from _a . Though in general variable, it frequently suffices for practical purposes to use a constant shift factor (the constant being characteristic of the type of viscometer used). This does not apply to dilute solutions or any fluids with two plateau values for . For plastic fluids, it is shown that Casson or Bingham behavior can — if valid at all — only describe the high shear stress limit of _a .     

Pulsatile flow in a tube with wall injection and suction

This paper studies similarity solutions for pulsatile flow in a tube with wall injection and suction. The Navier-Stokes equations are reduced to a system of three ordinary differential equations. Two of the equations represent the effects of suction and injection on the steady flow while the third represents the effects of suction and injection on pulsatile flow. Since the equations for steady flow have been studied previously, the analysis centers on the third equation. This equation is solved numerically and by the method of matched asymptotic expansions. The exact numerical solutions compare well with the asymptotic solutions.The effects of suction and injection on pulsatile flow are the following: a) Small values of suction can cause a resonance-like effect for low frequency pulsatile flow. b) The annular effect still occurs but for large injection or suction the frequency at which this effect becomes dominant depends on the cross-flow Reynolds number. c) The maximum shear stress at the wall is decreased by injection, but may be increased or decreased by suction.Nomenclature a radius of the tube - a ₀ ² i ² - A₀, B₀, C₀, D₀, E₀ constant coefficients appearing in the expression for pressure - b a non-dimensionalized length - b ₀ ² i ^{2 2} - b _k complex coefficients of a power series - B - C ₁, C ₂, D complex constants - d - D _1,2 - f() F(a ^1/2)/aV - f ₀,f ₁,... functions of order one used in asymptotic expansions of f() - F(r) rv _r - g() - G(r) a steady component of velocity in axial direction - h() 4/C₀ a ² H(a ^1/2) - h ₀,h ₁,h ₂,...;l ₀,l ₁,l ₂,... functions of order one used in asymptotic expansions for h() in outer regions - H(r) complex valued function giving unsteady component of velocity - H ₀, H ₁, H ₂, ... K ₀, K ₁, K ₂, ...; L ₀, L ₁, L ₂, ... functions of order one used in asymptotic expansions for h() in inner regions - i - J ₀, J ₁, Y ₀, Y ₁ Bessel functions of first and second kind - k - K Rk/2b ² - O order symbol - p pressure - p ₁(z, t) arbitrary function related to pressure - r radial coordinate - r ₀ (1+16 ^{4 4})^1/4 - R Va/, the crossflow Reynolds number - t time - u() G(r)/V - v _r radial velocity - v _z axial velocity - V constant velocity at which fluid is injected or extracted - z axial coordinate - ² a ²/4 - 4.196 - small parameter; =–2/R (Sect. 4); =–R/2 (Sect. 5); =2/R(Sect. 6) - r ²/a ² - * 0.262 - Arctan (4 ^{2 2}) - , inner variables - kinematic viscosity - b - * zero of g() - density - (r, t) arbitrary function related to axial velocity - frequency     

Concentration-dependent behaviour of the shear viscosity of coal-fuel oil suspensions

A function correlating the relative viscosity of a suspension of solid particles in liquids to their concentration is derived here theoretically using only general thermodynamic ideas, with out any consideration of microscopic hydrodynamic models. This function ( _r = exp (1/2B ^* C ²)) has a great advantage over the many different functions proposed in literature, for it depends on a single parameter,B ^*, and is therefore concise. To test the validity of this function, a least-squares regression analysis was undertaken of available data on the viscosity and concentration of suspensions of coal particles in fuel oil, which promise to be a useful alternative to fuel oil in the near future. The proposed function was found to accurately describe the concentration-dependent behaviour of the relative viscosity of these suspensions. Furthermore, an attempt was made to obtain information about the factors affecting the value ofB ^*, however the results were only qualitative because of, among other things, the inaccuracy of the viscosity measurements in such highly viscous fluids. shear viscosity of the suspension - ₀ shear viscosity of the Newtonian suspending medium - _r = /₀ relative viscosity - solid volume concentration - c solid weight concentration - _m maximum attainable volume concentration of solids - solid volume concentration at which the relative viscosity of the suspension becomes infinite - c _m maximum attainable solid weight concentration - _s density of the solid phase - _l density of the liquid phase - _m density of the suspension - k _n coefficients of theø-power series expansion of _r - { _j } sets of parameters specifying the thermodynamic state of the solid phase of a suspension - T absolute temperature (K) - f (c, T, _j) formal expression for the relative variation of the viscosity with concentration = [1 / (/c)] _T,j - d median size of the granulometric distribution - _B plastic or Bingham viscosity - K consistency factor - n flow index - g ([c _m –c],T, _j ) function including an asymptotic divergence asc tends toc _m , formally describing the concentration dependent behaviour of the shear viscosity of a suspension - A (T, _j) regression analysis parameters - B (T, _j) regression analysis parameters - B ^* (T, _j ) regression analysis parameters     

The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)