PIV measurements of turbulent flow in planar mixing layer

A turbulent mixing layer consists of two different flow types, i.e. shear layer (shear-flow turbulence) and free stream regions (nearly homogeneous turbulence). The inherent non-uniform seeding tracer distributions observed around the interfaces between the shear layer and two free stream regions usually lead to a difficulty in particle image velocimetry (PIV) measurements. A parametric study on the application of PIV to the measurement of velocity field in a planar mixing layer is made by means of six factors, including interrogation window size, aspect ratio of interrogation window, interrogation window offset, threshold of data validation, sharpening spatial filters (Prewitt and Sobel masks), and smoothing spatial filter (median mask). The objective of this study is to obtain accurate turbulent measurements in both mean and fluctuating velocities using PIV under an appropriate parametric setting. The optimal levels, which are trade-off in between the accuracy and fine spatial resolution of velocity field measurements, are determined with the aid of the Taguchi method. It is shown that the PIV measurements made with this optimal set of parameters are in good agreement with the measurements made by a two-component hot-wire anemometer. Case independency of the proposed optimal set of parameters on the flow condition of the mixing layer is validated through the applications to two additional tests under the different experimental conditions in changing solely either velocity ratio of high-speed to low-speed free stream velocities or Reynolds number.     

Momentum transport in a turbulent mixing layer

The turbulent momentum transport phenomena in a two-dimensional mixing layer are investigated numerically by a discrete vortex method. The numerical model and calculations are verified through a comparison with existing numerical simulations and experimental measurements. The main emphasis is placed on the exploration of the detailed time-dependent instantaneous local momentum fluctuations and on the comparison of numerical results with available experimental measurements. The current simulations confirm qualitatively the various trends in the turbulent momentum flux and fluctuating components of the velocity in the mixing layer found with several experimental results. The study shows that similarity exists in turbulent momentum quantities along the axial direction of the mixing layer. The calculations also show a definite correlation between the passage of a large-scale structure and a burst in the turbulent momentum flux. The probability density functions of the fluctuating quantities are shown to be mostly Gaussian-like, with only a few exceptions.     

Reconstruction of large coherent structures from SPIV measurements in a forced turbulent mixing layer

A procedure is suggested here for reconstructing the time variation of a three-dimensional (3D) coherent velocity field, based on applying Least Square Method to a very limited number of phase-locked measurements. The measurements are performed in a spanwise plane of a forced turbulent mixing layer by employing the stereo particle image velocimetry system. The forcing is produced by oscillating two- and three-dimensional (3D) flappers placed at the edge of the splitter plate. The feasibility and validity of the procedure for velocity field reconstruction are checked by using Hot Wire Measurements. Very clear patterns are observed of two types of coherent structures: spanwise K-H billows (rolls) and streamwise vortices. These are due to primary and secondary instabilities and their time evolution is presented.     

Three-component LDA measurements in a turbulent boundary layer

An experimental study of a three-dimensional, pressuredriven, attached turbulent boundary-layer flow was made at Mach 0.4. Both the mean velocities and the full Reynolds stress tensor were measured simultaneously by a three-component LDA system. Value of the resultant shear stress to turbulent kinetic energy ratio varied between 0.1 and 0.2 and did not remain constant across the boundary-layer. Slopes of the streamwise and azimuthal mixing length distributions in the wall region were around 0.4 and 1.2, respectively. Skew angle of the turbulent shear stress was larger than skew angle of the velocity gradient.     

RMS spanwise vorticity measurements in a turbulent boundary layer

Turbulent boundary layer measurements of the rms spanwise vorticity t_z with a four hot-wire probe are in reasonable agreement with direct numerical simulations and other published measurements at comparable Reynolds numbers. It is shown that a reasonable accurate approximation for _z can be obtained with only two parallel hot wires.     

A turbulent boundary layer over a two-dimensional rough wall

Particle image velocimetry (PIV) measurements and planar laser induced fluorescence (PLIF) visualizations have been made in a turbulent boundary layer over a rough wall. The wall roughness consisted of square bars placed transversely to the flow at a pitch to height ratio of λ/k = 11 for the PLIF experiments and λ/k = 8 and 16 for the PIV measurements. The ratio between the boundary layer thickness and the roughness height k/δ was about 20 for the PLIF and 38 for the PIV. Both the PLIF and PIV data showed that the near-wall region of the flow was populated by unstable quasi-coherent structures which could be associated to shear layers originating at the trailing edge of the roughness elements. The streamwise mean velocity profile presented a downward shift which varied marginally between the two cases of λ/k, in agreement with previous measurements and DNS results. The data indicated that the Reynolds stresses normalized by the wall units are higher for the case λ/k = 16 than those for λ/k = 8 in the outer region of the flow, suggesting that the roughness density effects could be felt well beyond the near-wall region of the flow. As expected the roughness disturbed dramatically the sublayer which in turn altered the turbulence production mechanism. The turbulence production is maximum at a distance of about 0.5k above the roughness elements. When normalized by the wall units, the turbulence production is found to be smaller than that of a smooth wall. It is argued that the production of turbulence is correlated with the form drag.     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

