Further contributions on the flow past a stationary and confined cylinder: Creeping and slowly moving flow of Power law fluids

The steady flow of generalized Newtonian fluid around a stationary cylinder placed between two parallel plates was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of the Reynolds number Re and the Power law index n values. The values of the Reynolds number, based on physical and rheological properties, cylinder radius and bulk velocity, were varied between 0.0001≤Re≤10, while the Power law index values mapped the 0.50≤n≤1.50 range, allowing for the investigation of both shear-thinning and shear-thickening effects at the creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on the most important characteristics of fluid flow past circular cylinder. It is shown that for the creeping flow regime there exist finite sized redevelopment length, drag and loss coefficient. Last but not least, the present numerical results indicate that the shear-thinning viscous behaviour decreases the onset of flow separation.     

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.     

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.     

GPU accelerated simulations of three-dimensional flow of power-law fluids in a driven cube

Newtonian fluid flow in two- and three-dimensional cavities with a moving wall has been studied extensively in a number of previous works. However, relatively a fewer number of studies have considered the motion of non-Newtonian fluids such as shear thinning and shear thickening power law fluids. In this paper, we have simulated the three-dimensional, non-Newtonian flow of a power law fluid in a cubic cavity driven by shear from the top wall. We have used an in-house developed fractional step code, implemented on a Graphics Processor Unit. Three Reynolds numbers have been studied with power law index set to 0.5, 1.0 and 1.5. The flow patterns, viscosity distributions and velocity profiles are presented for Reynolds numbers of 100, 400 and 1000. All three Reynolds numbers are found to yield steady state flows. Tabulated values of velocity are given for the nine cases studied, including the Newtonian cases.     

Turbulent characteristics of shear-thinning fluids in recirculating flows

 A miniaturised fibre optic Laser-Doppler anemometer was used to carry out a detailed hydrodynamic investigation of the flow downstream of a sudden expansion with 0.1–0.2% by weight shear-thinning aqueous solutions of xanthan gum. Upstream of the sudden expansion the pipe flow was fully-developed and the xanthan gum solutions exhibited drag reduction with corresponding lower radial and tangential normal Reynolds stresses, but higher axial Reynolds stress near the wall and a flatter axial mean velocity profile in comparison with Newtonian flow. The recirculation bubble length was reduced by more than 20% relative to the high Reynolds number Newtonian flow, and this was attributed to the occurrence further upstream of high turbulence for the non-Newtonian solutions, because of advection of turbulence and earlier high turbulence production in the shear layer. Comparisons with the measurements of Escudier and Smith (1999) with similar fluids emphasized the dominating role of inlet turbulence. The present downstream turbulence field was less anisotropic, and had lower maximum axial Reynolds stresses (by 16%) but higher radial turbulence (20%) than theirs. They reported considerably longer recirculating bubble lengths than we do for similar non-Newtonian fluids and Reynolds numbers. Received: 23 February 1999/Accepted: 28 April 1999     

Laminar velocity profile development in straight pipes of circular cross section

Velocity profile development has been studied experimentally in Newtonian and some non-Newtonian fluids. The entry length for the development of 99% of the terminal axial velocity from an initially flat profile has been found to be given byZ _e = 1.1–0.112N(Re) for laminar flow Reynolds numbers between 1 and 1500 with Newtonian fluids. There were substantial increases in this length for weakly visco-elastic aqueous solutions of polyethylene oxide and polyacrylamide in the Reynolds number range (between 1 and 10) where these have been studied.     

A numerical investigation of the flow over a pair of identical square cylinders in a tandem arrangement

