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.     

Visualization of a locally-forced separated flow over a backward-facing step

A laboratory water channel experiment was made of the separated flow over a backward-facing step. The flow was excited by a sinusoidally oscillating jet issuing from a separation line. The slit was connected to a cavity in which water was forced through a rigid pipe by a scotch-yoke system. The Reynolds number based on the step height (H) was fixed at Re _H =1200. The forcing frequency was varied in the range 0.305?St _H ?0.955 at the forcing amplitude A ₀=0.3. Time-averaged flow measurements were made by a LDV system, especially in the recirculating region behind the backward-facing step. To characterize the large-scale vortex evolution due to the local forcing, flow visualizations were performed by a dye tracer method with fluorescent ink. The vortex amalgamation process was captured at the effective forcing frequency (St _H =0.477) for laminar separation. This vortex merging process enhances flow mixing, which leads to the shortening of the reattachment length.     

A numerical study of flow over a confined backward-facing step

Numerical solutions using the SIMPLE algorithms for laminar flow over a backward-facing step are presented. Five differencing schemes were used: hybrid; quadratic upwind (QUICK); second-order upwind (SOUD); central-differencing and a novel scheme named second-order upwind biased (SOUBD). The SOUBD scheme is shown to be part of a family of schemes which include the central-differencing, SOUD and QUICK schemes for uniform grids. The results of the backward-facing step problem are presented and are compared with other numerical solutions and experimental data to evaluate the accuracy of the differencing schemes. The accuracy of the differencing schemes was ascertained by using uniform grids of various grid densities. The QUICK, SOUBD and SOUD schemes gave very similar accurate results. The hybrid scheme suffered from excessive diffusion except for the finest grids and the central-differencing scheme only converged for the finest grids.     

The entrance effect of laminar flow over a backward-facing step geometry

The study investigates the entrance effect for flow over a backward-facing step by comparing predictions that set the inlet boundary at various locations upstream of the sudden expansion. Differences are most significant in the sudden expansion region. If the geometry has an inlet channel, then shorter reattachment and separation lengths are predicted. Comparisons with experimental data indicate that better agreement is found using a long inlet channel, but only for low Reynolds numbers where the experimental error is less significant. For certain cases, predictions with a high expansion number are perturbed by the entrance effect more than low-expansion-number predictions; however, the effect is localized in the sudden expansion region. Channels with low expansion numbers always experience a greater entrance effect after some distance upstream and downstream of the sudden expansion. The boundary layer growth in the inlet channel was examined using a uniform inlet velocity profile. © 1997 John Wiley & Sons, Ltd.     

Finite element computation of a turbulent flow over a two-dimensional backward-facing step

This paper is devoted to the computation of turbulent flows by a Galerkin finite element method. Effects of turbulence on the mean field are taken into account by means of a k-? turbulence model. The wall region is treated through wall laws and, more specifically, Reichardt's law. An inlet profile for ? is proposed as a numerical treatment for physically meaningless values of k and ?. Results obtained for a recirculating flow in a two-dimensional channel with a sudden expansion in width are presented and compared with experimental values.     

Numerical simulations of laminar flow over a 3D backward-facing step

A numerical investigation of laminar flow over a three-dimensional backward-facing step is presented with comparisons with detailed experimental data, available in the literature, serving to validate the numerical results. The continuity constraint method, implemented via a finite element weak statement, was employed to solve the unsteady three-dimensional Navier–Stokes equations for incompressible laminar isothermal flow. Two-dimensional numerical simulations of this step geometry underestimate the experimentally determined extent of the primary separation region for Reynolds numbers Re greater than 400. It has been postulated that this disagreement between physical and computational experiments is due to the onset of three-dimensional flow near Re ≈ 400. This paper presents a full three-dimensional simulation of the step geometry for 100⩽ Re⩽ 800 and correctly predicts the primary reattachment lengths, thus confirming the influence of three-dimensionality. Previous numerical studies have discussed possible instability modes which could induce a sudden onset of three-dimensional flow at certain critical Reynolds numbers. The current study explores the influence of the sidewall on the development of three-dimensional flow for Re greater than 400. Of particular interest is the characterization of three-dimensional vortices in the primary separation region immediately downstream of the step. The complex interaction of a wall jet, located at the step plane near the sidewall, with the mainstream flow reveals a mechanism for the increasing penetration (with increasing Reynolds number) of three-dimensional flow structures into a region of essentially two-dimensional flow near the midplane of the channel. The character and extent of the sidewall-induced flow are investigated for 100⩽Re⩽ 800. © 1997 John Wiley & Sons, Ltd.     

Finite element computation of a turbulent flow over a two-dimensional backward-facing step

