Steady shear banding in complex fluids

Shear banding (SB) is manifested by the abrupt “demixing” of the flow into regions of high and low shear rate. In this paper, we first relate analytically the rheological parameters of the fluid with the range of shear rates and stresses of SB occurrence. For this, we accept that the origin of shear banding is constitutive, and adopt a non-linear viscoelastic expression able to accommodate the double-valuedness of the stress with flow intensity, under certain conditions. We then implement the model for the case of pressure driven flow through a cylindrical pipe; we derive approximate expressions for the velocity profile in the two-banded regions (core and outer annular), the overall throughput in the presence or absence of “spurt”, and the radial location limits of the shear rate discontinuity.     

A study on coherent structures and drag-reduction in the wall turbulence with polymer additives by TRPIV

An experimental measurement was performed using time-resolved particle image velocimetry (TRPIV) to investigate the spatial topological character of coherent structures in wall-bounded turbulence of polymer additive solution. The fully developed near-wall turbulent flow fields with and without polymer additives at the same Reynolds number were measured by TRPIV in a water channel. The comparisons of turbulent statistics confirm that due to viscoelastic structure of long-chain polymers, the wall-normal velocity fluctuation and Reynolds shear stress in the near-wall region are suppressed significantly. Furthermore, it is noted that such a behavior of polymers is closely related to the decease of the motion of the second and forth quadrants, i.e., the ejection and sweep events, in the near-wall region. The spatial topological mode of coherent structures during bursts has been extracted by the new mu-level criteria based on locally averaged velocity structure function. Although the general shapes of coherent structures are unchanged by polymer additives, the fluctuating velocity, velocity gradient, velocity strain rate and vorticity of coherent structures during burst events are suppressed in the polymer additive solution compared with that in water. The results show that due to the polymer additives the occurrence and intensity of coherent structures are suppressed, leading to drag reduction.     

Effect of pressure-dependent slip on flow curve multiplicity

Various microstructural pictures for slip at polymer/solid interfaces lead to relations which have a region where multiple values of slip velocity are predicted for the same shear stress. This leads to the expectation of multivalued flow curves, which has been verified in specific cases by numerous researchers. We study the effect of pressure dependence on flow curve multiplicity using a simple multivalued slip relation to model the phenomena of hysteresis and spurt flow in polymer extrusion. A continuation technique is used to trace out the boundaries of the region of flow curve multiplicity as pressure drop and die length to diameter (L/D) ratio are changed. Results for Newtonian, shear thinning and viscoelastic constitutive equations show that, despite the multivalued nature of the slip model, multiplicity (and thus hysteresis) is absent at high L/D.  For the sake of completeness, we also carry out time-dependent simulations at constant piston speed taking fluid compressibility into account. These simulations show that oscillations in the pressure drop and exit volumetric flow rate result only if the system is operated in the multiplicity region of the steady state flow curve, in agreement with the results of similar simulations by researchers using various multivalued slip models without pressure dependence. The results demonstrate that a multivalued slip model does not guarantee multiplicity in the flow curve for the constant pressure drop operation, nor oscillations for constant piston speed operation. Received: 18 August 1997 Accepted: 30 March 1998     

Theoretical and experimental investigation of the flow of a viscoelastic fluid in the gap between two rotating disks

A theory of the nonlinear viscoelastic behavior of polymer fluids has been constructed in [1]. The theory was used in [2] to investigate the motion of a nonlinear viscoelastic medium under steady and unsteady deformation rates in simple shear flow, and a comparison was made with experiment. The experiments in [2], which were performed on a cone-plate Weissenberg rheogoniometer, indicate that this arrangement is unsuitable for measurements of normal stresses under unsteady conditions in fluids with a fairly high viscosity. Below, we will show the suitability of using a disk-disk Weissenberg rheogoniometer to measure normal stresses in this case for unsteady conditions (transition to steady flow and stress relaxation). In this regard, a theoretical study of the flow of a viscoelastic fluid in the gap between rotating disks is needed. Note that in this case new information will be obtained from a comparison with simple uniform shear flow, since in the flow of a polymer between two disks all three normal stress components contribute to the axial force, while in the gap between a cone and a plate only the first normal stress difference contributes to the normal force.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 25–30, March–April, 1976.     

A rheological model for materials which support coexistent shear rates

