Showing papers on "Vortex shedding published in 2009"

Three-dimensional flows around low-aspect-ratio flat-plate wings at low Reynolds numbers

Abstract: Three-dimensional flows over impulsively translated low-aspect-ratio flat plates are investigated for Reynolds numbers of 300 and 500, with a focus on the unsteady vortex dynamics at post-stall angles of attack. Numerical simulations, validated by an oil tow-tank experiment, are performed to study the influence of aspect ratio, angle of attack and planform geometry on the wake vortices and the resulting forces on the plate. Immediately following the impulsive start, the separated flows create wake vortices that share the same topology for all aspect ratios. At large time, the tip vortices significantly influence the vortex dynamics and the corresponding forces on the wings. Depending on the aspect ratio, angle of attack and Reynolds number, the flow at large time reaches a stable steady state, a periodic cycle or aperiodic shedding. For cases of high angles of attack, an asymmetric wake develops in the spanwise direction at large time. The present results are compared to higher Reynolds number flows. Some non-rectangular planforms are also considered to examine the difference in the wakes and forces. After the impulsive start, the time at which maximum lift occurs is fairly constant for a wide range of flow conditions during the initial transient. Due to the influence of the tip vortices, the three-dimensional dynamics of the wake vortices are found to be quite different from the two-dimensional von Karman vortex street in terms of stability and shedding frequency.

On vortex shedding from an airfoil in low-Reynolds-number flows

Abstract: Development of coherent structures in the separated shear layer and wake of an airfoil in low-Reynolds-number flows was studied experimentally for a range of airfoil chord Reynolds numbers, 55 × 10 3 ≤ Re c ≤ 210 × 10 3 , and three angles of attack, α = 0°, 5° and 10°. To illustrate the effect of separated shear layer development on the characteristics of coherent structures, experiments were conducted for two flow regimes common to airfoil operation at low Reynolds numbers: (i) boundary layer separation without reattachment and (ii) separation bubble formation. The results demonstrate that roll-up vortices form in the separated shear layer due to the amplification of natural disturbances, and these structures play a key role in flow transition to turbulence. The final stage of transition in the separated shear layer, associated with the growth of a sub-harmonic component of fundamental disturbances, is linked to the merging of the roll-up vortices. Turbulent wake vortex shedding is shown to occur for both flow regimes investigated. Each of the two flow regimes produces distinctly different characteristics of the roll-up and wake vortices. The study focuses on frequency scaling of the investigated coherent structures and the effect of flow regime on the frequency scaling. Analysis of the results and available data from previous experiments shows that the fundamental frequency of the shear layer vortices exhibits a power law dependency on the Reynolds number for both flow regimes. In contrast, the wake vortex shedding frequency is shown to vary linearly with the Reynolds number. An alternative frequency scaling is proposed, which results in a good collapse of experimental data across the investigated range of Reynolds numbers.

Vortex wakes of a flapping foil

Abstract: We present an experimental study of a symmetric foil performing pitching oscillations in a vertically flowing soap film. By varying the frequency and amplitude of the oscillation we visualize a variety of wakes with up to 16 vortices per oscillation period, including von Karman vortex street, inverted von Karman vortex street, 2P wake, 2P+2S wake and novel wakes ranging from 4P to 8P. We map out the wake types in a phase diagram spanned by the width-based Strouhal number and the dimensionless amplitude. We follow the time evolution of the vortex formation near the round leading edge and the shedding process at the sharp trailing edge in detail. This allows us to identify the origins of the vortices in the 2P wake, to understand that two distinct 2P regions are present in the phase diagram due to the timing of the vortex shedding at the leading edge and the trailing edge and to propose a simple model for the vorticity generation. We use the model to describe the transition from 2P wake to 2S wake with increasing oscillation frequency and the transition from the von Karman wake, typically associated with drag, to the inverted von Karman wake, typically associated with thrust generation.

Vortex-induced vibrations of two cylinders in tandem arrangement in the proximity-wake interference region.

