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.     

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     

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     

Rheological models for gypsum plaster pastes

Shear stress and shear rate data obtained for gypsum plaster pastes were correlated by means of different rheological models. The pastes were prepared from a commercial calcium sulfate hemihydrate at various water/plaster ratios ranging from 100/150 to 100/190. The tests were performed at 25°C using a rotating coaxial cylinder viscosimeter. The measurements were accomplished by applying a step-wise decreasing shear rate sequence. Discrimination among the models was made: (1) on the basis of the fitting goodness; (2) by checking the physical meaning of the calculated parameters; (3) on the basis of the stability of the parameters and of their prediction capacity beyond the limits of the experimental data. In the light of above, the Casson model seemed to be most effective for application to gypsum plaster pastes. K Consistency - n Power-law index - N Number of experimental data - P Number of parameters - Shear rate (s^–1) - ₀ Viscosity (Pa · s) - _d Dispersing medium viscosity (Pa · s) - _p Plastic viscosity (Pa · s) - Viscosity at zero shear rate (Pa · s) - Viscosity at infinite shear rate (Pa · s) - [] Intrinsic viscosity - ² Variance - Shear stress (Pa) - ₀ Yield stress (Pa) - Solid volume fraction - _m Maximum solid volume fraction     

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     

Driven torsion pendulum for measuring the complex shear modulus in a steady shear flow

To investigate the viscoelastic behavior of fluid dispersions under steady shear flow conditions, an apparatus for parallel superimposed oscillations has been constructed which consists of a rotating cup containing the liquid under investigation in which a torsional pendulum is immersed. By measuring the resonance frequency and bandwidth of the resonator in both liquid and in air, the frequency and steady-shear-rate-dependent complex shear modulus can be obtained. By exchange of the resonator lumps it is possible to use the instrument at four different frequencies: 85, 284, 740, and 2440 Hz while the steady shear rate can be varied from 1 to 55 s^–1. After treatment of the theoretical background, design, and measuring procedure, the calibration with a number of Newtonian liquids is described and the accuracy of the instrument is discussed.Notation a radius of the lump - A geometrical constant - b inner radius of the sample holder - c constant - C ₁, C ₂ apparatus constants - D damping of the pendulum - e _x , e _y , e _z Cartesian basis - e _r , e , e _z orthonormal cylindrical basis - E geometrical constant - E _t , ⁰ E _t , _t relative strain tensor - f function of shear rate - F _t relative deformation tensor - G (t) memory function - G ^* complex shear modulus - G Re(G ^* ) - G Im(G ^* ) - h distance between plates - H ^* transfer function - , functional - i imaginary unit: i ²= – 1 - I moment of inertia - J _exc excitation current - J ₀ amplitude of J _exc - k ^* = k – ik complex wave number - K torsional constant - K fourth order tensor - l length of the lump - L mutual inductance - M _dr driving torque - M _liq torque exerted by the liquid - ⁰ M _liq, _liq steady state and dynamic part of Mliq - n power of the shear rate - p isotropic pressure - Q quality factor - r radial position - R,R ₀, R _c Re(Z ^*, Z ₀ ^* , Z _c ^* ) - s time - t, t time - T temperature - T, ⁰ T, stress tensor - u velocity - U lock-in output - ₀ velocity - V _det detector output voltage - V _sig, V _cr signal and cross-talk part of V _det - x Cartesian coordinate - X , X ₀, X _c Im(Z ^*, Z ₀ ^* , Z _c ^* ) - y Cartesian coordinate - z Cartesian coordinate, axial position     

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     

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     

Zum Mischen rheologisch inhomogener Stoffsysteme auf Einschneckenmaschinen

Zusammenfassung Das Mischen von Stoffen mit unterschiedlichen rheologischen Eigenschaften in Schneckenmaschinen ist in der Kunststoffauf- und -verarbeitung eine Standardaufgabe. Trotzdem gibt es hierfür kein zufriedenstellendes mathematisch-physikalisches Modell. Daher werden zunächst einfache Mischmodelle diskutiert. Auf der Basis dieser Modelle wird dann unter Berücksichtigung der Besonderheiten des Plastifizierextruderprozesses eine Mischgütebeziehung mathematisch formuliert. Die experimentelle Überprüfung erfolgt mit Hilfe der Grauwertanalyse extrudierter Zweistoffsysteme, bei denen ein Stoff mit Ruß eingefärbt war. Da der Mischprozeß hochgradig stochastisch ist, streuen die Meßergebnisse. Unter Berücksichtigung dieses Tatbestandes ist der theoretische Ansatz zufriedenstellend.

---

