Modification of static-wing tip vortex via a slender half-delta wing

The modification of the tip vortex generated by a rectangular NACA 0012 wing via a tip-mounted slender half-delta wing (HDW) was attempted experimentally at Re=2.81×10⁵. In addition to the increase in lift with increasing HDW deflection, compared to the baseline wing, the roll-up process of the tip vortex was also found to be significantly modified, as a result of the breakdown of the HDW vortex. The addition of the HDW also caused an increased total drag. Fortunately, the lift-induced drag was found to be reduced compared to its baseline counterpart for 0° and 5° HDW deflections. The change in the lift-induced drag also translates into a virtually unchanged profile drag, regardless of HDW deflection. In short, the largest lift-induced drag reduction achieved by the zero-deflection HDW resulted in an improved lift-to-drag ratio, at high angle-of-attack range, compared to the baseline wing.     

Near-field tip vortex behind a swept wing model

The near-field flow structure of a tip vortex behind a sweptback and tapered NACA 0015 wing was investigated and compared with a rectangular wing at the same lift force and Re=1.81×10⁵. The tangential velocity decreased with the downstream distance while increased with the airfoil incidence. The core radius was about 3% of the root chord c _r, regardless of the downstream distance and α for α<8°. The core axial velocity was always wake-like. The core Γ_c and total Γ_o circulation of the tip vortex remained nearly constant up to x/c _r=3.5 and had a Γ_c/Γ_o ratio of 0.63. The total circulation of the tip vortex accounted for only about 40% of the bound root circulation Γ_b. For a rectangular wing, the axial flow exhibited islands of wake- and jet-like velocity distributions with Γ_c/Γ_o=0.75 and Γ_o/Γ_b=0.70. For the sweptback and tapered wing tested, the inner region of the tip vortex flow exhibited a self-similar behavior for x/c _r≥1.0. The lift force computed from the spanwise circulation distributions agreed well with the force-balance data. A large difference in the lift-induced drag was, however, observed between the wake integral method and the inviscid lifting-line theory.     

Lift-induced drag of a cambered wing for Re ≤ 1 × 106

An investigation of the dependence of the lift-induced drag coefficient C _Di of a square-tipped, cambered wing model on Reynolds number for Re ≤ 1 × 10⁶ was conducted. Computed based on the vorticity distribution inferred from the near-field cross-flow velocity measurements of the tip vortex, different C _Di prediction schemes were used. The effect of measurement plane size and grid resolution on the C _Di calculations was also identified. The C _Di estimated by the integral method was found to increase with increasing Re and was below the C _Di = C _l²/πeAR prediction. Limits on the measurement plane size and grid resolution were determined to be at least 40% larger than the vortex outside diameter and no larger than 0.63% chord, respectively, in order to provide a good estimate of the induced drag.     

Experimental investigation on tip vortices and aerodynamics

The tip vortices and aerodynamics of a NACA0012 wing in the vicinity of the ground were studied in a wind tunnel. The wing tip vortex structures and lift/drag forces were measured by a seven-hole probe and a force balance, respectively. The evolution of the flow structures and aerodynamics with a ground height were analyzed. The vorticity of tip vortices was found to reduce with the decreasing of the ground height, and the position of vortex-core moved gradually to the outboard of the wing tip. Therefore, the down-wash flow induced by the tip vortices was weakened. However, vortex breakdown occurred as the wing lowered to the ground. From the experimental results of aerodynamics, the maximum lift-to-drag ratio was observed when the angle of attack was 2.5° and the ground clearance was 0.2.     

Large aerodynamic forces on a sweeping wing at low Reynolds number

The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically. After an initial start from rest, the wing is made to execute an azimuthal rotation (sweeping) at a large angle of attack and constant angular velocity. The Reynolds number (Re) considered in the present note is 480 (Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root). During the constant-speed sweeping motion, the stall is absent and large and approximately constant lift and drag coefficients can be maintained. The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows. Soon after the initial start, a vortex ring, which consists of the leading-edge vortex (LEV), the starting vortex, and the two wing-tip vortices, is formed in the wake of the wing. During the subsequent motion of the wing, a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength. This prevents the LEV from shedding. As a result, the size of the vortex ring increases approximately linearly with time, resulting in an approximately constant time rate of the first moment of vorticity, or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wingspan. The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force. The project supported by the National Natural Science Foundation of China (10232010)     

Control of vortex shedding in an axisymmetric bluff body wake