Abstract: We investigate numerically vortex-induced vibrations (VIV) of two identical two-dimensional elastically mounted cylinders in tandem in the proximity–wake interference regime at Reynolds number Re = 200 for systems having both one (transverse vibrations) and two (transverse and in-line) degrees of freedom (1-DOF and 2-DOF, respectively). For the 1-DOF system the computed results are in good qualitative agreement with available experiments at higher Reynolds numbers. Similar to these experiments our simulations reveal: (1) larger amplitudes of motion and a wider lock-in region for the tandem arrangement when compared with an isolated cylinder; (2) that at low reduced velocities the vibration amplitude of the front cylinder exceeds that of the rear cylinder; and (3) that above a threshold reduced velocity, large-amplitude VIV are excited for the rear cylinder with amplitudes significantly larger than those of the front cylinder. By analysing the simulated flow patterns we identify the VIV excitation mechanisms that lead to such complex responses and elucidate the near-wake vorticity dynamics and vortex-shedding modes excited in each case. We show that at low reduced velocities vortex shedding provides the initial excitation mechanism, which gives rise to a vertical separation between the two cylinders. When this vertical separation exceeds one cylinder diameter, however, a significant portion of the incoming flow is able to pass through the gap between the two cylinders and the gap-flow mechanism starts to dominate the VIV dynamics. The gap flow is able to periodically force either the top or the bottom shear layer of the front cylinder into the gap region, setting off a series of very complex vortex-to-vortex and vortex-to-cylinder interactions, which induces pressure gradients that result in a large oscillatory force in phase with the vortex shedding and lead to the experimentally observed larger vibration amplitudes. When the vortex shedding is the dominant mechanism the front cylinder vibration amplitude is larger than that of the rear cylinder. The reversing of this trend above a threshold reduced velocity is associated with the onset of the gap flow. The important role of the gap flow is further illustrated via a series of simulations for the 2-DOF system. We show that when the gap-flow mechanism is triggered, the 2-DOF system can develop and sustain large VIV amplitudes comparable to those observed in the corresponding (same reduced velocity) 1-DOF system. For sufficiently high reduced velocities, however, the two cylinders in the 2-DOF system approach each other, thus significantly reducing the size of the gap region. In such cases the gap flow is entirely eliminated, and the two cylinders vibrate together as a single body with vibration amplitudes up to 50% lower than the amplitudes of the corresponding 1-DOF in which the gap flow is active. Three-dimensional simulations are also carried out to examine the adequacy of two-dimensional simulations for describing the dynamic response of the tandem system at Re = 200. It is shown that even though the wake transitions to a weakly three-dimensional state when the gap flow is active, the three-dimensional modes are too weak to affect the dynamic response of the system, which is found to be identical to that obtained from the two-dimensional computations.

Resonance and propulsion performance of a heaving flexible wing

Abstract: The influence of the bending rigidity of a flexible heaving wing on its propulsive performance in a two-dimensional imposed parallel flow is investigated in the inviscid limit. Potential flow theory is used to describe the flow over the flapping wing. The vortical wake of the wing is accounted for by the shedding of point vortices with unsteady intensity from the wing’s trailing edge. The trailing-edge flapping amplitude is shown to be maximal for a discrete set of values of the rigidity, at which a resonance occurs between the forcing frequency and a natural frequency of the system. A quantitative comparison of the position of these resonances with linear stability analysis results is presented. Such resonances induce maximum values of the mean developed thrust and power input. The flapping efficiency is also shown to be greatly enhanced by flexibility.

Shallow and deep dynamic stall for flapping low Reynolds number airfoils

