Shear rheology of polymer solutions near the critical condition for elastic instability

The use of constant viscosity, highly elastic polymer solutions, so called Boger fluids, has been remarkably successful in elucidating the behavior of polymeric materials under flowing conditions. However, the behavior of these fluids is still complicated by many different physical processes occurring within a narrow window of observation time and applied shear rate. In this study, we investigate the long-time shear behavior of an ideal Boger fluid: a well characterized, athermal, dilute, binary solution of high molecular weight polystyrene in oligomeric polystyrene. Rheological measurements show that under an applied steady shear flow, this family of polymer solutions undergoes a transient decay of normal stresses on a timescale much longer than the polymer molecule's relaxation time. Rheological and flow visualization results demonstrate that the observed phenomenon is not caused by polymer degradation, phase separation, viscous heating, or secondary flows from elastic instabilities. Although the timescale is much shorter than that associated with polymer migration in the same solutions (MacDonald and Muller, 1996), the appearance of this phenomenon only at the rates where migration has been observed suggests that it may be a prerequisite for observing migration. In addition, we note that through sufficient preshearing of the sample, the normal stress decrease suppresses the elastic instability. These results show that there is considerable uncertainty in choosing the appropriate measure of the fluid relaxation time for consistently modeling the critical condition for the elastic instability, the decay of normal stresses, and the migration of polymer species.     

The importance of downstream events in microfluidic viscoelastic entry flows: Consequences of increasing the constriction length

Polymer solutions and melts can exhibit large upstream corner and lip vortices, unstable and diverging flow and an enhanced pressure drop when flowing through a geometry containing a constriction. In the present work, we use a planar microfluidic device to show that the length of the downstream constriction plays an important role in the upstream kinematics and the extra pressure drop. That is, the elastic flow phenomena observed upstream of a constriction during entry flows of polymer solutions are not exclusively a result of the stretching dynamics induced by the converging flow—the downstream relaxation events are, at least, equally important. Flow visualization experiments with semi-dilute solutions of a high molecular weight polymer showed that large stable symmetric vortices could be reduced to highly chaotic asymmetric flow, merely by increasing the length of the constriction—the Reynolds number and elasticity number were both held constant. This was accompanied by a higher extra pressure. These results support the hypothesis that elastic flow instabilities originate downstream of the constriction (at the expansion) and move progressively upstream with time and/or flowrate. These findings may also partly explain the discrepancies commonly observed between the results of entry flow experiments and numerical simulations, in which the downstream geometry is very rarely considered. Lastly, we illustrate how to minimize the occurrence of unstable flow upstream of a constriction, which is a necessary condition for closed microrheometry devices used to characterize low viscosity elastic fluids.     

Displacement of yield-stress fluids in a fracture

We consider a displacement of several yield-stress fluids in a Hele-Shaw cell. The topic is relevant to the development of a model for the flow of multiple phases inside a narrow fracture with application to hydraulically fracturing a hydrocarbon-bearing underground formation. Existing models for fracturing flows include only pure power-law models without yield stress, and the present work is aimed at filling this gap. The fluids are assumed to be immiscible and incompressible. We consider fluid advection in a plane channel in the presence of density gradients. Gravity is taken into account, so that there can be slumping and gravitational convection. We use the lubrication approximation so that governing equations are reduced to a 2D width-averaged system formed by the quasi-linear elliptic equation for pressure and transport equations for volume concentrations of fluids. The numerical solution is obtained using a finite-difference method. The pressure equation is solved using an iterative algorithm and the Multigrid method, while the transport equations are solved using a second-order TVD flux-limiting scheme with the superbee limiter. This numerical model is validated against three different sets of experiments: (i) gravitational slumping of fluids in a closed Hele-Shaw cell, (ii) viscous fingering of fluids with a high viscosity contrast due to the Saffman–Taylor (S–T) instability in a Hele-Shaw cell at microgravity conditions, (iii) displacement of Bingham fluids in a Hele-Shaw cell with the development of fingers due to the S–T instability. Good agreement is observed between simulations and laboratory data. The model is then used to investigate the joint effect of fingering and slumping. Numerical simulations show that the slumping rate of yield-stress fluid is significantly less pronounced than that of a Newtonian fluid with the same density and viscosity. If a low-viscosity Newtonian fluid is injected after a yield-stress one, the S–T instability at the interface leads to the development of fingers. As a result, fingers penetrating into a fluid with a finite yield stress locally decrease the pressure gradient and unyielded zones develop as a consequence.     

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.     

Purely elastic interfacial instabilities in superposed flow of polymeric fluids

