Effect of spanwise-varying local forcing on turbulent separated flow over a backward-facing step

An experimental study was made of turbulent separated flows over a backward-facing step. A local forcing was given to the separated flow by means of a sinusoidally oscillating jet issuing from a thin slit near the separation line. To produce a spanwise-varying local forcing at the separation edge, a banded thin tape covered the slit. Effects of the spanwise-varying local forcings on the separated flow were scrutinized by altering the spatially banded blocking width (w) and the open slit distance (g). An optimal value of w/g was sought, which led to the minimum reattachment length (x _R ). The effect of spanwise-varying local forcing on x _R was found to be slight compared to the case of two-dimensional forcing (w=0). The experiment was made at Re _H =33000 and A ₀=0.018 by changing the forcing frequency (0?St _H ?1.0).     

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.     

The manipulated transitional backward-facing step flow: an experimental and direct numerical simulation investigation

Results from a joint experimental and direct numerical simulation (DNS) investigation are presented for the flow over a backward-facing step manipulated by low-amplitude time-periodic (harmonic) blowing/suction excitation through a narrow slot at the edge of the step. For a Reynolds number of Re_h=3000 (based on step height, h, and inflow velocity, U_o) and for laminar inflow, a 33% reduction of the mean recirculation length (in comparison to the non-manipulated reference case) could be obtained with a forcing amplitude of the order of one per cent of U_o. Based on the momentum thickness, θ, of the incoming laminar boundary layer (at the edge of the step), the corresponding optimum Strouhal number is St_θ=f_optθ/U_o=0.012. From the experimental data it can be concluded that, in our flow case, the optimum frequency, f_opt=50 Hz , was the most amplified frequency in the transition-to-turbulence process of the separated laminar shear layer. Detailed comparison of the experimental data with data from the numerical simulation shows that DNS and experimental data agree up to second-order statistics. The joint experimental and numerical investigations exhibit a complementary nature in the sense that, on the one hand, the main advantage of the experiment was the relative ease with which a wide range of forcing parameters could be tested and, on the other hand, DNS could provide spatio-temporal details of the flow which could not be so easily obtained in the experiment.     

Control of turbulent separated flow over a backward-facing step by local forcing

An experimental study was made of the flow over a backward-facing step. Excitations were given to separated flow by means of a sinusoidally oscillating jet issuing from a thin slit near the separation line. The Reynolds number based on the step height (H) varied 13000 Re _H 33000. Effect of local forcing on the flow structure was scrutinized by altering the forcing amplitude (0 A ₀ 0.07) and forcing frequency (0 St _H 5.0). Small localized forcing near the separation edge enhanced the shear-layer growth rate and produced a large roll-up vortex at the separation edge. A large vortex in the shear layer gave rise to a higher rate of entrainment, which lead to a reduction in reattachment length as compared to the unforced flow. The normalized minimum reattachment length (x _r )_min/x _x0 was obtained at St 0.01. The most effective forcing frequency was found to be comparable to the shedding frequency of the separated shear layer.List of symbols a ₀ forcing amplitude=(Q _forced–Q _unforced)/U ₀ - AR aspect ratio=W/H - C _p wall-pressure coefficient=(P-P ₀)/(l/2) U ₀ ² - ER expansion ratio=(2H+H)/2H - f _f forcing frequency, Hz - f _s shedding frequency, Hz - g slit width = 1.0 ± 0.1 mm - H step height = 50 mm - P wall-static pressure, Pa - P ₀ wall-static pressure at x/H= -2.0, Pa - Q _forced total velocity measured at reference position for forced flow, m/s - Q _unforced total velocity measured at reference position for unforced flow, m/s - Re _H Reynolds number based on H and U ₀,= U ₀ H/v - St _H Reduced forcing frequency, Strouhal number = f _f H/U ₀ - St Reduced forcing frequency based on the momentum thickness = f _f /U ₀ - U, V streamwise and vertical time-mean velocity, m/s - u streamwise fluctuation velocity, m/s - U ₀ free-stream velocity, m/s - r.m.s. intensity of streamwise velocity fluctuation, m/s - x _r reattachment length, m - X _{r 0} reattachment length for A ₀ = 0, m - x, y, z distance of streamwise, vertical and spanwise respectively, m - W width of test section = 625 mm Greek symbols boundary-layer thickness, cm - ^* displacement thickness, cm - _p forward-flow time fraction - density of air for measurement, kg/m³ - v kinematic viscosity of air for measurement, m²/s - momentum thickness, cm     