Abstract: We consider a combined experimental (based on flow visualization, direct force measurement and phaseaveraged 2D particle image velocimetry in a water tunnel), computational (2D Reynolds-averaged Navier-Stokes) and theoretical (Theodorsen’s formula) approach to study the fluid physics of rigid-airfoil pitch-plunge in nominally two-dimensional conditions. Shallow-stall (combined pitch-plunge) and deep-stall (pure-plunge) are compared at a reduced frequency commensurate with flapping-flight in cruise in nature. Objectives include assessment of how well attached-flow theory can predict lift coefficient even in the presence of significant separation, and how well 2D velocimetry and 2D computation can mutually validate one another. The shallow-stall case shows promising agreement between computation and experiment, while in the deepstall case, the computation’s prediction of flow separation lags that of the experiment, but eventually evinces qualitatively similar leading edge vortex size. Dye injection was found to give good qualitative match with particle image velocimetry in describing leading edge vortex formation and return to flow reattachment, and also gave evidence of strong spanwise growth of flow separation after leadingedge vortex formation. Reynolds number effects, in the range of 10,000-60,000, were found to influence the size of laminar separation in those phases of motion where instantaneous angle of attack was well below stall, but have limited effect on post-stall flowfield behavior. Discrepancy in lift coefficient time history between experiment, theory and computation was mutually comparable, with no clear failure of Theodorsen’s formula. This is surprising and encouraging, especially for the deep-stall case, because the theory’s assumptions are clearly violated, while its prediction of lift coefficient remains useful for capturing general trends.

On the stability and extension of reduced-order Galerkin models in incompressible flows

Abstract: Proper orthogonal decomposition (POD) has been used to develop a reduced-order model of the hydrodynamic forces acting on a circular cylinder Direct numerical simulations of the incompressible Navier–Stokes equations have been performed using a parallel computational fluid dynamics (CFD) code to simulate the flow past a circular cylinder Snapshots of the velocity and pressure fields are used to calculate the divergence-free velocity and pressure modes, respectively We use the dominant of these velocity POD modes (a small number of eigenfunctions or modes) in a Galerkin procedure to project the Navier–Stokes equations onto a low-dimensional space, thereby reducing the distributed-parameter problem into a finite-dimensional nonlinear dynamical system in time The solution of the reduced dynamical system is a limit cycle corresponding to vortex shedding We investigate the stability of the limit cycle by using long-time integration and propose to use a shooting technique to home on the system limit cycle We obtain the pressure-Poisson equation by taking the divergence of the Navier–Stokes equation and then projecting it onto the pressure POD modes The pressure is then decomposed into lift and drag components and compared with the CFD results

Dynamics of laminar separation bubbles at low-Reynolds-number aerofoils

Abstract: The laminar separation bubble on an SD7003 aerofoil at a Reynolds number Re = 66000 was investigated to determine the dominant frequencies of the transition process and the flapping of the bubble. The measurements were performed with a high-resolution time-resolved particle image velocimetry (TR-PIV) system. Contrary to typical measurements performed through conventional PIV, the different modes can be identified by applying TR-PIV. The interaction between the shed vortices is analysed, and their significance for the production of turbulence is presented. In the shear layer above the bubble the generation and amplification of vortices due to Kelvin–Helmholtz instabilities is observed. It is found that these instabilities have a weak coherence in the spanwise direction. In a later stage of transition these vortices lead to a three-dimensional breakdown to turbulence.

MTV measurements of the vortical field in the wake of an airfoil oscillating at high reduced frequency

Abstract: We present an experimental investigation of the flow structure and vorticity field in the wake of a NACA-0012 airfoil pitching sinusoidally at small amplitude and high reduced frequencies. Molecular tagging velocimetry is used to quantify the characteristics of the vortex array (circulation, peak vorticity, core size, spatial arrangement) and its downstream evolution over the first chord length as a function of reduced frequency. The measured mean and fluctuating velocity fields are used to estimate the mean force on the airfoil and explore the connection between flow structure and thrust generation.Results show that strong concentrated vortices form very rapidly within the first wavelength of oscillation and exhibit interesting dynamics that depend on oscillation frequency. With increasing reduced frequency the transverse alignment of the vortex array changes from an orientation corresponding to velocity deficit (wake profile) to one with velocity excess (reverse Karman street with jet profile). It is found, however, that the switch in the vortex array orientation does not coincide with the condition for crossover from drag to thrust. The mean force is estimated from a more complete control volume analysis, which takes into account the streamwise velocity fluctuations and the pressure term. Results clearly show that neglecting these terms can lead to a large overestimation of the mean force in strongly fluctuating velocity fields that are characteristic of airfoils executing highly unsteady motions. Our measurements show a decrease in the peak vorticity, as the vortices convect downstream, by an amount that is more than can be attributed to viscous diffusion. It is found that the presence of small levels of axial velocity gradients within the vortex cores, levels that can be difficult to measure experimentally, can lead to a measurable decrease in the peak vorticity even at the centre of the flow facility in a flow that is expected to be primarily two-dimensional.

