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.     

Injection and suction of an upper-convected maxwell fluid through a porous-walled tube

The steady-state, similarity solutions of the flow of an upper-convected Maxwell fluid through a tube with a porous wall are constructed by asymptotic and numerical analyses as functions of the direction of flow through the tube, the amount of elasticity in the fluid, as measured by the Deborah number De, and the degree of fluid slip along the tube wall. Fluid slip is assumed to be proportional to the local shear stress and is measured by a slip parameter β that ranges between no-slip (β = 1) and perfect slip (β = 0). The most interesting results are for fluid injection into the tube. For β = 1, the family of flows emanating from the Newtonian limit (De = 0) has a limit point where it turns back to lower values of De. These solutions become asymptotic to De = 0) and develop an O(De) boundary layer near the tube wall with singularly high stresses matched to homogeneous elongational flow in the core. This solution structure persists for all nonzero values of the slip parameter. For β ≠ 1, a family of exact solutions is found with extensional kinematics, but nonzero shear stress convected into the tube through the wall. These flows differ for low De from the Newtonian asymptote only by the absence of the boundary layer at the tube wall. Finite difference calculations evolve smoothly between the Newtonian-like and extensional solutions because of approximation error due to under-resolution of the boundary layer. The radial gradient of the axial normal stress of the extensional flow is infinite at the centerline of the tube for De > 1; this singularity causes failure of the finite difference approximations for these Deborah numbers unless the variables are rescaled to take the asymptotic behavior into account.     

Conjugate heat transfer in the unbounded flow of a viscoelastic fluid past a sphere

This work addresses the conjugate heat transfer of a simplified PTT fluid flowing past an unbounded sphere in the Stokes regime (Re = 0.01). The problem is numerically solved with the finite-volume method assuming axisymmetry, absence of natural convection and constant physical properties. The sphere generates heat at a constant and uniform rate, and the analysis is conducted for a range of Deborah (0 ≤ De ≤ 100), Prandtl (10⁰ ≤ Pr ≤ 10⁵) and Brinkman (0 ≤ Br ≤ 100) numbers, in the presence or absence of thermal contact resistance at the solid–fluid interface and for different conductivity ratios (0.1 ≤ κ ≤ 10). The drag coefficient shows a monotonic decrease with De, whereas the normalized stresses on the sphere surface and in the wake first increase and then decrease with De. A negative wake was observed for the two solvent viscosity ratios tested (β = 0.1 and 0.5), being more intense for the more elastic fluid. In the absence of viscous dissipation, the average Nusselt number starts to decrease with De after an initial increase. Heat transfer enhancement relative to an equivalent Newtonian fluid was observed for the whole range of conditions tested. The dimensionless temperature of the sphere decreases and becomes more homogeneous when its thermal conductivity increases in relation to the conductivity of the fluid, although small changes are observed in the Nusselt number. The thermal contact resistance at the interface increases the average temperature of the sphere, without affecting significantly the shape of the temperature profiles inside the sphere. When viscous dissipation is considered, significant changes are observed in the heat transfer process as Br increases. Overall, a simplified PTT fluid can moderately enhance heat transfer compared to a Newtonian fluid, but increasing De does not necessarily improve heat exchange.     

Dynamics of hairpin vortices generated by a mixing tab in a channel flow

To better understand mixing by hairpin vortices, time-series particle image velocimetry (PIV) was applied to the wake of a trapezoidal-shaped passive mixing tab mounted at the bottom of a square turbulent channel (Re _h =2,080 based on the tab height). Instantaneous velocity/vorticity fields were obtained in sequences of 10 Hz in the tab wake in the center plane (x–y) and in a plane (x–z) parallel to the wall. Periodically-shed hairpin vortices were clearly identified and seen to rise as they advected downstream. Experimental evidence shows that the vortex-induced ejection of the near-wall viscous fluid to the immediate upstream is important to the dynamics of hairpin vortices. It can increase the strength of the hairpin vortices in the near tab region and cause generation of secondary hairpin vortices further downstream when the hairpin heads are farther away from the wall. Measurements also reveal the existence of a type of new secondary vortice with the opposite-sign spanwise vorticity. The distribution of vortex loci in the x–y plane shows that the hairpin vortices and the reverse vortices are spatially segregated in distinct layers. Turbulence statistics, including mean velocity profiles, Reynolds stresses, and turbulent kinetic energy dissipation rate distributions, were obtained from the PIV data. These statistical quantities clearly reveal imprints of the identified vortex structures and provide insight into mixing effectiveness. Received: 24 February 2000/Accepted: 24 October 2000     

