Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

Comparison between uniaxial and biaxial elongational flow behavior of viscoelastic fluids as predicted by differential constitutive equations

Behavior of polymer melts in biaxial as well as uniaxial elongational flow is studied based on the predictions of three constitutive models (Leonov, Giesekus, and Larson) with single relaxation mode. Transient elongational viscosities in both flows are calculated for three constitutive models, and steady-state elongational viscosities are obtained as functions of strain rates for the Giesekus and the Larson models.Change of elongational flow behavior with adjustable parameter is investigated in each model. Steady-state viscosities _E and _B are obtained for the Leonov model only when the strain-hardening parameter is smaller than the critical value _cr determined in each flow. In this model, uniaxial elongational viscosity _E increases with increasing strain rate , while biaxial elongational viscosity _B decreases with increasing biaxial strain rate _B . The Giesekus model predictions depend on the anisotropy parameter . _E and _B increase with strain rates for small _B while they decrease for large . When is 0.5, _E in increasing, but _B is decreasing. The Larson model predicts strain-softening behavior for both flows when the chain-contraction parameter > 0.5. On the other hand, when is small, the steady-state viscosities of this model show distinct maximum around = _B = 1.0 with relaxation time . The maximum is more prominent in _E than in _B .     

Fractal-time stochastic processes and dynamic functions of polymeric systems

Dynamic material functions of polymeric systems are calculated via a defect-diffusion model. The random motion of defects is modelled by a fractaltime stochastic process. It is shown that the dynamic functions of polymeric solutions can be approximated by the defect-diffusion process of the mixed type. The relaxation modulus of Kohlrausch type is obtained for a fractal-time defect-diffusion process, and it is shown that this modulus is capable of portraying the dynamic behavior of typical viscoelastic solutions.The Fourier transforms of the Kohlrausch function are calculated to obtain and. A three-parameter model for and is compared with the previous calculations. Experimental measurements for five polymer solutions are compared with model predictions. D rate of deformation tensor - G(t) mechanical relaxation modulus - H relaxation spectrum - I(t) flux of defects - P _n (s) probability of finding a walker ats aftern-steps - P generating function ofP _n (s) - s(t) fraction of surviving defects - , () gamma function (incomplete) - ₀ zero shear viscosity - ^* () complex viscosity - frequency - t _n n-th moment - F[] Fourier transform - f ^* (u) Laplace transform off(t) - , components of ^* - G _f, _f ^* fractional model - G ₃, ₃ ^* three parameter model - complex conjugate ofz - material time derivative ofD     

Fragile networks and rheology of concentrated suspensions

Suspensions consisting of particles of colloidal dimensions have been reported to form connected structures. When attractive forces act between particles in suspension they may flocculate and, depending on particle concentration, shear history and other parameters, flocs may build-up in a three-dimensional network which spans the suspension sample. In this paper a floc network model is introduced to interpret the elastic behavior of flocculated suspensions at small deformations. Elastic percolation concepts are used to explain the variation of the elastic modulus with concentration. Data taken from the suspension rheology literature, and new results with suspensions of magnetic -Fe₂O₃ and non-magnetic -Fe₂O₃ particles in mineral oil are interpreted with the model proposed.Non-zero elastic modulus appeared at threshold particle concentrations of about 0.7 vol.% and 0.4 vol.% of the magnetic and non-magnetic suspensions, respectively. The difference is attributed to the denser flocs formed by magnetic suspensions. The volume fraction of particles in the flocs was estimated from the threshold particle concentration by transforming this concentration into a critical volume concentration of flocs, and identifying this critical concentration with the theoretical percolation threshold of three-dimensional networks of different coordination numbers. The results obtained indicate that the flocs are low-density structures, in agreement with cryo-scanning electron micrographs. Above the critical concentration the dynamic elastic modulus G was found to follow a scaling law of the type G ( _f - _f ^c ) ^f , where _f is the volume fraction of flocs in suspension, and _f ^c is its threshold value. For magnetic suspensions the exponent f was found to rise from a low value of about 1.0 to a value of 2.26 as particle concentration was increased. For the non-magnetic a similar change in f was observed; f changed from 0.95 to 3.6. Two other flocculated suspension systems taken from the literature showed a similar change in exponent. This suggests the possibility of a change in the mechanism of stress transport in the suspension as concentration increases, i.e., from a floc-floc bond-bending force mechanism to a rigidity percolation mechanism.     