Purely elastic interfacial stability of superposed plane Poiseuille flow of polymeric liquids has been investigated utilizing both asymptotic and numerical techniques. It is shown that these instabilities are caused by an unfavorable jump in the first normal stress difference across the fluid interface. To determine the significance of these instabilities in finite experimental geometries, a comparison between the maximum growth rates of purely elastic instabilities with instabilities driven primarily by a viscosity or a combined viscosity and elasticity difference is made. Based on this comparison, it is shown that purely elastic interfacial instabilities can play a major role in superposed flow of polymeric liquids in finite experimental geometries.     

An experimental investigation of initial oscillations in a radial Hele-Shaw cell

Hele-Shaw cell is a laboratory device consisting of two parallel plates of glass separated by a thin gap. In this cell, in the flow of two immiscible fluids, when a fluid of higher viscosity is displaced by a fluid of lower viscosity, the less viscous fluid is observed to form “fingers” into the more viscous one due to the unstable interface. The Saffman-Taylor or viscous finger instability has been examined and modeled for over forty years for the rectilinear Hele-Shaw cell and about half as long for the radial Hele-Shaw cell. In this paper, we study, in detail, the early development of viscous instabilities in a radial Hele-Shaw cell. This source flow configuration has been chosen so that the instability can be monitored precisely. The objective of this study is to examine the onset of fingering, i.e. initial number of fingers that form, and the evolution of interface instability. Our experiments suggest that there may be some order in this formation process and one can model this aspect by considering the unsteady velocity components and predicting temporal changes in wavenumber responsible for the initial number of fingers and may be later accounting for the fingertip oscillations and splitting. We injected a water-based fluid into an oil in a radial Hele-Shaw cell at constant flow rate and recorded the movement of the less viscous droplet as it evolved. The relative curvature changes on the expanding droplet boundary was plotted with the angular positions about the interface and subtracting out the average radius, resulting in a plot of the change in amplitude with respect to time for the interface configuration. Three unstable configured tests at kinematic viscosity contrast (v ^O) of 0.34, 0.68, and 0.94 were run at approximately the same flow rate (2π cm²/s). The droplet exhibited oscillatory movement for these unstable configuration. The amplitude and the rate of oscillations were measured from digitized data. The smaller the viscosity difference, the smaller was the amplitude growth rate and resulted in a longer time to form visible finger initiation. This work was supported by National Science Foundation, grant number EID-9017555. We also like to thank Dr. Len Schwartz, Professor of Mechanical Engineering at the University of Delware for his insight and helpful suggestions.     

Variability of a typical flow in a Hele-Shaw cell

The influence of different factors on a four-eddy convective flow, typical of a Hele-Shaw cell, with reconnection of the diagonal eddies is studied experimentally and theoretically. It is shown that this self-oscillating flow develops in the Hele-Shaw cell over a wide range of values of the super-criticality in very different situations and may be both regular and stochastic. This flow is observed in cavities of different dimensions, for homogeneous fluids and mixtures, for heat-conducting and thermally insulated wide sides of the cavity, and also in the presence of vibration and other external actions.     

Shear viscosity of suspensions of aligned non-Brownian fibres

We report on the steady-state shear viscosity of suspensions of fibres dispersed in Newtonian fluids, in a wide range of volume fractions throughout the dilute and semi-dilute regimes. We show that the apparent shear-thinning behaviour, which is sometimes observed in the semi-dilute regime at intermediate shear rates, is an experimental artefact due to the presence of transient clusters of entangled fibres in the suspensions. At high shear rates, the fibres are aligned and the suspensions exhibit Newtonian behaviour. In this regime, the viscosity is a function of volume fraction and fibre aspect ratio only. The data can be rescaled onto a universal curve using a variable that accounts for the average contribution of the particles to the bulk stress. All these results are discussed in relation to recent theories. Received: 19 January 1999 Accepted: 17 June 1999     

Large Deborah number flows around confined microfluidic cylinders

Viscoelastic flow around a confined cylinder at high Deborah numbers is studied using microfluidic channels. By varying fluid properties and flow rates, a systematic study of the roles of elasticity and inertia is accomplished. Two new elastic flow instabilities that occur at high Deborah numbers are identified. A downstream instability of disordered and temporally varying streamlines is observed at a Deborah number above 10. This instability is a precursor to an unsteady vortex that develops upstream of the cylinder at higher Deborah numbers. Both instabilities occur at moderate Reynolds numbers but are fundamentally elastic. The size and steadiness of the upstream vortex are primarily controlled by the Deborah and the elasticity number.     

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     

Flow instability and shear localization in a drilling mud