Structure of the turbulent boundary layer over a rod-roughened wall

Turbulent coherent structures near a rod-roughened wall are scrutinized by analyzing instantaneous flow fields obtained from direct numerical simulations (DNSs) of a turbulent boundary layer (TBL). The roughness elements used are periodically arranged two-dimensional spanwise rods, and the roughness height is k/δ = 0.05 where δ is the boundary layer thickness. The Reynolds number based on the momentum thickness is varied in the range Re_θ = 300–1400. The effect of surface roughness is examined by comparing the characteristics of the TBLs over smooth and rough walls. Although introduction of roughness elements onto the smooth wall affects the Reynolds stresses throughout the entire boundary layer when scaled by the friction velocity, the roughness has little effect on the vorticity fluctuations in the outer layer. Pressure-strain tensors of the transport equation for the Reynolds stresses and quadrant analysis disclose that the redistribution of turbulent kinetic energy of the rough wall is similar to that of the smooth wall, and that the roughness has little effect on the relative contributions of ejection and sweep motions in the outer layer. To elucidate the modifications of the near-wall vortical structure induced by surface roughness, we used two-point correlations, joint weighted probability density function, and linear stochastic estimation. Finally, we demonstrate the existence of coherent structures in the instantaneous flow field over the rod-roughened surface.     

Instantaneous velocity and pressure measurements in turbulent mixing layers

This paper presents a comprehensive comparison of the mean velocity and turbulence measurements from a four-hole pressure probe, also known as the Cobra probe, and an X-probe in plane mixing layers. The objective is to validate the measurement accuracy of the Cobra probe in a flow where the turbulence reaches high levels, but whose properties are well known. The comparison is made for the mean velocities, Reynolds stresses, triple products, and spectra, and demonstrates that the Cobra probe has reasonable accuracy for some of these quantities, such as the mean streamwise velocity and primary shear stress, but not for others, such as the mean normal velocity. The correlation of the pressure and the streamwise velocity, measured by the Cobra probe, behaves correctly in the potential flow. However, the correlation of the pressure and the cross-stream velocity, which appears in the transport equation for the turbulent kinetic energy, and the pressure redistribution term in the corresponding equation for the streamwise normal stress, are poorly measured.     

Flow field analysis of a turbulent boundary layer over a riblet surface

The near-wall flow structures of a turbulent boundary layer over a riblet surface with semi-circular grooves were investigated experimentally for the cases of drag decreasing (s ⁺=25.2) and drag increasing (s ⁺=40.6). One thousand instantaneous velocity fields over riblets were measured using the velocity field measurement technique and compared with those above a smooth flat plate. The field of view was 6.75 × 6.75 mm² in physical dimension, containing two grooves. Those instantaneous velocity fields were ensemble averaged to get turbulent statistics including turbulent intensities and turbulent kinetic energy. To see the global flow structure qualitatively, flow visualization was also carried out using the synchronized smoke-wire technique under the same experimental conditions. For the case of drag decreasing (s ⁺=25.2), most of the streamwise vortices stay above the riblets, interacting with the riblet tips frequently. The riblet tips impede the spanwise movement of the streamwise vortices and induce secondary vortices. The normalized rms velocity fluctuations and turbulent kinetic energy are small near the riblet surface, compared with those over a smooth flat plate. Inside the riblet valleys, these are sufficiently small that the increased wetted surface area of the riblets can be compensated. In addition, in the outer region (y ⁺ > 30), these values are almost equal to or slightly smaller than those for the smooth plate. For the case of drag increasing (s ⁺=40.6), however, most of the streamwise vortices stay inside the riblet valleys and contact directly with the riblet surface. The high-speed down-wash flow penetrating into the riblet valley interacts actively with the wetted riblet surface and increases the skin friction. The rms velocity fluctuations and turbulent kinetic energy have larger values compared with those over a smooth flat plate. Received: 24 March 1999/Accepted: 10 March 2000     

Spanwise effects on instabilities of compressible flow over a long rectangular cavity

