Influence of rounding corners on the wake of a finite-length cylinder: An experimental study

This study aims to investigate experimentally the influence of rounding corners (r) as well as aspect ratio (AR) on the flow structures of a surface-mounted finite cylinder. The cylinders with sharp (r* = r/D = 0) and rounded corners (r*=0.167, 0.25 and 0.5) and aspect ratio or height-to-width/diameter ratio (AR = H/D) between 2 and 7 are utilized. The experiments are based on the five-hole probe and hot-wire measurements as well as the oil flow visualization. Wake measurements are made in an open return wind tunnel at the Reynolds number, Re = 1.6 × 10⁴, where Re is defined based on the side width/diameter (D) of the cylinder cross-section and the freestream velocity. It is found that r* and AR have significant effects on the flow structure from the perspective of wake topology, strength of streamwise vortices, and vortex shedding frequency. For all r* considered, the wake is characterized by a quadrupole type (both the tip and base vortices are present) at AR = 7, while a dipole type occurs for AR = 2 and 4 (the base vortices are absent). The strength (circulation) of the streamwise vortex structures is affected by r*. For all AR examined in the present study, the strengths of tip and base vortex structures decrease with increasing r*. The oil flow visualization demonstrates that the features of the horseshoe vortex are sensitive to r* and AR. With increasing r*, the location of the separation line moves downstream and the distance between horseshoe vortex legs decreases. Velocity measurements reveal that the downwash flow enhances with increasing r*. It is also found that the Strouhal number increases progressively by 60% as r* increases from 0 to 0.5, regardless of AR.     

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.     

Flow structures and force characteristics for flat plate in oscillatory flows withKc number from 2 to 40 and in combined flows

