Turbulent expansion flow of low molecular weight shear-thinning solutions

A Laser-Doppler anemometer and a pressure transducer were used to carry out detailed measurements of the mean and root mean square of the velocity and wall-pressure in an axisymmetric sudden expansion flow, with 0.4 and 0.5% by weight shear-thinning aqueous solutions of a low molecular weight polymer (6,000), after appropriate rheological characterisation. In spite of their very low molecular weight, these solutions still exhibited elongational elastic effects through drag reduction of up to 35% relative to Newtonian turbulent pipe flow, as shown by Pereira and Pinho (1994). The results showed small variations of the recirculation bubble length with polymer concentration and Reynolds number and reductions of the normal Reynolds stresses of up to 30%, especially in the tangential and radial directions. The reduction in normal Reynolds stresses within the shear layer is an elongational elastic effect, but this elasticity needs to be considerably more intense, such as with high molecular weight polymers, in order to strongly affect the mean flow characteristics. The observed mean flow patterns with these low molecular weight polymer solutions were indeed similar to those exhibited by Newtonian and inelastic fluids.     

Performance modeling and analysis of blood flow in elastic arteries

Nomenclature　　τ　　wallshearstressγshearrateτy yieldstressηc Cassonviscosityktheconsistencyindexnnon_Newtonianindexτp shearstressofthepthelementωangularvelocityRvessel’sradiusCwavespeedM　　magneticparameter (Hartmannnumber)u,w　velocitycomponentinther_andz_directions,respectivelyP　　pressureα　　unsteadinessparameter k , R meanparametersTp relaxationtimeofthepthelementρ densityIntroductionTheimportancetoatherogenesisofarterialflowphenomenasuchasflowseparation ,recirculationands…     

Extensional flow behavior of polymer solutions and particle suspensions in a spinning motion

Mechanical spinning of fluid filaments was used to generate an extensional flow, in which rheological measurements were obtained for a Newtonian fluid, two aqueous polymer solutions, and two fluid suspensions of rod-shaped particles. The tensile stress was determined by measuring the tensile force of the fluid filament while the kinematics were determined from photographic measurement of the filament profile and the assumption of a flat velocity profile. The measured tensile stresses for the Newtonian fluid matched predicted stresses, thereby confirming the validity of the experimental technique.The spinning behavior of each polymer solution could be correlated as stress versus extension rate. The apparent “spinning viscosity” increased with increasing rate of extension, in contrast to shear-thinning behavior in viscometric flow. For the fluid suspensions, the presence of rod-shaped particles increased the apparent viscosity far more in extensional flow than in shear. Tensile stresses calculated from a theoretical formula for suspensions proposed by Batchelor agreed rather well with experiment. Some general criteria for the interpretation of the spinning experiment are proposed, and some microrheological implications of the present findings are discussed.     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

Extensional behavior influence on viscoelastic turbulent channel flow

In this work we use in the simulation of a viscoelastic turbulent channel flow a modification of the finitely extensible of non-linear elastic dumbbells with the Peterlin approximation (FENE-P) constitutive model for dilute polymer solutions, applicable to high extensional deformations. The new feature introduced by this modification is that the free energy of the polymer (since it is assumed to be entirely entropically driven) remains always bounded (FENE-PB). The characteristics of the model under steady shear flow, pure elongational flow and transient extensional behavior are presented. It is found that the FENE-PB model is more shear thinning than FENE-P. Most importantly, it also shows a higher extensional viscosity than the FENE-P model. Although the steady-state Trouton ratio asymptotically reaches at high extensional rates the same limit as the FENE-P model, the transition from the Newtonian value is sharper and faster. We use the FENE-PB model in direct numerical simulations (DNS) of viscoelastic turbulent channel flow using spectral approximations. The results for various statistics of the flow and the polymer conformation, when compared against those obtained with the original FENE-P model and the same rheological parameters, show an enhanced polymer-induced drag reduction effect and enhanced deformation of the polymer molecules. This indicates that it is not only the asymptotic but also details from the extensional rheological behavior that matter in quantitatively specifying turbulent viscoelastic flow behavior.     

Asymmetry in transitional pipe flow of drag-reducing polymer solutions

New experimental data are presented and discussed for fully developed pipe flow of shear-thinning, viscoelastic polymer solutions in the transitional regime between laminar and turbulent flow. The data confirm that such transitional flows exhibit significant departures from axisymmetry in contrast to the fully developed pipe flow of Newtonian fluids or both laminar and turbulent flows of such drag-reducing liquids. The azimuthal structure of the asymmetry is investigated together with its axial development and also the velocity fluctuation levels. These data do not lead to an explanation for the asymmetry but do suggest that the influence of the flow geometry both upstream and downstream can be ruled out.     