Neck propagation as a shock wave

Neck propagation in the stretching of elastic solid filaments having a yield point was analyzed using the space one-dimensional thin filament governing equations developed previously by the authors and other researchers. Constitutive model for the filament was assumed to be expressible as engineering tensile stress(X) (tensile force) given as a function of elongational strain with the(X) curve having a yield point maxima followed by a minima and a breaking point greater than the yield point maxima. Also incorporated into the model is the hysteresis of irreversible plastic deformation. When inertia is taken into consideration, the thin filament equations were found to reduce to the nonlinear wave equation ² (X)/ ² =C ₁ ² X/ ² where is Lagrangean space coordinate, is time, andC ₁ is inertia coefficient. The above nonlinear wave equation yields a solutionX(, ) having a stepwise discontinuity inX which propagates along the axis. The zero speed limit of the step wave solution was found to describe the above neck propagation occurring in solid filaments. Furthermore, it was recognized that the nonlinear wave equation was known for many years to also govern the plastic shock wave which propagates axially within a metal rod subjected to a very strong impact on its end. The one-dimensional atmospheric shock wave also was known to be governed by the nonlinear wave equation upon making certain simplifying assumptions. The above and other evidences lead to the conclusion that neck propagation occurring in the extension of solid filament obeying the above(X) function can be formally described as a shock wave.     

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     

Analyse und Berechnung der Strömungsvorgänge im zylindrischen Scherspalt von Plastifizieraggregaten

Zusammenfassung Dieser Aufsatz zeigt eine Möglichkeit auf, zylindrische Scherteile einer Plastifiziereinheit, auf der strukturviskose Materialien verarbeitet werden, approximativ zu berechnen. Es ist möglich, den Volumenstrom und Druckabfall, die mittlere Schergeschwindigkeit, Scherdeformation und Schubspannung im Scherspalt zu approximieren. Durch diese Gleichungen wird eine Abschätzung der Verteil- und Zerteilvorgänge im Scherelement möglich.

---

A method is described for approximating the flow in cylindrical shearing gaps of plasticating extruder, which is applicable to shear thinning materials. It is possible to calculate the through-put and pressure drop as well as the shear rate, strain and shear stress in the gap. With these equations the distribution and separation process in shearing gaps can be evaluated.

---

D Zylinderdurchmesser - d ₁ Schnecken-Kerndurchmesser der Meteringzone - d _s Durchmesser des zylindrischen Scherteils - K Konstante im Potenzfließgesetz - K _0T Koeffizient des Potenzfließgesetzes - L ₁ Länge der Anlaufschräge - L _s Länge des zylindrischen Scherteils - n Fließindex - n ₀ Drehzahl - p Druckabfall über der Scherteillänge - s Scherspalthöhe - T _M Massetemperatur - ₀ Umfangsgeschwindigkeit - _0x Geschwindigkeitskomponente inx-Richtung - _x, _z Geschwindigkeit inx- bzw.z-Richtung als Funktion der Koordinatey - Volumenstrom - x, z Ortskoordinaten - Exponent des Potenzfließgesetzes - Schergeschwindigkeit - mittlere Schergeschwindigkeit - Viskosität - dimensionslose Höhe - Dichte der Schmelze - Schubspannung - _yx, _yz Schubspannungskomponenten - _xx, _zz Normalspannungskomponenten - _ps dimensionsloser Druckgradient - dimensionsloser Volumenstrom - _x, _z dimensionslose Geschwindigkeit inx- bzw.z-Richtung     

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

Some aspects of dilute suspension theory

A theory proposed by the author as representative of the flow of a general suspension contains three interaction forces, f, S and N. For a quasi-concentrated suspension and for a dilute suspension, N and S, N are omitted, respectively. For the latter special case, we treat diffusion of a fluid through an elastic solid. For a quasi-concentrated suspension, we show that F and S depend on the gradient of the motion gradient. We demonstrate the existence of interesting phenomena: non-simple behavior, dissipative effects, generalized lift and drag forces.Presented at the second conference Recent Developments in Structured Continua, May 23 – 25, 1990, in Sherbrooke, Québec, Canada.     

An experimental study on effects of solid concentration and ζ-potential on pseudoplastic flow of suspensions

