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.     

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.     

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.     

Effects of viscosity ratio on deformation of a viscoelastic drop in a Newtonian matrix under steady shear

Deformation of an Oldroyd B drop in a Newtonian matrix under steady shear is simulated using a front tracking finite difference method for varying viscosity ratio. For drop viscosity lower than that of the matrix, the long-time steady deformation behavior is similar to that of the viscosity matched system—the drop shows reduced deformation with increasing Deborah number due to the increased inhibiting viscoelastic normal stress inside the drop. However for higher viscosity ratio systems, the drop response is non-monotonic—the steady drop deformation first decreases with increasing Deborah number but above a critical Deborah number, it increases with further increase in Deborah number, reaching higher than the viscous case value for some viscosity ratios. We explain the increase in deformation with Deborah number by noting that at higher viscosity ratios, strain rate inside the drop is reduced, thereby reducing the inhibiting viscoelastic stress. Furthermore, similar to the viscosity matched system, the drop inclination angle increases with increasing Deborah number. A drop aligned more with the maximum stretching axis at 45 degree of the imposed shear, experiences increased viscous stretching. With increased ratio of polymeric viscosity to total drop viscosity, the drop deformation decreases and the inclination angle increases. Our simulation results compare favorably with a number of experimental and computational results from other researchers.     

Flow of linear molecules through a 4:1:4 contraction–expansion using non-equilibrium molecular dynamics: Extensional rheology and pressure drop

In this work, non-equilibrium molecular dynamics simulations are used to generate the flow of linear polymer chains (monomer-springs with FENE potential) and a Lennard–Jones fluid (Newtonian fluid) through a contraction–expansion (4:1:4) geometry. An external force field simulating a constant pressure gradient upstream the contraction region induces the flow, where the confining action of the walls is represented by a Lennard–Jones potential. The equations of motion are solved through a multiple-step integration algorithm coupled to a Nosé-Hoover dynamics [S. Nose, A unified formulation of the constant temperature molecular dynamics methods, J. Chem. Phys. 81 (1984) 511–519], i.e., to simulate a thermostat, which maintains a constant temperature. In this investigation, we assume that the energy removed by the thermostat is related to the viscous dissipation along the contraction–expansion geometry. A non-linear increasing function between the pressure drop and the mean velocity along the contraction for the linear molecules is found, being an order of magnitude larger than that predicted for the Lennard–Jones fluid. The pressure drop of both systems (the linear molecules and Lennard–Jones fluid) is related to the dissipated energy at the contraction entry. The large deformation that the linear molecules experience and the evolution of the normal stress at the contraction entry follow a different trajectory in the relaxation process past the contraction, generating large hysteresis loops. The area enclosed by these cycles is related to the dissipated energy. Large shear stresses developed near the re-entrant corners as well as the vortex formation, dependent on the Deborah number, are also predicted at the exit of the contraction. To our knowledge, for the first time, the excessive pressure losses found in experimental contraction flows can be explained theoretically.     

Chaotic streamlines in the flow of knotted and unknotted vortices

This paper describes the motion and the flow induced by a thin tubular vortex coiled on a torus. The vortex is defined by the number of turns, p, that it makes round the torus symmetry axis and the number of turns, q, that it makes round the torus centerline. All toroidal filamentary vortices are found to progress along and to rotate round the torus symmetry axis in an almost steady manner while approximately preserving their shape. The flow, observed in a frame moving with the vortex, possesses two stagnation points. The stream tube emanating from the forward stagnation point and the stream tube ending at the backward stagnation point transversely intersect along a finite number of streamlines. This produces a three-dimensional chaotic tangle whose geometry depends primarily on the value of p. Inside this chaotic shell there are two major stability tubes: the first one envelopes the vortex whereas the second one runs parallel to it and possesses the same topology. When p > 2 there is an additional stability tube enveloping the torus centerline.     

Viscoelastic fluid behavior in annulus using Giesekus model