This paper describes a numerical study of the two‐dimensional and three‐dimensional unsteady flow over two square cylinders arranged in an in‐line configuration for Reynolds numbers from 40 to 1000 and a gap spacing of 4D, where D is the cross‐sectional dimension of the cylinders. The effect of the cylinder spacing, in the range G = 0.3D to 12D, was also studied for selected Reynolds numbers, that is, Re = 130, 150 and 500. An incompressible finite volume code with a collocated grid arrangement was employed to carry out the flow simulations. Instantaneous and time‐averaged and spanwise‐averaged vorticity, pressure, and streamlines are computed and compared for different Reynolds numbers and gap spacings. The time averaged global quantities such as the Strouhal number, the mean and the RMS values of the drag force, the base suction pressure, the lift force and the pressure coefficient are also calculated and compared with the results of a single cylinder. Three major regimes are distinguished according to the normalized gap spacing between cylinders, that is, the single slender‐body regime (G < 0.5), the reattach regime (G < 4) and co‐shedding or binary vortex regime (G ≥4). Hysteresis with different vortex patterns is observed in a certain range of the gap spacings and also for the onset of the vortex shedding. Copyright © 2011 John Wiley & Sons, Ltd.     

Forced‐convection heat transfer from tandem square cylinders in cross flow at low Reynolds numbers

This paper presents the results of a numerical study on the flow characteristics and heat transfer over two equal square cylinders in a tandem arrangement. Spacing between the cylinders is five widths of the cylinder and the Reynolds number ranges from 1 to 200, Pr=0.71. Both steady and unsteady incompressible laminar flow in the 2D regime are performed with a finite volume code based on the SIMPLEC algorithm and non‐staggered grid. A study of the effects of spatial resolution and blockage on the results is provided. In this study, the instantaneous and mean streamlines, vorticity and isotherm patterns for different Reynolds numbers are presented and discussed. In addition, the global quantities such as pressure and viscous drag coefficients, RMS lift and drag coefficients, recirculation length, Strouhal number and Nusselt number are determined and discussed for various Reynolds numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulation of turbulent flows around a circular cylinder using nonlinear eddy-viscosity modelling: steady and oscillatory ambient flows

Steady incident flow past a circular cylinder for sub- to supercritical Reynolds number has been simulated as an unsteady Reynolds-averaged Navier–Stokes (RANS) equation problem using nonlinear eddy-viscosity modelling assuming two-dimensional flow. The model of Craft et al. (Int. J. Heat Fluid Flow 17 (1996) 108), with adjustment of the coefficients of the ‘cubic’ terms, predicts the drag crisis at a Reynolds number of about 2×10⁵ due to the onset of turbulence upstream of separation and associated changes in Strouhal number and separation positions. Slightly above this value, at critical Reynolds numbers, drag is overestimated because attached separation bubbles are not simulated. These do not occur at supercritical Reynolds numbers and drag coefficient, Strouhal number and separation positions are in approximate agreement with experimental measurements (which show considerable scatter). Fluctuating lift predictions are similar to sectional values measured experimentally for subcritical Reynolds numbers but corresponding measurements have not been made at supercritical Reynolds numbers. For oscillatory ambient flow, in-line forces, as defined by drag and inertia coefficients, have been compared with the experimental values of Sarpkaya (J. Fluid Mech. 165 (1986) 61) for values of the frequency parameter, β=D²/νT, equal to 1035 and 11240 and Keulegan–Carpenter numbers, KC=U₀T/D, between 0.2 and 15 (D is cylinder diameter, ν is kinematic viscosity, T is oscillation period, and U₀ is the amplitude of oscillating velocity). Variations with KC are qualitatively reproduced and magnitudes show best agreement when there is separation with a large-scale wake, for which the turbulence model is intended. Lift coefficients, frequency and transverse vortex shedding patterns for β=1035 are consistent with available experimental information for β≈250−500. For β=11240, it is predicted that separation is delayed due to more prominent turbulence effects, reducing drag and lift coefficients and causing the wake to be more in line with the flow direction than transverse to it. While these oscillatory flows are highly complex, attached separation bubbles are unlikely and the flows probably two dimensional.     

Turbulent pipe flow of a drag-reducing rigid “rod-like” polymer solution