The pseudoplastic flow of suspensions, alumina or styrene-acrylamide copolymer particles in water or an aqueous solution of glycerin has been studied by the step-shear-rate method. The relation between the shear rate,D, and the shear stress,, in the step-shear-rate measurements, where the state of dispersion was considered to be constant, was expressed as = AD ^1/2 +CD. The effective solid volume fraction,ø _F, andA were dependent on the shear rate and expressed byø _F =aD ^b andA = D . Combining the above relations, the steady flow curve was expressed by = D ^{1/2 +} + ₀ (1 – a D ^b/0.74)^–1.85 D, where ₀ is the viscosity of the medium.With an increase in solid volume fraction and a decreases in the absolute value of the-potential, the flow behavior of the suspensions changed from Newtonian ( = = b = 0), slightly pseudoplastic ( = b = 0), pseudoplastic ( = 0) to a Bingham-like behavior.The change in viscosity of the medium had an effect on the change in the effective volume fraction.     

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     

Transient stress responses predicted by the internal viscosity model in elongational flow

Predictions are made for the elongational-flow transient rheological properties of the dilute-solution internal viscosity (IV) model developed earlier by Bazua and Williams. Specifically, the elongational viscosity growth function _e ⁺ (t) is presented for abrupt changes in the elongational strain rate . For calculating _e ⁺, a novel treatment of the initial rotation of chain submolecules is required; such rotation occurs in response to the macroscopic step change of at t = 0. Representative are results presented for N = 100 (N = number of submolecules) and = 1000 f and 10000 f (where is the IV coefficient and f is the bead friction coefficient), using h ^* = 0 (as in the original Rouse model) for the hydrodynamic interaction. The major role of IV is to cause the following changes relative to the Rouse model: 1) abrupt stress jump at t = 0 for _e ⁺; 2) general time-retardance of response. There is no qualitative change from the Rouse-model prediction of unbounded il growth when exceeds a critical value ( ), and calculations of submolecule strains at various show that the unbounded- _e ^– behavior arises from unlimited submolecule strains when . However, the time-retardance could delay such growth beyond the timescale of most experiments and spinning processes, so that the instability might not be detected. Finally, _e ⁺ (t) and _e ( ) in the limit are presented for N = 1 and compared with exact predictions for the analogous rigid-rod molecule; close agreement lends support for the new physical approximation introduced for solving the transient dynamics for any N.     

The stability of the helical flow of pseudoplastic liquids in a narrow annular gap with a rotating inner cylinder

The stability of a laminar helical flow of pseudoplastic liquids in an annular gap with a rotating inner cylinder is investigated theoretically. The analysis is carried out under the assumption of a torroidal form of the secondary flow (torroidal Taylor vortices) for the narrow gap geometry. The power law model has been applied to describe the pseudoplasticity of liquids. The problem of the stability has been formulated with the aid of the method of small disturbances, and solved using the Galerkin method. In order to describe the stability limit the Reynolds and Taylor numbers defined with the aid of the mean viscosity value have been introduced. It has been found that pseudoplasticity has a considerably destabilizing influence on the Couette motion as well as on the helical flow in the initial range of the Reynolds number values (Re<30). A decrease of the flow index value,n, is accompanied by a decrease of the critical value of the Taylor number. This destabilizing effect of pseudoplasticity vanishes in the range of the larger values of the Reynolds number. In the rangeRe>30, the stability limit of the flow of pseudoplastic liquids can be described by a general dependence of the critical valueTa _c onRe, which is consistent with results obtained for the case of Newtonian fluids. a frequency number (Eq. (27)), 1/s - b wave number (Eq. (27)), 1/m - B = M/N parameter - d = R ₂ –R ₁ gap width, m - f(y, B, k) function of viscosity distribution (Eq. (7)) - f ₀ (x) function of viscosity distribution (narrow gap Eq. (35)) - F(x) = V(x)/V _m dimensionless distribution of axial flow velocity - G(x) = U(x) _i dimensionless distribution of angular flow velocity - K consistency coefficient, N sⁿ/m² - M = (P/L)R ₂ parameter of the stress field (Eq. (1)), N/m² - M ₀ torque per unit length, N - n flow index - N = M ₀/(2R ₂ ² ) parameter of the stress field (Eq. (1)), N/m² - p = 1/2n–1/2 parameter - pressure disturbance amplitude, N/m² - p pressure disturbance, N/m² - (P/L) pressure drop per unit length of the gap, N/m² - r radial coordinate, m - r _m location of the maximum value of the axial velocity, m - R ₁,R ₂ inner, outer radius of the annulus, m - Re = V _m 2d/ _m Reynolds number - S = (P/L · d/N) parameteer of the stress field (narrow gap) - t time, s - Ta = _i d ^3/2 R ₁ ^1/2 / _m Taylor number - U tangential velocity, m/s - U _i tangential velocity at the surface of the inner cylinder, m/s - V axial velocity, m/s - V _m mean axial velocity, m/s - V disturbance vector of velocity field, m/s - amplitude of theV _k -disturbance, m/s - X, Y, Z functions in Eqs. (36–38) - y = r/R ₂ dimensionless radial coordinate - x = (r—(R ₁+R ₂)/2)d radial coordinate (narrow gap) - L ₁ L ₄ linear operators in Eqs. (36–38) - = ad/V _m dimensionless frequency number - = b·d dimensionless wave number - component of the rate of strain tensor, 1/s - component of the rate of strain tensor corresponding to the disturbance, 1/s - = R ₁/R ₂ radius ratio - apparent viscosity, Ns/m² - ₀ apparent viscosity in the main flow, Ns/m² - µ disturbance of the apparent viscosity, Ns/m² - µ _m mean apparent viscosity, Ns/m² - density, kg/m³ - _ij component of the stress tensor, N/m² - angular velocity, rad/s - _i angular velocity of the inner cylinder, rad/s     

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     