Squeeze film flow of ideal elastic liquids

An exact solution is presented for the squeeze film flow of an Oldroyd B. fluid. The solution demonstrates that the flow kinematics is similar to the Newtonian (or Maxwellian) one. Theoretical predictions for constant velocity squeezing are compared to experimental observation for well characterized non-shear thinning elastic fluids. It is shown both theoretically and experimentally that the effect of elasticity in a constant velocity squeeze film flow is to always reduce the load relative to the inelastic (Newtonian) prediction and that this load reduction falls between the upper and lower asymptote prediction by the exact solution for the Oldroyd B fluid. The upper load asymptote is given by the Stefan solution for the viscosity of the polymer solution and the lower asymptote is given by the Stefan solution for the viscosity of the solvent. Experimental observations agree with the theoretical prediction for the Oldroyd B fluid at low shear rates where it is shown that the steady and dynamic flow properties of the test fluids used in the experimental program are well represented by the Oldroyd B constitutive equation. With the exception of the work of Lee et al. [6] for constant load squeezing of a Maxwell fluid, this work represents one of the few cases where experimental observation of large effects due to elasticity are indeed predicted with a constitutive equation which actually describes the steady and dynamic shear properties of the fluids used in the experimental program.     

A pseudo‐viscoelastic approach for transient polymeric flow exiting a channel

The influence of elasticity of a fluid exiting a channel is examined on transient coating downstream. A hybrid spectral/boundary element approach is proposed to solve the problem. The flow inside the channel is assumed to be fully developed. A viscoelastic instability of one‐dimensional plane Couette flow is first determined for a large class of Oldroyd fluids with added viscosity, which typically represent polymer solutions composed of a Newtonian solvent and a polymeric solute. The Johnson–Segalman equation is used as the constitutive model. The velocity profile inside the channel is taken as the exit profile for the emerging free‐surface flow. The flow is assumed to be Newtonian as it emerges from the channel. An estimate of the magnitude of the rate‐of‐strain tensor components in the free‐surface region reveals that they are generally smaller than the shear rate inside the channel. The evolution of the flow front is simulated using the boundary element method. For the channel flow, the problem is reduced to a nonlinear dynamical system using the Galerkin projection method. Stability analysis indicates that the channel velocity may be linear or non‐linear depending on the range of the Weissenberg number. The evolution of the coating flow at the exit is examined for steady as well as transient (monotonic and oscillatory) channel flow. It is found that adverse flow can exist as a result of fluid elasticity, which can hinder the process of blade coating. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Anomalous predictions of pressure drop and heat transfer in ducts of arbitrary cross-section with modified power law fluids

 The apparent viscosities of purely viscous non-Newtonian fluids are shear rate dependent. At low shear rates, many of such fluids exhibit Newtonian behaviour while at higher shear rates non-Newtonian, power law characteristics exist. Between these two ranges, the fluid's viscous properties are neither Newtonian or power law. Utilizing an apparent viscosity constitutive equation called the “Modified Power Law” which accounts for the above behavior, solutions have been obtained for forced convection flows. A shear rate similarity parameter is identified which specifies both the shear rate range for a given fluid and set of operating conditions and the appropriate solution for that range. The results of numerical solutions for the friction factor–Reynolds number product and for the Nusselt number as a function of a dimensionless shear rate parameter have been presented for forced fully developed laminer duct flows of different cross-sections with modified power law fluids. Experimental data is also presented showing the suitability of the “Modified Power Law” constitutive equation to represent the apparent viscosity of various polymer solutions. Received on 21 August 2000     

Viscosity enhancement in the flow of hydrolysed poly(acrylamide) saline solutions around spheres: implications for enhanced oil recovery

We have studied dilute aqueous solutions of hydrolysed poly(acrylamide), in various ionic environments, in flow around single spheres and around two spheres aligned on the axis of flow. The spheres are held on flexible cantilevers, while the polymer solutions, or solvent, are drawn past at controlled flow rates. We estimate the specific viscosities of the various solutions as a function of the strain rate over strain rates encompassing both the shear thinning and extension thickening regimes. For flow of solutions without added salts around a single sphere, we observe shear thinning followed by a significant increase in the non-Newtonian viscosity with increasing strain rate. The shear thinning reduces the maximal extensional viscosities of the solutions, which has important implications regarding the effectiveness of hydrolysed poly(acrylamide) in oil field applications. For flow of polymer solutions around two axially aligned spheres, we observe a significant reduction in the non-Newtonian forces experienced by the downstream sphere in comparison to the upstream sphere. We consider that this is salient to the understanding of non-Newtonian viscosification in porous media flow.     