In the present experimental study the effect of a control disc mounted at the rear of an axisymmetric blunt-based body of revolution, first studied by Mair, is investigated in the Reynolds number range 3×10³≤Re_D≤5×10⁴ . As the distance of the control disc from the blunt base is increased, four vortex shedding regimes are identified: at small distances there is no effect, then a sharp increase of vortex shedding activity and total drag is observed, followed by an interval with reduced activity and drag and finally at large distances a regime where the flow around the main body and disc become essentially independent, i.e. where the drag forces of the two elements become additive. The near and far wake velocity fields corresponding to the different regimes are documented with time- and phase-averaged hot-wire and LDA measurements, with spectral analysis of the data and with flow visualizations of the near wake. The results are used to develop an improved understanding of the instability mechanism leading to high vortex shedding activity.     

A Numerical Study of the Flow Around a Model Winged Seed in Auto-Rotation

In this study the flow around a winged-seed in auto-rotation is characterized using direct numerical simulations (DNS) at Reynolds number in the range 80–240, based on the descent speed and a characteristic chord length. In this range, the flow is approximately steady when observed from a reference frame fixed to the seed. For all cases, the flow structure consists of a wing tip vortex which describes a helical path, a vortex shed behind the nut of the seed and a stable leading edge vortex above the wing surface which merges with the tip vortex. With increasing Reynolds number, the leading edge vortex becomes more intense and gets closer to the wing surface. The simulation results also show the formation of a spanwise flow on the upper surface of the wing, moving fluid towards the wing tip in a region downstream and beneath the leading edge vortex. This spanwise flow is rather weak inside the core of the leading edge vortex, and the analysis of the streamlines show a very weak transport of vorticity along the vortex for the cases under consideration. The analysis of the flow suggests that the stabilization of the leading edge vortex is mainly due to non-inertial accelerations, although viscous effects may contribute, specially at lower Re. Furthermore, the leading edge vortex has been characterized by analysing the flow variables averaged along cross-sections of the vortex. While some quantities, like the spanwise velocity or the pressure inside the vortex, are rather insensitive to the threshold used to define the leading edge vortex, the same is not true for the circulation of the vortex or its averaged spanwise vorticity, due to the viscous nature of the vortex. Finally, it is observed that the spanwise vorticity scales with the angular rotation of the seed for the different Re.     

Effects of a trapped vortex cell on a thick wing airfoil

The effects of a trapped vortex cell (TVC) on the aerodynamic performance of a NACA0024 wing model were investigated experimentally at Re = 10⁶ and 6.67×10⁵6.67\times 10^{5}. The static pressure distributions around the model and the wake velocity profiles were measured to obtain lift and drag coefficients, for both the clean airfoil and the controlled configurations. Suction was applied in the cavity region to stabilize the trapped vortex. For comparison, a classical boundary layer suction configuration was also tested. The drag coefficient curve of the TVC-controlled airfoil showed sharp discontinuities and bifurcative behavior, generating two drag modes. A strong influence of the angle of attack, the suction rate and the Reynolds number on the drag coefficient was observed. With respect to the clean airfoil, the control led to a drag reduction only if the suction was high enough. Compared to the classical boundary layer suction configuration, the drag reduction was higher for the same amount of suction only in a specific range of incidence, i.e., α = −2° to α = 6° and only for the higher Reynolds number. For all the other conditions, the classical boundary layer suction configuration gave better drag performances. Moderate increments of lift were observed for the TVC-controlled airfoil at low incidence, while a 20% lift enhancement was observed in the stall region with respect to the baseline. However, the same lift increments were also observed for the classical boundary layer suction configuration. Pressure fluctuation measurements in the cavity region suggested a very complex interaction of several flow features. The two drag modes were characterized by typical unsteady phenomena observed in rectangular cavity flows, namely the shear layer mode and the wake mode.     

Effect of apex strake incidence-angle on the vortex development and interaction of a double-delta wing

The vortex flow characteristics of a sharp-edged delta wing with an apex strake was investigated through the visualization and particle image velocimetry (PIV) measurement of the wing-leeward flow region, and the wing-surface pressure measurement. The wing model was a flat-plate, and 65°-sweep cropped-delta wing with sharp leading edges. The apex strake was also a flat-plate wing with a cropped-delta shape of 65°/90° sweep, and it can change its incidence angle. The flow Reynolds number was 2.2 × 10⁵ for the flow visualization and 8.2 × 10⁵ for the PIV and wing-surface pressure measurements. The physics of the vortex flow in the wing-leeward flow region and the suction-pressure distribution on the wing upper-surface were interrelated and analyzed. The effect of a positive (negative) strake incidence-angle was the upward movement of the strake and wing vortices away from (downward movement of the strake and wing vortices toward) the wing-upper surface and the delayed (enhanced) coiling interaction between them. This change of vortex flow characteristics projected directly on the suction pressure distribution on the wing upper-surface.     

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.     