In this work we study a version of the three constant differential-type Oldroyd constitutive relation which allows distinct objective time derivatives for the extra stress and the stretching. We integrate the constitutive equation and determine an equivalent history integral representation for this model for the general class of viscometric motions. For certain choices of the material parameters and initial conditions, we find that this model allows for the development of shear rate discontinuities in the flow domain as a steady viscometric flow is achieved. Correspondingly, we also give evidence that intense shear rate oscillations may occur during the transient period as an impulsively started viscometric flow in a channel tends to a steady state under a constant critical shear stress. This critical shear stress lies in an interval of values for which the material experiences the phenomenon of “flow yielding”. A qualitative comparison with experimental data is made for certain creams and greases. The material instabilities inherent in this constitutive theory for viscometric motions are suggestive of the instabilities that occur in many viscoelastic fluids such as sharkskin patterns, wavy fracture, and spurt flow.     

The transmission of turbulent boundary layer unsteady pressure and shear stress through a viscoelastic layer

The transmission of unsteady pressure and shear stress, generated by a turbulent boundary layer in water, through a viscoelastic layer backed by a rigid plate is investigated. Analytical models are used to estimate the unsteady pressure and shear stress from 10 to 1000 Hz for a flat plate boundary layer with zero pressure gradient. Additionally, models for the transfer of the unsteady pressures and shear stress through the viscoelastic layer are developed. The models are used to predict the unsteady pressure fluctuations, or flow noise, which would be seen by a finite size sensor embedded under the elastomer layer. The unsteady pressure levels are found to be 20 dB greater than the unsteady shear stress levels across all frequency ranges computed, in agreement with recent measurements. The unsteady pressure transfer functions have a peak at the shear wavenumber and are larger than the shear stress transfer magnitudes from 10 to 50 Hz. The unsteady shear stress transfer functions have a peak at the acoustic wavenumber and are larger than the pressure transfer magnitudes from 50 to 1000 Hz. Over the frequency range examined, the unsteady pressures were found to be the dominant contributor to the sensor flow noise due to the considerably larger magnitude of the unsteady pressures on the top of the viscoelastic layer.     

Viscoelastic plastic rheological model for particle filled polymer melts

A viscoelastic plastic model for suspension of small particles in polymer melts has been developed. In this model, the total stress is assumed to be the sum of stress in the polymer matrix and the filler network. A nonlinear viscoelastic model along with a yield criterion were used to represent the stresses in the polymer matrix and the filler network, respectively. The yield function is defined in terms of differential equations with an internal parameter. The internal parameter models the evolution of structure changes during floc rupture and restoration. The theoretical results were obtained for steady and oscillatory shear flow and compared with experimental data for particle filled thermoplastic melt. The experimental data included the steady state shear strress over a wide range of shear rates, the transient stress in a start up shear flow, stress relaxation after cessation of a steady state shear flow, the step shear and the oscillatory shear flow at various amplitudes.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

Asymptotic analysis of an end-loaded, transversely isotropic, elastic, semi-infinite strip weak in shear

We use linear elasticity to study a transversely isotropic (or specially orthotropic), semiinfinite slab in plane strain, free of traction on its faces and at infinity and subject to edge loads or displacements that produce stresses and displacements that decay in the axial direction. The governing equations (which are identical to those for a strip in plane stress, free of traction on its long sides and at infinity, and subject to tractions or displacements on its short side) are reduced, in the standard way. to a fourth-order partial differential equation with boundary conditions for a dimensionless Airy stress function ƒ. We study the asymptotic solutions to this equation for four sets of end conditions—traction, mixed (two), displacement—as g3, the ratio of the shear modulus to the geometric mean of the axial and transverse extensional moduli, approaches zero. In all cases, the solutions for ƒ consist of a “wide” boundary layer that decays slowly in the axial direction (over a distance that is long compared to the width of the strip) plus a “narrow” boundary layer that decays rapidly in the axial direction (over a distance that is short compared to the width of the strip). Moreover, we find that the narrow boundary layer has a “sinuous” part that varies rapidly in the transverse direction, but which, to lowest order, does not enter the boundary conditions nor affect the transverse normal stress or the displacements. Because the exact biorthogonality condition for the cigenfunctions associated with ƒ can be replaced by simpler orthogonality conditions in the limit as →b 0, we are able to obtain, to lowest order, explicit formulae for the coeflicients in the eigenfunction expansions of ƒ for the four different end conditions.     