A procedure to determine the viscosity function from experimental data of capillary flow

The accurate calculation of the viscosity η as function of the shear rate &γdot; from capillary viscometry is still a matter of debate in the literature. In fact, this problem involves the inversion of an integral equation, which leads to multiple solutions due to the unavoidable noise present in the experimental data. The purpose of this work is to develop an efficient procedure to determine the viscosity function from experimental data of capillary flow without presenting the difficulties inherent in other methods discussed previously in the literature. The system identification procedure is used here to estimate the parameters of a viscosity model, which is appropriately selected for the fluid under study through preliminary calculations involving the apparent shear rate – shear stress data. Once the model is chosen by satisfying criteria for the fit goodness and its parameters are evaluated, a smooth and continuous function η(γdot;) is obtained in the range of experimental shear rates. The procedure proposed is also applicable to fluids in shear flow that present two Newtonian plateaus, as it is typically found in macromolecular dilute solutions. The mean value theorem of continuous functions is used to reduce significantly the computational time. Received: 15 November 1999 Accepted: 7 November 2000     

A non‐linear dynamical system approach to finite amplitude Taylor‐Vortex flow of shear‐thinning fluids

The effect of shear thinning on the stability of the Taylor–Couette flow is explored for a Carreau–Bird fluid in the narrow‐gap limit. The Galerkin projection method is used to derive a low‐order dynamical system from the conservation of mass and momentum equations. In comparison with the Newtonian system, the present equations include additional non‐linear coupling in the velocity components through the viscosity. It is found that the critical Taylor number, corresponding to the loss of stability of the circular Couette flow, becomes lower as the shear‐thinning effect increases. That is, shear thinning tends to precipitate the onset of Taylor vortex flow, which coincides with the onset of a supercritical bifurcation. Comparison with existing measurements of the effect of shear thinning on the critical Taylor and wave numbers show good agreement. The Taylor vortex cellular structure loses its stability in turn, as the Taylor number reaches a critical value. At this point, an inverse Hopf bifurcation emerges. In contrast to Newtonian flow, the bifurcation diagrams exhibit a turning point that sharpens with shear‐thinning effect. Copyright © 2004 John Wiley & Sons, Ltd.     

Exponential shear flow of branched polyethylenes in rotational parallel-plate geometry

Exponential shear flow, as a strong flow with the potential to generate a high degree of molecular stretching, has attracted considerable interest in recent years. So far, exponential shear flow has been realized by either sliding-plate or cone-and-plate (CP) geometry. Both geometries guarantee homogeneous shear flow. Here, we present experimental data on exponential shear flow of several long-chain branched polyethylene melts with different degrees of strain hardening obtained by using parallel-plate (PP) geometry in a rotational rheometer. This type of geometry, which is standard in linear-viscoelastic characterization of polymer materials, produces inhomogeneous shear flow. A comparison of exponential shear flow data obtained by PP and CP geometry is made. Additionally, the experimental data are compared to predictions of the rubber-like liquid (RLL) and the molecular stress function (MSF) theories. For this purpose, the relaxation spectra of the polymer melts considered were obtained by standard linear-viscoelastic characterization. In addition, two irrotational parameters and one rotational parameter are required by the MSF theory. While the irrotational parameters were obtained from fitting to elongational viscosity data, the value of the rotational parameter was used as given in the literature. It can be concluded that viable experimental data in exponential shear flow can be obtained by PP geometry. For finite linear-viscoelasticity (RLL theory), predictions of reduced shear stress for CP and PP geometry coincide, but nonlinear material behavior (as modeled by the MSF theory) leads to small differences between both geometries. Furthermore, it is shown that the MSF predictions are in excellent agreement with the experimental data in exponential shear flow and that this type of flow leads to much less chain stretching than elongational flow.Dedicated to the memory of Prof. Arthur S. Lodge (1922–2005).     

Evaluation of flush mounted hot-film sensors for skin friction reduction measurements in viscoelastic polymer solutions

