Instantaneous flow field above the free end of finite-height cylinders and prisms

The flow above the free ends of surface-mounted finite-height circular cylinders and square prisms was studied experimentally using particle image velocimetry (PIV). Cylinders and prisms with aspect ratios of AR = 9, 7, 5, and 3 were tested at a Reynolds number of Re = 4.2 × 10⁴. The bodies were mounted normal to a ground plane and were partially immersed in a turbulent zero-pressure-gradient boundary layer, where the boundary layer thickness relative to the body width was δ/D = 1.6. PIV measurements were made above the free ends of the bodies in a vertical plane aligned with the flow centreline. The present PIV results provide insight into the effects of aspect ratio and body shape on the instantaneous flow field. The recirculation zone under the separated shear layer is larger for the square prism of AR = 3 compared to the more slender prism of AR = 9. Also, for a square prism with low aspect ratio (AR = 3), the influence of the reverse flow over the free end surface becomes more significant compared to that for a higher aspect ratio (AR = 9). For the circular cylinder, a cross-stream vortex forms within the recirculation zone. As the aspect ratio of the cylinder decreases, the reattachment point of the separated flow on the free end surface moves closer to the trailing edge. For both the square prism and circular cylinder cases, the instantaneous velocity vector field and associated in-plane vorticity field revealed small-scale structures mostly generated by the separated shear layer.     

Influence of aspect ratio on the mean flow field of a surface-mounted finite-height square prism

The flow around surface-mounted, finite-height square prisms at a Reynolds number of Re = 4.2 × 10⁴ was investigated experimentally in a low-speed wind tunnel using particle image velocimetry. The thickness of the boundary layer on the ground plane relative to the width of the prism was δ/D = 1.5. Four prism aspect ratios were tested, AR = 9, 7, 5, and 3, to study how the aspect ratio influences the flow field close to the prism. Upstream of the prism, lowering the aspect ratio from AR = 9 to AR = 3 causes the stagnation point on the upstream face to move closer to the free end, but there is no influence on the location and strength of the horseshoe vortex. Lowering the aspect ratio from AR = 9 to AR = 3 causes the cross-stream vortices in the upper and lower halves of the wake to move downstream and upstream, respectively; the latter vortex is absent for AR = 3, suggesting this prism sits below the critical aspect ratio. Above the free end of the prism, within the region of separated flow, lowering the aspect ratio from AR = 9 to AR = 3 shifts the location of the cross-stream vortex farther downstream. For the prism of AR = 3, reverse flow above the free end is stronger yet more unsteady compared to the more slender prisms, while the streamwise edge vortices are smaller and weaker.     

The effect of a splitter plate on the flow around a finite prism

The effect of a wake-mounted splitter plate on the flow around a surface-mounted finite-height square prism was investigated experimentally in a low-speed wind tunnel. Measurements of the mean drag force and vortex shedding frequency were made at Re=7.4×10⁴ for square prisms of aspect ratios AR=9, 7, 5 and 3. Measurements of the mean wake velocity field were made with a seven-hole pressure probe at Re=3.7×10⁴ for square prisms of AR=9 and 5. The relative thickness of the boundary layer on the ground plane was δ/D=1.5–1.6 (where D is the side length of the prism). The splitter plates were mounted vertically from the ground plane on the wake centreline, with a negligible gap between the leading edge of the plate and rear of the prism. The splitter plate heights were always the same as the heights of prisms, while the splitter plate lengths ranged from L/D=1 to 7. Compared to previously published results for an “infinite” square prism, a splitter plate is less effective at drag reduction, but more effective at vortex shedding suppression, when used with a finite-height square prism. Significant reduction in drag was realized only for short prisms (of AR≤5) when long splitter plates (of L/D≥5) were used. In contrast, a splitter plate of length L/D=3 was sufficient to suppress vortex shedding for all aspect ratios tested. Compared to previous results for finite-height circular cylinders, finite-height square prisms typically need longer splitter plates for vortex shedding suppression. The effect of the splitter plate on the mean wake was to narrow the wake width close to the ground plane, stretch and weaken the streamwise vortex structures, and increase the lateral entrainment of ambient fluid towards the wake centreline. The splitter plate has little effect on the mean downwash. Long splitter plates resulted in the formation of additional streamwise vortex structures in the upper part of the wake.     