PIV measurements of the wake behind a rotationally oscillating circular cylinder

The near-wake behind a circular cylinder undergoing rotational oscillatory motion with a relatively high forcing frequency has been investigated experimentally. Experiments were carried out varying the ratio of the forcing frequency f_f to the natural vortex shedding frequency f_n in the range of 0.0 (stationary) to 1.6 at an oscillation amplitude of θ_A=30° and Reynolds number of Re=4.14×10³. Depending on the frequency ratio (F_R=f_f /f_n), the near-wake flow could be divided into three regimes—non-lock-on (F_R=0.4), transition (F_R=0.8, 1.6) and lock-on (F_R=1.0) regimes—with markedly different flow structures. When the frequency ratio was less than 1.0 (F_R⩽1.0), the rotational oscillatory motion of the cylinder decreased the length of the vortex formation region and enhanced the mutual interaction between large-scale vortices across the wake centerline. The entrainment of ambient fluid seemed to play an important role in controlling the near-wake flow and shear-layer instability. In addition, strong vortex motion was observed throughout the near-wake region. The flow characteristics changed markedly beyond the lock-on flow regime (F_R=1.0) due to the high frequency forcing. At F_R=1.6, the high frequency forcing decreased the size of the large-scale vortices by suppressing the lateral extent of the wake. In addition, the interactions between the vortices shed from both sides of the cylinder were not so strong at this forcing frequency. As a consequence, the flow entrainment and momentum transfer into the wake center region were reduced. The turbulent kinetic energy was large in the region near the edge of the recirculation region, where the vortices shed from both sides of the cylinder cross the wake centerline for all frequency ratios except for the case of F_R=1.6. The temporally resolved quantitative flow information extracted in the present work is useful for understanding the effects of open-loop active flow control on the near-wake flow structure.     

Numerical solutions of pulsating flow and heat transfer characteristics in a channel with a backward-facing step

The incompressible laminar flow of air and heat transfer in a channel with a backward-facing step is studied for steady cases and for pulsatile inlet conditions. For steady flows the influence of the inlet velocity profile, the height of the step and the Reynolds number on the reattachment length is investigated. A parabolic entrance profile was used for pulsatile flow. It was found with amplitude of oscillation of one by Re=100 that the primary vortex breakdown through one pulsatile cycle. The wall shear rate in the separation zone varied markedly with pulsatile flows and the wall heat transfer remained relatively constant. The time-average pulsatile heat transfer at the walls was greater as with steady flow with the same mean Reynolds number.     

Multi-resolution analysis of the large-scale coherent structure in a turbulent separation bubble affected by an unsteady wake

Multi-resolution analysis (MRA) was applied to the large-scale coherent structure in a turbulent separation bubble affected by an unsteady wake. The unsteady wake was generated using a spoked-wheel type wake generator, which was installed in front of the separation bubble. The wake generator was rotated either clockwise (CW) or counter-clockwise (CCW) with a normalized passing frequency of St_H=0.2. The Reynolds number based on the half-thickness of the blunt body was Re_H=5600. To show the unsteady dynamic flow structures between the ‘cutting’ and ‘wrapping’ regimes, a MRA using the maximal overlap discrete wavelet transform (MODWT) was performed. This method enabled delineation of the coherent structure of the turbulent separation bubble through a scale-resolved analysis. Reconstruction of the flow field in combination with conditional averaging was attempted. Flapping motions as well as sawtooth movements of the unsteady separation bubble were analyzed using the MODWT. The unsteady wakes decayed faster in the system with CCW rotation than in that with CW rotation.     

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.     

Unsteady near wake of a flat disk normal to a wall