Rheology and flow-birefringence from viscoelastic polymer-clay solutions

 The shear orientation of viscoelastic clay-polymer solutions was investigated by means of rheology and flow birefringence (Δn). The polymer chains are in dynamic adsorption/desorption equilibrium with the clay particles to form a “network”. The elastic behavior of the network was characterized by constant stress, oscillatory shear, and stress relaxation experiments. Constant stress experiments indicated a yield stress upon which shear flow started and no strain recovery could be observed. Oscillatory shear experiments showed a broad elastic region followed by flow when a critical strain was reached. Stress relaxation experiments showed several relaxation times when the same critical strain was reached. Experiments under steady flow characterized the transient behavior of the network. With increasing steady shear rate a pronounced minimum in birefringence was observed at a critical shear rate. The shear rate dependent viscosity showed near power law behavior and no corresponding critical feature. While birefringence detects orientational effects on a microscopic length scale, rheology averages over macroscopic changes in the sample. The same degree of orientation could be achieved under constant shear rate or constant stress conditions. Received: 25 January 2001 Accepted: 22 May 2001     

Constitutive relations for the shear band evolution in granular matter under large strain

A so-called ＂split-bottom ring shear cell＂ leads to wide shear bands under slow, quasi-static deformation. Unlike normal cylindrical Couette shear cells or rheometers, the bottom plate is split such that the outer part of it can move with the outer wall, while the other part （inner disk） is immobile. From discrete element simulations （DEM）, several continuum fields like the density, velocity, deformation gradient and stress are computed and evaluated with the goal to formulate objective constitutive relations for the powder flow behavior. From a single simulation, by applying time- and （local） space-averaging, a non-linear yield surface is obtained with peculiar stress dependence. The anisotropy is always smaller than the macroscopic friction coefficient. However, the lower bound of anisotropy increases with the strain rate, approaching the maximum according to a stretched exponential with a specific rate that is consistent with a shear path of about one particle diameter.     

Modeling of blood rheological behavior in unsteady-state shearing flows

Issues of blood flow modeling under unsteady-state conditions at moderate shear rates are considered using a blood rheological model accounting both for the viscoelastic properties and the thixotropy caused by erythrocyte aggregation. The resulting shear stress versus time relations for single shear rate steps and the dependence of the complex viscosity components on the shear rate amplitude in oscillating shear flow show good qualitative correspondence with the experimental data reported in the literature.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 26–30, November–December, 1995.     

Structure and rheology of wormlike micelles

In a semi-dilute aqueous solution under certain conditions, surfactant molecules will self assemble to form wormlike micelles. The micelles are dynamic in structure since they can break and reform, providing an additional mode of relaxation. The viscoelastic properties of the wormlike micelles can be predicted using simple theological models. For many surfactant solutions the mechanical data can be related to the optical data by the stress-optical rule. From the viscoelastic data it is possible to estimate the breaking time of the micelle. The techniques of birefringence and small angle light scattering are used to study the microstructure of a surfactant solution under simple shear and extensional flow. The sample under investigation is a solution of cetyltrimethylammonium bromide and sodium salicylate in water, with a salt to surfactant ratio of 7.7. Below a critical shear rate, the birefringence increases linearly with shear rate and the stress-optical rule is valid. The SALS patterns reveal distinctive butterfly patterns indicating that scattering is a result of concentration fluctuations that moderately couple to the flow. However, above a critical shear rate the birefringence plateaus and the stress-optical rule is no longer valid. SALS patterns show both a bright streak and a butterfly pattern. The bright streak is caused by elongated structures aligned in the direction of the flow. The oriented structures occur when the characteristic time of flow is faster than the breaking time of the micelles.Dedicated to Prof. Dr. J. Meissner on the occasion of his retirement from the chair of Polymer Physics at the Eidgenössische Technische Hochschule (ETH) Zürich, Switzerland     

On the Generation of Discontinuous Shearing Motions of a Non-Newtonian Fluid

