An experimental investigation of negative wakes behind spheres settling in a shear-thinning viscoelastic fluid

We present detailed experimental results examining “negative wakes” behind spheres settling along the centerline of a tube containing a viscoelastic aqueous polyacrylamide solution. Negative wakes are found for all Deborah numbers (2.43≤De(˙γ)≤8.75) and sphere-to-tube aspect ratios (0.060≤a/R≤0.396) examined. The wake structures are investigated using laser-Doppler velocimetry (LDV) to examine the centerline fluid velocity around the sphere and digital particle image velocimetry (DPIV) for full-field velocity profiles. For a fixed aspect ratio, the magnitude of the most negative velocity, U _min , in the wake is seen to increase with increasing De. Additionally, as the Deborah number becomes larger, the location of this minimum velocity shifts farther downstream. When normalized with the sphere radius and the steady state velocity of the sphere, the axial velocity profiles become self-similar to the point of the minimum velocity. Beyond this point, the wake structure varies weakly with aspect ratio and De, and it extends more than 20 radii downstream. Inertial effects at high Reynolds numbers are observed to shift the entire negative wake farther downstream. Using DPIV to investigate the transient kinematic response of the fluid to the initial acceleration of the sphere from rest, it is seen that the wake develops from the nonlinear fluid response at large strains. Measurements of the transient uniaxial extensional viscosity of this weakly strain-hardening fluid using a filament stretching rheometer show that the existence of a negative wake is consistent with theoretical arguments based on the opposing roles of extensional stresses and shearing stresses in the wake of the sphere. Received: 10 November 1997 Accepted: 1 May 1998     

Uniform flow of viscoelastic fluids past a confined falling cylinder

Uniform steady flow of viscoelastic fluids past a cylinder placed between two moving parallel plates is investigated numerically with a finite-volume method. This configuration is equivalent to the steady settling of a cylinder in a viscoelastic fluid, and here, a 50% blockage ratio is considered. Five constitutive models are employed (UCM, Oldroyd-B, FENE-CR, PTT and Giesekus) to assess the effect of rheological properties on the flow kinematics and wake patterns. Simulations were carried out under creeping flow conditions, using very fine meshes, especially in the wake of the cylinder where large normal stresses are observed at high Deborah numbers. Some of the results are compared with numerical data from the literature, mainly in terms of a drag coefficient, and significant discrepancies are found, especially for the constant-viscosity constitutive models. Accurate solutions could be obtained up to maximum Deborah numbers clearly in excess of those reported in the literature, especially with the PTT and FENE-CR models. The existence or not of a negative wake is identified for each set of model parameters.     

Structure of stratified flow around a cylinder at low internal froude number

The flow pattern around a horizontal cylinder towed at constant velocity in a continuously stratified fluid is visualized by the shadow method. The velocities in the leading flow disturbance, i. e., in the flow-blocking region ahead of the cylinder, are presented. In the body wake, a new class of small-size structures in the density gradient field is revealed against the background of a smooth velocity profile. The evolution of the flow pattern with variation of the parameters of body motion is studied. Institute of Mechanics Problems, Moscow 117526. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 80–88, January–February, 1999.     

Flow of wormlike micelle solutions past a confined circular cylinder

In this paper, we present the results of an investigation into the flow of a series of viscoelastic wormlike micelle solutions past a confined circular cylinder. Although this benchmark flow has been studied in great detail for polymer solutions, this paper reports the first experiments to use a viscoelastic wormlike micelle solution as the test fluid. The flow kinematics, stability and pressure drop were examined for two different wormlike micelle solutions over a wide range of Deborah numbers and cylinder to channel aspect ratios. A combination of particle image velocimetry and pressure drop measurements were used to characterize the flow kinematics, while flow-induced birefringence measurements were used to measure the micelle deformation and alignment in the flow. The pressure drop was found to decrease initially due to the shear thinning of the test fluid before increasing at higher flow rates as elastic effects begin to dominate the flow. Above a critical Deborah number, an elastic instability was observed for just one of the test fluids studied, the other remained stable for all Deborah number tested. Flow-induced birefringence and velocimetry measurements showed that observed instability originates in the extensional flow in the wake of the cylinder and appears not as periodic counter-rotating vortices as has been observed in the flow of polymer solutions past circular cylinders, but as a chaotic rupture event in the wake of the cylinder that propagates axially along the cylinder. Reducing the cylinder to channel aspect ratio and the degree of shearing introduced by the channel walls had a weak impact on the stability of the flow. These measurements, when taken in conjunction with previous work on flow of wormlike micelle solutions through a periodic array of cylinders, definitively show that the instability can be attributed to a breakdown of the wormlike micelle solutions in the extensional flow in the wake of the cylinder.     