The effect of incidence angle on the mean wake of surface-mounted finite-height square prisms

The mean wake of a surface-mounted finite-height square prism was studied experimentally in a low-speed wind tunnel to explore the combined effects of incidence angle (α) and aspect ratio (AR). Measurements of the mean wake velocity field were made with a seven-hole pressure probe for finite square prisms of AR = 9, 7, 5 and 3, at a Reynolds number of Re = 3.7 × 10⁴, for incidence angles from α = 0° to 45°. The relative thickness of the boundary layer on the ground plane, compared to the prism width, was δ/D = 1.5. As the incidence angle increases from α = 0° to 15°, the mean recirculation zone shortens and the mean wake shifts in the direction opposite to that of the mean lift force. The downwash is also deflected to this side of the wake and the mean streamwise vortex structures in the upper part of the wake become strongly asymmetric. The shortest mean recirculation zone, and the greatest asymmetry in the mean wake, is found at the critical incidence angle of α_critical ≈ 15°. As the incidence angle increases from α = 15° to 45°, the mean recirculation zone lengthens and the mean streamwise vortex structures regain their symmetry. These vortices also elongate in the wall-normal direction and become contiguous with the horseshoe vortex trailing arms. The mean wake of the prism of AR = 3 has some differences, such as an absence of induced streamwise vorticity near the ground plane, which support its classification as lying below the critical aspect ratio for the present flow conditions.     

Early stages of an impulsively started unsteady laminar flow past tapered trapezoidal cylinders

Characteristics of the developing recirculation region behind a tapered trapezoidal cylinder and its interaction with the separating shear layer from the leading edges were studied numerically for an impulsively started laminar flow. An unsteady stream function–vorticity formulation was used. The Reynolds numbers considered range from 25 to 1000. Pressure contours, surface pressure coefficient, wake length and drag coefficient were studied through the streamline flow field. Main flow and subflow regimes were identified by an analysis of the evolution of the flow characteristics. It was found that typically, for a given trapezoidal cylinder, flow starts with no separation. As time advances, the symmetrical standing zone of recirculation develops aft of the trapezoidal cylinder. The rate of growth in width, length and structure of the aft end eddies depends on the Reynolds number. In time, separated flow from the leading edges of the trapezoidal cylinder also develops and forms growing separation bubbles on the upper and lower inclined surfaces of the trapezoidal cylinder. As time advances, the separation bubbles on the upper and lower inclined surfaces of the cylinder grow towards the downstream regions and eventually merge with the swelling symmetrical eddies aft of the cylinder. This merging of the flows creates a complex flow regime with a disturbed tertiary flow zone near the merging junction. Eventually, depending on the Reynolds number and the tapered angle of the trapezoidal cylinder, the flow develops into a specific category of symmetrical standing recirculatory flow with its own distinct characteristics. Comparisons with the available results of other investigators showed very good agreement. © 1998 John Wiley & Sons, Ltd.     

EFFECTS OF LID-DRIVEN CAVITY SHAPE ON THE FLOW ESTABLISHMENT PHASE

Experiments are carried out to study the flow establishment phase inside closed cavities submitted to the impulsive translation from rest, of one of their walls at a Reynolds number of 1000. Three standard industrially machined or molded cylindrical cavity shapes are studied and are compared with respect to the efficiency of mixing process: square, rectangular and semi-circular of length-to-width ratio of 2:1. The flow structures in the mid-cross-section are analysed by means of fine topological and kinematic visualization series using two complementary techniques: continuous dye filament and discrete solid tracers both coupled with a laser sheet illumination. Particular attention is given to vorticity propagation and primary/secondary eddy formations. Although a roughly similar vortex generation is observed in all examined cavities, important differences appear with time. The semi-circular cavity flow results in a much more homogeneous and uniform recirculation with no secondary flow recirculation zone. On the contrary, the square and rectangular cavity flows develop a better flow mass dispersion and, respectively, one and two secondary eddies. At the final time of observation (t^*=12), both semi-circular and rectangular cavity flows seem to reach their steady state whereas the square one continues to evolve. Comparisons with 2-D computational results of other authors illustrate the three-dimensional flow aspect present in experiments.     

TANDEM CYLINDERS IN IMPULSIVELY STARTED FLOW