The system under study models unsteady, one-dimensional shear flow of a highly elastic and viscous incompressible non-Newtonian fluid with fading memory under isothermal conditions. The flow, in a channel, is driven by a constant pressure gradient, is symmetric about the center line, and satisfies a no-slip boundary condition at the wall. The non-Newtonian contribution to the stress is assumed to obey a differential constitutive law (due to Oldroyd, Johnson & Segalman), the key feature of which is a non-monotone relation between the total steady shear stress and strain rate. In a regime in which the Reynolds number is much smaller than the Deborah (or Weissenberg) number, one obtains a degenerate, singularly perturbed system of nonlinear reaction-diffusion equations. It is shown that if the driving pressure gradient exceeds a critical value (the local shear stress maximum of the steady stress vs. strain rate relation), then the solution to the governing system, starting from rest at , tends as to a particular discontinuous steady state solution (the “top-jumping” steady state), except in a small neighborhood of the discontinuity. This discontinuous steady state is shown to be nonlinearly stable in a precise sense with respect to perturbations yielding smooth initial data. Such discontinuous steady states have been proposed to explain “spurting” flows, which exhibit a large increase in mean flow rate when the driving pressure is raised above a critical value. (Accepted April 22, 1996)     

Elongational flow behavior of polymeric fluids according to the Leonov model

The viscoelastic behavior of polymeric systems based upon the Leonov model has been examined for (i) the stress growth at constant strain rate, (ii) the stress growth at constant speed and (iii) the elastic recovery in elongational flow. The model parameters have been determined from the available rheological data obtained either in steady shear flow (shear viscosity and first normal-stress difference as a function of shear rate) or oscillatory flow (storage and loss moduli as a function of frequency in the linear region) or from extensional flow at very small strain rates (time-dependent elongation viscosity in the linear viscoelastic limit). In addition, the effect of the parameter characterizing the strain-hardening of the material during elongation has also been studied. The estimation of this parameter has been based upon the structural characteristics of the polymer chain which include the critical molecular weight and molecular weight of an independent segment. Five different polymer melts have been considered with varying number of modes (maximum four modes). Resulting predictions are in fair agreement with corresponding experimental data in the literature.     

Characterization of micellar structure dynamics for a drag-reducing surfactant solution under shear: normal stress studies and flow geometry effects

Some surfactant solutions have been observed to exhibit a strong drag reduction behavior in turbulent flow. This effect is generally believed to result from the formation of large cylindrical micelles or micellar structures. To characterize and understand better these fluids, we have studied the transient rheological properties of an efficient drag-reducing aqueous solution: tris (2-hydroxyethyl) tallowalkyl ammonium acetate (TTAA) with added sodium salicylate (NaSal) as counter ion. For a 5/5 mM equimolar TTAA/NaSal solution, there is no measurable first normal stress difference (N ₁) immediately after the inception of shear, but N ₁ begins to increase after a well-defined induction time — presumably as shear-induced structures (SIS) are formed — and it finally reaches a fluctuating plateau region where its average value is two orders of magnitude larger than that of the shear stress. The SIS buildup times obtained by first normal stress measurements were approximately inversely proportional to the shear rate, which is consistent with a kinetic process during which individual micelles are incorporated through shear into large micellar structures. The SIS buildup after a strong preshear and the relaxation processes after flow cessation were also studied and quantified with first normal stress difference measurements. The SIS buildup times and final state were also found to be highly dependent on flow geometry. With an increase in gap between parallel plates, for example, the SIS buildup times decreased, whereas the plateau viscosity increased.     

Some experimental results on the development of Couette flow for non-Newtonian fluids

This paper reports some experimental results on the time development of a Couette flow following the start-up of shear flow using the technique of two-color flow birefringence. Measurements obtained on collagen solutions are consistent with two theoretical studies which predict that for some viscoelastic liquids, momentum is transferred from the moving Couette cell boundary to the interior of the fluid through a velocity wave propagating and reflecting between the cell boundaries. This non-Newtonian phenomenon, exhibited as an oscillatory response in the measured birefringence and orientation angle, is observed at shear rates above a critical value when the response time of the polymer solution approaches the flow development time in the Couette flow cell.     

Numerical and experimental investigation of turbulent characteristics in a drag-reducing flow with surfactant additives

Cetyltrimethyl ammonium chloride (CTAC) surfactant additives, because of their long-life characteristics, can be used as promising drag-reducers in district heating and cooling systems. In the present study we performed both numerical and experimental tests for a 75 ppm CTAC surfactant drag-reducing channel flow. A two-component PIV system was used to measure the instantaneous streamwise and wall-normal velocity components. A Giesekus constitutive equation was adopted to model the extra stress due to the surfactant additives, with the constitutive parameters being determined by well-fitting apparent shear viscosities, as measured by an Advanced Rheometric Expansion System (ARES) rheometer. In the numerical study, we connected the realistic rheological properties with the drag-reduction rate. This is different from previous numerical studies in which the model parameters were set artificially. By performing consistent comparisons between numerical and experimental results, we have obtained an insight into the mechanism of the additive-induced drag-reduction phenomena.

