Experimental Investigation of Three-Dimensional Vortex Structures Downstream of Vortex Generators Over a Backward-Facing Step

An experimental investigation of vortex generators has been carried out in turbulent backward-facing step (BFS) flow. The Reynolds number, based on a freestream velocity U₀ = 10 m/s and a step height h = 30 mm, was Re_h = 2.0 × 10⁴. Low-profile wedge-type vortex generators (VGs) were implemented on the horizontal surface upstream of the step. High-resolution planar particle image velocimetry (2D-2C PIV) was used to measure the separated shear layer, recirculation region and reattachment area downstream of the BFS in a single field of view. Besides, time-resolved tomographic particle image velocimetry (TR-Tomo-PIV) was also employed to measure the flow flied of the turbulent shear layer downstream of the BFS within a three-dimensional volume of 50 × 50 × 10 mm³ at a sampling frequency of 1 kHz. The flow control result shows that time-averaged reattachment length downstream of the BFS is reduced by 29.1 % due to the application of the VGs. Meanwhile, the Reynolds shear stress downstream of the VGs is considerably increased. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) have been applied to the 3D velocity vector fields to analyze the complex vortex structures in the spatial and temporal approaches, respectively. A coherent bandwidth of Strouhal number 0.3 < St_h < 0.6 is found in the VG-induced vortices, and moreover, Λ-shaped three-dimensional vortex structures at St_h = 0.37 are revealed in the energy and dynamic approaches complementarily.     

PIV measurement of separated flow over a blunt plate with different chord-to-thickness ratios

The influence of the chord-to-thickness ratio (c/t) on the spatial characteristics of the separated shear layer over a blunt plate and the leading-edge vortices embedded in the separated shear layer was studied extensively using planar particle image velocimetry (PIV). Three systems corresponding to different shedding modes were chosen for the comparative study: c/t=3, 6 and 9. The Reynolds number based on the plate's thickness (t) was Re_t=1×10³. A gigapixel CCD camera was used to acquire images with a spatial resolution of 0.06t×0.06t in the measurement range of 9.5t×4.5t. Distributions of statistical quantities, such as the streamline pattern, streamwise velocity fluctuation intensity, shear stress and reverse flow intermittency, showed that the separated shear layer in the system with c/t=3 did not reattach to the plate's surface, while the near‐wake behind the trailing edge was highly unstable because the energetic leading-edge vortices were shed into the wake. The separated shear layer of the system with c/t=6 periodically reattached to the plate's surface, which resulted in intensified fluctuations of the near wake behind the trailing edge. In the longest system (c/t=9), the separated shear layer always reattached to the plate's surface far upstream from the trailing edge, which did not induce large fluctuations of the near wake. Furthermore, the proper orthogonal decomposition (POD) was extensively employed to filter the original velocity fields spatially to identify the large-scale vortices immersed in the separated shear layer easily. The distribution of the v-v correlation coefficients of the spatially filtered flow fields reflected the organized large-scale vortices in the three systems. The number of alternations of the positive and negative correlation coefficients across the flow field were determined to be 1, 2 and 3 for the systems with c/t=3, 6 and 9, respectively; this is in agreement with the shedding mode of each system. The distribution of the swirling strength of the separated shear layer accurately determined the positions and structures of the large-scale vortices formed above the plate surface.     

Flow field investigation in rotating rib-roughened channel by means of particle image velocimetry

The turbulent velocity field over the rib-roughened wall of an orthogonally rotating channel is investigated by means of two-dimensional particle image velocimetry (PIV). The flow direction is outward, with a bulk Reynolds number of 1.5 × 10⁴ and a rotation number ranging from 0.3 to 0.38. The measurements are obtained along the wall-normal/streamwise plane at mid-span. The PIV system rotates with the channel, allowing to measure directly the relative flow velocity with high spatial resolution. Coriolis forces affect the stability of the boundary layer and free shear layer. Due to the different levels of shear layer entrainment, the reattachment point is moved downstream (upstream) under stabilizing (destabilizing) rotation, with respect to the stationary case. Further increase in rotation number pushes further the reattachment point in stabilizing rotation, but does not change the recirculation length in destabilizing rotation. Turbulent activity is inhibited along the leading wall, both in the boundary layer and in the separated shear layer; the opposite is true along the trailing wall. Coriolis forces affect indirectly the production of turbulent kinetic energy via the Reynolds shear stresses and the mean shear. Two-point correlation is used to characterize the coherent motion of the separated shear layer. Destabilizing rotation is found to promote large-scale coherent motions and accordingly leads to larger integral length scales; on the other hand, the spanwise vortices created in the separating shear layer downstream of the rib are less organized and tend to be disrupted by the three-dimensional turbulence promoted by the rotation. The latter observation is consistent with the distributions of span-wise vortices detected in instantaneous flow realizations.     