The impulsively started flow field for circular cylinders of equal diameter arranged in tandem was investigated using flow visualization and particle image velocimetry (PIV), over a longitudinal pitch ratio range ofL /D=1·0–3·0, and for Reynolds numbers from Re=1200–3800. The PIV technique was used to obtain a time history of the instantaneous in-plane vorticity field from the moment of impulsive start, from which the spatial and temporal development of the flow was studied. Measurements of vortex strength and vortex position relative to the cylinders were obtained from these data. Three types of fluid behaviour were identified based on L/D: single bluff-body behaviour when the cylinders are in contact, constrained streamwise growth and lateral expansion of the gap recirculation zones at small and intermediate L/D, and independent formation of recirculation zones similar to a single impulsively started circular cylinder at larger L/D.     

NUMERICAL PREDICTION OF FORCE ON RECTANGULAR CYLINDERS IN OSCILLATING VISCOUS FLOW

Numerical results are presented for an oscillating viscous flow past a square cylinder with square and rounded corners and a diamond cylinder with square corners at Keulegan–Carpenter numbers up to 5. This unsteady flow problem is formulated by the two-dimensional Navier–Stokes equations in vorticity and stream-function form on body-fitted coordinates and solved by a finite-difference method. Second-order Adams-Bashforth and central-difference schemes are used to discretize the vorticity transport equation while a third-order upwinding scheme is incorporated to represent the nonlinear convective terms. Since the vorticity distribution has a mathematical singularity at a sharp corner and since the force coefficients are found in experiments to be sensitive to the corner radius of rectangular cylinders, a grid-generation technique is applied to provide an efficient mesh system for this complex flow. Local grid concentration near the sharp corners, instead of any artificial treatment of the sharp corners being introduced, is used in order to obtain high numerical resolution. The elliptic partial differential equation for stream function and vorticity in the transformed plane is solved by a multigrid iteration method. For an oscillating flow past a rectangular cylinder, vortex detachment occurs at irregular high frequency modes at KC numbers larger than 3 for a square cylinder, larger than 1 for a diamond cylinder and larger than 3 for a square cylinder with rounded corners. The calculated drag and inertia coefficients are in very good agreement with the experimental data. The calculated vortex patterns are used to explain some of the force coefficient behavior.     

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.     

Particle response in a planar sudden expansion flow

Particle velocity and concentration statistics were measured in a vertically downward planar sudden expansion flow for large-eddy particle Stokes numbers (τ_pU_o/5H) ranging from 0.5 to 7.4. Particles with Stokes numbers greater than 3 did not enter the recirculation zone, exhibited substantial attenuation of cross-stream velocity fluctuations, and had large streamwise velocity fluctuations in regions of strong velocity gradient. The smallest particles filled the recirculation zone and showed strong response to the large eddies in the flow. Phase-locked particle concentration measurements showed that these particles were centrifuged away from vortex cores and concentrated between vortices. Intermediate-size particles with Stokes numbers of 1.4 were injected intermittently into the recirculation zone as tongues of particles moving down between vortices. Particle Reynolds number was found to have negligible effect on the particle velocity statistics.     

Time-averaged topological flow patterns and their influence on vortex shedding of a square cylinder in crossflow at incidence

Flow characteristics around the square cylinder and their influence on the wake properties are studied. Time-averaged flow patterns on the surfaces of square cylinder in a cross-stream at incidence are experimentally probed by surface-oil flow technique and analyzed by flow topology for Reynolds numbers between 3.9×10⁴ and 9.4×10⁴ as the incidence angle changes from 0° to 45°. Vortex shedding characteristics are measured by a single-wire hot-wire anemometer for Reynolds numbers between 5×10³ and 1.2×10⁵. The effects of topological flow patterns on the wake properties then are revealed and discussed. Flows around the square cylinder are identified as three categories: the subcritical, supercritical, and wedge flows according to the prominently different features of the topological flow patterns. The Strouhal number of vortex shedding, turbulence in the wake, and wake width present drastically different behaviors in different characteristic flow regimes. A critical incidence angle of 15° separates the subcritical and supercritical regimes. At the critical incidence angle the wake width and shear-layer turbulence present minimum values. The minimum wake width appearing at the critical incidence angle, which leads to the maximum Strouhal number, is due to the reattachment of one of the separated boundary layer to the lateral face of the square cylinder. If the Strouhal numbers are calculated based on the wake width instead of the cross-stream projection width of cylinder, the data in the subcritical and supercritical regimes are well correlated into two groups, which would approach constants at high Reynolds numbers.     