Flow and blockage of highly concentrated granular suspensions in non-Newtonian fluid

Studying the flow of highly concentrated granular suspensions represents a great challenge since they are characterized by a rather complex rheological behavior. In addition, macroscopic heterogeneities may be induced by the flow during rheological measurements due to the eventual relative motion between the liquid and the granular phases that may take place under certain conditions. Solid–liquid separation may ultimately lead to flow blockage. In the present investigation we consider experimentally the influence of the rheological properties of the suspending fluid on the transition between the flow and blockage of a concentrated suspension in a squeeze set-up geometry. The suspending fluid consists of an aqueous Xanthan solution for which rheological properties can be tuned by changing the polymer concentration. For each polymer concentration, it is shown that there exist flow parameters (squeeze velocity and gap thickness) for which one has a transition between homogeneous flow of the suspension and its blockage. Blockage diagrams, delimiting flowability and blockage zones, are then determined. Physical arguments are given to relate the evolution of the blockage diagrams to the flow parameters and rheological properties of suspending fluid.     

Influence of polymer additives in a flow on the hydrodynamic characteristics of two-dimensional profiles

The problem of separationless flow of homogeneous dilute polymer solutions over two-dimensional profiles is considered. The complete flow is divided by the outer edge of the boundary layer and the wake into two regions: a region of irrotational flow and a region of viscous flow — the boundary layer and wake. The characteristics of the two regions are matched at their boundary. The problem is solved by successive approximation with allowance for the mutual influence of the two regions on each other. The influence of the irrotational region on the viscous region is taken into account through the distribution of the pressure on the boundary of the wake and the boundary layer. The influence of the viscous part of the flow is taken into account by the introduction of an associated vortex whose intensity is equal to the integral of the vorticity in the complete viscous region, and also by the introduction of additional velocities on the boundary of the wake and the boundary layer. These deform the streamlines in the irrotational part of the flow and ensure that they match the flow pattern in the real fluid. The results of the calculations of the hydrodynamic characteristics of a Zhukovskii profile are compared with experimental data. The influence of the introduction into the flow of polymer additives on the distributed and total characteristics of the flow at a number of Reynolds numbers is analyzed for the example of the modified profile NACA66.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 35–41, October–December, 1981.     

Dynamics of horizontal viscoelastic fluid filaments

The tension force of a thinning high-molecular polymer solution filament is measured using the filament itself as a force sensor. The axial filament stresses and the effects of fluid flow from the filament into adjacent drops are estimated. It is shown that these effects are insignificant for polymer solutions in a low-viscosity solvent (water) but substantial for solutions in a high-viscosity fluid (glycerine). A modification of the standard rheological capillary filament method is proposed. This modification makes it possible to exclude any hypotheses concerning the stress distribution pattern within the filament. Periodic transverse oscillations of the filament axis are revealed and analyzed.     

Flow around individual Taylor bubbles rising in stagnant polyacrylamide (PAA) solutions

The flow around single Taylor bubbles rising in stagnant non-Newtonian solutions of polyacrylamide (PAA) polymer was studied using a technique employing simultaneous particle image velocimetry (PIV) and shadowgraphy. Solutions with different weight percentages of polymer, varying from 0.01 to 0.80 wt.%, were used to cover a wide range of flow regimes. The rheological fluid properties and pipe dimension yielded Reynolds numbers between 2 and 1160, and Deborah numbers up to 115. The shape of the bubbles rising in the different solutions was compared and quantified by fitting correlations. The flow around the nose of the bubbles was found to be similar for all conditions studied. Velocity profiles were measured and analysed in the liquid film around the bubbles. A comparison of bubble wake flow patterns was made. For the 0.10 and 0.20 wt.% PAA solutions, long wakes with a recirculation region were observed. Below the wakes, a flow of stretched liquid was found. Negative wakes were also observed for the more concentrated solutions.     

Nitrogen/non-Newtonian fluid two-phase upward flow in non-circular microchannels

This paper presents experimental investigations on nitrogen/non-Newtonian fluid two-phase flow in vertical noncircular microchannels, which have square or triangular cross-section with the hydraulic diameters being D_h = 2.5, 2.886 and 0.866 mm, respectively, by visualization method. Three non-Newtonian aqueous solutions with typical rheological properties, i.e., 0.4% carboxymethyl cellulose (CMC), 0.2% polyacrylamide (PAM) and 0.2% xanthan gum (XG) are chosen as the working fluids. The common flow patterns are identified as slug flow, churn flow and annular flow. The dispersed bubble flow is only found in the case with nitrogen/CMC solution two-phase flow in the largest channel. A new flow pattern of nitrogen/PAM solution two-phase flow, named chained bubble/slug flow, is observed in all the test channels. The flow regime maps are also developed and the results show that the rheological properties of the non-Newtonian fluid have remarkable influence on the flow pattern transitions. The geometrical factors of the microchannel such as the cross-section shape and hydraulic diameter of the channel can also affect the flow regime map. Finally, the results obtained in this work are compared with the available flow pattern transitions.     