The unsteady wake of a flat disk (diameter D) located at a distance of H from a flat plate has been experimentally investigated at a Reynolds number Re _D  = 1.3 × 10⁵. Tests have been performed for a range of gap ratio (H/D), spanning from 0.3 to 1.75. The leading edge of the flat plate is either streamlined (elliptical) or blunt (square). These configurations have been studied with PIV, high speed PIV and multi-arrayed off-set fluctuating pressure measurements. The results show a progressive increase of the complexity of the flow and of the interaction as the gap ratio decreases. For large values of H/D (1.75), the interaction is weak and the power spectral densities (PSD) exhibit a strong peak associated with the vortex shedding events (St = 0.131) – St = fD/U _∞ is the Strouhal number. For lower values of H/D (0.75), the magnitude of the wall fluctuating pressure increases significantly. A large band contribution is associated with the unsteady wake structure and turbulence. A slight increase of the shedding frequency (St = 0.145) is observed. A critical value of the gap ratio (about 0.35) has been determined. Below this critical value, a three-dimensional separated region is observed and the natural vortex shedding process is very strongly altered. These changes induce a great modification of the fluctuating pressure at the wall. Each interaction reacts in a different way to perturbed upstream conditions. In particular, the disk is an overwhelming perturbation for the lowest H/D value studied here and the relative influence of the upstream turbulence on the wall fluctuating pressure below the near wake region is moderate.     

Vortex pairing in an axisymmetric jet using two-frequency acoustic forcing at low to moderate strouhal numbers

 Hot-wire measurement and multi-smoke wire flow visualization method are employed to study vortex pairing in the jet column mode under two-frequency forcing with controlled initial phase differences. For the range of 0.3<St _D <0.6, vortex pairing can be easily controlled by means of the fundamental and its subharmonic forcing with varying initial phase differences. As stable vortex pairing dominates, the variation of the subharmonic component with the initial phase difference changes from a sine shape to a cusp-like shape. The harmonics of the subharmonic also show similar trends. The detuning induces the amplitude and phase modulations of the u-signal in the time trace and the sideband growth in the spectra. The u-signal reflects the subharmonic variation with the initial phase difference in its envelope. For 0.6<St _D <0.9, non-pairing advection of vortices due to improper phase difference is sometimes observed under single-frequency forcing. In this case, vortex pairing can be made to occur by the addition of a subharmonic with very small amplitude. As the initial level of this subharmonic is increased, the onset position of vortex pairing moves upstream. In this range, the initial phase difference is not an effective parameter in controlling vortex pairing. Received: 22 May 1997 / Accepted: 16 October 1997     

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

Wall modeling for implicit large-eddy simulation and immersed-interface methods

We propose and analyze a wall model based on the turbulent boundary layer equations (TBLE) for implicit large-eddy simulation (LES) of high Reynolds number wall-bounded flows in conjunction with a conservative immersed-interface method for mapping complex boundaries onto Cartesian meshes. Both implicit subgrid-scale model and immersed-interface treatment of boundaries offer high computational efficiency for complex flow configurations. The wall model operates directly on the Cartesian computational mesh without the need for a dual boundary-conforming mesh. The combination of wall model and implicit LES is investigated in detail for turbulent channel flow at friction Reynolds numbers from Re _τ  = 395 up to Re _τ =100,000 on very coarse meshes. The TBLE wall model with implicit LES gives results of better quality than current explicit LES based on eddy viscosity subgrid-scale models with similar wall models. A straightforward formulation of the wall model performs well at moderately large Reynolds numbers. A logarithmic-layer mismatch, observed only at very large Reynolds numbers, is removed by introducing a new structure-based damping function. The performance of the overall approach is assessed for two generic configurations with flow separation: the backward-facing step at Re _h = 5,000 and the periodic hill at Re _H = 10,595 and Re _H = 37,000 on very coarse meshes. The results confirm the observations made for the channel flow with respect to the good prediction quality and indicate that the combination of implicit LES, immersed-interface method, and TBLE-based wall modeling is a viable approach for simulating complex aerodynamic flows at high Reynolds numbers. They also reflect the limitations of TBLE-based wall models.     

Array measurements of the surface pressure beneath a forced axi-symmetric separation bubble