Fully developed turbulent pipe flow of an aqueous solution of a rigid “rod-like” polymer, scleroglucan, at concentrations of 0.005% (w/w) and 0.01% (w/w) has been investigated experimentally. Fanning friction factors were determined from pressure-drop measurements for the Newtonian solvent (water) and the polymer solutions and so levels of drag reduction for the latter. Mean axial velocity u and complete Reynolds normal stress data, i.e. u′, v′ and w′, were measured by means of a laser Doppler anemometer at three different Reynolds numbers for each fluid. The measurements indicate that the effectiveness of scleroglucan as a drag-reducing agent is only mildly dependent on Reynolds number. The turbulence structure essentially resembles that of flexible polymer solutions which also lead to low levels of drag reduction.     

Effects of shear-thinning/thickening nature on the stability of Couette-like flow with uniform crossflow

The linear stability of wall-injected pressure- driven Couette-like flow in power-law fluids is studied. Previous study on this kind of flow for Newtonian fluids by Nicoud and Angilella [Phys. Rev. E 56, 3000 （1997）] was extended to power-law fluids to understand the effects of shear-thinning/thickening nature on the flow stability. A related expression between the critical crossflow Reynolds number for Newtonian fluids and that for power-law fluids is obtained as the streamwise Reynolds number is large enough based on numerical computations, and verified theoretically in the case of a limiting condition with the power-law index.     

Drag of a sphere in dilute polymer solutions in high Reynolds number range

We have measured by means of four ultrasonic transducers the fall velocity of a sphere at high Reynolds number range in dilute polyacrylamide solutions which have viscoelastic effects. The polymer solutions were 5, 20 and 50ppm in the concentration. Basset-Bousinessq-Oseen equation for the falling sphere was analyzed numerically on Newtonian fluids in order to compare with the fall velocity of a sphere in the polymer solutions, and the experimental data of the fall velocity in tap water is in agreement with the range of no effect of the test tank wall. In polymer solutions, it was shown that the fall velocity is larger than that in Newtonian fluids within the critical Reynolds number range such that the drag reduction occurs and is smaller than that of Newtonian fluids over the range. The experimental data for the drag reduction ratio of polymer solutions is arranged by Weissenberg number calculating the experimental data of the first normal stress differences. It was shown that the maximum drag reduction ratio in the polymer solutions lies in the range of We=3∼10. Received: 15 October 1997 Accepted: 12 May 1998     

Experimental and numerical investigations of pressure drop in a rectangular duct with modified power law fluids

Numerical solutions are presented for fully developed laminar flow for a modified power law fluid (MPL) in a rectangular duct. The solutions are applicable to pseudoplastic fluids over a wide shear rate range from Newtonian behavior at low shear rates, through a transition region, to power law behavior at higher shear rates. The analysis identified a dimensionless shear rate parameter which, for a given set of operating conditions, specifies where in the shear rate range a particular system is operating, i.e. in the Newtonian, transition, or power law regions. The numerical results of the friction factor times Reynolds number for the Newtonian and power law region are compared with previously published results showing agreement within 0.05% in the Newtonian region, and 0.9% and 5.1% in the power law region. Rheological flow curves were measured for three CMC-7H4 solutions and were found to be well represented by the MPL constitutive equation. The friction factor times Reynolds number values were measured in the transition region for which previous measurements were unavailable. Good agreement was found between experiment and calculation thus confirming the validity of the analysis.     

Hydrodynamic forces acting on a rigid fixed sphere in early transitional regimes

A spectral – spectral-element code is used to investigate the hydrodynamic forces acting on a fixed sphere placed in a uniform flow in the Reynolds number interval [10–320] covering the early stages of transition, i.e. the steady axisymmetric regime with detached flow, the steady non-axisymmetric and the unsteady periodic regimes of the sphere wake. The mentioned changes of regimes, shown by several authors to be related to a regular and a Hopf bifurcations in the wake, result in significant changes of hydrodynamic action of the flow on the sphere. In the present paper, we show that the loss of axisymmetry is accompanied not only by an onset of lift but also of a torque and we give accurate values of drag, lift and torque in the whole interval of investigated Reynolds numbers. Among other results show, moreover, that each bifurcation is accompanied also by a change of the trend of the drag versus Reynolds number dependence, the overall qualitative effect of instabilities being an increase of drag.     