Aerodynamic forces and three-dimensional flow structures in the mean wake of a surface-mounted finite-height square prism

Different flow models have been proposed for the flow around surface-mounted finite-height square prisms, but there is still a lack of consensus about the origin and connection of the streamwise tip vortices with the other elements of the wake. This numerical study was performed to address this gap, in addition to clarifying the relationship of the near-wake structures with the far wake and the near-wall flow, which is associated with the fluid forces. A large-eddy simulation approach was adopted to solve the flow around a surface-mounted finite-height square prism with an aspect ratio of AR = 3 and a Reynolds number Re = 500. The mean drag and normal forces and the bending moment for the prism were quantitatively compared in terms of skin-friction and pressure contributions, and related to the near-wall flow. Both three-dimensional visualizations and planar projections of the time-averaged flow field were used to identify, qualitatively, the main structures of the wake, including the horseshoe vortex, corner vortices and regions of high streamwise vorticity in the upper part of the wake. These features showed the same qualitative behavior as reported in high Reynolds number studies. It was found that some regions of high streamwise vorticity magnitude, like the tip vortices, are associated with the three-dimensional bending of the flow, and the tip vortices did not continuously extend to the free end of the prism. The three-dimensional flow analysis, which integrated different observations of the flow field around surface-mounted finite-height square prisms, also revealed that the mean near-wake structure is composed of two sections of different origin and location of dominance.     

Numerical study of vortex shedding from a rotating cylinder immersed in a uniform flow field

A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=10³ and α≤3; the other with Re=10⁴ and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd.     

Influence of aspect ratio on the flow above the free end of a surface-mounted finite cylinder

The flow above the free end of a surface-mounted finite-height cylinder was studied in a low-speed wind tunnel using particle image velocimetry (PIV). Velocity measurements were made in vertical and horizontal measurement planes above the free end of finite cylinders of aspect ratios AR = 9, 7, 5 and 3, at a Reynolds number of Re = 4.2 × 10⁴. The relative thickness of the boundary layer on the ground plane was δ/D = 1.7. Flow separating from the leading edge formed a prominent recirculation zone on the free-end surface. The legs of the mean arch vortex contained within the recirculation zone terminate on the free-end surface on either side of the centreline. Separated flow from the leading edge attaches onto the upper surface of the cylinder along a prominent attachment line. Local separation downstream of the leading edge is also induced by the reverse flow and arch vortex circulation within the recirculation zone. As the cylinder aspect ratio is lowered from AR = 9 to AR = 3, the thickness of the recirculation zone increases, the arch vortex centre moves downstream and higher above the free-end surface, the attachment position moves downstream, and the termination points of the arch vortex move upstream. A lowering of the aspect ratio therefore results in accentuated curvature of the arch vortex line. Changes in aspect ratio also influence the vorticity generation in the near-wake region and the shape of the attachment line.     

A quantitative study of the flow around an impulsively started circular cylinder

The detailed flow structure behind an impulsively started circular cylinder has been investigated experimentally. The Reynolds number based on the steady state velocity and the diameter of the cylinder was 500 to 3,000. This work is unique in that unsteady spatial velocities were measured simultaneously by a quantitative visualization technique — Laser Induced Photochemical Anemometry (LIPA). The surface vorticity at g/q = π/2 and vorticity distribution behind the cylinder in the Lagrangian coordinates (i.e. coordinates fixed on the cylinder) were calculated from the measured velocities. The surface vorticity shows in the early stage of flow development a close agreement with the previous results obtained by analytical and numerical approaches. The large-field velocity and vorticity information provides an insight into the formation process of the vortices downstream of the cylinder. In addition to the quantitative information, the results of visualized flow pattern obtained by LIPA technique are also presented. A preliminary version of this paper was presented at the Twelfth Symposium on Turbulence, University of Missouri-Rolla, Sept. 24–26, 1990     

A mixed algorithm for the calculation of rapidly varying fluid flows: The impulsively started circular cylinder