Three-dimensional vorticity patterns of cylinder wakes

Abstract: The vortex organization of cylinder wakes is experimentally studied by time-resolved tomographic Particle Image Velocimetry at Reynolds numbers ranging from 180 to 5,540. Time resolved measurements are performed at Re = 180, 360 and 540, whereas the transitional (Re = 1,080) and turbulent regimes (Re = 5,540) are investigated by snapshots separated in phase by more than π/4. The vortex structure evolution is visualized by the 3D vorticity field, revealing a regular shedding at the lowest Reynolds, whereas at Re > 500 the Benard-Karman vortex street exhibits counter-rotating stream-wise vortex pairs (characteristic of Mode B) dominating the 3D motion. The regime at Re = 360 produces a transitional pattern where the counter-rotating vortex pairs (Mode B), coexist with profoundly distorted shedding of oblique elements forming a chain of rhombus-like vortex cells. In the turbulent flow regime (Re = 5,540) a large increase in the range of flow scales is directly observed with the appearance of Kelvin-Helmholtz type vortices in the separated shear layer consistently with what is abundantly reported in literature. The statistical description of the secondary structures is inferred from a 3D autocorrelation analysis yielding two span-wise wavelengths for the counter-rotating pairs, an inner length given by (twice) the distance between counter-rotating elements and an outer one given by the distance between pairs. The uncertainty analysis of the present tomographic PIV experiments reveals that this approach is suited for the investigation of vortex wakes with a typical error of 2 and 10% on the velocity and vorticity vectors, respectively.

Prediction of vortex-induced vibration response by employing controlled motion

Abstract: In order to predict response and wake modes for elastically mounted circular cylinders in a fluid flow, we employ controlled-vibration experiments, comprised of prescribed transverse vibration of a cylinder in the flow, over a wide regime of amplitude and frequency. A key to this study is the compilation of high-resolution contour plots of fluid force, in the plane of normalized amplitude and wavelength. With such resolution, we are able to discover discontinuities in the force and phase contours, which enable us to clearly identify boundaries separating different fluid-forcing regimes. These appear remarkably similar to boundaries separating different vortex-formation modes in the map of regimes by Williamson & Roshko (J. Fluids Struct., vol. 2, 1988, pp. 355–381). Vorticity measurements exhibit the 2S, 2P and P + S vortex modes, as well as a regime in which the vortex formation is not synchronized with the body vibration. By employing such fine-resolution data, we discover a high-amplitude regime in which two vortex-formation modes overlap. Associated with this overlap regime, we identify a new distinct mode of vortex formation comprised of two pairs of vortices formed per cycle, where the secondary vortex in each pair is much weaker than the primary vortex. This vortex mode, which we define as the 2POVERLAP mode (2PO), is significant because it is responsible for generating the peak resonant response of the body. We find that the wake can switch intermittently between the 2P and 2PO modes, even as the cylinder is vibrating with constant amplitude and frequency. By examining the energy transfer from fluid to body motion, we predict a free-vibration response which agrees closely with measurements for an elastically mounted cylinder. In this work, we introduce the concept of an ‘energy portrait’, which is a plot of the energy transfer into the body motion and the energy dissipated by damping, as a function of normalized amplitude. Such a plot allows us to identify stable and unstable amplitude-response solutions, dependent on the rate of change of net energy transfer with amplitude (the sign of dE*/dA*). Our energy portraits show how the vibration system may exhibit a hysteretic mode transition or intermittent mode switching, both of which correspond with such phenomena measured from free vibration. Finally, we define the complete regime in the amplitude–wavelength plane in which free vibration may exist, which requires not only a periodic component of positive excitation but also stability of the equilibrium solutions.

Unsteady fluid-structure interactions of membrane airfoils at low Reynolds numbers