An analytical solution is derived for the steady state, laminar, axial, fully developed flow of a viscoelastic fluid obeying the Giesekus model without any retardation time in a concentric annulus.An approximation is used for the estimation of radial normal stress. The influence of Deborah number (De) and the mobility factor (α) on the velocity profile, axial pressure gradient are investigated and results show strong effects of mobility factor and Deborah number on above parameters.     

A new mixed finite element for calculating viscoelastic flow

A new mixed finite element has allowed us to calculate flows of Maxwell-B and Oldroyd-B fluids at very high values of the Deborah number, De. The element is divided into several bilinear sub-elements for the stresses, while streamline-upwinding is used for discretizing the constitutive equation. The method is applied to the stick-slip problem, the flow through a tapered contraction and the flow through four-to-one abrupt plane and circular contractions. Important corner vortices develop at high values of De in the circular contraction. We have not encountered upper limits for the Deborah number in our calculations with Oldroyd-B fluids.     

Linear stability and dynamics of viscoelastic flows using time-dependent numerical simulations

Ultimately, numerical simulation of viscoelastic flows will prove most useful if the calculations can predict the details of steady-state processing conditions as well as the linear stability and non-linear dynamics of these states. We use finite element spatial discretization coupled with a semi-implicit θ-method for time integration to explore the linear and non-linear dynamics of two, two-dimensional viscoelastic flows: plane Couette flow and pressure-driven flow past a linear, periodic array of cylinders in a channel. For the upper convected Maxwell (UCM) fluid, the linear stability analysis for the plane Couette flow can be performed in closed form and the two most dangerous, although always stable, eigenvalues and eigenfunctions are known in closed form. The eigenfunctions are non-orthogonal in the usual inner product and hence, the linear dynamics are expected to exhibit non-normal (non-exponential) behavior at intermediate times. This is demonstrated by numerical integration and by the definition of a suitable growth function based on the eigenvalues and the eigenvectors. Transient growth of the disturbances at intermediate times is predicted by the analysis for the UCM fluid and is demonstrated in linear dynamical simulations for the Oldroyd-B model. Simulations for the fully non-linear equations show the amplification of this transient growth that is caused by non-linear coupling between the non-orthogonal eigenvectors. The finite element analysis of linear stability to two-dimensional disturbances is extended to the two-dimensional flow past a linear, periodic array of cylinders in a channel, where the steady-state motion itself is known only from numerical calculations. For a single cylinder or widely separated cylinders, the flow is stable for the range of Deborah number (De) accessible in the calculations. Moreover, the dependence of the most dangerous eigenvalue on De≡λV/R resembles its behavior in simple shear flow, as does the spatial structure of the associated eigenfunction. However, for closely spaced cylinders, an instability is predicted with the critical Deborah number De_c scaling linearly with the dimensionless separation distance L between the cylinders, that is, the critical Deborah number De_Lc≡λV/L is shown to be an O(1) constant. The unstable eigenfunction appears as a family of two-dimensional vortices close to the channel wall which travel downstream. This instability is possibly caused by the interaction between a shear mode which approaches neutral stability for De ≫ 1 and the periodic modulation caused by the presence of the cylinders. Nonlinear time-dependent simulations show that this secondary flow eventually evolves into a stable limit cycle, indicative of a supercritical Hopf bifurcation from the steady base state.     

Instabilities in multi-hole converging flow of viscoelastic fluids