Dual-plane PIV technique to determine the complete velocity gradient tensor in a turbulent boundary layer

Simultaneous dual-plane PIV experiments, which utilized three cameras to measure velocity components in two differentially separated planes, were performed in streamwise-spanwise planes in the log region of a turbulent boundary layer at a moderate Reynolds number (Re 1100). Stereoscopic data were obtained in one plane with two cameras, and standard PIV data were obtained in the other with a single camera. The scattered light from the two planes was separated onto respective cameras by using orthogonal polarizations. The acquired datasets were used in tandem with continuity to compute all 9 velocity gradients, the complete vorticity vector and other invariant quantities. These derived quantities were employed to analyze and interpret the structural characteristics and features of the boundary layer. Sample results of the vorticity vector are consistent with the presence of hairpin-shaped vortices inclined downstream along the streamwise direction. These vortices envelop low speed zones and generate Reynolds shear stress that enhances turbulence production. Computation of inclination angles of individual eddy cores using the vorticity vector suggests that the most probable inclination angle is 35° to the streamwise-spanwise plane with a resulting projected eddy inclination of 43° in the streamwise-wall-normal plane.

---

Flow past a cylinder: shear layer instability and drag crisis

Flow past a circular cylinder for Re=100 to 10⁷ is studied numerically by solving the unsteady incompressible two‐dimensional Navier–Stokes equations via a stabilized finite element formulation. It is well known that beyond Re ～ 200 the flow develops significant three‐dimensional features. Therefore, two‐dimensional computations are expected to fall well short of predicting the flow accurately at high Re. It is fairly well accepted that the shear layer instability is primarily a two‐dimensional phenomenon. The frequency of the shear layer vortices, from the present computations, agree quite well with the Re^0.67 variation observed by other researchers from experimental measurements. The main objective of this paper is to investigate a possible relationship between the drag crisis (sudden loss of drag at Re ～ 2 × 10⁵) and the instability of the separated shear layer. As Re is increased the transition point of shear layer, beyond which it is unstable, moves upstream. At the critical Reynolds number the transition point is located very close to the point of flow separation. As a result, the shear layer eddies cause mixing of the flow in the boundary layer. This energizes the boundary layer and leads to its reattachment. The delay in flow separation is associated with narrowing of wake, increase in Reynolds shear stress near the shoulder of the cylinder and a significant reduction in the drag and base suction coefficients. The spatial and temporal power spectra for the kinetic energy of the Re=10⁶ flow are computed. As in two‐dimensional isotropic turbulence, E(k) varies as k^?5/3 for wavenumbers higher than energy injection scale and as k^?3 for lower wavenumbers. The present computations suggest that the shear layer vortices play a major role in the transition of boundary layer from laminar to turbulent state. Copyright © 2004 John Wiley & Sons, Ltd.     

Three-dimensional wake structure measurement using a modified PIV technique

Three-dimensional vortical structures have been measured in a circular-cylinder wake using particle imaging velocimetry (PIV) for the Reynolds number range of 2×10³ to 1×10⁴. The PIV was modified, compared with the conventional one, in terms of its light sheet arrangement to capture reliably streamwise vortices. While in agreement with previous reports, the presently measured spanwise structures complement the data in the literature in the streamwise evolution of the near-wake spanwise vortex in size, strength, streamwise and lateral convection velocities, shedding new light upon vigorous interactions between oppositely signed spanwise structures. The longitudinal vortices display mushroom patterns in the (x, z)-plane in the immediate proximity to the cylinder. Their most likely inclination in the (x, y)-plane is inferred from the measurements in different (x, z)-planes. The longitudinal vortices in the (y, z)-plane show alternate change in sign, though not discernible at x/d > 15. They decay in the maximum vorticity and circulation rapidly from x/d = 5 to 10 and slowly for x/d > 10, and are further compared with the spanwise vortices in size, strength and rate of decay.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.     

Flow around a NACA0018 airfoil with a cavity and its dynamical response to acoustic forcing