Mixing of polymer resins with different rheological properties is a usual demand in plastics processing using screw extruders. A mathematical model describing this processing problem sufficiently is not known, however. Therefore, simple mixing models will be discussed. Based on these, a concept for the calculation of mixing homogeneity will be presented, including the particular requirement of the plasticating screw process. An experimental investigation utilizes the grey-value analysis of extruded two-component materials, which in one phase is carbon-black filled. Considering the fact that the mixing process is highly random, the theoretical model leads to a good level of aggreement with the scattering measurement data.

---

b Schneckenkanalbreite - B Bandbreite der Grauwerte - c Konstante - mittlere Konzentration, bezogen auf die Grauwertbandbreite - h Höhe, Gangtiefe, Schneckenkanalhöhe - h ₀ Gangtiefe der Einzugszone - h ₁ Gangtiefe der Ausstoßzone - L Länge - gemittelte Schmelzebettlänge - n Exponent des Potenzfließgesetzes - s Standardabweichung der Grauwerte bezogen auf die Grauwertbandbreite - S Standardabweichung der Grauwerte - t Verweilzeit - t ₁ kürzeste Verweilzeit - mittlere Verweilzeit - ₀ Umfangsgeschwindigkeit - mittlere Geschwindigkeit - V Volumenstrom - w Dicke eines Kontrollelements - w Ausstreichdicke eines Kontrollelements - x Koordinate - Mittelwert der Grauwerte - y Koordinate - Scherdeformationswinkel - Scherdeformation - mittlere Scherdeformation - Schergeschwindigkeit - Viskosität - ₁ dimensionslose kürzeste Verweilzeit - dimensionsloser Volumenstrom - _LSM laminarer Schermischgrad - _{LSM, the} theoretischer laminarer Schermischgrad - _{LSM, exp} experimenteller laminarer Schermischgrad - ² Varianz der Verweilzeit im Schmelzebett - Schubspannung - Gangsteigungswinkel der Schnecke - ø Volumenanteil - dimensionslose Kennzahl     

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.     

Slit die viscometry at shear rates up to 5 × 106s−1: An analytical correction for small viscous heating errors

If the viscosity can be expressed in the form = (T)f(), the walls are at a constant temperatureT ₀, and the extra stress, velocity and temperature fields are fully developed, then the wall shear rate can be calculated by applying the Weissenberg-Rabinowitsch operator toF _c Q instead of to the flow rateQ, whereF _c is a correction factor which differs from 1 when the temperature field is non-uniform; the isothermal equation relating the wall shear stress and pressure gradient is still valid. For the case in whcih = ₀|| ⁿ /(1 +(T–T ₀)), wheren, ₀, and are independent of shear stress and temperatureT, an exact analytical expression forF _c in terms of the Nahme-Griffith numberNa andn is obtained. Use of this expression gives agreement with data obtained for degassed decalin ( = 2.5 mPa s) from a new miniature slit-die viscometer at shear rates up to 5 × 10⁶s^–1; here, the correction is only 7%,Na is 1.3, andGz, the Graetz number at the die exit, is 119. For a Cannon standard liquidS6 ( = 9 mPa s), agreement extends up to 5 × 10⁵s^–1; at 2×10⁶s^–1 (whereNa = 7.2 andGz = 231), the corrections are 11% (measured) and 36% (calculated).Notation x, y Cartesian coordinates - v _x ,v velocity inx-direction, dimensionless velocity - p _xx ,p _yy normal stress onx- andy-planes - N ₁ first normal stress difference - shear stress ony-planes acting inx-direction - _w value of shear stress at the wall - shear rate, shear rate at the wall - Q, Q flow rate (Eqs. (2.13), (2.15)) - T, T ₀ temperature, temperature at the wall - ø, dimensionless temperature (Eqs. (2.24), (2.25)) - h, w half of die height, width of die - R diameter of a tube - , ₀ viscosity, viscosity atT = T ₀ - viscosity-temperature coefficient - k thermal conductivity - c _p specific heat at constant pressure - n, m dimensionless parameters characterizing shear stress dependence of viscosity - Na Nahme Griffith number (Eq. (2.21)) - Gz Graetz number (Eq. (5.1)) - F _c viscous heating correction factor (Eq. (2.18)) - ( ) a function characterizing temperature dependence of viscosity (Eq. (2.8)) - J _k ( ) Bessel function of the first kind This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

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     

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      

Response of an elastic Bingham fluid to oscillatory shear