Studies of the onset of instabilities were conducted on single hole and multi-hole contractions using laser speckle visualization. A well characterized elastic fluid was used with constant viscosity of 13.1 Pa · s and elasticity characterized by a longest relaxation time constant of 2.233 s. The onset of instabilities was characterized in terms of the Deborah number and the contraction ratio. Three types of instabilities were observed: pulsing vortices, azimuthally rotating vortices, and swirling vortices. For the single hole contractions the critical Deborah number for instability increased from 4.4 to 5.07 to 5.25 as the contraction ratio increased from 4: 1 to 8: 1 to 12: 1. The magnitude of the instabilities was much greater for the 4: 1 contraction than for the other two contraction ratios. For the multi-hole contraction a square array of nine holes was used and the ratio of the hole diameter to hole spacing was varied. The height of the vortices is very similar for the single hole and multi-hole contractions at low Deborah numbers. At high Deborah numbers the effect of adjacent holes is to reduce the height of the vortices by a factor of three. For the 4: 1 spacing no secondary vortex was observed below a Deborah number of De = 3.7. Secondary vortices occurred for the 8:1 and 10:1 spacing at all Deborah numbers. Unstable pulsing vortices appeared for all spacings at a critical Deborah number around 5.5. Adjacent holes decreased the strength of the unsteady vortex motions. The centerline velocities were measured for the multi-hole contraction at shear rates of 5, 30, and 300 s^–1. The elongational strain rates are similar at a low shear rate of 5 s^–1. As shear rate is increased the onset of stretching occurs closer to the plane of the contraction for the smaller contraction ratios.     

Numerical solutions of pulsating flow and heat transfer characteristics in a channel with a backward-facing step

The incompressible laminar flow of air and heat transfer in a channel with a backward-facing step is studied for steady cases and for pulsatile inlet conditions. For steady flows the influence of the inlet velocity profile, the height of the step and the Reynolds number on the reattachment length is investigated. A parabolic entrance profile was used for pulsatile flow. It was found with amplitude of oscillation of one by Re=100 that the primary vortex breakdown through one pulsatile cycle. The wall shear rate in the separation zone varied markedly with pulsatile flows and the wall heat transfer remained relatively constant. The time-average pulsatile heat transfer at the walls was greater as with steady flow with the same mean Reynolds number.     

A numerical study of steady and unsteady viscoelastic flow past bounded cylinders

We consider two-dimensional, inertia-free, flow of a constant-viscosity viscoelastic fluid obeying the FENE-CR equation past a cylinder placed symmetrically in a channel, with a blockage ratio of 0.5. Through numerical simulations we show that the flow becomes unsteady when the Deborah number (using the usual definition) is greater than De ≈ 1.3, for an extensibility parameter of the model of L² = 144. The transition from steady to unsteady flow is characterised by a small pulsating recirculation zone of size approximately equal to 0.15 cylinder radius attached to the downstream face of the cylinder. There is also a rise in drag coefficient, which shows a sinusoidal variation with time. The results suggest a possible triggering mechanism leading to the steady three-dimensional Gortler-type vortical structures, which have been observed in experiments of the flow of a viscoelastic fluid around cylinders. The results reveal that the reason for failure of the search for steady numerical solutions at relatively high Deborah numbers is that the two-dimensional flow separates and eventually becomes unsteady. For a lower extensibility parameter, L² = 100, a similar recirculation is formed given rise to a small standing eddy behind the cylinder which becomes unsteady and pulsates in time for Deborah numbers larger than De ≈ 4.0–4.5.     

Hydroelastic instabilities in viscoelastic flow past a cylinder confined in a channel

Experimental investigation of viscoelastic flow past a cylinder is conducted for a polyisobutylene-based polymer solution. High-image-density particle image velocimetry is utilized to quantitatively determine the spatial features of elastic wake instabilities. The viscoelastic flow bifurcates from steady two-dimensional flow to steady three-dimensional flow for values of the Deborah number (dimensionless flow rate) greater than a critical value. These hydroelastic flow transitions are manifested in the form of three-dimensional cells spaced periodically along the axis of the cylinder. The elastic flow structures do not exist in the Newtonian counterpart of creeping flow past a cylinder. Received: 7 October 1998 / Accepted: 22 April 1999     

Experimental study on the flow regimes past a confined prism undergoing self-sustained oscillations