Stokes paradox for power-law flow around a cylinder

A simple analysis for power-law fluids shows that the Stokes paradox for creeping flow around a cylinder is removed for shear-thinning (n < 1) but not for shear-thickening (n 1) fluids. An approximate drag value is found for n < 1 and is compared with computed results.     

The flow past a sphere in a cylindrical tube: effects of intertia,shear-thinning and elasticity

Numerical solutions are presented for the flow past a sphere placed at the centreline of a cylindrical tube for Reynolds numbers ranging from 0 to 150, using a boundary element method. Fluids are modelled by a variety of constitutive equations including the Newtonian, the Carreau and the Phan-Thien-Tanner models. The influence of inertia, shear-thinning and fluid elasticity on the flow field, drag and the pressure drop force-drag ratio is examined. Some results are compared with available experimental data.     

The stratified flow of gas and non-newtonian liquid in horizontal pipes

Predictions of pressure drop and holdup are presented for the stratified flow of gas and non-Newtonian liquid obeying the Ostwald-de Waele power law model. The model of Taitel & Dukler (1976) for gas/Newtonian liquid flow is extended to liquids possessing either shear-thinning or shear-thickening laminar flow behaviour and computed results are given for flow behaviour indices in the range $\text{0.1 ≤ n ≤ 2}$. In particular, conditions are defined for drag reduction of the liquid flow by the presence of the gas. It is concluded that drag reduction occurs over the largest ranges of liquid and gas flow rates at the lowest $\text{n}$ values, provided that liquid flow remains laminar, but that maximum drag reduction may be expected for shear-thickening liquids with $\text{n}$ values of 2 or greater. Ratios of the liquid flow rate in the presence of gas to that for liquid flow alone under a constant pressure gradient are also presented. These ratios frequently exceed unity and are greatest for highly shear-thinning liquids.Although the Taitel & Dukler approach is consistent with experiments on gas/Newtonian liquid flow, and, in addition, appears to be valid for immiscible Newtonian liquid-liquid systems, provided that the viscosity ratio of the two phases is at least five, experiments are required to confirm its applicability for gas/non-Newtonian systems.     

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.     

New contributions on laminar flow of inelastic non-Newtonian fluid in the two-dimensional symmetric expansion: Creeping and slowly moving flow conditions

The steady flow of generalized Newtonian fluid in a two-dimensional 1:3 sudden expansion was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of Reynolds number and Power law index values. The values of generalized Reynolds number, based on physical and rheological properties, upstream channel height and bulk velocity, were varied between 0.0001 ≤ Re_gen ≤ 10, while the Power law index values mapped the 0.60 ≤ n ≤ 1.40 range, allowing for the investigation of both shear-thinning and shear-thickening effects at creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on most important characteristics of recirculating fluid flow in the downstream section of sudden expansion geometry. It is shown that for the creeping flow regime there exist finite sized redevelopment length, extra pressure drop (Couette correction) and recirculation zones (also called as Moffatt vortices) that are influenced by the non-Newtonian viscous behaviour.     

Numerical simulation of shear effect on vortex shedding behind a square cylinder

This paper is concerned with the numerical simulation of the flow structure around a square cylinder in a uniform shear flow. The calculations were conducted by solving the unsteady 2D Navier–Stokes equations with a finite difference method. The effect of the shear parameter K of the approaching flow on the vortex-shedding Strouhal number and the force coefficients acting on the square cylinder is investigated in the range K=0·0–0·25 at various Reynolds numbers from 500 to 1500. The computational results are compared with some existing experimental data and previous studies. The effect of shear rate on the Strouhal number and the force acting on the cylinder has a tendency to reduce the oscillation. The Strouhal number, mean drag and amplitude of the fluctuating force tend to decrease as the shear rate increases, but show no significant change at low shear rate. Increasing the Reynolds number decreases the Strouhal number and increases the force acting on the cylinder. At high shear rate the shedding frequencies of the fluctuating drag and lift coefficients are identical. © 1997 John Wiley & Sons, Ltd.