A mixed discrete Fourier transform-Finite difference algorithm is developed and used for the calculation of rapidly changing viscous fluid flows past a circular cylinder. The numerical approach has been designed to overcome certain difficulties arising for high Reynolds number simulations. The foremost advantage of the technique lies in its fast calculations of the convolution sums portraying the convective terms of the governing equations. Third-order spatial discretizations and fourth-order time marching are implemented. New schemes are proposed for the boundary conditions at the solid wall and at large distances. The techniques are tested on a case study with other schemes (summarized by Roache¹) in order to obtain an optimal choice. Definite indications on the stability and accuracy of boundary condition schemes are achieved. Support for the statement of dominant importance of boundary conditions is also given. A comparison of computational results with experimental data is presented for the case study of the flow past an impulsively started cylinder at Reynolds number 20. The time development of the symmetrical zone of recirculation, which is formed at an early stage of the flow, has been studied for 300 ≤ Re ≤ 9500 by means of the proposed algorithm. Computational results, comparisons with experimental data² and discussion of upper limits of validity of the procedure will be presented in a companion paper.     

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.     

Transition phenomena in the wake of an inclined square cylinder

The division of flow regimes in a square cylinder wake at various angles of attack (α) is studied. This study provides evidence of the existence of modes A and B instabilities in the wake of an inclined square cylinder. The critical Reynolds numbers for the inception of these instability modes were identified through the determination of discontinuities in the Strouhal number versus Reynolds number curves. The spectra and time traces of wake streamwise velocity were observed to display three distinct patterns in different flow regimes. Streamwise vortices with different wavelengths at various Reynolds numbers were visualized. A PIV technique was employed to quantitatively measure the parameters of wake vortices. The wavelengths of the streamwise vortices in the modes A and B regimes were measured by using the auto-correlation method. From the present investigation, the square cylinder wake at various angles of attack undergoes a similar transition path to that of a circular cylinder, although various quantitative parameters measured which include the critical Reynolds numbers, spanwise wavelength of secondary vortices, and the circulation and vorticity of wake vortices all show an α dependence.     

Cross-stream forces and velocities of fixed and freely suspended particles in viscoelastic Poiseuille flow: Perturbation and numerical analyses

The cross-stream migration of a circular particles (or infinitely long cylinder) in two dimensional, inertia-less viscoelastic pressure-driven flows is examined through complementary finite element simulations and second-order fluid perturbation analyses for small Deborah number (De), where De is defined as the fluid relaxation time divided by the characteristic flow time. A neutrally buoyant, freely suspended particle migrates toward the center of the channel for all particle sizes and cross-stream positions due to the coupled effects of the linear and quadratic variations of the imposed velocity. A particle that is held at a fixed position, in contrast, experiences a cross-stream force directed toward the wall as a result of the coupled effects of the local shear flow and the flow relative to the particle.     

NEAR-WAKE FLOW STRUCTURE OF A CYLINDER WITH A HELICAL SURFACE PERTURBATION

The near-wake of a circular cylinder having a helical wire pattern about its surface is characterized using a technique of high-image-density velocimetry. Patterns of vorticity in three orthogonal planes show substantial influence of a wire having a diameter an order of magnitude smaller than the cylinder diameter. The distinctive patterns of vorticity in these three planes are associated with lack of formation of large-scale Kármán-like clusters of vorticity (ω_z) in the near-wake region of the cylinder. The instantaneous structure of the separating spanwise vorticity (ω_z) layers on either side of the cylinder involve small-scale concentrations of vorticity analogous to the well-known Kelvin–Helmholtz vortices from a smooth cylinder. Moreover, a dual vorticity layer, i.e., two adjacent layers of like vorticity (ω_z), can form from one side of the cylinder. Along the span of the cylinder, distributions of instantaneous velocity and transverse vorticity (ω_y) show a spatially periodic sequence of wake-like patterns, each of which has features in common with the very near-wake of a two-dimensional bluff body, including a large velocity defect bounded by vorticity layers with embedded small-scale vorticity concentrations. In the cross-flow plane of the wake, patterns of streamwise vorticity (ω_x) show small-scale, counter-rotating pairs of vorticity concentrations (ω_x) emanating from the inclined helical perturbation, rather than isolated concentrations of vorticity of like sign, which would indicate single streamwise vortices. All of the aforementioned patterns of small-scale vorticity concentrations are scaled according to the local wake width/local pitch of the helical wire pattern in the respective plane of observation.