The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, wyenKc number varies from 2 to 40, the vortex pattern changes from a “harmonic wave” shaped (in a range of smallKc number) to a slight inclined “harmonic wave” shaped (in a range of moderateKc numbers), then to inclined vortex clusters with an angle of 50° to the oncoming flow direction (atKc=20), at last, asKc number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other authors. The existence of minimum point of inertia force coefficientC _m nearKc=20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a “longitudinal wave” structure, and a vortex cluster shape corresponding to the ratios ofU _m toU ₀ which are ofO (10⁻¹)O(1) andO(10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demonstrated qualitatively. In addition to this, the lock in phenomenon of vortex shedding has been checked. The project supported by National Natural Science Foundation of China & LNM, Institute of Mechanics, CAS     

On the formation of vortex rings and pairs

The axisymmetric vortex sheet model developed by Nitsche & Krasny (1994) has been extended to study the formation of vortex rings (pairs) at the edge of circular (2D) tube and opening. Computations based on this model are in good agreement with the experiments (Didden (1979) for circular tube and Auerbach (1987) for 2D tube and opening). Using this new model, evidences are provided to show that the main failure of the similarity theory (the false prediction of axial trajectory of vortex ring) is due to its ignorance of the self-induced ring velocity (mutual induction for vortex pair). We further reason why the similarity theory succeeds in its prediction of radial movement of vortex ring. The effects of various parameters such as turning angle α and piston speedU _p (t) on the formation of vortex ring are investigated. Numerical result shows that turning angle α has no effect on circulation shed τ. We also discuss Glezer (1988)'s summary on the influence ofU _p upon the shedding circulation, and finally give the variation of core distribution of vortex ring with α andU _p (t). The project is supported by National Natural Science Foundation of China and Doctoral Program of Institution of Higher Education     

A Lagrangian vortex method for unbounded flows

A technique is presented for velocity calculations on the highly distorted node distributions typical of those found in Lagrangian vortex methods. The method solves the partial differential equation for streamfunction directly on the nodes, via a sparse, symmetric system of equations that can be solved using standard iterative solvers. When implemented in a triangulated vortex method, the technique gives computation times which scale as N^1.23, where N is the number of nodes. The computation scheme is derived for two‐dimensional problems and applied to the prediction of the evolution of perturbed multipolar vortices. Due to the numerical performance of the method, it has been possible to examine such evolution at higher and lower Reynolds numbers than have been considered in published numerical studies. Copyright © 2008 John Wiley & Sons, Ltd.     

Rankine combined vortex interaction with a rectangular prism

Large eddy simulation is utilised to study the three-dimensional interaction between a travelling Rankine combined vortex and a rectangular prism. The study examines the strength and the topology of a vortex during the interaction with a prism that is much wider than the vortex core diameter. The physics of the interaction is revealed for the straight (β = 0°) and the oblique (β = 45°) impacts. For both cases, the low-level portion of the vortex undergoes displacements in the streamwise and the lateral directions. Also the vortex shape and the core vorticity are substantially disrupted. Behind the prism the full vortex circulation is recovered after a considerable distance. This created a low-velocity region. The sheltering effect of the prism is noticed for both straight and oblique impacts. The flow velocities in the sheltering region, right behind the prism, are reduced by more than 42% compared to the maximum flow speeds before the interaction.     

The effects of corrugation and wing planform on the aerodynamic force production of sweeping model insect wings

The effects of corrugation and wing planform (shape and aspect ratio) on the aerodynamic force production of model insect wings in sweeping (rotating after an initial start) motion at Reynolds number 200 and 3500 at angle of attack 40° are investigated, using the method of computational fluid dynamics. A representative wing corrugation is considered. Wing-shape and aspect ratio (AR) of ten representative insect wings are considered; they are the wings of fruit fly, cranefly, dronefly, hoverfly, ladybird, bumblebee, honeybee, lacewing (forewing), hawkmoth and dragonfly (forewing), respectively (AR of these wings varies greatly, from 2.84 to 5.45). The following facts are shown. (1) The corrugated and flat-plate wings produce approximately the same aerodynamic forces. This is because for a sweeping wing at large angle of attack, the length scale of the corrugation is much smaller than the size of the separated flow region or the size of the leading edge vortex (LEV). (2) The variation in wing shape can have considerable effects on the aerodynamic force; but it has only minor effects on the force coefficients when the velocity at r ₂ (the radius of the second moment of wing area) is used as the reference velocity; i.e. the force coefficients are almost unaffected by the variation in wing shape. (3) The effects of AR are remarkably small: when AR increases from 2.8 to 5.5, the force coefficients vary only slightly; flowfield results show that when AR is relatively large, the part of the LEV on the outer part of the wings sheds during the sweeping motion. As AR is increased, on one hand, the force coefficients will be increased due to the reduction of 3-dimensional flow effects; on the other hand, they will be decreased due to the shedding of part of the LEV; these two effects approximately cancel each other, resulting in only minor change of the force coefficients. The project supported by the National Natural Science Foundation of China (10232010 and 10472008) and Ph. D. Student Foundation of Chinese Ministry of Education (20030006022) The English text was polished by Keren Wang.     

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.     

Experimental study on the wake behind tapered circular cylinders

The flow around tapered cylinders can act as basic models for numerous bluff body flows with a spanwise variation of either the body shape or the inflow conditions. The well-known vortex street is influenced by strong three-dimensional effects from the spanwise variation of the shedding frequency, namely oblique vortex shedding and vortex dislocations. Stereo-PIV was chosen to study these phenomena, since it allows analyzing planes with the full three-component, instantaneous velocity fields and local, time-dependent variations in the same setting. Hence, detailed aspects of the vortex dislocation phenomenon are presented. Single vortex dislocation events are presented through the local variation of the three measured velocity components u, v and w. Longer time-series reveal both period and location of these dislocation events, as well as quantity and sizes of the cells of constant shedding velocity in between them. The influence of the Reynolds number and the cylinder aspect ratio on the vortex cells could be shown. The analysis of the vortex shedding behavior shows good agreement with previously published results. At the same time, the applied PIV technique provides more spatial information than point-based measurements and offers insight into a Reynolds number range that is currently out of reach of Direct Numerical Simulations.     

Suppression of force fluctuations in flow past an aerofoil

Force fluctuations on a solid body are associated with unsteadiness in the wake, e.g. vortex shedding. Therefore, the control of force fluctuations can be realised by suppressing the flow unsteadiness. A NACA0024 aerofoil closed with a round trailing edge is chosen to represent the solid body throughout this investigation, with the Reynolds number fixed at Re = 1000 and angle of attack α ≤ 15^o, at which the uncontrolled flow is two-dimensional. A linear optimal control is calculated by analysing the distribution of sensitivity of unsteadiness to control around the entire surface of the body. The nonlinear effects of the calculated control, which can be actuated through surface-normal suction and blowing across the surface of the aerofoil, are tested through two-dimensional direct numerical simulations. It is observed that a surface-normal velocity control with a maximum magnitude less than 8% of the free stream velocity completely suppresses unsteadiness at α = 10° with an overall drag reduction of 14% and a 138% increase of lift.     

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.     

Investigations of wake flows of a flat plate in steady,oscillatory and combined flows

Numerical study on near wake flows of a flat plate in three kinds of oncoming flows is made by using the discrete vortex model and improved vorticity creation method. For steady oncoming flow, both gross and detailed features of the wake flow are calculated and discussed. Then, in harmonic oscillatory oncoming flow two different wake flow patterns withK _c=2,4 and 10 are obtained respectively. Our results present a new wake flow pattern for lowKc numbers (Kc<5) describing vortex shedding, pairing and moving in a period of the oscillatory flow starting from rest. The calculated drag and inertia force coefficients are closer to experimental data from the U-tube than the previous results of vortex simulation. For in-line combined oncoming flow the vortex lock-in and dynamic characteristics are simulated. The results are shown to be in good agreement with experiments. The project supported by National Natural Science Fundation of China and LNM of Institute of Mechanics. CAS     

Shock–vortex interactions in a soap film

This work experimentally visualizes the interaction of a quasi-one-dimensional moving shock wave with a two-dimensional vortex in a soap film for the first time. A vertical soap film shock tube was used to generate a quasi-one-dimensional moving shock wave and a NACA-0012 airfoil intruded into the soap film was towed to shed the starting vortex. The interesting interaction phenomena were then visualized using a traditional high-speed flash photography. The concentration of sodium dodecyl sulfate (SDS) used was 0.5 CMC (critical micelle concentration) to keep the surfactant molecules behave as two-dimensional gases. A sequence of pictures shows that the shock is distorted non-symmetrically as it passes through the spiral vortex flow field and the vortex structure is compressed in the direction normal to the shock. These flow features observed in soap films are qualitatively similar to their counterparts in gases. In addition, the visualization of the interactions of a quasi-one-dimensional moving shock wave with a Kármán vortex street are presented.      

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

3D tangentially injected swirling recirculating flow in a nozzle with a slotted‐tube—Effects of groove parameters

A numerical prediction for 3D swirling recirculating flow in an air‐jet spinning nozzle with a slotted‐tube is carried out with the realizable k–ε turbulence model. The effects of the groove parameters on the flow and yarn properties are investigated. The simulation results show that some factors, such as reverse flow upstream of the injector, vortex breakdown downstream of the injector, corner recirculation zone (CRZ) behind the step and vortex ring in the groove caused by the groove geometric variation, are significantly related to fluid flow, and consequently to yarn properties. With increasing groove height, the length of the CRZ increases, while the initial vortex ring in the groove decreases and a same direction rotating vortex forms in the bottom of the groove. Similarly, as the groove width increases, the extent of both vortex breakdown in downstream of the injectors and the vortex ring in the groove increases slightly, whereas the CRZ lengths in stream‐wise direction decrease. Some factors, such as the negative tangential velocities, the size of the vortex rings in the grooves and the CRZ, are constant for nozzles with different groove lengths. Copyright © 2009 John Wiley & Sons, Ltd.     

Large eddy simulation of a vortex ring impinging on a three-dimensional circular cylinder

A vortex ring impacting a three-dimensional circular cylinder is studied using large eddy simulation (LES) for a Reynolds number Re = 4 × 10⁴ based on the initial translation speed and diameter of the vortex ring. We have investigated the evolution of vortical structures and identified three typical evolution phases. When the primary vortex closely approaches to the cylinder, a secondary vortex is generated and its segment parts move inward to the primary vortex ring. then two large-scale loop-like vortices are formed to evolve in opposite directions. Thirdly, the two loop-like vortices collide with each other to form complicated small-scale vortical structures. Moreover, a series of hair-pin vortices are generated due to the stretching and deformation of the tertiary vortex. The trajectories of vortical structures and the relevant evolution speeds are analyzed. The total kinetic energy and enstrophy are investigated to reveal their properties relevant to the three evolution phases.     

Study on the Fuel Air Mixing Induced by a Shock Wave Propagating into a H2-Air Interface

2nd-order upwind TVD scheme was used to solve the laminar, fully Navier-Stokes equations. The numerical simulations were done on the propagation of a shock wave with Ma _S = 2 and 4 into a hydrogen and air mixture in a duct and a duct with a rearward step. The results indicate that a swirling vortex may be generated in the lopsided interface behind the moving shock. Meanwhile, the complex shock system is also formed in this shear flow region. A large swirling vortex is produced and the fuel mixing can be enhanced by a shock wave at low Mach number. But in a duct with a rearward step, the shock almost disappears in hydrogen for Ma _S = 2. The shock in hydrogen will become strong if Ma _S is large. Similar to the condition of a shock moving in a duct full of hydrogen and air, a large vortex can be formed in the shear flow region. The large swirling vortex even gets through the reflected shock and impacts on the lower wall. Then, the distribution of hydrogen behind the rearward step is divided into two regions. The transition from regular reflection to Mach reflection was observed as well in case Ma _S = 4.     

Wakes of elliptical cylinders at low Reynolds number

The wakes of elliptical cylinders are numerically investigated at a Reynolds number Re_D = 150. ANSYS-Fluent, based on the finite volume method, is used to simulate two-dimensional Newtonian fluid flow. The cylinder cross-sectional aspect ratio (AR) is varied from 0.25 to 1.0 (circular cylinder), and the angle of attack (α) of the cylinder is changed as α = 0° – 90°. With the changes in AR and α, three distinct wake patterns (patterns I, II, III) are observed, associated with different characteristics of fluid forces. Steady wake (pattern I) is characterised by two steady bubbles forming behind the cylinder, occurring at AR < 0.37 and α < 2.5°. Time-mean drag and fluctuating lift coefficients are small. Pattern II refers to Karman wake followed by steady wake (AR ≥ 0.37 – 0.67, depending on α) with the Karman street transitioning to two steady shear layers downstream. An inflection angle α_i is identified where the time-mean drag of the elliptical cylinder is identical to that of a circular cylinder. Pattern III is the Karman wake followed by secondary wake (AR ≤ 0.67, α > 52°), where the Karman street forming behind the cylinder is modified to a secondary vortex street with a low frequency. The Time-mean drag coefficient is maximum for this pattern.     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Optimized incompressible smoothed particle hydrodynamics methods and validations

The solution of the Poisson's equation used by the incompressible smoothed particle hydrodynamics (ISPH) methods for estimating the pressure field is expensive in CPU time. The CPU time, consumed by the inversion of the operator ∇(1/ρ∇) and the estimation of the right hand side of the Poisson's equation, increases with the number N of particles used in a purely Lagrangian framework. In this work, this default of ISPH methods is overcome by solving the Poisson's equation on a Cartesian grid. This SPH-mesh coupling is equivalent to the particle in cell method. In a first step, in order to analyze its efficiency, the optimized version of two ISPH methods (divergence free and density invariant) is compared with the standard weakly compressible SPH method through two benchmarks of incompressible bidimensional flows characterized by the Reynolds number Re, Lamb-Oseen vortex (10 ≤Re≤ 100) and lid-driven cavity flow (100 ≤Re≤ 1000). In a second step, the numerical results obtained by the three SPH methods are compared to laboratory experimental data of a dam break flow in order to show the performance of the SPH-mesh coupling in a practical and complex flow problem. As in the configuration of the experimental setup, the numerical results are obtained for a Reynolds number Re = 3.8 10⁶.