Three-dimensional flow of Newtonian and Boger fluids in square–square contractions

The flow of a Newtonian fluid and a Boger fluid through sudden square–square contractions was investigated experimentally aiming to characterize the flow and provide quantitative data for benchmarking in a complex three-dimensional flow. Visualizations of the flow patterns were undertaken using streak-line photography, detailed velocity field measurements were conducted using particle image velocimetry (PIV) and pressure drop measurements were performed in various geometries with different contraction ratios. For the Newtonian fluid, the experimental results are compared with numerical simulations performed using a finite volume method, and excellent agreement is found for the range of Reynolds number tested (Re₂ ≤ 23). For the viscoelastic case, recirculations are still present upstream of the contraction but we also observe other complex flow patterns that are dependent on contraction ratio (CR) and Deborah number (De₂) for the range of conditions studied: CR = 2.4, 4, 8, 12 and De₂ ≤ 150. For low contraction ratios strong divergent flow is observed upstream of the contraction, whereas for high contraction ratios there is no upstream divergent flow, except in the vicinity of the re-entrant corner where a localized atypical divergent flow is observed. For all contraction ratios studied, at sufficiently high Deborah numbers, strong elastic vortex enhancement upstream of the contraction is observed, which leads to the onset of a periodic complex flow at higher flow rates. The vortices observed under steady flow are not closed, and fluid elasticity was found to modify the flow direction within the recirculations as compared to that found for Newtonian fluids. The entry pressure drop, quantified using a Couette correction, was found to increase with the Deborah number for the higher contraction ratios.     

Flow field characteristics in the near field region of particle-laden coaxial jets

The effects of solid particles on the flow structure in the near field region of a coaxial water jet are investigated non-intrusively using molecular tagging velocimetry. Glass beads of 240 μm and specific gravity SG of 2.46 are used at three volume loadings of γ_v=0.03, 0.06, and 0.09% in the central water jet with a Reynolds number of 4.1×10⁴. Measurements are acquired for four annular to central jet velocity ratios in the range 0.11≤ U _o/U _i≤1.15 at downstream distances up to six inner jet diameters and the results are analyzed for the effects of solid particles on the characteristics of flow. It is found that the addition of particles does not affect the mean fluid velocity profile in this region. The results also indicate a small and moderate enhancement of axial turbulent velocity and radial gradients of velocity fluctuations, respectively, due to the presence of particles.     

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.     

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.     

Synchronized analysis of an unsteady laminar flow downstream of a circular cylinder centred between two parallel walls using PIV and mass transfer probes

Experiments are carried out in the wake of a cylinder of d _c  = 10 mm diameter placed symmetrically between two parallel walls with a blockage ratio r = 1/3 and a Reynolds number varying between 75 ≤ Re ≤ 277. Particle image velocimetry is exerted to obtain the instantaneous velocity components in the cylinder wake. A snapshot proper orthogonal decomposition (POD) is also applied to these PIV results in order to extract the dominant modes through the implementation of an inhomogeneous filtering of these different snapshots, apart from an interpolation to estimate the wall shear rate at the lower wall downstream the cylinder. Mass transfer circular probes are placed at the lower wall downstream this obstacle so as to further determine the time evolution of the wall shear rate, by bringing the inverse method to bear on the convective-diffusion equation. Comparisons between the two synchronized techniques demonstrate that electrochemical method can give more accurate information about the coherent structures present in the flow and about the interaction of the von Kármán vortices with the walls of the tunnel as well. The comparison between the two measurement techniques in the flow regions concerns the spatiotemporal evolutions of the wall shear rate obtained from PIV measurements and the wall shear rate using mass transfer probes. Discrepancy between the PIV measurements and the electrochemical ones near the wall, where the secondary vortices P ₁′ are generated at wall, are caused by a PIV bias and a limitations of the singular mass transfer probes.     

The motion of long bubbles through viscoelastic fluids in capillary tubes

