Steady and modulated flow of an Ostwald fluid around a rotating disk

The solutions of the continuity equation and the equations of motion of the flow in the vicinity of a rotating disk have been established for an Ostwald fluid under steady-state conditions and in modulated flow around a mean value. Under steady-state conditions, the kinematics of the flow is scarcely dependent on the rheological parameters close to the disk, however, for n < 1 long-range effects have been put forward. For modulated flow, in the high-frequency range, a behaviour very different from that observed for a Newtonian fluid was found. In the low-frequency range an asymptotic solution has been proposed which is of special interest in mass transfer problems.Presented at the Second Conference of European Rheologists, Prague, June 17–20, 1986     

Viscometric flows of polymeric melts treated according to the rational theory of Eckart's continuum. I. General theory

The purely rational theory of Eckart continua (i.e. elastic bodies with a variable relaxed state) is applied to viscometric flows of polymeric melts. The main assumptions are thermodynamic non-interaction of inelastic behaviour and of non-elastic stress, as well as elastic isotropy. After establishing the time-dependent differential equations of viscometric flow, these equations are simplified to a set of algebraic equations describing steady-state flow. From this we deduce two general equations connecting the three elastic steady-state viscometric functions which do not depend upon the elastic behaviour. The law of rubber elasticity used in this paper is described in the Appendix.     

Response of blood flow through an artery under stenotic conditions

This paper presents an analytical study on the behavoiur of blood flow in an artery having a stenosis. This is basically formulated through the use of a suitable mathematical model. The arterial segment under consideration is simulated by an anisotropically elastic cylindrical tube filled with a viscous incompressible fluid representing blood. The analysis is carried out for an artery with mild local narrowing in its lumen forming a stenosis. Particular emphasis has been paid to the effect of the surrounding connective tissues on the motion of the arterial wall. Blood is treated as a Newtonian fluid. The analysis is restricted to propagation of small amplitude harmonic waves, generated due to the flow of blood whose wave length is large compared to the radius of the arterial segment. The effect of the shape of stenosis on the resistance to blood flow has been well illustrated quantitatively through numerical computations of the resulting expressions. A quantitative analysis is also made for the variation of the phase velocity, as well as the velocity of wave propagation and the flow rate, in order to illustrate the applicability of the model.     

The flow past a sphere in a cylindrical tube: effects of intertia,shear-thinning and elasticity

Numerical solutions are presented for the flow past a sphere placed at the centreline of a cylindrical tube for Reynolds numbers ranging from 0 to 150, using a boundary element method. Fluids are modelled by a variety of constitutive equations including the Newtonian, the Carreau and the Phan-Thien-Tanner models. The influence of inertia, shear-thinning and fluid elasticity on the flow field, drag and the pressure drop force-drag ratio is examined. Some results are compared with available experimental data.     

A model for the thermorheological behavior of viscoelastic fluids

Using a power-law ansatz for the temperature dependence of the shear modulus on the level of internal variables, the thermorheological behavior is modeled for viscoelastic fluids of a special group of rheological constitutive equations (rate-type models). The model parameter introduced characterizes thermoelastic contributions. The relation between the model parameter and the physical quantities appearing in deformation processes is discussed. Based on the chosen temperature dependence of the shear modulus, thermodynamically consistent equations like the nonlinear rheological constitutive equation and the temperature equation are derived. The special cases of entirely entropy and energy elastic fluids are also considered. The thermorheological behavior (exo-, - or endothermal processes) of a viscoelastic fluid in a stress-growth experiment followed by relaxation is analyzed with respect to the model parameter.     

A thermodynamic approach to non-linear viscoelasticity

We show that Extended Irreversible Thermodynamics (EIT) may be regarded as a unifying thermodynamic scheme which generates a variety of non-linear viscoelastic equations of state. In particular, we show how the so-called eight-constants Oldroyd model, with all its invariant properties, may be derived from the postulates of EIT. It is also shown that the structure of the thermodynamic based rheological equations of state is in complete agreement with that of the representation theorems formulation of constitutive relationships originally proposed by Rivlin and Ericksen. We discuss how more complex rheological models can be derived from higher order approximations of the theory. Finally, the connections between our approach and other related works are briefly discussed.     