To have a better knowledge of problems occurring with drilling fluids in complex wells, we carried out a detailed rheological analysis of a typical drilling mud at low shear rates using both conventional rheometry and MRI velocimetry. We show the existence of a viscosity bifurcation effect: Below a critical stress value, the mud tends to completely stop flowing, whereas beyond this critical stress, it reaches an apparent shear rate larger than a finite (critical) value, and no stable flows can be obtained between this critical shear rate value and zero. These results are confirmed by MRI velocity profiles, which exhibit a slope break at the interface between the solid and the liquid phases inside the Couette geometry. Moreover, this viscosity bifurcation is a transient phenomenon, the progressive development of which can be observed by MRI. A further examination of MRI data shows that, in the transient regime, the shear rate does not vary monotonously in the rheometer gap and is particularly large along the outer (rough) cylinder, which might be at the origin of the development of a region of constant shear rate in the apparent flow curve.     

Shear-banding structure orientated in the vorticity direction observed for equimolar micellar solution

 The non-monotonic shear flow of a viscoelastic equimolar aqueous surfactant solution (cetylpyridinium chloride-sodium salicylate) is investigated rheologically and optically in a transparent strain-controlled Taylor Couette flow cell. As reported before, this particular wormlike micellar solution exhibits first a shear thinning and then a pronounced shear-thickening behavior. Once this shear-thickening regime is reached, a transient phase separation/shear banding of the solution into turbid and clear ring-like patterns orientated perpendicular to the vorticity axis, i.e., stacked like pancakes, is observed (Wheeler et al. 1998; Fischer 2000). The solution exhibit several unique features as no induction period of the shear induced phase, no structural build-up at the inner rotating cylinder, jumping pancake structure of clear and turbid ringlike phases, and oscillating shear stresses appear once the pancake structure is present. According to our analysis this flow phenomenon is not purely a mechanical or rheological driven hydrodynamic instability but one has to take into account structural changes of the oriented micellar aggregates (flow induced non-equilibrium phase transition) as proposed by several authors. Although this particular flow behavior and the underlying mixture of shear induced phases and mechanical instabilities is not fully understood yet, some classification characteristics based on a recent theoretical approach by Schmitt et al. (1995) and Porte et al. (1997) where a strong coupling between the flow instability (non-homogeneous flow profile due to the bands) and the structural changes causes the observed transient phenomena can be derived. In reference to the presented model the observed orientation of the rings is typical for complex fluids that undergo a spinodal phase separation coupled with a thermodynamic flow instability. In contrast to other shear banding phenomena, this one is observed in parallel plate, cone-plate, and Couette flow cell as well as under controlled stress and controlled rate conditions. Therefore, it adds an additional aspect to the present discussion on shear banding phenomena, i.e., the coupling of hydrodynamics and phase transition of rheological complex fluids. Received: 8 January 2001 Accepted: 15 May 2001     

Anomalous motion of viscous fingers in surfactant solutions in a Hele–Shaw cell

Viscous fingering in surfactant solutions in a rectangular Hele–Shaw cell was investigated. Test fluids were aqueous solutions of cetyltrimethylammonium bromide (CTAB) with sodium salicylate (NaSal) as a counter ion, and the ratio of mole concentration of CTAB and that of NaSal was 1–7.7. Two fluids that had a mole concentration different from that of CTAB were used. Air was injected into the cell and the growth of the interface between air and a CTAB/NaSal solution was observed. The fingertip grew similar to the finger growth in shear-thinning fluids at low pressure gradients. It took a cuspidate shape at the intermediate pressure gradient, and a sudden protrusion at a critical shear rate occurred. In high shear rate regions, the finger behaved as in a less shear-thinning fluid. These phenomena relate to rheological properties of the test fluids. Comparison with flow curves for CTAB/NaSal systems showed that the critical shear rate related to the shear rate at which a bending point appeared in the flow curve.     

气流作用下同轴带电射流的不稳定性研究

通过对气体驱动同轴电流动聚焦的实验模型进行简化,开展了电场力和惯性力共同作用下同轴带电射流的不稳定性理论研究.在流动为无黏、不可压缩、无旋的假设下,建立了三层流体带电射流物理模型并得到了扰动在时间域内发展演化的解析形式色散关系,利用正则模方法求解色散方程发现了流动的不稳定模态,进而分析了主要控制参数对不稳定模态的影响.结果表明,在参考状态下轴对称模态的最不稳定增长率最大,因此轴对称扰动控制整个流场.外层气流速度越高,气体惯性力越大,射流的界面越容易失稳.内外层液-液同轴射流之间的速度差越大,射流越不稳定.表面张力对射流不稳定性起到促进作用.轴向电场对射流不稳定性具有双重影响:当加载电场强度较小时,射流不稳定性被抑制;当施加电压大于某一临界值时,轴向电场会促进射流失稳.临界电压的大小与界面上自由电荷密度和射流表面扰动发展关系密切.这些结果与已有的实验现象吻合,能够对实验的过程控制提供理论指导.     

Time dependent flow in equimolar micellar solutions: transient behaviour of the shear stress and first normal stress difference in shear induced structures coupled with flow instabilities