An array of microphones is used to study the space–time characteristics of the wall-pressure field beneath a forced separation bubble downstream of an axi-symmetric backward-facing step. To excite the flow, an externally driven Helmholtz resonator is employed. A unique aspect of the present study is the utilization of an amplitude-modulated forcing scheme in order to avoid contamination of the measured hydrodynamic pressure fluctuations by acoustic radiation from the forcing device. The results lead to the hypothesis that the optimal forcing frequency is achieved when the forced disturbance originates near the center of the unforced separation bubble in the limit of very low levels of forcing. Moreover, a frequency–wavenumber spectrum analysis highlights the possibility for achieving separation control while minimizing potential acoustic radiation due to coupling between the forced disturbance and resonant modes of the underlying surface.     

Turbulent flow over a rough backward-facing step

This work characterizes the impacts of the realistic roughness due to deposition of foreign materials on the turbulent flows at surface transition from elevated rough-wall to smooth-wall. High resolution PIV measurements were performed in the streamwise-wall-normal (x–y) planes at two different spanwise positions in both smooth and rough backward-facing step flows. The experiment conditions were set at a Reynolds number of 3450 based on the free stream velocity U_∞ and the mean step height h, expansion ratio of 1.01, and the ratio of incoming boundary layer thickness to the step height, δ/h, of 8. The mean flow structures are observed to be modified by the roughness and they illustrate three-dimensional features in rough backward-facing step flows. The mean reattachment length X_r is significantly reduced by the roughness at one PIV measurement position while is slightly increased by the different roughness topography at the other measurement position. The mean velocity profiles at the reattachment point indicate that the studied roughness weakens the perturbation of the step to the incoming turbulent flow. Comparisons of Reynolds normal and shear stresses, productions of normal stresses, quadrant analysis of the instantaneous shear-stress contributing events, and mean spanwise vorticity reveal that the turbulence in the separated shear layer is reduced by the studied roughness. The results also indicate an earlier separation of the turbulent boundary layer over the current rough step, probably due to the adverse pressure gradient produced by the roughness topography even before the step.     

Flow visualization of supersonic laminar flow over a backward-facing step via NPLS

An experimental study on a supersonic laminar flow over a backward-facing step of 5 mm height was undertaken in a low-noise indraft wind tunnel. To investigate the fine structures of Ma = 3.0 and 3.8 laminar flow over a backward-facing step, nanotracer planar laser scattering was adopted for flow visualization. Flow structures, including supersonic laminar boundary layer, separation, reattachment, redeveloping turbulent boundary layer, expansion wave fan and reattachment shock, were revealed in the transient flow fields. In the Ma = 3.0 BFS (backward-facing step) flow, by measuring four typical regions, it could be found that the emergence of weak shock waves was related to the K–H (Kelvin–Helmholtz) vortex which appeared in the free shear layer and that the convergence of these waves into a reattachment shock was distinct. Based on large numbers of measurements, the structure of time-averaging flow field could be gained. Reattachment occurred at the location downstream from the step, about 7–7.5 h distance. After reattachment, the recovery boundary layer developed into turbulence quickly and its thickness increased at an angle of 4.6°. At the location of X = 14h, the redeveloping boundary layer was about ten times thicker than its original thickness, but it still had not changed into fully developed turbulence. However, in the Ma = 3.8 flow, the emergence of weak shock waves could be seen seldom, due to the decrease of expansion. The reattachment point was thought to be near X = 15h according to the averaging result. The reattachment shock was not legible, which meant the expansion and compression effects were not intensive.     

Investigation of flow mechanism of a robotic fish swimming by using flow visualization synchronized with hydrodynamic force measurement

