On the development of Taylor vortices in a vertical annulus with a heated rotating inner cylinder

A numerical study has been conducted to determine the heat transfer characteristics and flow patterns which develop around a rotating, heated vertical cylinder enclosed within a stationary concentric cylinder. A tall annulus (aspect ratio of 10) with fixed, adiabatic horizontal end-plates and a radius ratio of 0·5 has been considered. Furthermore, the effect that the introduction of buoyancy forces by heating the inner cylinder has on the development of the Taylor vortex flow is examined. It is observed that the formation of the Taylor vortices is delayed until the rotational parameter σ = Gr/Re² has a value below unity for any given Reynolds number Re which is above the critical value Re_crit for the formation of Taylor vortices in an isothermal flow. Also, the Taylor cells first appear at the top of the annulus. As σ is gradually decreased below unity, bifurcations to other states are observed. The final structure of the secondary flow is noticeably distorted in the mixed-convection mode, with the size of the Taylor cells varying greatly along the height of the annulus. This distortion diminishes as σ is further decreased, until the isothermal flow pattern is nearly recovered below σ = 0·01.     

Characterization of fluid flow patterns and heat transfer in horizontal channel mixed convection

Two mechanisms of roll initiation are highlighted in a horizontal channel flow, uniformly heated from below, at constant heat flux (Γ = 10, Pr = 7, 50 ≤ Re ≤ 100, 0 ≤ Ra ≤ 10⁶). The first mechanism is the classical one, it occurs for low Rayleigh numbers and is initiated by the lateral wall effect. The second occurs for higher Rayleigh numbers and combines the previous effect with a supercritical vertical temperature gradient in the lower boundary layer, which simultaneously triggers pairs of rolls in the whole zone in between the two lateral rolls. We have found that in the present configuration, the transition between the two roll initiation mechanisms occurs for Ra/Re ² ≈ 18. Consequently, the heat transfer is significantly enhanced compared to the pure forced convection case owing to the flow pattern responsible of the continuous flooding the heated wall with cold fluid.     

Average and extremal properties of heat transfer and shear stress on a wall surface in Rayleigh–Bénard convection

In this study surface-averaged and extremal properties of heat transfer and shear stress on the upper wall surface of Rayleigh–Bénard convection are numerically examined. The Prandtl number was raised up to 10³, and the Rayleigh number was changed between 10⁴ and 10⁷. As a result, average Nusselt number Nu and shear rate τ/Pr depends on Pr, Ra, and the entire numerical results are distributed between two correlation equations corresponding to small and large Pr. The small and large Pr equations are closely related to steady and unsteady flow regimes, respectively. Nevertheless, a single relation τ/Pr ~ Nu ^3.0 exists to explain the entire results. Similarly the change of local maximal properties Nu _max and τ _max/Pr depends on Pr, Ra, and these values are also distributed between two correlation equations corresponding to small and large Pr cases. Despite such complicated dependence we can obtain a correlation equation as a form of τ _max/Pr ~ Nu _max^2.6, which has not been obtained theoretically.     

Large-scale flow structures and heat transport of turbulent forced and mixed convection in a closed rectangular cavity

Results of an experimental study of flow structure formation and heat transport in turbulent forced and mixed convection are presented. The experiments were conducted in a rectangular cavity with a square cross section, which has an aspect ratio between length and height of Γ_xz = 5. Air at atmospheric pressure was used as working fluid. The air inflow was supplied through a slot below the ceiling, while exhausting was provided by another slot, which is located directly above the floor. Both vents extend over the whole length of the cell. In order to induce thermal convection the bottom of the cell is heated while the ceiling is maintained at a constant temperature. This configuration allows to generate and study mixed convection under well defined conditions. Results of forced convection at Re = 1.07 × 10⁴ as well as mixed convection at 1.01 × 10⁴ ? Re ? 3.4 × 10⁴ and Ra = 2.4 × 10⁸ (3.3 ? Ar ? 0.3), which were obtained by means of Particle Image Velocimetry and local temperature measurements, are presented. For purely forced convection a 2D mean wind, which can be approximated by a solid body rotation, is found. With increasing Archimedes number this structure becomes unstable, leading to a transition of the solid body rotation into additional smaller convection rolls. Proper orthogonal decomposition of the instantaneous velocity fields has been performed for further analysis of these coherent large-scale structures. Their fingerprint is found in the spatial temperature distribution of the out flowing air at the end of the outlet channel, which reveals a temporally stable profile with two maxima over the length of the outlet. Moreover a maximum in the global enthalpy transport by the fluid is found at Ar ≈ 0.6.     