An experimental study based on Particle Image Velocimetry (PIV) is presented with the objective of studying the flow regimes that appear in the flow past a confined prism undergoing self-sustained oscillations at low Reynolds numbers (Re). The square-section prism, placed inside a 3D square cross-section vertical channel with a confinement ratio of 1/2.5, was tethered to the channel walls and, therefore, it was allowed to move freely transverse to the incoming flow. Re (based on the prism cross-section height) was varied in the range from 100 to 700. Three different prism to fluid density ratios (m¹) were considered: 0.56, 0.70, and 0.91. These two parameters, Re and m¹, were used to map the results obtained. In particular, it was found that five different regimes appear: (1) steady prism with steady recirculation bubble, (2) steady prism with unsteady vortex shedding wake, (3) large amplitude low frequency oscillating prism with unsteady vortex shedding wake, (4) small amplitude high frequency oscillating prism with unsteady vortex shedding wake, and (5) irregular/chaotic motion of both the prism and the wake. The PIV results and associated numerical simulations were used to analyze the different prism and wake states.     

Gas displacement of polymer melts in a cylinder: Experiments and viscoelastic simulations

Series of isothermal gas-displacement of a polystyrene melt in a circular cylinder were performed. The experiments show an increase in the steady fractional coverage above a Newtonian level (m = 0.6) at very low Deborah numbers, where the melt behaves like a diluted un-entangled system. At higher Deborah numbers, where the melt behaves as an entangled melt system, the steady fractional coverage decreases.     

Reynolds-number-effects in flow around a rectangular cylinder with aspect ratio 1:5

The paper reports on experiments carried out over a wide range of Reynolds numbers in a high pressure wind tunnel. The model was a sharp-edged rectangular cylinder with aspect ratio height/width 1:5 (width/span ratio 1:10.8), which was investigated in both basic orientations, lengthwise (4×10³<Re<4×10⁵) and perpendicular to the flow (2.7×10⁴<Re<6.4×10⁵). The Reynolds number is based on the height of the model normal to the flow. Steady and unsteady forces were measured with a piezoelectric balance. Thus along with steady (i.e. time averaged values) including the base pressure coefficient, also power spectra and probability density functions were measured yielding for example Strouhal numbers, higher statistical moments, etc. A response diagram for the vortex resonance phenomenon was taken for the natural bending motion of the slender model. If lift coefficient for constant angle of attack is plotted against Reynolds number, a significant Reynolds number effect is seen. For α=4°, the curve shows an inflection point and the lift varies between 0.3 and 0.6. For α=6° and 2° there are similar variations shifted to lower and higher values of Re, respectively. Probably the shapes of separation bubbles that depend on the Reynolds number are responsible for these effects. No Reynolds number effects were observed when the long side was normal to the flow, an orientation where reattachment at the side walls is not possible. Comparing both basic cases (α=0° and 90°), the interpretation of the probability distributions of lift force leads to the conclusion that the possibility of reattachment (α=0°) seems to enhance the degree of order in the vortex shedding process.     

Bifurcation phenomena in viscoelastic flows through a symmetric 1:4 expansion

In this work we present an investigation of viscoelastic flow in a planar sudden expansion with expansion ratio D/d = 4. We apply the modified FENE–CR constitutive model based on the non-linear finite extensibility dumbbells (FENE) model. The governing equations were solved using a finite volume method with the high-resolution CUBISTA scheme utilised for the discretisation of the convective terms in the stress and momentum equations. Our interest here is to investigate two-dimensional steady-state solutions where, above a critical Reynolds number, stable asymmetric flow states are known to occur. We report a systematic parametric investigation, clarifying the roles of Reynolds number (0.01 < Re < 100), Weissenberg number (0 < We < 100) and the solvent viscosity ratio (0.3 < β < 1). For most simulations the extensibility parameter of the FENE model was kept constant, at a value L² = 100, but some exploration of its effect in the range 100–500 shows a rather minor influence. The results given comprise flow patterns, streamlines and vortex sizes and intensities, and pressure and velocity distributions along the centreline (i.e. y = 0). For the Newtonian case, in agreement with previous studies, a bifurcation to asymmetric flow was observed for Reynolds numbers greater than about 36. In contrast viscoelasticity was found to stabilise the flow; setting β = 0.5 and We = 2 as typical values, resulted in symmetric flow up to a Reynolds number of about 46. We analyse these two cases in particular detail.     