The purpose of this investigation was to evaluate the performance of flush mounted hot-film sensors for mean wall shear stress measurement in turbulent flows of dilute drag reducing polymer solution. A series of pipe flow expriments were conducted over a range of Reynolds numbers and polymer solution concentrations to compare the level of skin friction drag reduction measured by hot-film sensors with values calculated from pipe pressure drop. It is shown that water calibrated hot-film sensors consistently underestimate the wall shear stress suggesting that Reynolds analogy is not valid in dilute polymer solutions. The Newtonian form of the relationship between the wall shear stress and the heat transfer remains valid in dilute polymer solutions. However, multiplicative and additive factors in the relationship are shown to increase linearly with the logarithm of the polymer concentration.     

Oscillatory squeezing flow of a biological material

Large-amplitude oscillatory squeezing flow data are reported for a complex biological material, which is highly shear-thinning in oscillatory shear flow. This soft tissue has a linear viscoelastic limit at a strain of approximately 0.2%. The oscillatory squeezing flow data at large strain are analyzed using two constitutive models: a bi-viscosity Newtonian model, and a non-linear Maxwell model. It is found that although both models may have the same response in shape, the later matches with our non-linear experimental data better. It is also concluded that the non-linear response of the material in large amplitude oscillatory flow is mainly due to the shear thinning of the material. Received: 9 February 2000/Accepted: 22 February 2000     

Unsteady flow of viscoelastic fluid between two cylinders using fractional Maxwell model

The unsteady flow of an incompressible fractional Maxwell fluid between two infinite coaxial cylinders is studied by means of integral transforms.The motion of the fluid is due to the inner cylinder that applies a time dependent torsional shear to the fluid.The exact solutions for velocity and shear stress are presented in series form in terms of some generalized functions.They can easily be particularized to give similar solutions for Maxwell and Newtonian fluids.Finally,the influence of pertinent parameters on the fluid motion,as well as a comparison between models,is highlighted by graphical illustrations.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

This is the second part of a study examining the mechanical properties and capillary flow of fiber suspensions in Newtonian fluids and in polymer solutions. In part I results for the viscous and elastic properties of the fiber suspensions were presented and it was shown that the fiber suspensions exhibited normal stresses in Newtonian as well as in viscoelastic suspending media. It was thus expected that circulating secondary flows would occur near the entrance to a capillary. Four types of fillers (glass, carbon, nylon and vinylon fibers) suspended in glycerin, HEC solutions and Separan solutions were investigated. The entrance flow patterns were visualized and the pressure fluctuations measured. The visualization enabled the eddies occurring in the fiber suspensions in Newtonian fluids to be analysed and classified into two tpyes. The results from the flow visualization were correlated with the pressure fluctuations. Empirical equations for the tube length correction factor due to the elasticity were obtained.     

Experimental investigations and numerical simulations for an open channel flow of a weak elastic polymer solution around a T-profile

The present paper is concerned with experimental and numerical investigations of planar complex flows of weak elastic polymer solutions (whose concentration are below the critical overlap concentration), characterised by small relaxation times (<0.1 s) and almost constant shear viscosities for small and medium shear rates. The main aim of the study is to detect to what extent a very small amount of elasticity present in a viscous fluid can influence its behaviour in complex flows, without introducing major modifications of classical rheological tests. The samples are polymer solutions of low PIB molecular weight dissolved in highly viscous Newtonian mineral oil. The analysed motion is steady, and takes place in an open channel around a T profile. Maximum values of the characteristic parameters for the experiments, the Reynolds and Weissenberg numbers, were 45 and 0.1, respectively. The experiments show a decrease of the wake length downstream the profile for weak elastic solutions in comparison to the Newtonian solvent. Actually, the same wake length as in the Newtonian case was obtained for tested polymer solutions, but at higher Re numbers. Numerical simulations using the Giesekus model predict the same behaviour and are consistent with experiments from both qualitative and quantitative point of views. The results of research conclude that, even in small amounts, the presence of elasticity in pure viscous liquids induces quantitative changes from Newtonian flow in complex dominant elongational flows, at elongational rates for which the sudden thickening of extensional viscosity is remarkable. The study is important, since it should enable better understanding and modelling of viscoelastic flows that involve dilute polymer solutions, or fluids with similar rheology; biofluid mechanics being one area of application of this research. Corroboration of experimental flow visualization with numerical simulation is currently a feasible method used to characterise weak elastic polymer solutions, since classical rheological techniques generally fail to obtain realistic values of relaxation time for these particular viscoelastic fluids. Corneliu Balan dedicates this paper to the anniversary of one hundred years from the birth of Academician Dumitru Dumitrescu (1904–1983), charismatic personality of the Romanian school of fluid mechanics.

---