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     

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     

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].     

Rheological characteristics of glass-fibre-filled polypropylene melts

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

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     

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     

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     

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     

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

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

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     

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.     

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.     

The rotation of a suspended axisymmetric ellipsoid in a magnetic field

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

Experimental studies of suspended ferromagnetic fibers in a magnetic field

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

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     

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     

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.     

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.     

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.