Steady vortex flows of a self-gravitating gas

An exact solution of equations of steady motion of a self-gravitating gas is found. This solution describes a vortex flow of the gas from the surface of a spherical source and is partly invariant with respect to the group of rotations (Ovsyannikov vortex, singular vortex). The factor-system of the solution is reduced to finite formulas and one ordinary differential equation of the third order. Various regimes of gas motion described by this solution are determined: unlimited spreading of the gas with swirling from the surface of a spherical source and gas exhaustion with formation of a sphere with elevated density at a finite distance from this source.     

Vortex shedding from two finite circular cylinders in a staggered configuration

Wind tunnel experiments were conducted to measure the vortex shedding frequencies for two circular cylinders of finite height arranged in a staggered configuration. The cylinders were mounted normal to a ground plane and were partially immersed in a flat-plate turbulent boundary layer. The Reynolds number based on the cylinder diameter was Re_D=2.4×10⁴, the cylinder aspect ratio was AR=9, the boundary layer thickness relative to the cylinder height was δ/H=0.4, the centre-to-centre pitch ratio was varied from P/D=1.125 to 5, and the incidence angle was incremented in small steps from α=0° to 90°. The Strouhal numbers were obtained behind the upstream and downstream cylinders using hot-wire anemometry. From the behaviour of the Strouhal number data obtained at the mid-height position, the staggered configuration could be broadly classified by the pitch ratio as closely spaced (P/D<1.5), moderately spaced (1.5?P/D?3), or widely spaced (P/D>3). The closely spaced staggered finite cylinders were characterized by the same Strouhal number measured behind both cylinders, an indication of single bluff-body behaviour. Moderately spaced staggered finite cylinders were characterized by two Strouhal numbers at most incidence angles. Widely spaced staggered cylinders were characterized by a single Strouhal number for both cylinders, indicative of synchronized vortex shedding from both cylinders at all incidence angles. For selected staggered configurations representative of closely spaced, moderately spaced, or widely spaced behaviour, Strouhal number measurements were also made along the vertical lengths of the cylinders, from the ground plane to the free end. The power spectra showed that for certain cylinder arrangements, because of the influences of the cylinder–wall junction and free-end flow fields, the Strouhal numbers and flow patterns change along the cylinder.     

Predictions of the excess pressure drop of Boger fluids through a 2:1:2 contraction–expansion geometry using non-equilibrium molecular dynamics

Non-equilibrium molecular dynamics are used to generate the flow of polymer solutions, specifically of Boger fluids, through a planar 2:1:2 contraction–expansion geometry. The solvent molecules are represented by Lennard–Jones particles, while linear molecules are described by spring-monomers with a finite extensible non-linear elastic spring potential. The equations for Poiseuille flow are solved using a multiple time-scale algorithm extended to non-equilibrium situations. Simulations are performed at constant temperature using Nose–Hoover dynamics. At simulation conditions, changes in concentration show no significant effect on molecular conformation, velocity profiles, and stress fields, while variations in the Deborah number have a strong influence on fluid response. Increasing the magnitude of the Deborah number (De), larger deformation rates are developed in the flow region. For a Deborah number of one, the non-dimensional pressure drop presents values lower than the correspondent Newtonian case. However, for large Deborah numbers, the pressure drop increases above the Newtonian reference. An effective excess pressure drop above the Newtonian value is predicted for Boger fluids along this geometry.