The stability properties of two-dimensional (2D) and three-dimensional (3D) compressible flows over a rectangular cavity with length-to-depth ratio of \(L/D=6\) are analyzed at a free-stream Mach number of \(M_\infty =0.6\) and depth-based Reynolds number of \(Re_D=502\). In this study, we closely examine the influence of three-dimensionality on the wake mode that has been reported to exhibit high-amplitude fluctuations from the formation and ejection of large-scale spanwise vortices. Direct numerical simulation (DNS) and bi-global stability analysis are utilized to study the stability characteristics of the wake mode. Using the bi-global stability analysis with the time-averaged flow as the base state, we capture the global stability properties of the wake mode at a spanwise wavenumber of \(\beta =0\). To uncover spanwise effects on the 2D wake mode, 3D DNS are performed with cavity width-to-depth ratio of \(W/D=1\) and 2. We find that the 2D wake mode is not present in the 3D cavity flow with \(W/D=2\), in which spanwise structures are observed near the rear region of the cavity. These 3D instabilities are further investigated via bi-global stability analysis for spanwise wavelengths of \(\lambda /D=0.5{-}2.0\) to reveal the eigenspectra of the 3D eigenmodes. Based on the findings of 2D and 3D global stability analysis, we conclude that the absence of the wake mode in 3D rectangular cavity flows is due to the release of kinetic energy from the spanwise vortices to the streamwise vortical structures that develops from the spanwise instabilities.     

Velocity measurements for a turbulent nonseparated flow over solid waves

The laser-Doppler velocimeter was used to obtain measurements of the streamwise velocity over solid sinusoidal waves of small enough amplitude that a nonseparated flow existed. The measurements provide a critical test for Reynolds stress closure models since they are particularly sensitive to happenings in the viscous wall region (y ⁺ < 40), for which present theories are of uncertain accuracy. The results are compared with calculations that use an eddy viscosity model that successfully describes measurements of the wall shear stress along waves of small enough amplitude that a linear response is obtained. These calculations are in approximate agreement with measurements because they exactly account for inertia and viscous effects. However, there are significant differences which point to the inadequacy of turbulence models. In particular, non-linear effects and the amplitudes of the wave-induced velocity variations are underpredicted.     

Large-scale structures of turbulent flows over an open cavity

Low Mach number turbulent flows over an open cavity were studied to investigate the quantitative characteristics of large-scale vortical structures responsible for self-sustained oscillations. Wind tunnel experiments with particle image velocimetry (PIV) were conducted in the range of the ratio of cavity length (L) to depth (D), 1<L/D<4, when the incoming boundary layer is turbulent at Re_θ=830 and 1810. Self-sustained oscillation modes were classified by varying the conditions of L/D and Re_θ. The oscillation modes were consistent with the number of vortical structures existing between the leading and trailing edges of the cavity. Proper orthogonal decomposition (POD) was employed to the spatial distributions of vertical velocity correlations on the lip line of cavity geometry. By examining the conditionally averaged distributions of the correlation coefficients of POD, the spatial characteristics of large-scale vortical structures for self-sustained oscillations were examined.     

Self-sustained oscillations of turbulent flows over an open cavity

The mechanism of self-sustained oscillations in laminar cavity flows has been well characterized; however, the occurrence of self-sustained oscillations in turbulent cavity flows has only previously been characterized by direct observation of flows. Here, the quantitative characteristics of vortical structures in turbulent flows over an open cavity were determined, and then statistical properties were examined for evidence of self-sustained oscillations. Specifically, instantaneous velocity fields were measured using PIV and wall pressure fluctuations were determined from microphone data. Cavity geometries of L/D = 1 and 2, where L and D are the length and depth of the cavity, respectively, were used under conditions where the incoming boundary layer was turbulent at Re _θ  = 830. Statistical analyses were applied based on the instantaneous velocity fields of PIV data. The spatial distributions of vertical velocity correlations (v–v) showed alternating patterns that reflect the organized nature of the large-scale vortical structures corresponding to the modes of N = 2 for L/D = 1 and N = 3 for L/D = 2. These values were consistent with the numbers of vortical structures obtained from a modified version of Rossiter’s equation. Furthermore the numbers of vortical structures determined in the statistical analyses were consistently observed in instantaneous distributions of the swirling strength (λ _ci). The incoming turbulent boundary layer can give rise to the formation of large-scale vortical structures responsible for self-sustained oscillations.     

Applicability of the boundary-layer approximation to the calculation of a plane turbulent mixing layer