When swimming in water by flapping its tail, a fish can overcome the drag from uniform flow and propel its body. The involved flow mechanism concerns 3-D and unsteady effects. This paper presents the investigation of the flow mechanism on the basis of a 3-D robotic fish model which has the typical geometry of body and tail with periodic flapping 2-freedom kinematical motion testing in the case of St = 0.78, Re = 6,600 and phase delay mode (φ = −75°), in which may have a greater or maximum propulsion (without consideration of the optimal efficiency). Using a special technique of dye visualization which can clearly show vortex sheet and vortices in detail and using the inner 3-component force balance and cable supporting system with the phase-lock technique, the 3-D flow structure visualized in the wake of fish and the hydrodynamic force measurement were synchronized and obtained. Under the mentioned flapping parameters, we found the key flow structure and its evolution, a pair of complex 3-D chain-shape vortex (S–H vortex-rings, S₁–H₁ and S₂–H₂, and their legs L₁ and L₂) flow structures, which attach the leading edge and the trailing edge, then shed, move downstream and outwards and distribute two anti-symmetric staggering arrays along with the wake of the fish model in different phase stages during the flapping period. It is different with in the case of St = 0.25–0.35. Its typical flow structure and evolution are described and the results prove that they are different from the viewpoints based on the investigation of 2-D cases. For precision of the dynamic force measurement, in this paper it was provided with the method and techniques by subtracting the inertial forces and the forces induced by buoyancy and gravity effect in water, etc. from original data measured. The evolution of the synchronized measuring forces directly matching with the flow structure was also described in this paper.     

DNS and LES of Passively Controlled Turbulent Backward-Facing Step Flow

A passive control approach (no external energy input) for an unsteady separated flow case was investigated numerically. A surface-mounted control fence was positioned upstream of a backward-facing step, and as an oncoming flow a thin and fully developed turbulent boundary layer with a thickness of δ/h = 0.8 was used. The objective of the passive control was to enhance the entrainment rate of the shear layer bounding the separation zone behind the step, thereby reducing the mean reattachment length,〈 X _r0 〉. Direct Numerical Simulations (DNS) and Large-Eddy Simulations (LES) at Re _h = 3000 (based on the step height, h, and the free stream velocity, U _∞) were carried out for the uncontrolled and the controlled flow case. The LES results were in good agreement with the DNS reference solutions. Adaptively controlled feedback simulations showed that a certain minimum distance between the step edge and the upstream position of the control fence is required to achieve a maximum reduction of the reattachment length.     

An experimental study of a spanwise structure around a reattachment region of a two-dimensional backward-facing step

 The flow field downstream of a two-dimensional backward-facing step is usually assumed to be independent of the direction along the span of the step. This assumption is made even though it is well known that the flow exhibits a three-dimensional vortex structure. This state of affairs is no doubt due to the lack of detailed information concerning the characteristics of the vortex structure. In this paper, we report our investigations of the flow structure around a reattachment region using an ultrasound velocity profiler to measure the spanwise velocity component as a function of the spanwise coordinate and time. The flow field is found to be very complex both in space and time. The low-frequency component of the spanwise velocity fluctuation becomes dominant in the near-wall region, with peaks in the power spectrum at frequencies fh/Uc=0.05 and fh/Uc=0.012. Using multiple ultrasound transducers, we also find that a streamwise vortex exists in the flow. Received: 20 March 2000 / Accepted: 15 January 2001 Published online: 29 November 2001     

Assessment of the organization of a turbulent separated and reattaching flow by measuring wall pressure fluctuations

The effect of local forcing on the organization of a turbulent separated and reattaching flow was assessed by measuring wall pressure fluctuations. Multi-arrayed microphones were installed on the surface to measure the simultaneous spatial and temporal wall pressure fluctuations. Local forcing at the separation edge was applied to the separated flow over a backward-facing step through a thin slit. The organization of the separated and reattaching flow was found to be greatest at the effective forcing frequency. The flow structure was diagnosed by analyzing several characteristics of the wall pressure fluctuations: the wall pressure fluctuation coefficients, wall pressure spectrum, wavenumber-frequency spectrum, coherence, cross-correlation, and multi-resolution autocorrelations of pressure fluctuations using the maximum overlap discrete wavelet transform and continuous wavelet transform. Features indicative of the amalgamation of vortices under the local forcing were observed; this amalgamation process accounted for the observed reduction of the reattachment length. Examination of the wall pressure fluctuations revealed that introduction of local forcing enhanced flapping motion as well as the streamwise and spanwise dispersions of vortical structures.