Recently we studied time dependent structural changes that are coupled with flow instabilities (Fischer 1998; Wheeler 1998; Fischer 2000). Within a stability analysis, a classification scheme for the feedback circuit of coupled shear-induced structure and flow instabilities was derived by Schmitt et al. (1995) and applied to our samples. Here, inhomogeneous flow layers of different concentration and viscosity are generated by shear-induced diffusion (spinodal demixing) and, as consequence, one no longer observes a homogeneous solution but a type of shear banding that is seen here for the first time. In this paper we present the behaviour of the first normal stress difference observed in the critical shear-rate regime where transient shear-induced structure is coupled with flow instability. Similar to the oscillations of the shear stresses (strain-controlled rheometer) one observes oscillations in the first normal stress difference. This behaviour indicates that elastic structures are built up and destroyed while the shear-induced structures occur and that the induced phase is more elastic than the initial one. Oscillations of shear stress and first normal stress difference are in phase and indicate that both phenomena are caused by the same mechanism. Received: 30 June 1999/Accepted: 14 December 1999     

The axisymmetric contraction–expansion: the role of extensional rheology on vortex growth dynamics and the enhanced pressure drop

The flow of a polystyrene Boger fluid through axisymmetric contraction–expansions having various contraction ratios (2≤β≤8) and varying degrees of re-entrant corner curvatures are studied experimentally over a large range of Deborah numbers. The ideal elastic fluid is dilute, monodisperse and well characterized in both shear and transient uniaxial extension. A large enhanced pressure drop above that of a Newtonian fluid is observed independent of contraction ratio and re-entrant corner curvature. Streak images, laser Doppler velocimetry (LDV) and digital particle image velocimetry (DPIV) are used to investigate the flow kinematics upstream of the contraction plane. LDV is used to measure velocity fluctuation in the mean flow field and to characterize a global elastic flow instability which occurs at large Deborah numbers. For a contraction ratio of β=2, a steady elastic lip vortex is observed while for contraction ratios of 4≤β≤8, no lip vortex is observed and a corner vortex is seen. Rounding the re-entrant corner leads to shifts in the onset of the flow transitions at larger Deborah numbers, but does not qualitatively change the overall structure of the flow field. We describe a simple rescaling of the deformation rate which incorporates the effects of lip curvature and allows measurements of vortex size, enhanced pressure drop and critical Deborah number for the onset of elastic instability to be collapsed onto master curves. Transient extensional rheology measurements are utilized to explain the significant differences in vortex growth pathways (i.e. elastic corner vortex versus lip vortex growth) observed between the polystyrene Boger fluids used in this research and polyisobutylene and polyacrylamide Boger fluids used in previous contraction flow experiments. We show that the role of contraction ratio on vortex growth dynamics can be rationalized by considering the dimensionless ratio of the elastic normal stress difference in steady shear flow to those in transient uniaxial extension. It appears that the differences in this normal stress ratio for different fluids at a given Deborah number arise from variations in solvent quality or excluded volume effects.     

Time-dependent plane Poiseuille flow of a Johnson–Segalman fluid

We numerically solve the time-dependent planar Poiseuille flow of a Johnson–Segalman fluid with added Newtonian viscosity. We consider the case where the shear stress/shear rate curve exhibits a maximum and a minimum at steady state. Beyond a critical volumetric flow rate, there exist infinite piecewise smooth solutions, in addition to the standard smooth one for the velocity. The corresponding stress components are characterized by jump discontinuities, the number of which may be more than one. Beyond a second critical volumetric flow rate, no smooth solutions exist. In agreement with linear stability analysis, the numerical calculations show that the steady-state solutions are unstable only if a part of the velocity profile corresponds to the negative-slope regime of the standard steady-state shear stress/shear rate curve. The time-dependent solutions are always bounded and converge to different stable steady states, depending on the initial perturbation. The asymptotic steady-state velocity solution obtained in start-up flow is smooth for volumetric flow rates less than the second critical value and piecewise smooth with only one kink otherwise. No selection mechanism was observed either for the final shear stress at the wall or for the location of the kink. No periodic solutions have been found for values of the dimensionless solvent viscosity as low as 0.01.     

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.     

Twisted Elastic Rings and the Rediscoveries of Michell's Instability

Elastic rings become unstable when sufficiently twisted. This fundamental instability plays an important role in the modeling of DNA mechanics and in cable engineering. In 1962, Zajac computed the value of the critical twist for the instability. This critical value was rediscovered in 1979 by Benham and independently by Le Bret in elastic models for DNA; unstable rings have since become an important example of elastic instabilities in rods both for the development of new methods and in applications. The purpose of this note is to show that the problem had been completely solved by John Henry Michell in 1889 in a rather elegant manner and to reflect on its history and modern developments.