Trapping of vortices in a cavity has been explored in recent years as a drag reduction measure for thick airfoils. If, however, trapping fails, then oscillation of the cavity flow may couple with elastic vibration modes of the airfoil. To examine this scenario, the effect of small amplitude vertical motion on the oscillation of the shear layer above the cavity is studied by acoustic forcing simulating a vertical translation of a modified NACA0018 profile. At low Reynolds numbers based on the chord (O(10⁴)), natural instability modes of this shear layer are observed for Strouhal numbers based on the cavity width of order unity. Acoustic forcing sufficiently close to the natural instability frequency induces a strong non-linear response due to lock-in of the shear layer. At higher Reynolds numbers (above 10⁵) for Strouhal number 0.6 or lower, no natural instabilities of the shear layer and only a linear response to forcing were observed. The dynamical pressure difference across the airfoil is then dominated by added mass effects, as was confirmed by numerical simulations.     

Near-wake flow structure downwind of a wind turbine in a turbulent boundary layer

Wind turbines operate in the surface layer of the atmospheric boundary layer, where they are subjected to strong wind shear and relatively high turbulence levels. These incoming boundary layer flow characteristics are expected to affect the structure of wind turbine wakes. The near-wake region is characterized by a complex coupled vortex system (including helicoidal tip vortices), unsteadiness and strong turbulence heterogeneity. Limited information about the spatial distribution of turbulence in the near wake, the vortex behavior and their influence on the downwind development of the far wake hinders our capability to predict wind turbine power production and fatigue loads in wind farms. This calls for a better understanding of the spatial distribution of the 3D flow and coherent turbulence structures in the near wake. Systematic wind-tunnel experiments were designed and carried out to characterize the structure of the near-wake flow downwind of a model wind turbine placed in a neutral boundary layer flow. A horizontal-axis, three-blade wind turbine model, with a rotor diameter of 13 cm and the hub height at 10.5 cm, occupied the lowest one-third of the boundary layer. High-resolution particle image velocimetry (PIV) was used to measure velocities in multiple vertical stream-wise planes (x–z) and vertical span-wise planes (y–z). In particular, we identified localized regions of strong vorticity and swirling strength, which are the signature of helicoidal tip vortices. These vortices are most pronounced at the top-tip level and persist up to a distance of two to three rotor diameters downwind. The measurements also reveal strong flow rotation and a highly non-axisymmetric distribution of the mean flow and turbulence structure in the near wake. The results provide new insight into the physical mechanisms that govern the development of the near wake of a wind turbine immersed in a neutral boundary layer. They also serve as important data for the development and validation of numerical models.     

Tomographic and time resolved PIV measurements on a finite cylinder mounted on a flat plate

Tomographic and time resolved PIV measurements were performed to examine the 3D flow topology and the flow dynamic above the upper surface of a low-aspect ratio cylinder at Re ≈ 1 ×  10⁵. This generic experiment is of fundamental interest because it represents flow features which are relevant to many applications such as laminar separation bubbles and turbulent reattachment. At Re  ≈ 1 × 10⁵, laminar separation bubbles arise on the side of the cylinder. Furthermore, on the top of the cylinder a separation with reattachment is of major interest. The tomographic PIV measurement, which allows to determine all three velocity components in a volume instantaneously, was applied to examine the flow topology and interaction between the boundary layer and wake structures on the top of the finite cylinder. In the instantaneous flow fields the tip vortices and the recirculation region becomes visible. However, it is also observed that the flow is quite unsteady due to the large separation occurring on the top of the cylinder. In order to study the temporal behaviour of the separation, time resolved PIV was applied. This technique allows capturing the dynamic processes in detail. The development of vortices in the separated shear layer is observed and in addition regions with different dominant frequencies are identified.     

Large eddy simulation of flows after a bluff body: Coherent structures and mixing properties

This paper performs large eddy simulations (LES) to investigate coherent structures in the flows after the Sydney bluff-body burner, a circular bluff body with an orifice at its center. The simulations are validated by comparison to existing experimental data. The Q function method is used to visualize the instantaneous vortex structures. Three kinds of structures are found, a cylindrical shell structure in the outer shear layer, a ring structure and some hairpin-like structures in the inner shear layer. An eduction scheme is employed to investigate the coherent structures in this flow. Some large streaks constituted by counter-rotating vortices are found in the outer shear layer and some well-organized strong structures are found in the inner shear layer. Finally, the influences of coherent structures on scalar mixing are studied and it is shown that scalar in the recirculation region is transported outward by coherent structures.     