The unraveling of a polymer chain in a strong extensional flow

By assuming that in a strong extensional flow a polymer molecule in dilute solution is quickly driven into a folded or kinked state in which drag and entropic elastic forces dominate over the Brownian force, we derive kink dynamics equations that describe the unraveling of the molecule in the extensional flow. Solving these equations numerically, we find that although the ends of the chain move, on average, affinely in the flow field until the chain is unfolded to about 1/3 of its fully extended length, large viscous stresses are produced because the solvent must flow around nonextending strands of polymer that lie between neighboring kinks. These predictions are compared with available experimental data and with other theoretical models.     

On the cone and plate deformation of a rubber-like material: Instability

In this report it is shown that the place boundary-value problem, for a small deformation superimposed on the large cone-and-plate deformation of a Mooney-Rivlin material, has no unique solution at some displacement angles of the cone. The kinematics of the small deformation are chosen such that the problem is reduced to a set of ordinary differential equations. With respect to this restricted class of kinematics the cone-and-plate deformation is unstable.     

A boundary element simulation of the unsteady motion of a sphere in a cylindrical tube containing a viscoelastic fluid

A boundary element method is used to simulate the unsteady motion of a sphere falling under gravity along the centreline of a cylindrical tube containing a viscoelastic fluid. The fluid is modelled by the upper-convected Maxwell constitutive equation. Results show that the viscoelasticity of the liquid leads to a damped oscillation in sphere velocity about its terminal value. The maximum sphere velocity, which occurs in the first overshoot, is approximately proportional to the square root of the Weissenberg number when the ratio of the sphere radius to the tube radius is sufficiently small. Particular attention is also paid to the wall effects. It is shown that a closer wall reduces the oscillatory amplitude of the sphere velocity but increases its frequency. The results suggest that the falling-ball technique, which is now widely used for viscosity measurement, might also be used for the determination of a relaxation time for a viscoelastic fluid.     

Extensional flow resistance of dilute polyacrylamide and surfactant solutions

The extensional flow behaviour of dilute aqueous solutions of a partiallyhy-drolyzed polyacrylamide and a surfactant were investigated in an extensional flow cell. The cell was designed such that fluids were subjected to steady shear before undergoing extensional motion in a converging channel. Extensional resistance was monitored by measuring the pressure drop through the channel. Such measurements were made over a range of extensional rates at fixed values of the upstream shear rate. Solutions of different concentrations were tested — up to 40 ppm of polyacrylamide and 450 ppm of surfactant — at various temperatures in the case of surfactant and for different types and amounts of salt in the case of polyacrylamide. Of the results, the more notable are that the extensional resistance of polyacrylamide solutions is affected much more by CaCl₂ than by NaCl and that surfactant solutions do not exhibit extensional resistance unless they are pre-sheared.     

Statistical mechanics of temporary polymer networks

The statistical theory of temporary polymer networks developed in parts I and II of this work is closed by a discussion on dynamical effects, including comparison with experiments. After a short review of the theory developed in parts I and II we discuss the situation in which a prescribed velocity gradient is imposed on the physical system formed by the temporary network and the solvent. The dynamical equation for the 2nd moments (derived from a generalized bead-spring model) contains the configuration-dependent transition probability, which depends on the second moments in a complex way. An approximate solution of these equations is obtained from a computer program. It shows the experimentally observed behavior, in particular the stress overshoot maximum after a sudden start of the flow.The derivation of the above equations contains simplifying assumptions which, however, leave the essential physics intact. The most important assumptions concern the relaxation time approach and the decay and formation processes. These were supposed to be dilute enough to break up the many-junction processes into one-junction processes. Both assumptions are analogous to frequently used assumptions in Boltzmann's kinetic theory.     

The rheology of powders

This paper studies the slow flow of powders. It is argued that since powders can flow like liquids, there must be equations similar to those of liquids. The phenomenon of a variable density, dilatancy, is described by an analogue of temperature called the compactivity X. Whereas, in thermal physicsT = E/S, powders are controlled byX = V/S. The equations for, v, T of a liquid are replaced by, v, X. An analogy for free energy is described, and the solution to some simple problems of packing and mixing are offered. As an example of rheology, it is shown that the simplest flow equations produce a transition to plug flow in appropriate circumstances.Delivered as a Gold Medal Lecture at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.     