Effect of velocity ratio on the streamwise vortex structures in the wake of a stack

The time-averaged velocity and streamwise vorticity fields within the wake of a stack were investigated in a low-speed wind tunnel using a seven-hole pressure probe. The experiments were conducted at a Reynolds number, based on the stack external diameter, of Re_D=2.3×10⁴. The stack, of aspect ratio AR=9, was mounted normal to a ground plane and was partially immersed in a flat-plate turbulent boundary layer, where the ratio of the boundary layer thickness to the stack height was δ/H≈0.5. The jet-to-cross-flow velocity ratio was varied from R=0 to 3, which covered the downwash, crosswind-dominated and jet-dominated flow regimes. In the downwash and crosswind-dominated flow regimes, two pairs of counter-rotating streamwise vortex structures were identified within the stack wake. The tip vortex pair located close to the free end of the stack, and the base vortex pair located close to the ground plane within the flat-plate boundary layer, were similar to those found in the wake of a finite circular cylinder, and were associated with the upwash and downwash flow fields within the stack wake, respectively. In the jet-dominated flow regime, a third pair of streamwise vortex structures was observed, referred to as the jet-wake vortex pair, which occurred within the jet-wake region above the free end of the stack. The jet-wake vortex pair had the same orientation as the base vortex pair and was associated with the jet rise. The peak vorticity and strength of the streamwise vortex structures were functions of the jet-to-cross-flow velocity ratio. For the tip vortex structures, their peak vorticity and strength reduced as the jet-to-cross-flow velocity ratio increased.     

PIV study of near-field tip vortex behind perforated Gurney flaps

The impact of Gurney flaps, of different heights and perforations, on the growth and development of a tip vortex, both along the tip and in the near field of a finite NACA 0012 wing, at Re = 1.05 × 10⁵ was investigated by using particle image velocimetry (PIV). Wind-tunnel force balance measurements were also made to supplement the PIV results. This study is a continuation of the work of Lee and Ko (Exp Fluids 46(6):1005–1019, 2009) on the near-wake measurements behind perforated Gurney flaps. The present results show that along the tip, the overall behavior of the secondary vortices and their interaction with the primary, or tip, vortex remained basically unchanged, regardless of flap height and perforation. The peak vorticity of the tip vortex, however, increased with flap height and always exhibited a local maximum at x/c = 0.8 (from the leading edge). In the near field, the strength and structure of the near-field tip vortex were found to vary greatly with the flap height and perforation. The small flaps produced a more concentrated tip vortex with an increased circulation, while the large Gurney flaps caused a disruption of the tip vortex. The disrupted vortex can, however, be re-established by the addition of flap perforation. The larger the flap perforation the more organized the tip vortex. The Gurney flaps have the potential to serve as an alternative off-design wake vortex control device.     

Effect on the flow and heat transfer characteristics for sinusoidal pulsating laminar flow in a heated square cylinder

Two-dimensional numerical simulation is performed to understand the effect of flow pulsation on the flow and heat transfer from a heated square cylinder at Re = 100. Numerical calculations are carried out by using a finite volume method based on the pressure-implicit with splitting of operators algorithm in a collocated grid. The effects of flow pulsation amplitude (0.2 ≤ A ≤ 0.8) and frequency (0 ≤ f _p  ≤ 20 Hz) on the detailed kinematics of flow (streamlines, vorticity patterns), the macroscopic parameters (drag coefficient, vortex shedding frequency) and heat transfer enhancement are presented in detail. The Strouhal number of vortices shedding, drag coefficient for non-pulsating flow are compared with the previously published data, and good agreement is found. The lock-on phenomenon is observed for a square cylinder in the present flow pulsation. When the pulsating frequency is within the lock-on regime, time averaged drag coefficient and heat transfer from the square cylinder is substantially augmented, and when the pulsating frequency in about the natural vortex shedding frequency, the heat transfer is also substantially enhanced. In addition, the influence of the pulsating amplitude on the time averaged drag coefficient, heat transfer enhancement and lock-on occurrence is discussed in detail.     

Calculation of all drag components for 3-D bodies in wind tunnel flows. A boundary element-wake approach