Abstract: Membrane wings are used both in nature and small aircraft as lifting surfaces. Separated flows are common at low Reynolds numbers and are the main sources of unsteadiness. Yet, the unsteady aspects of the fluid-structure interactions of membrane airfoils are largely unknown. An experimental study of unsteady aerodynamics of two-dimensional membrane airfoils at low Reynolds numbers has been conducted. Measurements of membrane shape with a high-speed camera were complemented with the simultaneous measurements of unsteady velocity field with a high frame-rate particle image velocimetry system and flow visualization. Vibrations of the membrane and mode shapes were investigated as a function of angle of attack and free stream velocity. While the mean membrane shape is not very sensitive to angle of attack, the amplitude and mode of the vibrations of the membrane depend on the relative location and the magnitude of the unsteadiness of the separated shear layer. The results indicate strong coupling of unsteady flow with the membrane oscillations. There is evidence of coupling of membrane oscillations with the vortex shedding in the wake, in particular, for the post-stall incidences. Comparison of rigid (but cambered) and flexible membrane airfoils shows that the flexibility might delay the stall. Hence this is a potential passive flow control method using flexibility in nature and engineering applications.

An unsteady point vortex method for coupled fluid–solid problems

Abstract: A method is proposed for the study of the two-dimensional coupled motion of a general sharp-edged solid body and a surrounding inviscid flow. The formation of vorticity at the body’s edges is accounted for by the shedding at each corner of point vortices whose intensity is adjusted at each time step to satisfy the regularity condition on the flow at the generating corner. The irreversible nature of vortex shedding is included in the model by requiring the vortices’ intensity to vary monotonically in time. A conservation of linear momentum argument is provided for the equation of motion of these point vortices (Brown–Michael equation). The forces and torques applied on the solid body are computed as explicit functions of the solid body velocity and the vortices’ position and intensity, thereby providing an explicit formulation of the vortex–solid coupled problem as a set of non-linear ordinary differential equations. The example of a falling card in a fluid initially at rest is then studied using this method. The stability of broadside-on fall is analysed and the shedding of vorticity from both plate edges is shown to destabilize this position, consistent with experimental studies and numerical simulations of this problem. The reduced-order representation of the fluid motion in terms of point vortices is used to understand the physical origin of this destabilization.

Flow control around a circular cylinder using pulsed dielectric barrier discharge surface plasma

Abstract: Dielectric barrier discharge (DBD) plasma actuators have been used to control the flow around a circular cylinder at Re=15 000, where the near-wake structure was studied using time-resolved particle image velocimetry with simultaneous measurements of the dynamic lift and drag forces. It was shown that the vortex shedding was suppressed when the surface plasma placed near the natural separation point was activated in a pulsed mode at nondimensional frequency, fp+, above 0.6 with a force coefficient, Cp, greater than 0.05%. Plasma actuator performance on flow control was summarized by mapping the changes in drag and lift fluctuations as a function of the forcing frequency and the force coefficient. They showed that more than 70% reduction in lift fluctuations was obtained with up to 32% drag reduction at fp+=2.0 and Cp=0.32%. Here, narrowing of the wake was observed as the plasma promoted shear-layer roll-ups at the forcing frequency. This, however, did not affect the shear layer on the opposite side of the...

Time-resolved wake structure and kinematics of bat flight

Abstract: We present synchronized time-resolved measurements of the wing kinematics and wake velocities for a medium sized bat, Cynopterus brachyotis, flying at low-medium speed in a closed-return wind tunnel. Measurements of the motion of the body and wing joints, as well as the resultant wake velocities in the Trefftz plane are recorded at 200 Hz (approximately 28–31 measurements per wing beat). Circulation profiles are found to be quite repeatable although variations in the flight profile are visible in the wake vortex structures. The circulation has almost constant strength over the middle half of the wing beat (defined according the vertical motion of the wrist, beginning with the downstroke). A strong streamwise vortex is observed to be shed from the wingtip, growing in strength during the downstroke, and persisting during much of the upstroke. At relatively low flight speeds (4.3 m/s), a closed vortex structure behind the bat is postulated.