Calculation of a mixing layer-a classical problem in the theory of jet flows-is usually performed with a boundary-layer approximation. If the velocity of one of the flows is small, then even in the case of an incompressible fluid the angle of expansion of the mixing layer reaches 20?, and in the case of flows with different densities this angle is even larger. Therefore, there is some doubt about the smallness of the neglected terms which refer to differentiation in the flow direction. In the present paper the problem of the mixing of two semiinfinite flows is solved both with a boundary-layer approximation and by using a complete system of Reynolds equations (the molecular viscosity is assumed small compared to the turbulent viscosity; Re→∞).     

Simultaneous concentration and velocity field measurements in a shock-accelerated mixing layer

A novel technique to obtain simultaneous velocity and concentration measurements is applied to the Richtmyer–Meshkov instability. After acceleration by a Mach 2.2 shock wave, the interface between the two gases develops into a turbulent mixing layer. A time-separated pair of acetone planar laser-induced fluorescence images are processed to yield concentration and, through application of the Advection-Corrected Correlation Image Velocimetry technique, velocity fields. This is the first application of this technique to shock-accelerated flows. We show that when applied to numerical simulations, this technique reproduces the velocity field to a similar quality as particle image velocimetry. When applied to the turbulent mixing layer of the experiments, information about the Reynolds number and anisotropy of the flow is obtained.     

Quantitative schlieren measurements of coherent structures in a cavity shear layer

Quantitative flow-field data were obtained in a planar shear layer spanning an open cavity with an extension of the schlieren method. The technique is based on the measurement of light-intensity fluctuations in a real-time schlieren image. Data were collected using a fiber-optic sensor embedded in the imaging screen coupled to a photodetector. Time-resolved measurements of the instantaneous density gradient at a point in the two dimensional flow cross section were thus obtained. Detailed surveys were carried out with both the optical instrument as well as a hot wire at a Mach number of 0.25 and with the optical instrument alone at a Mach number of 0.6. A comparison of the results shows that the non-intrusive technique can accurately measure the growth rates of instability waves in the initial “linear” region of the shear layer. The density-gradient fluctuations measured at different locations (and times) were synchronized by using a microphone inside the cavity as a reference and integrated to yield profiles of the density fluctuations associated with the dominant large-scale structures in the shear layer. Such quantitative visualization is expected to clarify the mechanism of sound generation by shear-layer impingement at the cavity trailing edge and elucidate the nature of this sound source. Received: 28 December 1999/Accepted: 15 March 2000     

Simultaneous velocity and concentration field measurements of passive-scalar mixing in a confined rectangular jet

Simultaneous velocity and concentration fields in a confined liquid-phase rectangular jet with a Reynolds number based on the hydraulic diameter of 50,000 (or 10,000 based on the velocity difference between streams and the jet exit dimension) and a Schmidt number of 1,250 were obtained by means of a combined particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) system. Data were collected at the jet exit and six further downstream locations. The velocity and concentration field data were analyzed for flow statistics such as turbulent fluxes, turbulent viscosity and diffusivity, and turbulent Schmidt number (Sc _T ). The streamwise turbulent flux was found to be larger than the transverse turbulent flux, and the mean concentration gradient was not aligned with the turbulent flux vector. The average Sc _T was found to vary both in streamwise and in cross stream directions and had a mean value around 0.8, a value consistent with the literature. Spatial correlation fields of turbulent fluxes and concentration were then determined. The R _u′ϕ′ correlation was elliptical in shape with a major axis tilted downward with respect to the streamwise axis, whereas the R _v′ϕ′ correlation was an ellipse with a major axis aligned with the cross-stream direction. Negative regions of R _u′ϕ′ were observed in the outer streams, and these negatively correlated regions decayed with downstream distance and finally disappeared altogether. The R _ϕ′ϕ′ correlation field was found to be an ellipse with the major axis inclined at about 45° with respect to the streamwise direction. Linear stochastic estimation was used to interpret spatial correlation data and to determine conditional flow structures. It is believed that a vortex street formed near the splitter plate is responsible for the negatively correlated region observed in the R _u′ϕ′ spatial correlations of turbulent fluxes. A positive concentration fluctuation event was observed to correspond to a finger of nearly uniform concentration fluid reaching out into the outer stream, whereas a negative event corresponds to a pocket of nearly uniform fluid being entrained from the outer stream into the center jet region. Large-scale vortical structures were observed in the conditional velocity fields with an elliptical shape and a streamwise major axis. The growth of the structure size increased linearly initially but then grew more slowly as the flow transitioned toward channel flow. Support of this work was provided by the National Science Foundation through grants CTS-9985678 and CTS-0336435 and by the Dow Chemical Company. The author greatly acknowledge Charles Lipp at Dow Chemical and Ken Junk at Emerson Fisher for their valuable assistance in the design and construction of the flow system.