A boundary element-wake numerical approach is developed and used to determine all drag components of a three dimensional body in a wind tunnel flow. The approach decomposes the total drag into three components; the profile drag, the cross flow drag (induced drag), and the tunnel-wall effect component, each with its own physical significance. Additionally, the cross flow drag component is divided into two components, the vortex component and the source (dilatation) component. In the present approach, the transverse kinematics relations are expressed as integral representations of the axial vorticity and the transverse dilatation (source strength). This advantage permits the vortex and the source drag computations to be performed only in the vortical area of the transverse wake and hence avoids excessive computations. Also, the procedures distinguish the contribution of the transverse dilatation to the cross flow drag. The validity of the present procedure is examined by comparing the present results against the experimental data of reference [1] for a car and wing models. The comparison shows that the present computed total drag, for the wing and the car models, agrees very well with the experimental data, provided that the wake data are measured at survey planes moderately distant from the body.     

Flow around two tandem square cylinders near a plane wall

An experimental study has been conducted to investigate the flow around two identical square cylinders in tandem arrangement and placed near a plane wall at a Reynolds number of 6,300. The inter-cylinder spacing ratio was varied from S ^* = 0.5 to 6, and the cylinder-to-wall gap ratio from G ^* = 0.25 to 2. Totally, 42 cases were considered to systematically examine the effects of wall proximity and the mutual interference between the two cylinders in the normalized gap–spacing (G ^*–S ^*) plane. The flow fields were captured using digital particle image velocimetry, in conjunction with measurements of the fluid forces (drag and lift) acting on the downstream cylinder using a piezoelectric load cell. The results show that the flow is highly dependent on the combined values of G ^* and S ^*. Categories relating to G ^* could be broadly classified as small-gap regime (G ^* < 0.5) at which periodic vortex shedding from the cylinders is suppressed, intermediate-gap regime (0.5 < G ^* < 1) where vortex shedding occurs but is under the influence of the wall proximity, and large-gap regime (G ^* > 1) where the wall effects become negligible. Similarly, the flow interference between the two cylinders can be divided into three basic categories as a function of S ^*, namely, shielding regime at S ^* < 1, reattachment regime at 1 < S ^* < 3, and impinging regime at S ^* > 3. Variations of force coefficients, amplitude spectra, Strouhal numbers, and Reynolds shear stress with G ^* and S ^* are presented to characterize the different flow regimes.     

The effect of a wake-mounted splitter plate on the flow around a surface-mounted finite-height circular cylinder

The influence of a wake-mounted splitter plate on the flow around a surface-mounted circular cylinder of finite height was investigated experimentally using a low-speed wind tunnel. The experiments were conducted at a Reynolds number of Re=7.4×10⁴ for cylinder aspect ratios of AR=9, 7, 5 and 3. The thickness of the boundary layer on the ground plane relative to the cylinder diameter was δ/D=1.5. The splitter plates were mounted on the wake centreline with negligible gap between the base of the cylinder and the leading edge of the plate. The lengths of the splitter plates, relative to the cylinder diameter, ranged from L/D=1 to 7, and the plate height was always equal to the cylinder height. Measurements of the mean drag force coefficient were obtained with a force balance, and measurements of the vortex shedding frequency were obtained with a single-component hot-wire probe situated in the wake of the cylinder–plate combination. Compared to the well-studied case involving an infinite circular cylinder, the splitter plate was found to be a less effective drag-reduction device for finite circular cylinders. Significant reduction in the mean drag coefficient was realized only for the finite circular cylinder of AR=9 with intermediate-length splitter plates of L/D=1–3. The mean drag coefficients of the other cylinders were almost unchanged. In terms of its effect on vortex shedding, a splitter plate of sufficient length was able to suppress Kármán vortex shedding for all of the finite circular cylinders tested. For AR=9, vortex shedding suppression occurred for L/D≥5, which is similar to the case of the infinite circular cylinder. For the smaller-aspect-ratio cylinders, however, the splitter plate was more effective than what occurs for the infinite circular cylinder: for AR=3, vortex shedding suppression occurred for all of the splitter plates tested (L/D≥1); for AR=5 and 7, vortex shedding suppression occurred for L/D≥1.5.     

Flow-induced vibrations of low aspect ratio rectangular membrane wings