Mixed convection boundary layer flow over a horizontal circular cylinder with Newtonian heating

The steady mixed convection boundary layer flow over a horizontal circular cylinder, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered in this study. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically using a numerical scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient Re ^1/2 C _f and the local wall temperature θ _w (x) as well as the velocity and temperature profiles with two parameters, namely the mixed convection parameter λ and the Prandtl number Pr.     

Thermal convection heat transfer of air/water layers enclosed in horizontal annuli with mixed boundary conditions

This paper presents a numerical study for thermal convection of air/water layers enclosed in a horizontal concentric/eccentric cylindrical annulus with the inner wall subjected to a constant heat flux and an isothermal condition at the outer wall. Results are generated for an annulus of radius ratio of 2.6 with three different vertically eccentric positions of the inner cylinder, the modified Rayleigh number up to 10⁷, and the modified Marangoni number varying from 0.0 to 10⁴. The heat transfer characteristics and fluid flow patterns driven by buoyancy and/or thermocapillary effects in air/water-filled annulus are illustrated by means of contour maps of healtline and streamline respectively. Under the mixed thermal boundary conditions considered, significant heat exchange across the air-water interface is found and appears to be unaffected by the thermocapillary convection.Es wird eine numerische Untersuchung der thermischen Konvektion von Luft/Wasser-Schichten, die von horizontalen konzentrischen/exzentrischen zylindrischen Ringkanälen umschlossen sind, deren innere Wand einer konstanten Wärmestromdichte und deren äußere Wand isothermen Bedingungen ausgesetzt sind, dargestellt. Ergebnisse werden für einen Ringkanal mit einem Radiusverhältnis von 2,6 mit drei verschiedenen senkrechten exzentrischen Positionen des inneren Zylinders, für die modifizierte Rayleigh-Zahl bis 10⁷ und für die modifizierte Marangoni-Zahl zwischen 0,0 und 10⁴ erhalten. Die Charakteristiken der Wärmeübertragung und des Fluid-Stromlinienbildes, hervorgerufen durch Auftriebskräfte und/oder thermokapillare Effekte in mit Luft/Wasser gefüllten Ringkanälen wurden mittels Schichtlinienplänen der Isothermen und der Stromlinien dargestellt. Bei den untersuchten gemischten thermischen Randbedingungen wurde ein signifikanter Wärmeaustausch über der Luft/Wasser-Grenzfläche gefunden, der scheinbar nicht durch die thermokapillare Konvektion gestört wird.     

Prediction of natural convection flow using network model and numerical simulations inside enclosure with distributed solid blocks

Steady state natural convection of a fluid with Pr ≈ 1 within a square enclosure containing uniformly distributed, conducting square solid blocks is investigated. The side walls are subjected to differential heating, while the top and bottom ones are kept adiabatic. The natural convection flow is predicted employing the nondimensional volumetric flow rate (Q_max^* Q_{\max }^{*} ) by using a network model and also using numerical simulations. For identical solid and fluid thermal conductivities (i.e. k _s  = k _f ), a parametric study of the effect of number of blocks (N ²), gap size (δ) and enclosure Rayleigh number (Ra) on Q_max^* Q_{\max }^{*} is performed using the two approaches. Network model predictions are observed to agree well with that from the simulations until Raδ³ ~ 12. Considering the enclosure with blocks as a porous medium, for a fixed enclosure Ra number, increasing the number of blocks for a fixed volumetric porosity leads to a decrease in enclosure permeability, which in turn reduces the flow rate. When the number of blocks is fixed, and for a given Ra number, the flow rate increases as the porosity increases by widening the gap between the blocks.     

Non-Darcy mixed convection from a vertical plate in saturated porous media-variable surface heat flux

Nonsimilarity solutions for non-Darcy mixed convection from a vertical impermeable surface embedded in a saturated porous medium are presented for variable surface heat flux (VHF) of the power-law form. The entire mixed convection region is divided into two regimes. One region covers the forced convection dominated regime and the other one covers the natural convection dominated regime. The governing equations are first transformed into a dimensionless form by the nonsimilar transformation and then solved by a finite-difference scheme. Computations are based on Keller Box method and a tolerance of iteration of 10⁻⁵ as a criterion for convergence. Three physical aspects are introduced. One measures the strength of mixed convection where the dimensionless parameter Ra^* _x /Pe^3/2 _x characterizes the effect of buoyancy forces on the forced convection; while the parameter Pe _x /Ra^*2/3 _x characterizes the effect of forced flow on the natural convection. The second aspect represents the effect of the inertial resistance where the parameter K′U _∞/ν is found to characterize the effect of inertial force in the forced convection dominated regime, while the parameter (K′U _∞/ν)(Ra^*2/3 _x /Pe _x ) characterizes the effect of inertial force in the natural convection dominated regime. The third aspect is the effect of the heating condition at the wall on the mixed convection, which is presented by m, the power index of the power-law form heating condition. Numerical results for both heating conditions are carried out. Distributions of dimensionless temperature and velocity profiles for both Darcy and non-Darcy models are presented. Received on 26 May 1997     

Taguchi design of three dimensional simulations for optimization of turbulent mixed convection in a cavity

This study discusses the application of Taguchi method in assessing maximum heat transfer rate for the turbulent mixed convection in an enclosure embedded with rotating isothermal cylinder. The simulations were planned based on Taguchi’s L16 orthogonal array with each trial performed under different conditions of position of the cylinder, Reynolds number (Re) and Rayleigh number (Ra). The thermal lattice Boltzmann based on D3Q19 methods without any turbulent submodels was purposed to simulate the flow and thermal fields. A relaxation time method with the stability constants is introduced to solve turbulent natural convection problems. Signal-to-noise ratios (S/N) analysis were carried out in order to determine the effects of process parameters and optimal factor settings. Finally, confirmation tests verified that Taguchi method achieved optimization of heat transfer rate with sufficient accuracy.     

Mixed convection on vertical plate for fluids of any prandtl number

A mixed convection parameter=(Ra) ^1/4/(Re)^1/2, with=Pr/(1+Pr) and=Pr/(1 +Pr)^1/2, is proposed to replace the conventional Richardson number, Gr/Re², for combined forced and free convection flow on an isothermal vertical plate. This parameter can readily be reduced to the controlling parameters for the relative importance of the forced and the free convection,Ra ^1/4/(Re ^1/2 Pr ^1/3) forPr 1, and (RaPr)^1/2/(RePr ^1/2 forPr 1. Furthermore, new coordinates and dependent variables are properly defined in terms of, so that the transformed nonsimilar boundary-layer equations give numerical solutions that are uniformly valid over the entire range of mixed convection intensity from forced convection limit to free convection limit for fluids of any Prandtl number from 0.001 to 10,000. The effects of mixed convection intensity and the Prandtl number on the velocity profiles, the temperature profiles, the wall friction, and the heat transfer rate are illustrated for both cases of buoyancy assisting and opposing flow conditions.

---

Mischkonvektion an einer vertikalen Platte für Fluide beliebiger Prandtl-Zahl

Zusammenfassung Für die kombinierte Zwangs- und freie Konvektion an einer isothermen senkrechten Platte wird ein Mischkonvektions-Parameter=( Ra) ^1/4 (Re)^1/2, mit=Pr/(1 +Pr) und=Pr/(1 +Pr)^1/2 vorgeschlagen, den die gebräuchliche Richardson-Zahl, Gr/Re², ersetzen soll. Dieser Parameter kann ohne weiteres auf die maßgebenden Kennzahlen für den relativen Einfluß der erzwungenen und der freien Konvektion reduziert werden,Ra ^1/4/(Re ^1/2 Pr ^1/3) fürPr 1 und (RaPr)^1/4/(RePr)^1/2 fürPr 1. Weiterhin werden neue Koordinaten und abhängige Variablen als Funktion von definiert, so daß für die transformierten Grenzschichtgleichungen numerische Lösungen erstellt werden können, die über den gesamten Bereich der Mischkonvektion, von der freien Konvektion bis zur Zwangskonvektion, für Fluide jeglicher Prandtl-Zahl von 0.001 bis 10.000 gleichmäßig gültig sind. Der Einfluß der Intensität der Mischkonvektion und der Prandtl-Zahl auf die Geschwindigkeitsprofile, die Temperaturprofile, die Wandreibung und den Wärmeübergangskoeffizienten werden für die beiden Fälle der Strömung in und entgegengesetzt zur Schwerkraftrichtung dargestellt.

---

Nomenclature C _f local friction coefficient - C _p specific heat capacity - f reduced stream function - g gravitational acceleration - Gr local Grashoff number,g T _w –T )x³/v² - Nu local Nusselt number - Pr Prandtl number,v/ - Ra local Rayleigh number,g T _w –T x ³/( v) - Re local Reynolds number,u x/v - Ri Richardson number,Gr/Re ² - T fluid temperature - T _w wall temperature - T free stream temperature - u velocity component in thex direction - u free stream velocity - v velocity component in they direction - x vertical coordinate measuring from the leading edge - y horizontal coordinate Greek symbols thermal diffusivity - thermal expansion coefficient - mixed convection parameter (Ra)1/4/Re)^1/2 - pseudo-similarity variable,(y/x) - ₀ conventional similarity variable,(y/x)Re ^1/2 - dimensionless temperature, (T–T T _W –T - unified mixed-flow parameter, [(Re) ^1/2 + (Ra)^1/4] - dynamic viscosity - kinematic viscosity - stretched streamwise coordinate or mixed convection parameter, [1 + (Re)^1/2/(Ra) ^1/4]^–1=/(1 +) - density - Pr/(1 + Pr) _w wall shear stress - stream function - Pr/(l+Pr)^1/3 This research was supported by a grand from the National Science Council of ROC     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Heat transfer to non-Newtonian flows over a cylinder in cross flow

Over a range of 102<Re^*<5800, 6.5<Pr^*<79, and 0.6<n<1, circumferential wall temperatures for water and aqueous polymer (purely viscous) solution flows over a smooth cylinder were measured experimentally. The cylinder was heated by passing direct electric current through it. Aqueous solutions of Carbopol 934 and EZ1 were used as power-law non-Newtonian fluids. The peripherally averaged heat transfer coefficient for purely viscous non-Newtonian fluids, at any fixed flow rate, decreases with increasing polymer concentration. A new correlation is proposed for predicting the peripherally averaged Nusselt number for power-law fluid flows over a heated cylinder in cross flow.     

A numerical study of the steady forced convection heat transfer from an unconfined circular cylinder

Forced convection heat transfer from an unconfined circular cylinder in the steady cross-flow regime has been studied using a finite volume method (FVM) implemented on a Cartesian grid system in the range as 10 ≤ Re ≤ 45 and 0.7 ≤ Pr ≤ 400. The numerical results are used to develop simple correlations for Nusselt number as a function of the pertinent dimensionless variables. In addition to average Nusselt number, the effects of Re, Pr and thermal boundary conditions on the temperature field near the cylinder and on the local Nusselt number distributions have also been presented to provide further physical insights into the nature of the flow. The rate of heat transfer increases with an increase in the Reynolds and/or Prandtl numbers. The uniform heat flux condition always shows higher value of heat transfer coefficient than the constant wall temperature at the surface of the cylinder for the same Reynolds and Prandtl numbers. The maximum difference between the two values is around 15–20%.     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

Experimental investigation of natural convection heat transfer in horizontal and inclined annular fluid layers

In the present study, an experimental investigation of heat transfer and fluid flow characteristics of buoyancy-driven flow in horizontal and inclined annuli bounded by concentric tubes has been carried out. The annulus inner surface is maintained at high temperature by applying heat flux to the inner tube while the annulus outer surface is maintained at low temperature by circulating cooling water at high mass flow rate around the outer tube. The experiments were carried out at a wide range of Rayleigh number (5 × 10⁴ < Ra < 5 × 10⁵) for different annulus gap widths (L/D _o = 0.23, 0.3, and 0.37) and different inclination of the annulus (α = 0°, 30° and 60°). The results showed that: (1) increasing the annulus gap width strongly increases the heat transfer rate, (2) the heat transfer rate slightly decreases with increasing the inclination of the annulus from the horizontal, and (3) increasing Ra increases the heat transfer rate for any L/D _o and at any inclination. Correlations of the heat transfer enhancement due to buoyancy driven flow in an annulus has been developed in terms of Ra, L/D _o and α. The prediction of the correlation has been compared with the present and previous data and fair agreement was found.     

Natural convection from a vertical cylinder in a thermally stratified porous medium

The problem of non-Darcy natural convection adjacent to a vertical cylinder embedded in a thermally stratified porous medium has been analyzed. Nonsimilarity solutions are obtained for the case that the ambient temperature increases linearly with height of the cylinder. A generalized flow model was used in the present study to include the effects of the macroscopic viscous term and the microscopic inertial force. Also, the thermal dispersion effect is considered in the energy equation. Thus, the main aim of this work is to examine the effects of thermal stratification and non-Darcy flow phenomena on the free convection flow and heat transfer characteristics. It was found that the present problem depends on six parameters, namely, the local thermal stratification parameter ξ, the boundary effect parameter Bp, the modified Grashof number Gr^*, wall temperature exponent m, the curvature parameter ω, and the modified Rayleigh number based on pore diameter Ra _d . The impacts of these governing parameters on the local heat transfer parameter are discussed in great detail. Also, representative velocity and temperature profiles are presented at selected values of the thermal stratification parameter. In general, the local heat transfer parameter is increased with increasing the values of m, ω, and Ra _d ; while it is decreased with increasing the values of ξ, Bp, and Gr^*. Received on 19 May 1998     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

Natural convection flow in a square cavity revisited: Laminar and turbulent models with wall functions

Numerical simulations have been undertaken for the benchmark problem of natural convection flow in a square cavity. The control volume method is used to solve the conservation equations for laminar and turbulent flows for a series of Rayleigh numbers (Ra) reaching values up to 10¹⁰. The k-? model has been used for turbulence modelling with and without logarithmic wall functions. Uniform and non-uniform (stretched) grids have been employed with increasing density to guarantee accurate solutions, especially near the walls for high Ra-values. ADI and SIP solvers are implemented to accelerate convergence. Excellent agreement is obtained with previous numerical solutions, while some discrepancies with others for high Ra-values may be due to a possibly different implementation of the wall functions. Comparisons with experimental data for heat transfer (Nusselt number) clearly demonstrates the limitations of the standard k-? model with logarithmic wall functions, which gives significant overpredictions.     

Mixed boundary elements for laminar flows

This paper presents a mixed boundary element formulation of the boundary domain integral method (BDIM) for solving diffusion–convective transport problems. The basic idea of mixed elements is the use of a continuous interpolation polynomial for conservative field function approximation and a discontinuous interpolation polynomial for its normal derivative along the boundary element. In this way, the advantages of continuous field function approximation are retained and its conservation is preserved while the normal flux values are approximated by interpolation nodal points with a uniquely defined normal direction. Due to the use of mixed boundary elements, the final discretized matrix system is overdetermined and a special solver based on the least squares method is applied. Driven cavity, natural and forced convection in a closed cavity are studied. Driven cavity results at Re=100, 400 and 1000 agree better with the benchmark solution than Finite Element Method or Finite Volume Method results for the same grid density with 21×21 degrees of freedom. The average Nusselt number values for natural convection 10³≤Ra≤10⁶ agree better than 0.1% with benchmark solutions for maximal calculated grid densities 61×61 degrees of freedom. Copyright © 1999 John Wiley & Sons, Ltd.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.