Time resolved investigations on flow field and quasi wall shear stress of an impingement configuration with pulsating jets by means of high speed PIV and a surface hot wire array

The effects of jet pulsation on flow field and quasi wall shear stress of an impingement configuration were investigated experimentally. The excitation Strouhal number and amplitude were varied as the most influential parameters. A line-array with three submerged air jets, and a confining plate were used. The flow field analysis by means of time resolved particle image velocimetry shows that the controlled excitation can considerably affect the near-field flow of an impinging jet array. These effects are visualized as organization of the coherent flow structures. Augmentation of the Kelvin–Helmholtz vortices in the jet shear layer depends on the Strouhal number and pulsation magnitude and can be associated with pairing of small scale vortices in the jet. A total maximum of vortex strength was observed when exciting with Sr = 0.82 and coincident high amplitudes.Time resolved interaction between impinging vortices and impingement plate boundary layer due to jet excitation was verified by using an array of 5 μm surface hot wires. Corresponding to the global flow field modification due to periodic jet pulsation, the impact of the vortex rings on the wall boundary layer is highly influenced by the above mentioned excitation parameters and reaches a maximum at Sr = 0.82.     

Time-resolved and volumetric PIV measurements of a transitional separation bubble on an SD7003 airfoil

To comprehensively understand the effects of Kelvin–Helmholtz instabilities on a transitional separation bubble on the suction side of an airfoil regarding as to flapping of the bubble and its impact on the airfoil performance, the temporal and spatial structure of the vortices occurring at the downstream end of the separation bubble is investigated. Since the bubble variation leads to a change of the pressure distribution, the investigation of the instantaneous velocity field is essential to understand the details of the overall airfoil performance. This vortex formation in the reattachment region on the upper surface of an SD7003 airfoil is analyzed in detail at different angles of attack. At a Reynolds number Re _c < 100,000 the laminar boundary layer separates at angles of attack >4°. Due to transition processes, turbulent reattachment of the separated shear layer occurs enclosing a locally confined recirculation region. To identify the location of the separation bubble and to describe the dynamics of the reattachment, a time-resolved PIV measurement in a single light-sheet is performed. To elucidate the spatial structure of the flow patterns in the reattachment region in time and space, a stereo scanning PIV set-up is applied. The flow field is recorded in at least ten successive light-sheet planes with two high-speed cameras enclosing a viewing angle of 65° to detect all three velocity components within a light-sheet leading to a time-resolved volumetric measurement due to a high scanning speed. The measurements evidence the development of quasi-periodic vortex structures. The temporal dynamics of the vortex roll-up, initialized by the Kelvin–Helmholtz (KH) instability, is shown as well as the spatial development of the vortex roll-up process. Based on these measurements a model for the evolving vortex structure consisting of the formation of c-shape vortices and their transformation into screwdriver vortices is introduced.     

Wavy structures in compressible mixing layers