Molecular orientation in non-Newtonian flow of dilute polymer solutions around spheres

A novel approach is presented to study the benchmark problem of flow around spheres in model dilute solutions of monodisperse samples of atactic polystyrene in di-octyl phthalate. Spheres are held stationary on flexible cantilevers of known spring-constant, k, while the polymer solutions are pumped past at controlled flow rates, allowing access to a wide range of Deborah number. In this way the non-Newtonian forces experienced by the spheres can be measured as a function of Deborah number, while detailed observations and measurements of birefringence are made, enabling assessment of macromolecular strain and orientation. In addition the flow field around a sphere has been measured in an a-PS solution. Experiments have been performed on a single sphere and on two spheres axially aligned in the direction of flow. The extensional flow around the downstream stagnation point of the single sphere is found to play a pivotal role in the development of molecular strain and stress, resulting in flow modification and subsequent non-Newtonian behaviour. The flow birefringence in the wake is found to modify severely the flow around a second, downstream, sphere, affecting the non-Newtonian forces encountered by the second sphere. This provides an explanation for the time interval dependent terminal velocity often observed when two spheres follow the same path through viscoelastic liquids.     

Investigation of the wake of a hypersonic sphere in air by means of an open microwave resonator

The use of an open microwave resonator with plane-parallel mirrors for plasma diagnostics was first proposed in [1]. A resonator with spherical mirrors, which provides better spatial resolution in addition to high sensitivity, was used later [2, 3] to investigate the wake flow of models moving through air at hypersonic velocities. The presence of free electrons in the flow field is caused by ionization processes behind the bow shock and in the model boundary layer in this case. However, only the results of measurements of the density of electrons are presented in [2, 3], and no information is given on another important plasma parameter: the effective collision frequency of electrons with heavy particles. In the present study we use a microwave (8-mm range) resonator for an experimental study of the flow of gas in the wake of a polymer (Kaprolon) sphere traveling through air at hypersonic velocity. The flow is visualized by the schlieren technique. The electron densities and effective collision frequency of electrons with heavy particles are determined as a function of the distance behind the sphere.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 156–160, July–August, 1978.     

Negative wake in the uniform flow past a cylinder

The upstream/downstream streamline shift and the associated negative wake generation (streamwise velocity overshoot in the wake) in a viscoelastic flow past a cylinder are studied in this paper, for the Oldroyd-B, UCM, PTT, and FENE-CR fluids, using the Discrete Elastic Viscous Split Stress Vorticity (DEVSS-ω) scheme (Dou HS, Phan-Thien N (1999). The flow of an Oldroyd-B fluid past a cylinder in a channel: adaptive viscosity vorticity (DAVSS-ω) formulation. J Non-Newtonian Fluid Mech 87:47–73). The numerical algorithm is a parallelized unstructured Finite Volume Method (FVM), running under a distributed computing environment through the Parallel Virtual Machine (PVM) library. It is demonstrated that both the normal stress and its gradient are responsible for the negative wake generation and streamline shifting. Fluid extensional rheology plays an important role in the generation of the negative wake. The negative wake can occur in flows where the fluid extensional viscosity does not increase rapidly with strain rate. The formation of the negative wake does not depend on whether the streamlines undergo an upstream or a downstream shift. Shear-thinning viscosity weakens the velocity overshoot and while shear-thinning first normal stress coefficient enhances the velocity overshoot. Wall proximity is not necessary for the velocity overshoot; however, it enhances the strength of the negative wake. For the Oldroyd-B fluid, the ratio of the solvent viscosity to the zero-shear viscosity plays an important role in the streamline shift. In addition, mesh dependent behaviour of normal stresses along the centreline at high De in most cylinder/sphere simulations is due to the convection of normal stress from the cylinder to the wake, which results in the maximum of the normal stress being located off the centreline by a short distance at high De.     

Hydrodynamics of a sphere in a stratified fluid

The flow pattern around a sphere moving at constant velocity in a fluid with an exponential density distribution is investigated by optical methods. The thin density boundary layer forming the high-gradient envelope of the wake is distinguished as one of the elements of the structure. The symmetry properties of the flow are investigated. The limits of applicability of the traditional approximation of weak stratification in the problem of excitation of attached internal waves are noted.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 3–9, January–February, 1989.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