Our simulation showed that the addition of surfactant additives introduces several changes in turbulent flow characteristics: (1) In the viscous sublayer, the mean velocity gradient becomes gentler due to the viscoelastic forces introduced by the additives. The buffer layer becomes expanded and the slope of the velocity profile in the logarithmic layer increases. (2) The locations where the streamwise velocity fluctuation and Reynolds shear stress attain their maximum value shifted from the wall region to the bulk flow region. (3) The root-mean-square velocity fluctuations in the wall-normal direction decrease for the drag-reducing flow. (4) The Reynolds shear stress decreases dramatically and the deficit of the Reynolds shear stress is mainly compensated by the viscoelastic shear stress. (5) The turbulent production becomes much smaller and its peak-value position moves toward the bulk flow region. All of these findings agree qualitatively with experimental measurements.

Regarding flow visualization, the violent streamwise vortices in the near wall region become dramatically suppressed, indicating that the additives weaken the ejection and sweeping motion, and thereby inhibit the generation of turbulence. The reduction in turbulence is accomplished by additive-introduced viscoelastic stress. Surfactant additives have dual effects on frictional drag: (1) introduce viscoelastic shear stress, which increases frictional drag; and (2) dampen the turbulent vortical structures, decrease the turbulent shear stress, and then decrease the frictional drag. Since the second effect is greater than the first one, drag-reduction occurs.     

Channel, tube, and Taylor–Couette flow of complex viscoelastic fluid models

We show how to formulate two-point boundary value problems to compute laminar channel, tube, and Taylor–Couette flow profiles for some complex viscoelastic fluid models of differential type. The models examined herein are the Pom-Pom Model [McLeish and Larson 42:81–110, (1998)] the Pompon Model [Öttinger 40:317–321, (2001)] and the Two Coupled Maxwell Modes Model (Beris and Edwards 1994). For the two-mode Upper-Convected Maxwell Model, we calculate analytical solutions for the three flow geometries and use the solutions to validate the numerical methodology. We illustrate how to calculate the velocity, pressure, conformation tensor, backbone orientation tensor, backbone stretch, and extra stress profiles for various models. For the Pom-Pom Model, we find that the two-point boundary value problem is numerically unstable, which is due to the aphysical non-monotonic shear stress vs shear rate prediction of the model. For the other two models, we compute laminar flow profiles over a wide range of pressure drops and inner cylinder velocities. The volumetric flow rate and the nonlinear viscoelastic material properties on the boundaries of the flow geometries are determined as functions of the applied pressure drop, allowing easy analysis of experimentally measurable quantities.     

Characteristic turbulent structure of a modified drag-reduced surfactant solution flow via dosing water from channel wall

The modification of the near-wall structure is very important for the control of wall turbulence. To ascertain the effect of near-wall modulation on the viscoelastic drag-reduced flow, the modified characteristics of a surfactant solution channel flow were investigated experimentally. The modulation was conducted on the boundary of the channel flow by injecting water from the whole surface of one side of the channel wall. The diffusion process of the injected water was observed by using the planar laser-induced fluorescence technique. The velocity statistics and characteristic structure including the spatial distributions of instantaneous streamwise velocity, swirling strength, and Reynolds shear stress were analyzed based on the velocity vectors acquired in the streamwise wall-normal plane by using the particle imaging velocimetry technique. The results indicated that the disturbance of the injected water was constricted within a finite range very near the dosing wall, and the Reynolds shear stress was increased in this region. However, the eventual drag reduction rate was found to be increased due to a relatively large decrement of viscoelastic shear stress in this near-wall region. Moreover, the flow structure under this modulation presented obvious regional characteristics. In the unstable disturbed region, the mixing of high-speed and low-speed fluids and the motions of ejection and sweep occurred actively. Many clockwise vortex cores were also found to be generated. This characteristic structure was similar to that in the ordinary turbulence of Newtonian fluid. Nevertheless, outside this disturbed region, the structure still maintained the characteristics of the drag-reduced flow with non-Newtonian viscoelastic additives. These results proved that the injected Newtonian fluid associated with the modified stress distribution creates a diverse characteristic structure and subsequent enhanced drag reduction. This investigation can provide the experimental basis for further study of turbulence control.