Semi-periodic structures namely inclined wavy structures(IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers(Mc=0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering(PLMS) technique. Two methods are developed to investigate the spatial distribution and geometry of IWS:(1) the dominant mode extraction(DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and(2) the phase tracking(PT) method, to identify the shape of IWS. The results suggest that pressure perturbations account for the formation of IWS in the initial mixing region and the joint efect of dilatation and coherent vortices enhances IWS in the welldeveloped region. The large transverse(cross-flow) scale of the IWS and their relation to coherent vortices(CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the efective tools to study the geometrical features of wavy structures in compressible shear flows.     

Experimental Investigation of Separated Shear Flow under Subharmonic Perturbations over a Backward-Facing Step

Subharmonic-perturbed shear flow downstream of a two-dimensional backward-facing step was experimentally investigated. The Reynolds number was Re_h = 2.0 ×10⁴, based on free-stream velocity and step height. Planar 2D-2C particle image velocimetry was employed to measure the separating and reattaching flow in the horizontal-vertical plane in the center position. The subharmonic perturbations were generated by an oscillating flap which was implemented over the step edge and driven by periodic Ampere force. The subharmonic frequency was 55 Hz as the half of the fundamental frequency of the turbulent shear layer. As a result of the subharmonic perturbations, the size of recirculation region behind the backward-facing step is reduced and the time-averaged reattachment length is 31.0% shorter than that of the natural flow. The evolution of vortices, including vortex roll-up, growth and breakdown process, is analyzed by using phase-averaging, cross-correlation function and proper orthogonal decomposition. It is found that Reynolds shear stress is considerably increased in which the vortices roll up and then break down further downstream. In particular, rapid growth of vortices based on the “step mode” occurs at approximate half of the recirculation region, caused by in interaction between the shear layer and the recirculation region. Furthermore, the coherent structures, which are represented by a phase-correlated POD mode pair, are reconstructed in phases in order to show regular patterns of the subharmonic-perturbed coherent structures.     

Turbulent fuel-air mixing study of jet in crossflow at different velocity ratios using LES

In order to study the mixing mechanism of fuel and air in gas turbine, large eddy simulation has been used to investigate the methane jet-in-crossflow with the velocity ratio (R) of 1.5 and 4. This study aims to explore the formation mechanism of vortices such as the hairpin vortices, hovering vortices and horseshoe vortices, the relationship between the fuel–air mixing and flow characteristics at different velocity ratios. The numerical methods in the present work are firstly validated with the experimental data in terms of mean and root mean square values of velocity. For R = 4, the shear layer vortices, horseshoe vortices, counter-rotating vortices pairs (CVP) and wake vortices can be observed, while the jet shear layer cannot be observed for R = 1.5. The hairpin vortices originating from the vortice-ring are lifted and shed from the downstream of the jet-outlet due to Kutta-Joukowski lift. The hairpin vortices are similar to CVP. The horseshoe vortices in R = 1.5 and 4 are formed due to the blockage of the jet (CH₄) and the crossflow (air) respectively, and its evolution is associated with the hovering vortices which only exist for R = 1.5. The uniform index and pr-obability density function are used for quantitative analysis of the mixing performance. The uniform index at X/D = 0 (fuel-inlet) and at X/D = 25 (outlet) are 0.033 and 0.335 for R = 1.5 and 0.130 and 0.047 for R = 4. For R = 4, the jet penetration is higher and the deflection angle of jet is smaller than that in case of R = 1.5. Higher R will provide more region for mixing, therefore uniform index is higher and the mixing is more uniform in the downstream.     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

PIV analysis of near-wake behind a sphere at a subcritical Reynolds number

The vortical structure of near-wake behind a sphere is investigated using a PIV technique in a circulating water channel at Re = 11,000. The measured velocity fields show a detailed vortical structure in the recirculation region such as recirculation vortices, reversed velocity zone, and out-of-plane vorticity distribution. The vorticity distribution of the sphere wake shows waviness in cross-sectional planes. The time-averaged turbulent structures are consistent with the visualized flow showing the onset of shear layer instability. The spatial distributions of turbulent intensities provide turbulent statistics for validating numerical predictions.     

Upward and downward inclination oil–water flows

The effect of upward (+5°, +10°) and downward (−5°) pipe inclinations on the flow patterns, hold up and pressure gradient during two-liquid phase flows was investigated experimentally for mixture velocities between 0.7 and 2.5 m/s and phase fractions between 10% and 90%. The investigations were performed in a 38 mm ID stainless steel test pipe with water and oil as test fluids. High-speed video recording and local impedance and conductivity probes were used to precisely identify the different flow patterns. In both positive and negative inclinations the dispersed oil-in-water regime extended to lower mixture velocities and higher oil fractions compared to horizontal flow. A new flow pattern, oil plug flow, appeared at both +5° and +10° inclination while the stratified wavy pattern disappeared at −5° inclination. The oil to water velocity ratio was higher for the upward than for the downward flows but in the majority of cases and all inclinations oil was flowing faster than water. At low mixture velocities the velocity ratio increased with oil fraction while it decreased at high velocities. The increase became more significant as the degree of inclination increased. The frictional pressure gradient in both upward and downward flows was in general lower than in horizontal flows while a minimum occurred at all inclinations at high mixture velocities during the transition from dispersed water-in-oil to dual continuous flow.     

Visualizations of large-scale vortices in flow about a blunt-faced flat plate

A new type of flow visualization method utilizing a smoke-wire, a high-speed camera with high framing rates and a laser light sheet was employed to delineate the unsteady processes of large-scale vortices in the separated shear layer about a blunt-faced flat plate at Re _H  = 560. The sequential images showed that the unsteady behavior of large-scale vortices in the separated shear layer varies as the shedding phases of large-scale vortices alter. Particularly, at a certain phase, a vortex-merging process between the two neighboring large-scale vortices took place. Received: 17 November 1998/Accepted: 1 November 1999