A model for the symmetry breaking of the reverse Bénard–von Kármán vortex street produced by a flapping foil

Abstract: The vortex streets produced by a flapping foil of span to chord aspect ratio of 4:1 are studied in a hydrodynamic tunnel experiment. In particular, the mechanisms giving rise to the symmetry breaking of the reverse Benard–von Karman (BvK) vortex street that characterizes fishlike swimming and forward flapping flight are examined. Two-dimensional particle image velocimetry (PIV) measurements in the midplane perpendicular to the span axis of the foil are used to characterize the different flow regimes. The deflection angle of the mean jet flow with respect to the horizontal observed in the average velocity field is used as a measure of the asymmetry of the vortex street. Time series of the vorticity field are used to calculate the advection velocity of the vortices with respect to the free stream, defined as the phase velocity Uphase, as well as the circulation Γ of each vortex and the spacing ξ between consecutive vortices in the near wake. The observation that the symmetry-breaking results from the formation of a dipolar structure from each couple of counter-rotating vortices shed on each flapping period serves as the starting point to build a model for the symmetry-breaking threshold. A symmetry-breaking criterion based on the relation between the phase velocity of the vortex street and an idealized self-advection velocity of two consecutive counter-rotating vortices in the near wake is established. The predicted threshold for symmetry breaking accounts well for the deflected wake regimes observed in the present experiments and may be useful to explain other experimental and numerical observations of similar deflected propulsive vortex streets reported in the literature.

Swimming performance of a biomimetic compliant fish-like robot

Abstract: Digital particle image velocimetry and fluorescent dye visualization are used to characterize the performance of fish-like swimming robots. During nominal swimming, these robots produce a ‘V’-shaped double wake, with two reverse-Karman streets in the far wake. The Reynolds number based on swimming speed and body length is approximately 7500, and the Strouhal number based on flapping frequency, flapping amplitude, and swimming speed is 0.86. It is found that swimming speed scales with the strength and geometry of a composite wake, which is constructed by freezing each vortex at the location of its centroid at the time of shedding. Specifically, we find that swimming speed scales linearly with vortex circulation. Also, swimming speed scales linearly with flapping frequency and the width of the composite wake. The thrust produced by the swimming robot is estimated using a simple vortex dynamics model, and we find satisfactory agreement between this estimate and measurements made during static load tests.

A generalized mean-field model of the natural and high-frequency actuated flow around a high-lift configuration

Abstract: A low-dimensional Galerkin model is proposed for the flow around a high-lift configuration, describing natural vortex shedding, the high-frequency actuated flow with increased lift and transients between both states. The form of the dynamical system has been derived from a generalized mean-field consideration. Steady state and transient URANS (unsteady Reynolds-averaged Navier–Stokes) simulation data are employed to derive the expansion modes and to calibrate the system parameters. The model identifies the mean field as the mediator between the high-frequency actuation and the low-frequency natural shedding instability.

Acoustics, Vortex Shedding, and Low-Frequency Dynamics Interaction in an Unstable Hybrid Rocket

Abstract: This paper deals with an experimental investigation into the stability behavior of a hybrid rocket where gaseous oxygen is fed with either an axial conical subsonic nozzle or a radial injector. The influence of the oxidizer-injection configurations on the motor stability is thoroughly examined. These distinct oxidizer-injection techniques allowed unveiling key and so far unreported features of the hybrid rocket combustion stability, especially emphasizing the role of vortex shedding which occurs in both the pre- and postcombustion chamber. Axial and radial injectors caused completely stable and unstable combustor operations, respectively, and this fact has been attributed to the fluid dynamics and unsteady heat release at the entrance of the fuel grain port. In particular, the unstable combustion in the radial-flow injector motor was dominated by low-frequency pressure oscillations, around 10-20 Hz. These low-frequency pressure oscillations were always accompanied by longitudinal acoustic modes. In some cases, the pressure oscillations abruptly increased, reaching peak-to-peak amplitude close to 70% of the mean chamber pressure, which is somewhat unusual for hybrid engines. Vortex shedding in the aft-mixing chamber is considered as the main driving mechanism of this latter behavior.