The present paper is devoted to the computation of turbulent flows by a Galerkin finite element method. Effects of turbulence on the mean field are taken into account by means of a (k-ε) turbulence model. The wall region is treated through wall laws and, more specifically, Reichardt's law. An inlet profile for ε is proposed as a numerical treatment for physically meaningless values of k and ε. Results obtained for a recirculating flow in a two-dimensional channel with a sudden expansion in width are presented and compared with experimental values.     

Three-dimensional coherent structure in a separated and reattaching flow over a backward-facing step

An experimental study was carried out to elucidate the large-scale vortical structure in a separated and reattaching flow over a backward-facing step. The Reynolds number based on the step height (H) was Re _H =33,000. The large-scale vortical structure was probed by means of three-dimensional velocity measurements performed at the recirculation zone (x/H=4.0) and the reattachment zone (x/H=7.5). A 32-channel microphone array extending in the streamwise and spanwise directions was used for sensing the wall pressure fluctuations. The relationship between the flow field and the relevant spatial mode of the pressure field was determined by examining the spatial box filtering. From the relevant spatial mode of the wall pressure fluctuations, a conditional averaging technique was employed to characterize the coherent structure. In addition, the cross-correlation between velocity and wall pressure fluctuations was calculated to identify the structure and the length scale of the large-scale vortex. The cross-correlation results revealed that the large-scale hairpin vortices have a three-dimensional structure, in agreement with previous findings. The present results clearly show the growth and downstream elongation of the hairpin vortices.List of symbols H step height, m - k turbulent kinetic energy, m²/s² - q freestream dynamic pressure, Pa - Re _H Reynolds number based on U ₀ and H,U _oH/ - U ₀ freestream velocity, m/s - U _c convection velocity, m/s - X ₀ streamwise coordinate of the measurement origin, m - x _R time mean reattachment length, m Greek symbols _p forward flow time fraction - cross-correlation coefficient - time delay, s - vorticity, m²/s     

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     

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.     

Effects of a turbulent wall jet on heat transfer over a non-confined backward-facing step

Experimental and numerical analysis of a turbulent wall jet on the heat transfer downstream of a non-confined backward-facing step are presented. Several configurations are studied to analyse the influence of the upstream flow and the height of the step. An infrared camera and a hot wire were used to visualize a temperature map and measure the instantaneous velocity, respectively. The main objective was to visualize and compare both the fluid flow and the heat transfer, by studying the skin friction coefficient C_f and the Nusselt number Nu_d, respectively. The latter is obtained by the calculation of the heat transfer coefficient, evaluated by inverse method. Both experimental data and numerical approach provide good agreement regarding the flow structure and thermal data for measuring the position and the value of characteristics scales in the recirculation zone. A correlation between the maximum heat transfer Nu_max and the maximum Reynolds number Re_max is presented. Similarities and differences are highlighted in the paper compared to confined configurations.     

Near-wall investigation of backward-facing step flows

The electrodiffusion technique has been used to investigate reattaching and recirculating flows behind a backward-facing step. The instantaneous wall shear rate vectors were determined using the current signals provided by a three-segment electrodiffusion probe. The near-wall extents of two counter-rotating recirculation zones located behind the step were determined under turbulent flow conditions in a water channel. The near-wall flow inside these recirculation zones was found to be very unsteady, with strong low-frequency fluctuations. The streamwise profiles of the wall shear stress were measured at several values of the Reynolds number and a high level of skin friction was obtained in the reverse-flow region. The strong dependence of the peak value of skin friction on the Reynolds number confirms the viscous-dominated character of the reverse flow appearing inside the recirculation zone. Received: 22 May 2000/Accepted: 28 March 2001     

Particle dispersion in a single-sided backward-facing step flow

The paper describes the particle dispersion in a single-sided backward-facing step flow. Particles of well-known sizes in the diameter range from 1 to 70 μm were suspended in an air flow and the particle motion over a step was measured by mean of a laser-Doppler anemometer. Thus, the local and integral flow quantities, i.e. the mean and turbulent velocity data could be measured precisely. In the experiments, monodispersed particle size distributions were used to exclude particle size related information ambiguity, known as triggering effects or size bias. The results of this study show qualitatively and quantitatively the difference in time-averaged particle dynamics for selected particle sizes in a backward-facing step flow. The experiments show, for different sizes, the changes in the particle velocity field in comparison with the velocity field of the continuous phase deduced from the 1 μm particles, and also imply the strong influences which different particle sizes have on flow data evaluation when size effects are not taken into account with particle-related optical measuring techniques.     

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).     

Large eddy simulation of separation control over a backward-facing step flow by suction