The response of an elastic Bingham fluid to oscillatory strain has been modeled and compared with experiments on an oil-in-water emulsion. The newly developed model includes elastic solid deformation below the yield stress (or strain), and Newtonian flow above the yield stress. In sinusoidal oscillatory deformations at low strain amplitudes the stress response is sinusoidal and in phase with the strain. At large strain amplitudes, above the yield stress, the stress response is non-linear and is out of phase with strain because of the storage and release of elastic recoverable strain. In oscillatory deformation between parallel disks the non-uniform strain in the radial direction causes the location of the yield surface to move in-and-out during each oscillation. The radial location of the yield surface is calculated and the resulting torque on the stationary disk is determined. Torque waveforms are calculated for various strains and frequencies and compared to experiments on a model oil-in-water emulsion. Model parameters are evaluated independently: the elastic modulus of the emulsion is determined from data at low strains, the yield strain is determined from the phase shift between torque and strain, and the Bingham viscosity is determined from the frequency dependence of the torque at high strains. Using these parameters the torque waveforms are predicted quantitatively for all strains and frequencies. In accord with the model predictions the phase shift is found to depend on strain but to be independent of frequency.Notation A plate strain amplitude (parallel plates) - A _R plate strain amplitude at disk edge (parallel disks) - G elastic modulus - m torque (parallel disks) - M normalized torque (parallel disks) = 2m/R ³₀ - N ratio of viscous to elastic stresses (parallel plates) =µ A/ ₀ ratio of viscous to elastic stresses (parallel disks) =µ A _R/₀ - r normalized radial position (parallel disks) =r/R - r radial position (parallel disks) - R disk radius (parallel disks) - t normalized time = t — /2 - t time - _E elastic strain - _P plate strain (displacement of top plate or disk divided by distance between plates or disks) - _PR plate strain at disk edge (parallel disks) - ₀ yield strain - _E normalized elastic strain = _E/₀ - _P normalized plate strain = _P/₀ - _PR normalized plate strain at disk edge (parallel disks) = _PR/₀ - ₀ normalized plate strain amplitude (parallel plates) =A/ ₀ — normalized plate strain amplitude at disk edge (parallel disks) =A _R/₀ - phase shift between _P andT (parallel plates) — phase shift between _PR andM (parallel disks) - µ Bingham viscosity - stress - ₀ yield stress - T normalized stress =/ ₀ - frequency     

Rheological properties of skim milk gels at various temperatures; interrelation between the dynamic moduli and the relaxation modulus

The rheological properties of rennet-induced skim milk gels were determined by two methods, i.e., via stress relaxation and dynamic tests. The stress relaxation modulusG _c (t) was calculated from the dynamic moduliG andG by using a simple approximation formula and by means of a more complex procedure, via calculation of the relaxation spectrum. Either calculation method gave the same results forG _c (t). The magnitude of the relaxation modulus obtained from the stress relaxation experiments was 10% to 20% lower than that calculated from the dynamic tests.Rennet-induced skim milk gels did not show an equilibrium modulus. An increase in temperature in the range from 20° to 35 °C resulted in lower moduli at a given time scale and faster relaxation. Dynamic measurements were also performed on acid-induced skim milk gels at various temperatures andG _c (t) was calculated. The moduli of the acid-induced gels were higher than those of the rennet-induced gels and a kind of permanent network seemed to exist, also at higher temperatures. G storage shear modulus,N·m^–2; - G loss shear modulus,N·m^–2; - G _c calculated storage shear modulus,N·m^–2; - G _c calculated loss shear modulus,N·m^–2; - G _e equilibrium shear modulus,N·m^–2; - G _ec calculated equilibrium shear modulus,N·m^–2; - G(t) relaxation shear modulus,N·m^–2; - G _c (t) calculated relaxation shear modulus,N·m^–2; - G ^*(t) pseudo relaxation shear modulus,N·m^–2; - H relaxation spectrum,N·m^–2; - t time,s; - relaxation time,s; - angular frequency, rad·s^–1. Partly presented at the Conference on Rheology of Food, Pharmaceutical and Biological Materials, Warwick, UK, September 13–15, 1989 [33].     

Apparent thickening behavior of dilute polystyrene solutions in extensional flows

An experimental investigation was undertaken to study the apparent thickening behavior of dilute polystyrene solutions in extensional flow. Among the parameters investigated were molecular weight, molecular weight distribution, concentration, thermodynamic solvent quality, and solvent viscosity. Apparent relative viscosity was measured as a function of wall shear rate for solutions flowing from a reservoir through a 0.1 mm I.D. tube. As increased, slight shear thinning behavior was observed up until a critical wall shear rate was exceeded, whereupon either a large increase in or small-scale thickening was observed depending on the particular solution under study. As molecular weight or concentration increased, decreased and, the jump in above , increased. increased as thermodynamic solvent quality improved. These results have been interpreted in terms of the polymer chains undergoing a coil-stretch transition at . The observation of a drop-off in at high (above ) was shown to be associated with inertial effects and not with chain fracture due to high extensional rates.     

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     

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     

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     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase