Experimental study of flow past a square cylinder at high Reynolds numbers

 Measurements of two-components of velocity in the wake of a square cylinder using a hot-wire anemometer are reported. Two Reynolds numbers, namely 8700 and 17,625, have been considered. The measurements were carried out in a low-speed, low-turbulence wind tunnel. Benchmark experiments at much lower Reynolds numbers show good agreement between the present experiments and published results. At higher Reynolds numbers, the experimental data reveal anticipated trends in terms of wake recovery and turbulence decay. Both velocity and velocity fluctuations show symmetry about the wake axis. The experimental data have been compared with the large eddy simulation (LES) calculation reported by Wang et al. [University of Illinois at Urbana – Champaign (1996) Report CFD 96-03] and LDV measurements of Lyn et al. [J Fluid Mech (1995) 304: 285–319]. The agreement among the three sets is generally acceptable in terms of the time-averaged velocity components, but not the velocity fluctuations. The turbulence fluctuations in the present experiments are seen to be lower than in the referred work. The differences have been traced to factors such as the aspect ratio, blockage ratio and upstream turbulence. Experiments with increased upstream turbulence did show a reduction in the discrepancy between the present experiments and the published data. An assessment of the experimental data in terms of physical mechanisms revealed that (a) streamwise normal stresses were correlated with the vortex centers, and (b) the turbulence kinetic energy profiles are similar to the turbulence shear stress. Spectral analysis of the velocity signals was carried out in the present work. Energy transfer from the mean flow to the streamwise velocity fluctuation was confirmed in the near wake. A redistribution of the kinetic energy between the streamwise and transverse components of velocity over a longer distance downstream was subsequently observed. Received: 17 May 1999/Accepted: 29 December 1999     

Numerical modelling of confined flow past a cylinder of square cross-section at various orientations

Results are presented for the unsteady, two-dimensional flow and heat transfer due to a square obstruction of diameter d located asymmetrically between the parallel sliding walls of a channel with length-to-height ratio W/H = 6·44. Analysis is based on the numerical solution of spatially and temporally second-order accurate finite difference approximations of the transport equations expressed in curvilinear co-ordinates. Laminar, constant property flow is assumed for obstruction configurations in which the blockage ratio is d/H = 0·192, the nearest-wall distances are g/d = 0·2, 0·5 and 1, the orientation angles are α=0°, 10° and 20° and the Reynolds numbers are Re=100, 500, and 1000. Preparatory testing of the numerical procedure was performed for a variety of documented flows to verify its physiconumerical accuracy and obtain estimates of the residual grid-dependent uncertainties in the variables calculated. Heat transfer, drag and lift coefficients and Strouhal numbers for the present flow were finally calculated to within 4%–7% of their grid-dependent values using non-uniformly spaced grids consisting of (x=99, y=55) nodes. Above a critical value of the Reynolds number, which depends on the geometrical parameters, the flow is characterized by alternate vortex shedding from the obstruction top and bottom surfaces. Streamline, vorticity and particle streakline plots provide qualitative impressions of the unsteady vortical flow. Especially noteworthy are the extremes in the lift coefficient which ranges from large positive values for an obstruction with g/d=0·2 and α=10° to negative values for one with g/d=0·5 and α=0°. Both the drag and lift coefficients as well as the Strouhal number exhibit non-monotonic variations with respect to the parameters explored. Asymmetries in the obstruction location and orientation account for relatively large vortex-induced periodic variations in heat transfer, especially along the wall nearest the obstruction. Notable differences are also predicted for the heat transfer coefficients of the individual obstruction surfaces as a function of the orientation angle.     

Control of mean and fluctuating forces on a circular cylinder at high Reynolds numbers

A narrow strip is used to control mean and fluctuating forces on a circular cylinder at Reynolds numbers from 2.0 × 10⁴ to 1.0 × 10⁵. The axes of the strip and cylinder are parallel. The control parameters are strip width ratio and strip position characterized by angle of attack and distance from the cylinder. Wind tunnel tests show that the vortex shedding from both sides of the cylinder can be suppressed, and mean drag and fluctuating lift on the cylinder can be reduced if the strip is installed in an effective zone downstream of the cylinder. A phenomenon of mono-side vortex shedding is found. The strip-induced local changes of velocity profiles in the near wake of the cylinder are measured, and the relation between base suction and peak value in the power spectrum of fluctuating lift is studied. The control mechanism is then discussed from different points of view. The project supported by the National Natural Science Foundation of China (10172087 and 10472124). The English text was polished by Yunming Chen.     

Nonstationary ideal incompressible fluid flows: Conditions of existence and uniqueness of solutions

The problem of formulating minimal conditions on input data that can guarantee the existence and uniqueness of solutions of the boundary value problems describing non-one-dimensional ideal incompressible fluid flow is considered using as an example the initial boundary value problem in a space-time cylinder constructed on a bounded flow domain with the nonpenetration condition on its boundary (which corresponds to fluid flow in a closed vessel). The existence problems are considered only for plane flows, and the uniqueness issues for three-dimensional flows as well. The required conditions are obtained in the form of conditions specifying that the vorticity belongs to definite functional Orlicz spaces. The results are compared with well-known results. Examples are given of admissible types of singularities for which the obtained results are valid, which is a physical interpretation of these results. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 130–145, July–August, 2008.     

Formation of the free surface of a viscous fluid volume inside a rotating horizontal cylinder

Two-dimensional viscous flow with a free surface in a horizontal cylinder rotating at a constant speed is investigated numerically using the boundary element method. It is shown that in the initial stage of rotation of the cylinder four different variants of the behavior of the free surface can be realized in the stage of transition from horizontal to steady-state form.     

Experimental and high-order LES analysis of the flow in near-wall region of a square cylinder

A coupled experimental/numerical analysis of turbulent flow past a square cylinder is performed at the ERCOFTAC Reynolds number Re = U_∞D/ν = 21,400, where U_∞ is the inflow velocity and D the cylinder height. Complementary Laser Doppler Velocimetry (LDV) and high-order large-eddy simulations (LES) approaches, based on a spectral vanishing technique (SVV-LES), provide a comprehensive data base including both instantaneous data and post-processed statistics. Beyond these results, an achievement of the paper is to investigate the coherent structures developing on the sides and in the wake of the cylinder with a special focus on the flow features in the near-wall region. The flow is found to separate at the leading edge of the cylinder with the occurence of three-dimensional Kelvin-Helmholtz (KH) pairings localized in the separating shear layer. The interaction between these KH vortical structures and Von Kármán vortex shedding (VK) in the near wake is discussed based on both visualisations and frequency analysis. In particular, signatures of VK and KH vortical structures are found on velocity time samples.     

Induced parallel vortex shedding from a circular cylinder at Re of O(10) by using the cylinder end suction technique

In the present work, the objective is to attempt to induce parallel vortex shedding at a moderately high Reynolds number (=1.578 × 10⁴) by using the cylinder end suction method, and measure the associated aerodynamic parameters.We first measured the aerodynamic parameters of a single circular cylinder without end suction, and showed that the quantities measured are in good agreement with equivalent data in the published literature. Next, by using different amount of end suction which resulted in increasing the cylinder end velocity by 1%, 2% and 2.5%, we were able to show that the above corresponded to the situation of under suction, optimal suction and over suction, respectively. With optimal suction, we demonstrated that the end suction method works at Re = 1.578 × 10⁴. The shape of the primary vortex shed became straighter than when there is no end suction, and parameters like cylinder surface pressure distribution, drag force per unit span, as well as vortex shedding frequency all showed negligible spanwise variation. Further careful analyses showed that when compared to the naturally existing curved vortex shedding, with parallel vortex shedding the mid-span drag per unit span became slightly smaller, but the drag averaged over the cylinder span became slightly larger. For cylinder surface pressure, it was found that cylinder end effects mainly influenced the surface pressure in the angular ranges −180°  β < −60° and 60° < β  180°. Without end suction, the cylinder surface pressure in the above ranges was found to increase (become less negative) slightly with |z/d|, but such increase disappeared when optimal end suction was applied, and the cylinder surface pressure distribution became spanwise location independent. As for the vortex shedding frequency (Strouhal number), although the Strouhal number showed spanwise variation when there is no end suction and negligible spanwise variation when optimal suction was applied, the difference between the spanwise averaged Strouhal number was quite negligible. With under suction, the spanwise dependence of various aerodynamic parameters existed, but was found to be not as significant as when no end suction was applied at all. With over suction, the flow situation was found to be practically no change from the optimal suction situation.     

An experimental study of the dynamic characteristics in boundary layer in the flow past a two-dimensional circular cylinder at subcritical reynolds numbers

An experimental study of the dynamic characteristics of flow past a two-dimensional circular cylinder is described. The fluctuationsoof wall shear stress, surface-pressure and velocity of the flow are measured with hot-film, hot-wire and pressure transducer. The frequency feature of fluctuations of wall shear stress is given. The cross-correlation functions of these fluctuations at any two points are calculated. The experimental results reveal that there is an overall syncronous fluctuation, at the shedding frequency, in boundary layer in the flow past a two-dimensional circular cylinder at subcritical Reynolds number.     

An experimental study on the effects of uniform injection through one perforated surface of a square cylinder on some aerodynamic parameters

It is well known that injection/suction (transpiration) through a perforated surface is an efficient way of influencing the characteristics of a turbulent boundary layer. Injection application creates a thicker boundary layer on a flat plate and it thus decreases drag. In aeronautical applications, suction is frequently used to delay boundary layer separation. This paper presents an experimental study on the effects of uniform injection through one perforated surface of a square cylinder on the pressure distribution and drag coefficient in a two-dimensional turbulent flow. For this purpose, surface pressure measurements around a square cylinder have been performed at three different Reynolds numbers in a wind tunnel. The parameters taken into account were injection rate, position of perforated surface (i.e., front, top, and rear), and pressure coefficient and drag coefficient. The results show that variation in pressure coefficient around the square cylinder and drag coefficient were influenced by the position of perforated surface and by injection rate.     

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

Numerical and experimental study of the fine structure of a stratified fluid flow over a circular cylinder

The fine structure of the flow field of a continuously stratified fluid around a circular cylinder for small values of the Froude number was investigated in laboratory and numerical experiments. The parameters of the leading perturbation, the internal-wave field, and the cylinder wake were calculated using a two-dimensional model. The existence of the previously experimentally observed high-gradient density layers in the wake that are parallel to the flow axis was for the first time confirmed by numerical calculations. Results of the numerical and experimental studies are in good agreement with each other and with analytical models for small values of the Froude number. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 43–54, January–February, 2007.     

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.     

Mixed convective heat and mass transfer in a non-Newtonian fluid at a peristaltic surface with temperature-dependent viscosity

We investigate the problem of the unsteady mixed convection peristaltic mechanism. The flow includes a temperature-dependent viscosity with thermal diffusion and diffusion-thermo effects. The peristaltic flow is between two vertical walls, one of which is deformed in the shape of traveling transversal waves exactly like peristaltic pumping and the other of which is a parallel flat plate wall. The equations of momentum, energy, and concentration are subject to a set of appropriate boundary conditions by assuming that the solution consists of two parts: a mean part and a perturbed part. The solution of the perturbed part has been obtained by using the long-wave approximation. The mean part has been solved and coincides with the approximation of Ostrach. The mean part (zeroth order), the first order, and the total solution of the problem have been evaluated numerically for several sets of values of the parameters entering the problem. The skin friction, and the rate of heat and mass transfer at the walls are obtained and illustrated graphically.     

Local flows of a deformable damageable medium under impact interaction with a cavitating fluid

Local viscoplastic-flow and damage processes in a deformable medium induced by the collapse of dispersed vapor-gas bubbles in the near-wall layer of a cavitating fluid in the presence of propagating shock waves are investigated. The study is based on a generalized model developed for describing nonlinear deformations and flow of damageable media and on the results obtained earlier for local fluid flows induced by limiting transitions of vapor-gas bubbles.     

Interaction of a heavy fluid flow in a channel with a transverse jet barrier

An experimental and numerical analysis of the interaction between a plane horizontal water flow in a rectangular channel (free water current) and a plane thin water jet (water jet curtain) is presented; the jet flows out vertically from either a slot nozzle in the bottom of the channel or the crest of a rigid spillway at a velocity appreciably (several times) greater than the water velocity in the channel. Numerical calculations were carried out using the STAR-CD software package preliminarily tested against the experimental data obtained. The dependence of the water level in the channel at a certain distance ahead of the jet barrier on the main jet parameters and the water flow rate in the horizontal channel is studied. It is found that in the region of the interface between the flows both steady and unsteady (self-oscillatory) flow patterns can be realized. Steady stream/jet interaction patterns of the “ejection” and “ejection-spillway” types are distinguished and a criterion separating these regimes is obtained. The notion of a rigid spillway equivalent to a jet curtain is introduced and an approximate dependence of its height on the relevant parameters of the problem is derived. The possibility of effectively controlling the water level ahead of a rigid spillway with a sharp edge by means of a plane water jet flowing from its crest is investigated. The boundary of transition to self-oscillation interaction patterns in the region of the flow interface is determined. The structure of these flows and a possible mechanism of their generation are described. Within the framework of the inviscid incompressible fluid model in the approximate formulation for a “thin” jet, an analytical dependence of the greatest possible depth of a reservoir filled with a heavy fluid at rest and screened by a vertical jet barrier on the jet parameters is obtained.     

On the influence of numerical schemes and subgrid–stress models on large eddy simulation of turbulent flow past a square cylinder

Influence of finite difference schemes and subgrid‐stress models on the large eddy simulation calculation of turbulent flow around a bluff body of square cylinder at a laboratory Reynolds number, has been examined. It is found that the type and the order of accuracy of finite‐difference schemes and the subgrid‐stress model for satisfactory results are dependent on each other, and the grid resolution and the Reynolds number. Using computational grids manageable by workstation‐level computers, with which the near‐wall region of the separating boundary layer cannot be resolved, central‐difference schemes of realistic orders of accuracy, either fully conservative or non‐conservative, suffer stability problems. The upwind‐biased schemes of third order and the Smagorinsky eddy‐viscosity subgrid model can give reasonable results resolving much of the energy‐containing turbulent eddies in the boundary layers and in the wake and representing the subgrid stresses in most parts of the flow. Noticeable improvements can be obtained by either using higher order difference schemes, increasing the grid resolution and/or by implementing a dynamic subgrid stress model, but each at a cost of increased computational time. For further improvements, the very small‐scale eddies near the upstream corners and in the laminar sublayers need to be resolved but would require a substantially larger number of grid points that are out of the range of easily accessible computers. Copyright © 2002 John Wiley & Sons, Ltd.     

A numerical study of the vortex motion in oscillating flow around a circular cylinder at low and middlieKc numbers

A new hybrid model, which is based on domain decomposition and proposed by the authors, is used for calculating the flow around a circular cylinder at low and middle Keulegan-Carpenter numbers (Kc=2−18) respectively. The vortex motion patterns in asymmetric regime, single pair (or transverse) regime and double pair (or diagonal) regime are successfully simulated. The calculated drag and inertial force coefficients are in better agreement with experimental data than other recent computational results. The project supported by the National Natural Science Foundations of China and the LNM, Institute of Mechanics, Academia Sinica