The penetration of long gas bubble through a viscoelastic fluid in a capillary tube has been studied in order to investigate the influence of viscoelastic material properties on the hydrodynamic coating thickness and local flow kinematics. Experiments are conducted for three tailored ideal elastic (Boger) fluids, designed to exhibit similar steady shear properties but substantially different elastic material functions. This allows for the isolation of elastic and extensional material effects on the bubble penetration process. The shear and extensional rheology of the fluid is characterized using rotational and filament stretching rheometers (FSR). The fluids are designed such that the steady-state extensional viscosity measured by the FSR at a Deborah number (De) greater than 1 differs over three orders of magnitude (Trouton ratio = 10³–10⁶). The experiment set up to measure the hydrodynamic coating thickness is designed to provide accurate data over a wide range of capillary numbers (0.01 < Ca < 100). The results indicate that the coating thickness in this process increases with an increase in the extensionally thickening nature of the fluid. Experiments are also conducted using several different capillary tube diameters (0.1 < D < 1 cm), in order to compare responses at similar Ca but different flow De. Suitable scaling methods and nonlinear viscoelastic constitutive equations are explored to characterize the displacement process for polymeric fluids. Bubble tip shapes at different De are recorded using a CCD camera, and measured using an edge detection algorithm. The influence of the mixed flow field on the bubble tip shape is examined. Particle tracking velocimetry experiments are conducted to compare the influence of viscoelastic properties on the velocity field in the vicinity of the bubble tip. Local shear and extension rates are calculated in the vicinity of the bubble tip from the velocity data. The results provide quantitative information on the influence of elastic and extensional properties on the bubble penetration process in gas-assisted injection molding. The bubble shape and velocity field information provides a basis for evaluating the performance of constitutive equations in mixed flow. Received: 19 January 1999 Accepted: 30 June 1999     

Interaction of circular cylinder wake with a short asymmetrically located downstream plate

This study reveals the interaction patterns of separated shear layers from a circular cylinder with a short downstream plate and their reflection on the frequency and the formation length of the vortices from the cylinder as a function of plate location relative to the cylinder. The effect of horizontal (G/D) and vertical (Z/D) distances between the cylinder and the plate on the near wake is studied via Digital Particle Image Velocimetry (DPIV) in a water channel for Reynolds numbers of 200, 400 and 750, based on the cylinder diameter D. It is shown that the interaction of wake with the plate of length D can be categorized depending on the horizontal and the vertical distances between the cylinder and the plate. For the vertical distance range of Z/D ≤ 0.7, there is a critical horizontal spacing before which the shear layers from the cylinder are inhibited to form vortices in front of the plate. Resulting elongated recirculation region between the plate and the cylinder suggests modification of the absolutely unstable near wake of free circular cylinder in favor of convective instability. Z/D = 0.9 provides a passage from Z/D ≤ 0.7 to ≥1.1 and is associated with a dominant effect on the near-wake characteristics of interaction of shear layers from the cylinder with those from the downstream plate. For Z/D ≥ 1.1, there is again, yet a smaller critical horizontal spacing after which vortices interact with decreased downstream plate interference. In this vertical separation distance range, a gap flow between the plate and the cylinder plays a determining role on the formation length and St number of vortices for small horizontal spacing values.     

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.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

Purely elastic flow asymmetries in flow-focusing devices

The flow of a viscoelastic fluid through a microfluidic flow-focusing device is investigated numerically with a finite-volume code using the upper-convected Maxwell (UCM) and Phan-Thien–Tanner (PTT) models. The conceived device is shaped much like a conventional planar “cross-slot” except for comprising three inlets and one exit arm. Strong viscoelastic effects are observed as a consequence of the high deformation rates. In fact, purely elastic instabilities that are entirely absent in the corresponding Newtonian fluid flow are seen to occur as the Deborah number (De) is increased above a critical threshold. From two-dimensional numerical simulations we are able to distinguish two types of instability, one in which the flow becomes asymmetric but remains steady, and a subsequent instability at higher De in which the flow becomes unsteady, oscillating in time. For the UCM model, the effects of the geometric parameters of the device (e.g. the relative width of the entrance branches, WR) and of the ratio of inlet average velocities (VR) on the onset of asymmetry are systematically examined. We observe that for high velocity ratios, the critical Deborah number is independent of VR (e.g. De_c ≈ 0.33 for WR = 1), but depends non-monotonically on the relative width of the entrance branches. Using the PTT model we are able to demonstrate that the extensional viscosity and the corresponding very large stresses are decisive for the onset of the steady-flow asymmetry.     