Separation control over a backward-facing step (BFS) flow by continuous suction was numerically investigated using the turbulence model of large eddy simulation (LES). The effect of suction control on the flow fields was scrutinised by altering the suction flow coefficient, and the results indicate that suction is not only very effective in shortening the reattachment length but also very influential in reducing the tangential velocity gradient and turbulence fluctuations of the reattached flows. With increasing increments of the absolute suction flow coefficient, the effect of suction control is more significant. Furthermore, the detailed flow fields (including the time-averaged stream and velocity fields) and turbulence characteristics (including the time-averaged resolved kinetic energy and RMS velocity) for the BFS models with or without suction are presented to discuss the mechanism of suction control. Comparisons of the time-averaged statistics between the numerical simulations and corresponding experiments are conducted, and it shows that the LES based on the dynamic kinetic energy subgrid-scale model (DKEM) can acquire exact results. Therefore, feasibility of the numerical methods to simulate suction-controlled models is validated.     

Open boundary condition symposium benchmark solution: Stratified flow over a backward-facing step

This paper presents a detailed numerical solution to a simplified version of two-dimensional stratified flow over a backward-facing step with a Froude number of 16/9, a Reynolds number of 800 and a Prandtl number of 1—one of the Open Boundary Condition Symposium test problems. The steady state solution was derived by integrating the time-dependent Boussinesq equations forward in time using a semi-implicit finite-element-based model on a 38400-element mesh. In addition to presenting the results derived on this grid, the paper also presents the results of a Richardson extrapolation calculation for a set of ‘key’ parameters. It is hoped that this solution can be used as a baseline to compare the performance of the various techniques discussed at the Open Boundary Condition Symposium.     

A test problem for outflow boundary conditions—flow over a backward-facing step

A numerical solution for steady incompressible flow over a two-dimensional backward-facing step is developed using a Galerkin-based finite element method. The Reynolds number for the simulations is 800. Computations are performed on an extended channel length to minimize the effect of the outflow boundary on the upstream recirculation zones. A thorough mesh refinement study is performed to validate the results. Extensive profile data at several channel locations are provided to allow future testing and evaluation of outflow boundary conditions.     

Control of backward-facing step flow using a flapping foil

The effect of oscillating a small foil in plunge on the reattachment of a separated shear layer in a two-dimensional backward-facing step flow has been studied using flow visualization and single component laser Doppler velocimetry (LDV) measurements. It has been shown that a jet instead of a wake is generated by the flapping action of the foil. Results indicate that this action induces strong mixing and entrainment when the foil is located within the recirculation flow region, thereby reducing the reattachment length by as much as 70%. Furthermore, it has been shown that the flapping foil is most effective in reducing the size of the separation zone when placed close to the wall and to the step. It is least effective when placed outside the separated shear layer or downstream of the reattachment zone. Received: 26 August 1999 / Accepted: 29 May 2001     

Particle image velocimetry measurements of a backward-facing step flow

Particle image velocimetry (PIV) measurements were carried out on a backward-facing step flow at a Reynolds number of Re_h=U_Xh/9=4,660 (based on step height and freestream velocity). In-plane velocity, out-of-plane vorticity, Reynolds stress and turbulent kinetic energy production measurements in the x-y and x-z planes of the flow are presented. Proper orthogonal decomposition was performed on both the fluctuating velocity and vorticity fields of the x-y plane PIV data using the method of snapshots. Low-order representations of the instantaneous velocity fields were reconstructed using the velocity modes. These reconstructions provided insight into the contribution that the various length scales make to the spatial distribution of mean and turbulent flow quantities such as Reynolds stress and turbulent kinetic energy production. Large scales are found to contribute to the Reynolds stresses and turbulent kinetic energy production downstream of reattachment, while small scales contribute to the intense Reynolds stresses in the vicinity of reattachment.     

Experimental study on the flow field behind a backward-facing step using a detonation-driven shock tunnel

The supersonic combustion RAM jet (SCRAM jet) engine is expected to be used in next-generation space planes and hypersonic airliners. To develop the engine, stabilized combustion in a supersonic flow field must be attained even though the residence time of flow is extremely short. A mixing process for breathed air and fuel injected into the supersonic flow field is therefore one of the most important design problems. Because the flow inside the SCRAM jet engine has high enthalpy, an experimental facility is required to produce the high-enthalpy flow field. In this study, a detonation-driven shock tunnel was built to produce a high-enthalpy flow, and a model SCRAM jet engine equipped with a backward-facing step was installed in the test section of the facility to visualize flow fields using a color schlieren technique and high-speed video camera. The fuel was injected perpendicularly to a Mach 3 flow behind the backward-facing step. The height of the step, the injection distance and injection pressure were varied to investigate the effects of the step on air/fuel mixing characteristics. The results show that the recirculation region increases as the fuel injection pressure increases. For injection behind the backward-facing step, mixing efficiency is much higher than with a flat plate. Also, the injection position has a significant influence on the size of the recirculation region generated behind the backward-facing step. The schlieren photograph and pressure histories measured on the bottom wall of the SCRAM jet engine model show that the fuel was ignited behind the step.Communicated by K. Takayama PACS 47.40.Ki