An experimental study of a low aspect ratio rectangular membrane wing in a wind tunnel was conducted for a Reynolds number range of 2.4×10⁴–4.8×10⁴. Time-accurate measurements of membrane deformation were combined with the flow field measurements. Analysis of the fluctuating deformation reveals chordwise and spanwise modes, which are due to the shedding of leading-edge vortices as well as tip vortices. At higher angles of attack, the second mode in the chordwise direction becomes dominant as the vortex shedding takes place. The dominant frequencies of the membrane vibrations are similar to those of two-dimensional membrane airfoils. Measured frequency of vortex shedding from the low aspect ratio rigid wing suggests that membrane vibrations occur at the natural frequencies close to the harmonics of the wake instabilities. Vortex shedding frequency from rigid wings shows remarkably small effect of aspect ratio even when it is as low as unity.     

Dynamic Mode Decomposition of a Wing-Body Junction Flow

Junction flows are subject to an intense adverse pressure gradient and three-dimensional separation when encountering a wall-mounted obstacle. A dynamically rich horseshoe vortex system is formed in this region. In this study the junction flow at the interaction of a wing and a flat plate is investigated. The numerical modelling is carried out using the three-dimensional large eddy simulation (LES) approach at the Reynolds number Re = 1.15×10⁵ based on the wing’s maximum thickness T and the free stream velocity U_ref. The comparison with the experimental results shows that the numerical simulations fairly accurately reproduce the phenomenon under study. The dynamic mode decomposition (DMD) of the resolved flow field is employed to obtain the coherent dynamics of the flow. To clearly demonstrate the oscillation characteristics and the horseshoe vortex structures of junction flow the velocity field in the plane of symmetry is decomposed with eduction of two dominant DMDmodes. These two DMDmodes are reconstituted and developed, together with the mean flow mode to explain the latent dynamics. Mode 1 reveals the merging of the horseshoe vortices and mode 2 is responsible for the process of fission and stretching.     

Flow around four cylinders arranged in a square configuration

This paper presents an experimental study of the flow around four circular cylinders arranged in a square configuration. The Reynolds number was fixed at Re=8000, the pitch-to-diameter ratio between adjacent cylinders was varied from P/D=2 to 5 and the incidence angle was changed from α=0° (in-line square configuration) to 45° (diamond configuration) at an interval of 7.5°. The flow field was measured using digital Particle Image Velocimetry (PIV) to examine the vortex shedding characteristics of the cylinders, together with direct measurement of fluid dynamic forces (lift and drag) on each cylinder using a piezoelectric load cell. Depending on the pitch ratio, the flow could be broadly classified as shielding regime (P/D≤2), shear layer reattachment regime (2.5≤P/D≤3.5) and vortex impinging regime (P/D≥4). However, this classification is valid only in the case that the cylinder array is arranged nearly in-line with the free stream (α≈0°), because the flow is also sensitive to α. As α increases from 0° to 45°, each cylinder experiences a transition of vortex shedding pattern from a one-frequency mode to a two-frequency mode. The flow interference among the cylinders is complicated, which could be non-synchronous, quasi-periodic or synchronized with a definite phase relationship with other cylinders depending on the combined value of α and P/D. The change in vortex pattern is also reflected by some integral parameters of the flow such as force coefficients, power spectra and Strouhal numbers.     

Dean instability and secondary flow structure in curved rectangular ducts

The pressure driven, fully developed turbulent flow of incompressible viscous fluid (water) in 120° curved ducts of rectangular cross-section is investigated experimentally and numerically. Three different types of curved duct (A-CL, B-SL and C-IL) with continuously varying curvature conform to blade profile as the inner and outer curvature walls to simplify and guide the impeller design of pumps. After validating the numerical method against Particle Image Velocimetry (PIV) measurements, the flow development in the ducts is analyzed in detail by Computational Fluid Dynamics (CFD) for a wide range of Reynolds numbers (Re = 2.4 × 10⁴–1.4 × 10⁵) and aspect ratios (Ar > 1.0, =1.0 and <1.0). The results clearly depict the existence of multiple Dean vortices along the duct: while the axial velocity profile is more related to an inner Dean vortex (called split base vortex), the wall pressure is more influenced by the Dean vortex attached to the inner curvature wall (called ICW Dean vortex). The induced multiple Dean vortices and the secondary flow patterns in the duct cannot be faithfully predicted by using traditional techniques. Therefore, a new criterion based on the vortex core velocities is proposed. With this approach, the effects of Re, Cr and Ar on the Dean instabilities in curved ducts are carefully studied. Decreasing Re promotes the generation of Dean vortices closer to the duct inlet, a trend that is as opposed to laminar flow. In addition, a new pair of vortices called entrainment Dean vortex occurs near the outlet of the curved duct with Ar = 1.0, which has not been previously reported in the literature.