Asymptotic analysis of the finite cone-and-plate flow of a non-newtonian fluid

Consider the shearing flow of a viscoelastic fluid trapped by surface tension between a cone and a plate. An asymptotic analysis of this problem in the limit of small gap angle has been done. This limit is realized in many practical situations. It is assumed that the Deborah number De, the Reynolds number Re, and the retardation parameter β are all order unity and that the shape of the free surface is very nearly spherical. Closed form analytic expressions are obtained for the leading terms of the primary and weak secondary motion of the fluid as well as the meniscus shape. It is found that the velocity field is bounded and continuous if and only if . There is a family of curves in the De-β plane on which the velocity field has a removable singularity at the origin. The secondary flow is made up of either one or two toroidal vortices. The meniscus has a bulge near the rotating cone and a trough near the stationary plate.     

Velocity fields and streamline patterns of miscible displacements in cylindrical tubes

Displacements of a viscous fluid by a miscible fluid of a lesser viscosity and density in cylindrical tubes were investigated experimentally. Details of velocity and Stokes streamline fields in vertical tubes were measured using a DPIV (digital particle image velocimetry) technique. In a reference frame moving with the fingertip, the streamline patterns around the fingertip obtained from the present measurements confirm the hypothesis of Taylor (1961) for the external patterns, and that of Petitjeans and Maxworthy (1996) for the internal patterns. As discussed in these papers, the dependent variable, m, a measure of the volume of viscous fluid left on the tube wall after the passage of the displacing finger, is a parameter that determines the flow pattern. When m>0.5 there is one stagnation point at the tip of the finger; when m<0.5 there are two stagnation points on the centerline, one at the tip and the other inside the fingertip, and a stagnation ring on the finger surface with a toroidal recirculation in the fingertip between the two stagnation points. The finger profile is obtained from the zero streamline of the streamline pattern.An erratum to this article can be found at     

Momentum and heat transport in a finite-length cylinder wake

This paper reports an experimental study of turbulent momentum and heat transport in the wake of a wall-mounted finite-length square cylinder, with its length-to-width ratio L/d = 3–7. The cylinder was slightly heated so that heat produced could be considered as a passive scalar. A moveable three-wire probe (a combination of an X-wire and a cold wire) was used to measure velocity and temperature fluctuations at a Reynolds number of 7,300 based on d and the free-stream velocity. Measurements were performed at 10 and 20d downstream of the cylinder at various spanwise locations. Results indicate that L/d has a pronounced effect on Reynolds stresses, temperature variance and heat fluxes. The downwash flow from the free end of the cylinder acts to suppress spanwise vortices and, along with the upwash flow from the cylinder base, makes the finite-length cylinder wake highly three-dimensional. Reynolds stresses, especially the lateral normal stress, are significantly reduced as a result of suppressed spanwise vortices at a small L/d. The downwash flow acts to separate the two rows of spanwise vortices further apart from the wake centerline, resulting in a twin-peak distribution in temperature variance. While the downwash flow entrains high-speed fluid into the wake, responsible for a small deficit in the time-averaged streamwise velocity near the free end, it does not alter appreciably the distribution of time-averaged temperature. It has been found that the downwash flow gives rise to a counter-gradient transport of momentum about the central region of the wake near the free end of the cylinder, though such a counter-gradient transport does not occur for heat transport.     

Digital particle image velocimetry/thermometry and application to the wake of a heated circular cylinder

 Digital particle image velocimetry/thermometry (DPIV/T) is a technique whereby the velocity and temperature fields are obtained using thermochromic liquid crystal (TLC) seeding particles in water. In this paper, the uncertainty levels associated with temperature and velocity measurements using DPIV/T are studied. The study shows that large uncertainties are encountered when the temperature is measured from individual TLC particles. Therefore, an averaging procedure is presented which can reduce the temperature uncertainties. The uncertainty is reduced by computing the average temperature of the particles within the common specified sampling window used for standard DPIV. Using this procedure, the velocity and temperature distributions of an unsteady wake behind a heated circular cylinder are measured experimentally at Re=610. The instantaneous DPIV/T measurements are shown to be useful for computing statistical flow quantities, such as mean and velocity-temperature correlations. Received: 3 January 2000/Accepted: 26 June 2000