Newtonian stratified flow through an abrupt expansion

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

Numerical simulation of delayed die swell

In a recent paper, Joseph et al. showed that, for a number of viscoelastic fluids, one can observe the phenomenon of delayed die swell beyond a critical extrusion velocity, or beyond a critical value of the viscoelastic Mach number. Giesekus had also observed that delayed die swell is a critical phenomenon.In the present paper, we find a set of material and flow parameters under which it is possible to simulate delayed die swell. For the viscoelastic flow calculation, we use the finite element algorithm with sub-elements for the stresses and streamline upwinding in the discretized constitutive equations. For the free surface, we use an implicit technique which allows us to implement Newton's method for solving the non-linear system of equations. The fluid is Oldroyd-B which, in the present problem, is a singular perturbation of the Maxwell fluid. The results show very little sensitivity to the size of the retardation time. We also show delayed die swell for a Giesekus fluid.This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

A linear analysis of instabilities in Couette flow of nematic liquid crystals

This paper presents a new method for determining the critical threshold values at which both homogeneous and roll-type instabilities may appear when thin layers of nematic liquid crystal are subjected to Couette flow. The anisotropic axis of the material is initially uniformly aligned parallel to the axial direction and the effect of a stabilising, axially oriented, uniform magnetic field is allowed for in calculations. Carrying out a linear stability analysis of the continuum equations, we employ a Fourier series method to derive an expression for determining the critical shear rate at which instability occurs. Thresholds for both homogeneous and roll instabilities are obtained and their variation with rotation rate and applied magnetic field strength is examined for the material MBBA.     

Optical and hydrodynamical regimes in dilute solutions of polyelectrolytes

A careful analysis of the equations governing dynamical disturbances in dilute solutions of polyelectrolytes (described on the basis of a conformational model with electromechanical interactions) allows one to place in evidence two extreme regimes, one called the optical regime, which justifies the use of optical methods for measuring the main instantaneous characteristics of the solutions, and the other that corresponds to electrohydrodynamical flows where a hierarchy can be exhibited for the various effects contributing to the flow equations.     

Velocity profiles in rectangular channel flow of liquid crystalline polymer solutions

The motion of aqueous solutions of hydroxypropylcellulose in a shallow channel was visualized by means of an electrochemical technique. Isotropic and anisotropic solutions of approximately the same viscosity were used. The velocity profiles of the liquid crystalline solution were found to be qualitatively very different from those of the isotropic one. In particular, anomalous maxima in the velocity near the side walls were observed.     

A simplified cooperative model of stress relaxation and other consolidation processes in solids

This paper explores the properties of a relaxation function derived from a differential equation mimicking the distribution mechanism of Bose-Einstein statistics in the time domain. Within a significant portion of the process, the relaxation quantity n decreases linearly with log time. The relation between dn/dt and n is an exponential one. In this respect, the present approach produces results largely equivalent to those obtained using the hypothesis of stress-dependent thermal activation or a box-like spectrum of relaxation times, . The spectrum of the model proposed here is discrete, with integer valued fractions of a characteristic -centering the equations.     

The motion of rodlike particles in the pressure-driven flow between two flat plates

The motion of freely suspended rodlike particles has been observed in the pressure-driven flow between the two flat plates of a Hele Shaw flow cell at low Reynolds numbers. Data are reported for rodlike particles with aspect ratios of 12.0 suspended in a Newtonian fluid for gap thickness to particle length ratios of 3, 6, and 20; and for rodlike particles with aspect ratios between 5 and 8 in a non-Newtonian fluid (79.25 wt.% water, 20.2 wt.% glycerine, and 0.55 wt.% polyacrylamide). For the Newtonian fluid, the time-dependent orientation of the particles near and far from walls was shown to be in quantitative agreement with Jeffery's theory for ellipsoids suspended in a simple shear flow if an effective aspect ratio is calculated from the experimental period of rotation. Particles aligned with the flow direction and less than a particle half-length from a wall interacted irreversibly with the wall. For the non-Newtonian fluid, the timedependent orientation far from a wall was shown to be in qualitative agreement with Leal's theory for a second-order fluid; however, particles that were aligned with the flow direction and were near walls did not rotate.