On the relationship between swelling of myofibrils and their suspension rheology

The swelling of myofibrils extracted from white bovine muscle was followed by measuring their suspension rheology. Swelling of the myofibril with increasing pH and ionic strength was accompanied by an increase in both the steady shear viscosity of the suspension and the dynamic viscoelastic properties. Swelling was continuously monitored by measuringG while the ionic strength of the suspension was being changed by dialysis. The relationship between the degree of swelling and the rheological parameters is complicated because myofibrils are rodshaped and swell radially and therefore swelling results in a change in shape. To allow for this an attempt was made to generalize the data by plotting viscosity andG againstcS/ø _m , wherec is the protein concentration in the suspension,S is the swollen volume of the myofibrils per weight of protein, and ø _m is the maximum packing fraction.The best fit to the data was represented by the equations _sp = 1.05 (cS/ _m – 0.84)^1.23 Pa · sG = 8.78 (cS/ _m – 0.67)^2.22 N m^–2. The scatter was greatest forG, possibly because at low degrees of shear the myofibrils were associated and this was confirmed by optical microscopy. Pronounced non-Newtonian behavior was observed and it was suggested that this was due to the disruption of aggregate structures, although at low concentration, orientation of the rods in the shear field may also be important.     

Temperature and concentration dependent viscosity and gelation temperature of ABA triblock copolymer solutions

In solutions of ABA-triblock copolymers in a poor solvent for A thermoreversible gelation can occur. A three-dimensional dynamic network may form and, given the polymer and the solvent, its structure will depend on temperature and polymer mass fraction. The zero-shear rate viscosity of solutions of the triblock-copolymer polystyrene-polyisoprene-polystyrene in n-tetradecane was measured as a function of temperature and polymer mass fraction, and analyzed; the polystyrene blocks contained about 100 monomers, the polyisoprene blocks about 2000 monomers. Empirically, in the viscosity at constant mass fraction plotted versus inverse temperature, two contributions could be discerned; one contribution dominating at high and the other one dominating at low temperatures. In a comparison with theory, the contribution dominating at low temperatures was identified with the Lodge transient network viscosity; some questions remain to be answered, however. An earlier proposal for defining the gelation temperature T _gel is specified for the systems considered, and leads to a gelation curve; T _gel as a function of polymer mass fraction.Mathematical symbols {} functional dependence; e.g., f{x} means f is a function of x - p _log logarithm to the base number p; e.g., ¹⁰log is the common logarithm - exp exponential function with base number e - sin trigonometric sine function - lim limit operation - – in integral sign: Cauchy Principal Value of integral, e.g., - derivative to x - partial derivative to x Latin symbols dimensionless constant - b constant with dimension of absolute temperature - constant with dimension of absolute temperature - B dimensionless constant - c mass fraction - dimensionless constant - constant with dimension of absolute temperature - d ^* dimensionless constant - D{0} constant with dimension of absolute temperature - e base number of natural (or Naperian) logarithm - g distribution function of inverse relaxation times - G relaxation strength relaxation function - h distribution function of relaxation times reaction constant enthalpy of a molecule - H Heaviside unit step function - i complex number defined by i ² = –1 - j{0} constant with dimension of viscosity - j index number - k Boltzmann's constant - k _H Huggins' coefficient - m mass of a molecule - n number - N number - p index number - s entropy of a molecule - t time - T absolute temperature Greek symbols as index: type of polymer molecule - as index: type of polymer molecule - shear as index: type of polymer molecule - shear rate - small variation; e.g. T is a small variation in T relative deviation Dirac delta distribution as index: type of polymer molecule - difference; e.g. is a difference in chemical potential - constant with dimension of absolute temperature - (complex) viscosity - constant with dimension of viscosity - [] intrinsic viscosity number - inverse of relaxation time - chemical potential - number pi; circle circumference divided by its diameter - mass per unit volume - relaxation time shear stress - angular frequency     