3D particle image velocimetry of the flow field around a sphere sedimenting near a wall: Part 1. Effects of Weissenberg number

The flow fields surrounding a sphere sedimenting through a liquid near a vertical wall are characterized using 3D stereoscopic particle-image velocimetry (PIV) experiments. Three different fluids, a Newtonian reference fluid, a constant (shear) viscosity Boger fluid, and a shear-thinning elastic fluid, are used to determine the effects of both elasticity and shear-thinning on the flow field. All three fluids have similar zero shear viscosities. The Weissenberg number is manipulated by varying the diameter and the composition of the ball. Significant differences are found for the different types of fluid, demonstrating both the influence of elasticity and shear-thinning on the velocity fields. In addition, the impact of the wall on the flow field is qualitatively different for each fluid. We find that the flow behind the sphere is strongly dependent on the fluid properties as well as the elasticity. Also, the presence of a negative wake is found for the shear-thinning fluid at high Weissenberg number (Wi > 1).     

Rheo-PIV analysis of the slip flow of a metallocene linear low-density polyethylene melt

The continuous extrusion of a metallocene linear low-density polyethylene through a transparent capillary die with and without slip was analyzed in this work by rheometrical measurements and particle image velocimetry (PIV). For this reason, a comparison was made between the rheological behaviors of the pure polymer and blended with a small amount of fluoropolymer polymer processing additive. Very good agreement was found between rheometrical and PIV measurements. The pure polymer exhibited stick-slip instabilities with nonhomogeneous slip at the die wall, whereas the blend showed stable flow. The slip velocity was measured directly from the velocity profiles and was negligible for the pure polymer before the stick-slip but increased monotonously as a function of the shear stress for the blend. The flow curves and the slip velocity as a function of the shear stress deviated from a power law and were well fitted by continuous “kink” functions. Comparison of PIV data with rheometrical ones permitted a direct proof of the basic assumption of the Mooney theory. Finally, the analysis of the velocity profiles showed that there is a maximum in the contribution of slip to the average fluid velocity, which is interpreted as the impossibility for the velocity profile to become plug like in the presence of shear thinning.     

Rheological studies and optimization of Herschel-Bulkley flow parameters of viscous karaya polymer suspensions using GA and PSO algorithms

Rheological characteristics of gum karaya suspensions which is proposed as a fracturing fluid were investigated with the main objective to determine the yield stress and other rheological parameters using various models. The flow hysteresis confirms the thixotropic behavior of fluid with increased structural breakdown at higher concentration and temperature. An empirical model developed for the studied samples accurately predicts the temperature and polymer concentration sensitivity of the apparent viscosity. The Herschel-Bulkley model showed the best fit to the experimental data; however, the yield stress obtained from some of the samples using nonlinear regression (NL) method reported physically insignificant, negative values. The proposed optimization technique, i.e., “Particle Swarm Optimization” offered the most realistic results with faster convergence over genetic algorithm making it a better choice. The oscillatory study provided more reliable yield stress values and revealed the thermogelation behavior of polymer suspensions making it suitable for high-temperature fracturing application.     

Incompressible Viscous Fluid Flows in a Thin Spherical Shell

Linearized stability of incompressible viscous fluid flows in a thin spherical shell is studied by using the two-dimensional Navier–Stokes equations on a sphere. The stationary flow on the sphere has two singularities (a sink and a source) at the North and South poles of the sphere. We prove analytically for the linearized Navier–Stokes equations that the stationary flow is asymptotically stable. When the spherical layer is truncated between two symmetrical rings, we study eigenvalues of the linearized equations numerically by using power series solutions and show that the stationary flow remains asymptotically stable for all Reynolds numbers.      

Investigation of the variability of the structure of a stratified wake behind a horizontal cylinder using optical and acoustic methods

The shadow flow pattern behind a horizontal cylinder uniformly towed in a stratified fluid with constant buoyancy frequency (in the imbedded vortex and turbulent wake regime) is recorded synchronously with acoustic echo sounding (basic frequency equal to 1 MHz) in a laboratory tank. Using computer processing, the illumination profiles in the schlieren pattern are constructed on scales comparable with the sounding acoustic ray width. Although the optical and acoustic profiles are not similar, nevertheless they enable the basic structural elements of the wake, including its high-gradient core to be identified, and their time variability traced. The features of the integral acoustic scattering characteristic, in particular, the volume scattering strength, allow this characteristic to be used, together with optical images of the flow pattern, to distinguish the flow regime identification criteria. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 5–17, May–June, 1998. The work was partly financed by the Russian Foundation for Basic Research (projects Nos. 96-05-64004 and 97-01-01013) and by INTAS Grant No. 93-1584.