Influence of the formation of associations of macromolecules on the rheo-optical properties of their dilute solutions. I. Theory

Assuming the formation of doublets in the flow according to a mass action law, the shear rate and the concentration dependence of the extinction angle, of the birefringence, and of the average coil expansion are calculated for dilute solutions of flexible macromolecules. It is shown that this reversible association process has a strong influence on the measurable parameters in a flow birefringence experiment. c concentration (g/cm³) - h ² mean square end-to-end distance at shear rate - h ₀ ² mean-square end-to-end distance at zero-shear rate - n refractive index of the solution (not very different from the solvent for a very dilute solution) - E mean coil expansion - K ₀,K constant of the mass action law - M molecular weight - R _G gas constant - T absolute temperature - ₁ – ₂ optical anisotropy of the segment - ₀ Deborah number: - Deborah number: - shear rate - ₀, reduced concentration - _s viscosity of the solvent - [] ₀ intrinsic viscosity at zero-shear rate - [] intrinsic viscosity at shear rate - extinction angle - N _a Avodagro's number - _n magnitude of the birefringence     

On the relaxation of shear and normal stresses of viscoelastic fluids following constant shear rate experiments

By means of a cone and plate rheometer the relaxation of the shear stress and the first normal stress difference in polymer liquids upon cessation of a constant shear rate were examined. The experiments were conducted mostly in a high shear rate region of relevance for the processing of these materials. The relaxation behavior at these shear rates can only be measured accurately under extremely precise specifications of the rheometer. To determine under which conditions the integral normal thrust is a convenient measure for the relaxing local first normal stress difference the radial distribution of the pressure in the shear gap was measured. The shape of relaxation of both the shear stress and the first normal stress difference could be closely approximated for the entire measured shear rate and time range by a two parameter statistical function. In the range of measured shear rates, one of the parameters, the standard deviationS, is equal for the shear and the normal stress, and is independent of the shear rate within the limit of experimental error. The second parameter, the mean relaxation timet _50, of the shear stress andt _50, of the first normal stress difference, can be calculated approximately from the viscosity function and only a single relaxation experiment.     

Determination of discrete relaxation and retardation time spectra from dynamic mechanical data

A powerful but still easy to use technique is proposed for the processing and analysis of dynamic mechanical data. The experimentally determined dynamic moduli,G() andG(), are converted into a discrete relaxation modulusG(t) and a discrete creep complianceJ(t). The discrete spectra are valid in a time window which corresponds to the frequency window of the input data. A nonlinear regression simultaneously adjust the parametersg _i , _i ,i = 1,2, N, of the discrete spectrum to obtain a best fit ofG, G, and it was found to be essential that bothg _i and _i are freely adjustable. The number of relaxation times,N, adjusts during the iterative calculations depending on the needs for avoiding ill-posedness and for improved fit. The solution is insensitive to the choice of initial valuesg _i,0, _i,0,N ₀. The numerical program was calibrated with the gel equation which gives analytical expressions both in the time and the frequency domain. The sensitivity of the solution was tested with model data which, by definition, are free of experimental error. From the relaxation time spectrum, a corresponding discrete set of parametersJ ₀,, J _d,i and _i of the creep complianceJ(t) can then readily be calculated using the Laplace transform.This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.